Changeset 21128


Ignore:
Timestamp:
08/10/16 13:45:09 (9 years ago)
Author:
glperez
Message:

CHG: Modified runmes in examples to use the any function in step conditional statements.
CHG: Implemented other examples in jenkins script. Still need to fill in steps for a couple of examples.

Location:
issm/trunk-jpl
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/examples/Greenland/runme.m

    r21057 r21128  
    210210
    211211if any(steps==8)
     212        disp('   Step 8: Plotting exercise');
     213        md = loadmodel('./Models/Greenland.HistoricTransient_200yr');
     214        figure
     215        surfmb=[]; for i=2:201; surfmb=[surfmb ...
     216                md.results.TransientSolution(i).SmbMassBalance]; end
     217        subplot(3,1,1); plot([1:200],mean(surfmb));
     218        title('Mean Surface mass balance');
     219        vel=[]; for i=2:201; vel=[vel md.results.TransientSolution(i).Vel]; end
     220        subplot(3,1,2); plot([1:200],mean(vel));
     221        title('Mean Velocity');
     222        volume=[]; for i=2:201; volume=[volume md.results.TransientSolution(i).IceVolume]; end
     223        subplot(3,1,3); plot([1:200],volume);
     224        title('Ice Volume'); xlabel('years');
    212225        disp('   Step 8: Plotting exercise');
    213226

  • issm/trunk-jpl/examples/Inversion/runme.m

    r21057 r21128  
    11step=1;
    2 if step==1
     2if any(step==1)
    33        %Generate observations
    44        md = model;
    … …  
    1313        save model1 md
    1414end
    15 if step==2
     15if any(step==2)
    1616        %Modify rheology, now constant
    1717        loadmodel('model1.mat');
    … …  
    2929        save model2 md
    3030end
    31 if step==3
     31if any(step==3)
    3232        %invert for ice rigidity
    3333        loadmodel('model2.mat');
    … …  
    4949                'data',md.results.StressbalanceSolution.Vel,'title','modeled velocities')
    5050end
    51 if step==4
     51if any(step==4)
    5252        %invert for ice rigidity
    5353        loadmodel('model2.mat');

  • issm/trunk-jpl/examples/Pig/runme.m

    r21057 r21128  
    11step=[1];
    22
    3 if step==1   %Mesh Generation #1
     3if any(step==1)   %Mesh Generation #1
    44        %Mesh parameters
    55        domain =['./DomainOutline.exp'];
    … …  
    5757end
    5858
    59 if step==2  %Masks #2
     59if any(step==2)  %Masks #2
    6060
    6161        md = loadmodel('./Models/PIG_Mesh_generation');
    … …  
    8686end
    8787
    88 if step==3  %Parameterization #3
     88if any(step==3)  %Parameterization #3
    8989
    9090        md = loadmodel('./Models/PIG_SetMask');
    … …  
    9898end
    9999
    100 if step==4  %Control Method #4
     100if any(step==4)  %Control Method #4
    101101
    102102        md = loadmodel('./Models/PIG_Parameterization');
    … …  
    141141end
    142142
    143 if step==5 %Plot #5
     143if any(step==5) %Plot #5
    144144
    145145        md = loadmodel('./Models/PIG_Control_drag');
    … …  
    158158end
    159159
    160 if step==6  %Higher-Order #6
     160if any(step==6)  %Higher-Order #6
    161161
    162162        % Load Model
    … …  
    174174end
    175175
    176 if step==7  %Plot #7
     176if any(step==7)  %Plot #7
    177177
    178178        mdHO = loadmodel('./Models/PIG_ModelHO');

  • issm/trunk-jpl/examples/PigSensitivity/runme.m

    r21057 r21128  
    11step=[1];
    22
    3 if step==1 %Transient Run #1
     3if any(step==1) %Transient Run #1
    44
    55        md = loadmodel('../Pig/Models/PIG_Control_drag');       
    … …  
    2626end
    2727
    28 if step==2 %High Melt #2
     28if any(step==2) %High Melt #2
    2929        md = loadmodel('./Models/PIG_Transient');       
    3030
    … …  
    4141end
    4242
    43 if step==3 %Ice Front retreat
     43if any(step==3) %Ice Front retreat
    4444        md = loadmodel('./Models/PIG_Transient');       
    4545
    … …  
    6060end
    6161
    62 if step==4 %High surface mass balance #3
     62if any(step==4) %High surface mass balance #3
    6363        %Load model
    6464

  • issm/trunk-jpl/jenkins/jenkins.sh

    r21127 r21128  
    211211
    212212#matlab tests
     213# {{{
    213214if [ $MATLAB_TEST -eq 1 ]; then
    214 #Launch all tests on different cpus {{{
     215#Launch all tests on different cpus
    215216for (( i=1;i<=$NUMCPUS_RUN;i++ ))
    216217do
    … …  
    286287#filter out windows characters:
    287288cat matlab_log.log | tr -cd '\11\12\40-\176' > matlab_log.log2 && mv matlab_log.log2 matlab_log.log
    288 #}}}
    289 fi
     289fi
     290# }}}
    290291
    291292#python tests
     293# {{{
    292294if [ $PYTHON_TEST -eq 1 ]; then
    293 #Launch all tests on different cpus {{{
     295#Launch all tests on different cpus
    294296PYTHON_START_TIME=$(timer);
    295297export PYTHONSTARTUP=$ISSM_DIR/src/m/dev/devpath.py
    … …  
    318320        cat python_log$i.log >> python_log.log
    319321done
    320 #}}}
    321 fi
    322 
     322fi
     323# }}}
     324
     325# Examples Test
     326# {{{
    323327# This test will allow us to check on the status of the examples.
    324328if [ $EXAMPLES_TEST -eq 1 ];
    … …  
    349353                                $SED -i.bak 's/steps=\[1\];/steps=\[1:8\];\n\ntry\n/' $FILE
    350354                                $SED -i.bak "s/if any(steps==8)/&${STEP_EIGHT}/" $FILE
     355                        elif [ "${dir}" == "./IceBridge" ];     
     356                        then
     357                                $SED -i.bak 's/steps=\[1\];/steps=\[1:5\];\n\ntry\n/' $FILE
     358                        elif [ "${dir}" == "./IceFlowModels" ];
     359                        then
     360                                # Almost nothing to this example
     361                                $SED -i.bak '1 s/^.*$/try\n\n&/' $FILE
     362                        elif [ "${dir}" == "./Inversion" ];     
     363                        then
     364                                $SED -i.bak 's/steps=\[1\];/steps=\[1:4\];\n\ntry\n/' $FILE
     365                        elif [ "${dir}" == "./Jakobshavn" ];   
     366                        then
     367                                $SED -i.bak 's/steps=\[1\];/steps=\[1:4\];\n\ntry\n/' $FILE
     368                        elif [ "${dir}" == "./Jakobshavn" ];   
     369                        then
     370                                $SED -i.bak 's/steps=\[1\];/steps=\[1:4\];\n\ntry\n/' $FILE
     371                        elif [ "${dir}" == "./Pig" ];   
     372                        then
     373                                # Step 6 is needed
     374                                STEP_SIX="\n disp('Needs work!')"
     375                                $SED -i.bak 's/steps=\[1\];/steps=\[1:7\];\n\ntry\n/' $FILE
     376                                $SED -i.bak "s/if any(steps==6)/&${STEP_SIX}/" $FILE
     377                        elif [ "${dir}" == "./PigSensitivity" ];       
     378                        then
     379                                # Step 4 is needed
     380                                STEP_FOUR="\n disp('Needs work!')"
     381                                $SED -i.bak 's/steps=\[1\];/steps=\[1:4\];\n\ntry\n/' $FILE
     382                                $SED -i.bak "s/if any(steps==6)/&${STEP_FOUR}/" $FILE
     383                        elif [ "${dir}" == "./SquareIceShelf" ];       
     384                        then
     385                                # Almost nothing to this example
     386                                $SED -i.bak '1 s/^.*$/try\n\n&/' $FILE
     387                        elif [ "${dir}" == "./UncertaintyQuantification" ];     
     388                        then
     389                                $SED -i.bak 's/steps=\[1\];/steps=\[1:7\];\n\ntry\n/' $FILE
    351390                        else
    352391                                echo "Not implemented yet!"
    … …  
    375414        done
    376415fi
     416# }}}
    377417
    378418#process logs to be junit compatible
Note: See TracChangeset for help on using the changeset viewer.