Ignore:
Timestamp:
08/01/18 14:40:07 (7 years ago)
Author:
caronlam
Message:

mask fix and missing new variable declaration

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/test/NightlyRun/test2002.py

    r22981 r23039  
    5555md.mask.groundedice_levelset=-np.ones((md.mesh.numberofvertices))
    5656
     57#make sure that the elements that have loads are fully grounded:
     58pos=np.nonzero(md.slr.deltathickness)[0]
     59md.mask.groundedice_levelset[md.mesh.elements[pos,:]-1]=1;
     60
    5761#make sure wherever there is an ice load, that the mask is set to ice:
    58 pos=np.nonzero(md.slr.deltathickness)[0]
    5962icemask[md.mesh.elements[pos,:]-1]=-1
     63md.mask.ice_levelset=icemask
    6064
     65md.slr.ocean_area_scaling=0
    6166
    6267#geometry
     
    7176md.materials.rheology_B=paterson(md.initialization.temperature)
    7277md.materials.rheology_n=3*np.ones((md.mesh.numberofelements))
     78
     79#New stuff
     80md.slr.spcthickness = np.nan((md.mesh.numberofvertices,));
     81md.slr.Ngia = np.zeros((md.mesh.numberofvertices,))
     82md.slr.Ugia = np.zeros((md.mesh.numberofvertices,))
    7383
    7484#Miscellaneous
Note: See TracChangeset for help on using the changeset viewer.