Changeset 19628


Ignore:
Timestamp:
10/13/15 14:32:31 (9 years ago)
Author:
schlegel
Message:

CHG: need to update smb

Location:
issm/trunk-jpl/examples/Greenland
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/examples/Greenland/Greenland.par

    r19163 r19628  
    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/Greenland/runme.m

    r19612 r19628  
    9292        smb = smb*md.materials.rho_freshwater/md.materials.rho_ice;
    9393        smb = [smb smb smb-1.0];
    94         md.surfaceforcings.mass_balance = [smb;1 10 20];
     94        md.smb.mass_balance = [smb;1 10 20];
    9595
    9696        %Set transient options, run for 20 years, saving every timestep
     
    102102        md.inversion.iscontrol=0;
    103103        md.transient.requested_outputs={'IceVolume','TotalSmb', ...
    104                                              'SurfaceforcingsMassBalance'};
     104                                             'SmbMassBalance'};
    105105        md.verbose=verbose('solution',true,'module',true,'convergence',true);
    106106
     
    124124        %Plot surface mass balance
    125125        surfmb=[]; for i=1:100; surfmb=[surfmb ...
    126                 md.results.TransientSolution(i).SurfaceforcingsMassBalance]; end
     126                md.results.TransientSolution(i).SmbMassBalance]; end
    127127        subplot(3,1,1); plot([0.2:0.2:20],mean(surfmb)); title('Mean Surface mass balance');
    128128
     
    181181        smb_mo = InterpFromMeshToMesh2d(index,x1(:),y1(:),smbmean(:),xi,yi);
    182182        smb = smb_mo*12/1000*md.materials.rho_freshwater/md.materials.rho_ice;
    183         md.surfaceforcings.mass_balance = [smb;1 ];
     183        md.smb.mass_balance = [smb;1 ];
    184184
    185185        %Set transient options, run for 20 years, saving every timestep
     
    191191        md.inversion.iscontrol=0;
    192192        md.transient.requested_outputs={'IceVolume','TotalSmb', ...
    193                 'SurfaceforcingsMassBalance'};
     193                'SmbMassBalance'};
    194194        md.verbose=verbose('solution',true,'module',true);
    195195
     
    272272                end
    273273        end
    274         md.surfaceforcings.mass_balance = ...
     274        md.smb.mass_balance = ...
    275275                [smb*12/1000*md.materials.rho_freshwater/md.materials.rho_ice; ...
    276276                [1/24:1/12:length(years_of_simulation)]];
     
    284284        md.inversion.iscontrol=0;
    285285        md.transient.requested_outputs={'IceVolume','TotalSmb', ...
    286                 'SurfaceforcingsMassBalance'};
     286                'SmbMassBalance'};
    287287        md.verbose=verbose('solution',true,'module',true);
    288288
Note: See TracChangeset for help on using the changeset viewer.