Changeset 14114
- Timestamp:
- 12/09/12 16:26:33 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/workshop/2012/Talks/09_PIG_Application/Exp_Par/Pig.par
r14045 r14114 36 36 md.geometry.bed = InterpFromGridToMesh(x1,y1,topg,md.mesh.x,md.mesh.y,0); 37 37 md.geometry.surface = InterpFromGridToMesh(x1,y1,usrf,md.mesh.x,md.mesh.y,0); 38 clear usrf, topg; 38 39 39 40 disp(' Constructing thickness'); … … 54 55 disp(' Interpolating temperatures'); 55 56 md.initialization.temperature=InterpFromGridToMesh(x1,y1,temp,md.mesh.x,md.mesh.y,0)+273.15+Temp_change; 57 clear temp; 56 58 57 59 disp(' Interpolating surface mass balance'); 58 60 md.surfaceforcings.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0); 59 61 md.surfaceforcings.mass_balance=md.surfaceforcings.mass_balance*md.materials.rho_water/md.materials.rho_ice; 62 clear smb; 60 63 61 64 disp(' Construct basal friction parameters'); … … 75 78 vx_obs=InterpFromGridToMesh(x2,y2,flipud(velx'),md.mesh.x,md.mesh.y,0); 76 79 vy_obs=InterpFromGridToMesh(x2,y2,flipud(vely'),md.mesh.x,md.mesh.y,0); 80 clear velx vely; 81 77 82 vel_obs=sqrt(vx_obs.^2+vy_obs.^2); 78 83 md.initialization.vx=vx_obs; … … 92 97 disp(' Set geothermal heat flux'); 93 98 md.basalforcings.geothermalflux=InterpFromGridToMesh(x1,y1,gflux,md.mesh.x,md.mesh.y,0); 99 clear gflux; 94 100 95 101 disp(' Set Pressure');
Note:
See TracChangeset
for help on using the changeset viewer.