Changeset 18202
- Timestamp:
- 06/30/14 18:19:09 (11 years ago)
- Location:
- issm/trunk-jpl/examples/Pig
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/examples/Pig/CheatSheet.m
r18198 r18202 1 1 if perform(org,'ModelHO') 2 2 % Load Model 3 md = loadmodel(org,'Control_drag');3 mdHO = loadmodel('./Models/PIG.Control_drag'); 4 4 md.inversion.iscontrol=0; 5 5 … … 15 15 16 16 % Save Model 17 save model(org,md);17 save ./Models/PIG.ModelHO md; 18 18 end -
issm/trunk-jpl/examples/Pig/Pig.par
r18198 r18202 129 129 disp(' Set boundary conditions'); 130 130 md=SetMarineIceSheetBC(md); 131 md.basalforcings. melting_rate = zeros(md.mesh.numberofvertices,1);131 md.basalforcings.floatingice_melting_rate = zeros(md.mesh.numberofvertices,1); 132 132 md.thermal.spctemperature = [md.initialization.temperature;1]; %impose observed temperature on surface 133 133 md.masstransport.spcthickness = NaN*ones(md.mesh.numberofvertices,1); -
issm/trunk-jpl/examples/Pig/runme.m
r18198 r18202 1 1 %Which steps to be performed 2 steps=[2] ; 3 4 name_prefix='PIG.'; 2 steps=[1:5] ; 3 5 4 %Run Steps 6 org=organizer('repository','./Models','prefix',name_prefix,'steps',steps);7 5 8 6 % {{{ Mesh Generation #1 9 if perform(org,'Mesh_generation'),10 11 md.miscellaneous.name= strcat(name_prefix,'Mesh_generation');7 if any(steps==1) 8 9 md.miscellaneous.name='PIG.Mesh_generation'; 12 10 13 11 %Mesh parameters … … 70 68 71 69 % Convert x,y coordinates (Polar stereo) to lat/lon 72 [md.mesh.lat,md.mesh.long]=xy2ll(md.mesh.x,md.mesh.y,-1 ,39,71);73 74 % Save model 75 save model(org,md);70 [md.mesh.lat,md.mesh.long]=xy2ll(md.mesh.x,md.mesh.y,-1); 71 72 % Save model 73 save ./Models/PIG.Mesh_generation md; 76 74 end 77 75 % }}} 78 76 79 77 % {{{ Masks #2 80 if perform(org,'SetMask'),81 82 md = loadmodel( org,'Mesh_generation');78 if any(steps==2) 79 80 md = loadmodel('./Models/PIG.Mesh_generation'); 83 81 84 82 % Load SeaRISe dataset for Antarctica … … 104 102 105 103 % Save model 106 save model(org,md);104 save ./Models/PIG.SetMask md; 107 105 end 108 106 % }}} 109 107 110 108 % {{{ Parameterization #3 111 if perform(org,'Parameterization')112 113 md = loadmodel( org,'SetMask');109 if any(steps==3) 110 111 md = loadmodel('./Models/PIG.SetMask'); 114 112 md = parameterize(md,'./Pig.par'); 115 113 … … 118 116 119 117 % Save model 120 save model(org,md);118 save ./Models/PIG.Parameterization md; 121 119 end 122 120 % }}} 123 121 124 122 % {{{ Control Method #4 125 if perform(org,'Control_drag')126 127 md = loadmodel( org,'Parameterization');123 if any(steps==4) 124 125 md = loadmodel('./Models/PIG.Parameterization'); 128 126 129 127 % Control general … … 163 161 164 162 % Save model 165 save model(org,md);163 save ./Models/PIG.Control_drag md; 166 164 end 167 165 % }}} 168 166 169 167 % {{{ Plot #5 170 if perform(org,'PlotSSA')171 172 md = loadmodel( org,'Control_drag');168 if any(steps==5) 169 170 md = loadmodel('./Models/PIG.Control_drag'); 173 171 174 172 plotmodel(md,'nlines',2,'ncols',2,'unit#all','km','axis#all','equal',... … … 187 185 188 186 % {{{ HO #6 189 if perform(org,'ModelHO')187 if any(steps==6) 190 188 191 189 % Load Model … … 205 203 206 204 % {{{ Plot #7 207 if perform(org,'PlotHO')208 209 mdHO = loadmodel( org,'ModelHO');210 mdSSA = loadmodel( org,'Control_drag');205 if any(steps==7) 206 207 mdHO = loadmodel('./Models/PIG.ModelHO'); 208 mdSSA = loadmodel('./Models/PIG.Control_drag'); 211 209 212 210 basal=find(mdHO.mesh.vertexonbase);
Note:
See TracChangeset
for help on using the changeset viewer.