Changeset 18202


Ignore:
Timestamp:
06/30/14 18:19:09 (11 years ago)
Author:
schlegel
Message:

CHG: update runmes for tutorials

Location:
issm/trunk-jpl/examples/Pig
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/examples/Pig/CheatSheet.m

    r18198 r18202  
    11if perform(org,'ModelHO')
    22  % Load Model
    3   md = loadmodel(org,'Control_drag');
     3  mdHO = loadmodel('./Models/PIG.Control_drag');
    44  md.inversion.iscontrol=0;
    55
     
    1515
    1616  % Save Model
    17   savemodel(org,md);
     17  save ./Models/PIG.ModelHO md;
    1818end
  • issm/trunk-jpl/examples/Pig/Pig.par

    r18198 r18202  
    129129disp('   Set boundary conditions');
    130130md=SetMarineIceSheetBC(md);
    131 md.basalforcings.melting_rate = zeros(md.mesh.numberofvertices,1);
     131md.basalforcings.floatingice_melting_rate = zeros(md.mesh.numberofvertices,1);
    132132md.thermal.spctemperature     = [md.initialization.temperature;1]; %impose observed temperature on surface
    133133md.masstransport.spcthickness    = NaN*ones(md.mesh.numberofvertices,1);
  • issm/trunk-jpl/examples/Pig/runme.m

    r18198 r18202  
    11%Which steps to be performed
    2 steps=[2] ;
    3 
    4 name_prefix='PIG.';
     2steps=[1:5] ;
     3
    54%Run Steps
    6 org=organizer('repository','./Models','prefix',name_prefix,'steps',steps);
    75
    86% {{{ Mesh Generation #1
    9 if perform(org,'Mesh_generation'),
    10 
    11         md.miscellaneous.name=strcat(name_prefix,'Mesh_generation');
     7if any(steps==1)
     8
     9        md.miscellaneous.name='PIG.Mesh_generation';
    1210
    1311        %Mesh parameters
     
    7068
    7169        % 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         savemodel(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;
    7674end
    7775% }}}
    7876
    7977% {{{ Masks #2
    80 if perform(org,'SetMask'),
    81        
    82         md = loadmodel(org,'Mesh_generation');
     78if any(steps==2)
     79
     80        md = loadmodel('./Models/PIG.Mesh_generation');
    8381
    8482        % Load SeaRISe dataset for Antarctica 
     
    104102       
    105103        % Save model
    106         savemodel(org,md);
     104        save ./Models/PIG.SetMask md;
    107105end
    108106% }}}
    109107
    110108% {{{ Parameterization #3
    111 if perform(org,'Parameterization')
    112 
    113         md = loadmodel(org,'SetMask');
     109if any(steps==3)
     110
     111        md = loadmodel('./Models/PIG.SetMask');
    114112        md = parameterize(md,'./Pig.par');
    115113
     
    118116       
    119117        % Save model
    120         savemodel(org,md);
     118        save ./Models/PIG.Parameterization md;
    121119end
    122120% }}}
    123121
    124122% {{{ Control Method #4
    125 if perform(org,'Control_drag')
    126 
    127         md = loadmodel(org,'Parameterization');
     123if any(steps==4)
     124
     125        md = loadmodel('./Models/PIG.Parameterization');
    128126
    129127        % Control general
     
    163161
    164162        % Save model
    165         savemodel(org,md);
     163        save ./Models/PIG.Control_drag md;
    166164end
    167165% }}}
    168166
    169167% {{{ Plot #5
    170 if perform(org,'PlotSSA')
    171 
    172         md = loadmodel(org,'Control_drag');
     168if any(steps==5)
     169
     170        md = loadmodel('./Models/PIG.Control_drag');
    173171
    174172        plotmodel(md,'nlines',2,'ncols',2,'unit#all','km','axis#all','equal',...
     
    187185
    188186% {{{ HO #6
    189 if perform(org,'ModelHO')
     187if any(steps==6)
    190188
    191189        % Load Model
     
    205203
    206204% {{{ Plot #7
    207 if perform(org,'PlotHO')
    208 
    209         mdHO = loadmodel(org,'ModelHO');
    210         mdSSA = loadmodel(org,'Control_drag');
     205if any(steps==7)
     206
     207        mdHO = loadmodel('./Models/PIG.ModelHO');
     208        mdSSA = loadmodel('./Models/PIG.Control_drag');
    211209
    212210        basal=find(mdHO.mesh.vertexonbase);
Note: See TracChangeset for help on using the changeset viewer.