Changeset 20532


Ignore:
Timestamp:
04/21/16 19:00:45 (9 years ago)
Author:
schlegel
Message:

CHG: update examples for new code

Location:
issm/trunk-jpl/examples
Files:
10 edited

Legend:

Unmodified
Added
Removed

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

    r19628 r20532  
    11clear all;
    2 steps=[1 2 3 4 5];
     2steps=[1];
    33
    44%Location of SeaRISE dataset
    … …  
    146146        lat  = ncread(ncbox,'lat');
    147147        lon  = ncread(ncbox,'lon');
    148         smbbox = ncread(ncbox,'SMB');
     148        smbbox = ncread(ncbox,'MassFlux');
    149149        [x1 y1]=ll2xy(lat,lon,+1,45,70);
    150150

  • issm/trunk-jpl/examples/ISMIP/CheatyRunme.m

    r19043 r20532  
    33%step 8 is specific to ISMIPF
    44
    5 steps=[1:7];
     5steps=[1];
    66
    77% parameter file to be used, choose between CheatyIsmipA.par or CheatyIsmipF.par

  • issm/trunk-jpl/examples/IceBridge/Greenland.par

    r19163 r20532  
    6262
    6363disp('   Interpolating surface mass balance');
    64 md.surfaceforcings.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0);
    65 md.surfaceforcings.mass_balance=md.surfaceforcings.mass_balance*md.materials.rho_water/md.materials.rho_ice;
     64md.smb.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0);
     65md.smb.mass_balance=md.smb.mass_balance*md.materials.rho_water/md.materials.rho_ice;
    6666
    6767disp('   Construct basal friction parameters');

  • issm/trunk-jpl/examples/IceBridge/Greenland_cheatsheet.par

    r19163 r20532  
    6666
    6767disp('   Interpolating surface mass balance');
    68 md.surfaceforcings.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0);
    69 md.surfaceforcings.mass_balance=md.surfaceforcings.mass_balance*md.materials.rho_water/md.materials.rho_ice;
     68md.smb.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0);
     69md.smb.mass_balance=md.smb.mass_balance*md.materials.rho_water/md.materials.rho_ice;
    7070
    7171disp('   Construct basal friction parameters');

  • issm/trunk-jpl/examples/IceBridge/Greenland_noOIB.par

    r19163 r20532  
    4040
    4141disp('   Interpolating surface mass balance');
    42 md.surfaceforcings.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0);
    43 md.surfaceforcings.mass_balance=md.surfaceforcings.mass_balance*md.materials.rho_water/md.materials.rho_ice;
     42md.smb.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0);
     43md.smb.mass_balance=md.smb.mass_balance*md.materials.rho_water/md.materials.rho_ice;
    4444
    4545disp('   Construct basal friction parameters');

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

    r18267 r20532  
    4646
    4747        md = setmask(md,'','');
    48         md = parameterize(md,'./Greenland_cheatsheet.par');
     48        md = parameterize(md,'./Greenland.par');
    4949        %md = parameterize(md,'./Greenland_noOIB.par');
    5050        md = setflowequation(md,'SSA','all');
    … …  
    108108        smb = InterpFromGridToMesh(x1,y1,smb',md.mesh.x,md.mesh.y,0)*1000/md.materials.rho_ice;
    109109        smb = [smb smb smb-1.0];
    110         md.surfaceforcings.mass_balance = [smb;1 10 20];
     110        md.smb.mass_balance = [smb;1 10 20];
    111111
    112112        %Set transient options, run for 20 years, saving every year
    … …  
    117117        %Additional options
    118118        md.inversion.iscontrol=0;
    119         md.transient.requested_outputs={'IceVolume','TotalSmb','SurfaceforcingsMassBalance'};
     119        md.transient.requested_outputs={'IceVolume','TotalSmb','SmbMassBalance'};
    120120        md.verbose=verbose('solution',true,'module',true,'convergence',true);
    121121
    … …  
    134134        plotmodel(md,'data',md.results.TransientSolution(end).Vel,'caxis',[1e-1 6000],...
    135135                'log', 10, 'title', 'Velocity (m/y)','gridded',1, ...
    136                 'data', md.results.TransientSolution(end).SurfaceforcingsMassBalance, ...
     136                'data', md.results.TransientSolution(end).SmbMassBalance, ...
    137137                'title', 'Surface mass balance (m/y)','gridded',1, ...
    138138                'data',md.results.TransientSolution(end).Thickness,...
    … …  
    146146        %Plot surface mass balance
    147147        surfmb=[]; for i=1:100; surfmb=[surfmb ...
    148                 md.results.TransientSolution(i).SurfaceforcingsMassBalance]; end
     148                md.results.TransientSolution(i).SmbMassBalance]; end
    149149        subplot(3,1,1); plot([0.2:0.2:20],mean(surfmb)); title('Mean Surface mass balance');
    150150

  • issm/trunk-jpl/examples/IceflowModels/EISMINT.par

    r18198 r20532  
    2929sb=10^-2/1000.; %m/yr/m
    3030rel=450.*1000.; %m
    31 md.surfaceforcings.mass_balance=min(smb_max,sb*(rel-radius));
     31md.smb.mass_balance=min(smb_max,sb*(rel-radius));
    3232
    3333disp('      creating velocities');

  • issm/trunk-jpl/examples/Jakobshavn/Jks.par

    r19163 r20532  
    3838
    3939disp('   Interpolating surface mass balance');
    40 md.surfaceforcings.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0);
    41 md.surfaceforcings.mass_balance=md.surfaceforcings.mass_balance*md.materials.rho_water/md.materials.rho_ice;
     40md.smb.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0);
     41md.smb.mass_balance=md.smb.mass_balance*md.materials.rho_water/md.materials.rho_ice;
    4242
    4343disp('   Construct basal friction parameters');

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

    r18201 r20532  
    1 steps=[4];
     1steps=[1];
    22
    33if any(steps==1)

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

    r18269 r20532  
    11%PIG Uncertainty Quantification Application
    2 steps=[1:3];
     2steps=[1];
    33
    44if any(steps==1)
    … …  
    140140
    141141        md.stressbalance.restol=10^-5; %tighten tolerances for UQ analyses
     142
     143        %Turn off verbosity
     144        md.verbose=verbose(0);
    142145       
    143146        %solve
    144147        md.qmu.isdakota=1; md.inversion.iscontrol=0;
    145148        md.cluster=generic('name',oshostname,'np',4);
     149
     150        %Dakota runs in parallel with a master/slave configuration.
     151        % At least 2 cpu's are needed to run the UQ
     152        md.qmu.params.evaluation_scheduling='master';
     153        md.qmu.params.processors_per_evaluation=md.cluster.np-1;
     154
    146155        md=solve(md,StressbalanceSolutionEnum,'overwrite','y');
    147156
    … …  
    217226        md.inversion.iscontrol=0;
    218227        md.cluster=generic('name',oshostname,'np',2);
     228
     229        %Dakota runs in parallel with a master/slave configuration.
     230        % At least 2 cpu's are needed to run the UQ
     231        md.qmu.params.evaluation_scheduling='master';
     232        md.qmu.params.processors_per_evaluation=md.cluster.np-1
     233
    219234        md.verbose=verbose('qmu',true);
    220235        md=solve(md,StressbalanceSolutionEnum,'overwrite','y');
Note: See TracChangeset for help on using the changeset viewer.