Changeset 24864


Ignore:
Timestamp:
05/15/20 20:28:17 (5 years ago)
Author:
Mathieu Morlighem
Message:

CHG: renaming groundedice_levelset -> ocean_levelset (negative if ocean present, positive outside

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

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/examples/AMR/mismip.par

    r23096 r24864  
    44md.geometry.surface=10*ones(md.mesh.numberofvertices,1);
    55md.geometry.thickness=md.geometry.surface-md.geometry.base;
    6 md.mask.groundedice_levelset=-1*ones(md.mesh.numberofvertices,1);
     6md.mask.ocean_levelset=-1*ones(md.mesh.numberofvertices,1);
    77
    88% creating basal drag

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

    r22818 r24864  
    5050        md.esa.love_l = love_numbers('l','CF'); md.esa.love_l(nlove+1:end)=[];
    5151
    52         md.mask.groundedice_levelset = ones(md.mesh.numberofvertices,1);
     52        md.mask.ocean_levelset = ones(md.mesh.numberofvertices,1);
    5353        md.mask.ice_levelset = ones(md.mesh.numberofvertices,1);
    5454        pos=find(md.esa.deltathickness~=0);

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

    r22818 r24864  
    6363
    6464        md.mask.ice_levelset = -ones(md.mesh.numberofvertices,1);
    65         md.mask.groundedice_levelset = ones(md.mesh.numberofvertices,1);
     65        md.mask.ocean_levelset = ones(md.mesh.numberofvertices,1);
    6666
    6767        di=md.materials.rho_ice/md.materials.rho_water;

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

    r20745 r24864  
    4545disp('   Construct basal friction parameters');
    4646md.friction.coefficient=30*ones(md.mesh.numberofvertices,1);
    47 pos=find(md.mask.groundedice_levelset<0);
     47pos=find(md.mask.ocean_levelset<0);
    4848md.friction.coefficient(pos)=0; %no friction applied on floating ice
    4949md.friction.p=ones(md.mesh.numberofelements,1);

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

    r21153 r24864  
    4545        % plot the given mask #md.mask to locate the field
    4646        %->
    47         plotmodel(md,'data',md.mask.groundedice_levelset);
     47        plotmodel(md,'data',md.mask.ocean_levelset);
    4848        % save the given model
    4949        %->

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

    r20746 r24864  
    6868disp('   Construct basal friction parameters');
    6969md.friction.coefficient=30*ones(md.mesh.numberofvertices,1);
    70 pos=find(md.mask.groundedice_levelset<0);
     70pos=find(md.mask.ocean_levelset<0);
    7171md.friction.coefficient(pos)=0; %no friction applied on floating ice
    7272md.friction.p=ones(md.mesh.numberofelements,1);

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

    r20532 r24864  
    4545disp('   Construct basal friction parameters');
    4646md.friction.coefficient=30*ones(md.mesh.numberofvertices,1);
    47 pos=find(md.mask.groundedice_levelset<0);
     47pos=find(md.mask.ocean_levelset<0);
    4848md.friction.coefficient(pos)=0; %no friction applied on floating ice
    4949md.friction.p=ones(md.mesh.numberofelements,1);

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

    r20746 r24864  
    7272disp('   Construct basal friction parameters');
    7373md.friction.coefficient=30*ones(md.mesh.numberofvertices,1);
    74 pos=find(md.mask.groundedice_levelset<0);
     74pos=find(md.mask.ocean_levelset<0);
    7575md.friction.coefficient(pos)=0; %no friction applied on floating ice
    7676md.friction.p=ones(md.mesh.numberofelements,1);

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

    r20532 r24864  
    1111disp('      creating drag');
    1212md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1);
    13 md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
     13md.friction.coefficient(find(md.mask.ocean_levelset<0.))=0.;
    1414md.friction.p=ones(md.mesh.numberofelements,1);
    1515md.friction.q=ones(md.mesh.numberofelements,1);

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

    r23137 r24864  
    3131%Friction
    3232md.friction.coefficient=20*ones(md.mesh.numberofvertices,1);
    33 md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
     33md.friction.coefficient(find(md.mask.ocean_levelset<0.))=0.;
    3434md.friction.p=ones(md.mesh.numberofelements,1);
    3535md.friction.q=ones(md.mesh.numberofelements,1);

  • issm/trunk-jpl/examples/Inversion/Square.par

    r18198 r24864  
    1212disp('      creating drag');
    1313md.friction.coefficient=200*ones(md.mesh.numberofvertices,1);
    14 md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
     14md.friction.coefficient(find(md.mask.ocean_levelset<0.))=0.;
    1515md.friction.p = ones(md.mesh.numberofelements,1);
    1616md.friction.q = ones(md.mesh.numberofelements,1);

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

    r20532 r24864  
    4343disp('   Construct basal friction parameters');
    4444md.friction.coefficient=30*ones(md.mesh.numberofvertices,1);
    45 md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
     45md.friction.coefficient(find(md.mask.ocean_levelset<0.))=0.;
    4646md.friction.p=ones(md.mesh.numberofelements,1);
    4747md.friction.q=ones(md.mesh.numberofelements,1);

  • issm/trunk-jpl/examples/LcurveAnalysis/Square.par

    r23752 r24864  
    1212disp('      creating drag');
    1313md.friction.coefficient=200*ones(md.mesh.numberofvertices,1);
    14 md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
     14md.friction.coefficient(find(md.mask.ocean_levelset<0.))=0.;
    1515md.friction.p = ones(md.mesh.numberofelements,1);
    1616md.friction.q = ones(md.mesh.numberofelements,1);

  • issm/trunk-jpl/examples/Pig/Pig.par

    r20764 r24864  
    3131
    3232%Get the node numbers of floating nodes
    33 pos=find(md.mask.groundedice_levelset<0);
     33pos=find(md.mask.ocean_levelset<0);
    3434
    3535%apply a flotation criterion on the precedingly defined nodes and
    … …  
    102102
    103103%no friction applied on floating ice
    104 pos=find(md.mask.groundedice_levelset<0);
     104pos=find(md.mask.ocean_levelset<0);
    105105md.friction.coefficient(pos)=0;
    106106md.groundingline.migration='SubelementMigration';

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

    r21607 r24864  
    7272
    7373        %fill in the md.mask structure
    74         md.mask.groundedice_levelset=groundedice; %ice is grounded for mask equal one
     74        md.mask.ocean_levelset=groundedice; %ice is grounded for mask equal one
    7575        md.mask.ice_levelset=-1*ones(md.mesh.numberofvertices,1);%ice is present when negatvie
    7676
    7777        %ploting
    78         plotmodel(md,'data',md.mask.groundedice_levelset,'title','grounded/floating','data',md.mask.ice_levelset,'title','ice/no-ice')
     78        plotmodel(md,'data',md.mask.ocean_levelset,'title','grounded/floating','data',md.mask.ice_levelset,'title','ice/no-ice')
    7979       
    8080        % Save model

  • issm/trunk-jpl/examples/Pig2/Pig.par

    r23388 r24864  
    3131
    3232%Get the node numbers of floating nodes
    33 pos=find(md.mask.groundedice_levelset<0);
     33pos=find(md.mask.ocean_levelset<0);
    3434
    3535%apply a flotation criterion on the precedingly defined nodes and
    … …  
    102102
    103103%no friction applied on floating ice
    104 pos=find(md.mask.groundedice_levelset<0);
     104pos=find(md.mask.ocean_levelset<0);
    105105md.friction.coefficient(pos)=0;
    106106md.groundingline.migration='SubelementMigration';

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

    r24311 r24864  
    6060
    6161        %fill in the md.mask structure
    62         md.mask.groundedice_levelset=groundedice; %ice is grounded for mask equal one
     62        md.mask.ocean_levelset=groundedice; %ice is grounded for mask equal one
    6363        md.mask.ice_levelset=-1*ones(md.mesh.numberofvertices,1);%ice is present when negatvie
    6464
    65         plotmodel(md,'data',md.mask.groundedice_levelset,'title','grounded/floating','data',md.mask.ice_levelset,'title','ice/no-ice')
     65        plotmodel(md,'data',md.mask.ocean_levelset,'title','grounded/floating','data',md.mask.ice_levelset,'title','ice/no-ice')
    6666       
    6767        save ./Models/PIG_SetMask md;
    … …  
    100100        md.inversion.min_parameters=md.materials.rheology_B;
    101101        md.inversion.max_parameters=md.materials.rheology_B;
    102         pos = find(md.mask.groundedice_levelset<0);
     102        pos = find(md.mask.ocean_levelset<0);
    103103        md.inversion.min_parameters(pos) = cuffey(273);
    104104        md.inversion.max_parameters(pos) = cuffey(200);
    … …  
    111111        % Solve
    112112        md.cluster=generic('name',oshostname,'np',2);
    113         mds=extract(md,md.mask.groundedice_levelset<0);
     113        mds=extract(md,md.mask.ocean_levelset<0);
    114114        mds=solve(mds,'Stressbalance');
    115115

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

    r21129 r24864  
    1212        md.transient.isthermal=0;
    1313       
    14         pos=find(md.mask.groundedice_levelset<0);
     14        pos=find(md.mask.ocean_levelset<0);
    1515        md.basalforcings.groundedice_melting_rate=zeros(md.mesh.numberofvertices,1);
    1616        md.basalforcings.floatingice_melting_rate=25*ones(md.mesh.numberofvertices,1);

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

    r22818 r24864  
    8484        md.mask.land_levelset = 1-md.mask.ocean_levelset;
    8585        md.mask.ice_levelset = ones(md.mesh.numberofvertices,1);
    86         md.mask.groundedice_levelset = -ones(md.mesh.numberofvertices,1);
     86        md.mask.ocean_levelset = -ones(md.mesh.numberofvertices,1);
    8787        pos=find(md.mesh.lat <-80);
    8888        md.mask.ice_levelset(pos(1))=-1; % ice yes! 
    89         md.mask.groundedice_levelset(pos(1))=1; % ice grounded! 
     89        md.mask.ocean_levelset(pos(1))=1; % ice grounded! 
    9090
    9191        di=md.materials.rho_ice/md.materials.rho_water;

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

    r22818 r24864  
    8585        md.slr.love_l = love_numbers('l','CM'); md.slr.love_l(nlove+1:end)=[];
    8686
    87         md.mask.groundedice_levelset = ones(md.mesh.numberofvertices,1);
     87        md.mask.ocean_levelset = ones(md.mesh.numberofvertices,1);
    8888        md.mask.ice_levelset = ones(md.mesh.numberofvertices,1);
    8989        pos=find(md.slr.deltathickness~=0);

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

    r23406 r24864  
    147147        md.slr.love_l = love_numbers('l','CM'); md.slr.love_l(nlove+1:end)=[];
    148148
    149         md.mask.groundedice_levelset = ones(md.mesh.numberofvertices,1);
     149        md.mask.ocean_levelset = ones(md.mesh.numberofvertices,1);
    150150        md.mask.ice_levelset = ones(md.mesh.numberofvertices,1);
    151151        pos=find(md.slr.deltathickness~=0);

  • issm/trunk-jpl/examples/SquareIceShelf/Square.par

    r18823 r24864  
    1212disp('      creating drag');
    1313md.friction.coefficient=200*ones(md.mesh.numberofvertices,1);
    14 md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
     14md.friction.coefficient(find(md.mask.ocean_levelset<0.))=0.;
    1515md.friction.p=ones(md.mesh.numberofelements,1);
    1616md.friction.q=ones(md.mesh.numberofelements,1);
Note: See TracChangeset for help on using the changeset viewer.