Changeset 21044
- Timestamp:
- 08/01/16 16:20:53 (9 years ago)
- Location:
- issm/trunk-jpl/test/NightlyRun
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test2002.m
r21001 r21044 15 15 pos=find(late <-80); 16 16 md.slr.deltathickness(pos)=-100; 17 %greenl nd17 %greenland 18 18 pos=find(late > 70 & late < 80 & longe>-60 & longe<-30); 19 19 md.slr.deltathickness(pos)=-100; -
issm/trunk-jpl/test/NightlyRun/test2002.py
r21001 r21044 20 20 #parameterize slr solution: 21 21 #slr loading: {{{ 22 md.slr.deltathickness=zeros((md.mesh.numberof vertices,1))22 md.slr.deltathickness=zeros((md.mesh.numberofelements,1)) 23 23 md.slr.sealevel=zeros((md.mesh.numberofvertices,1)) 24 24 #antarctica 25 pos=numpy.nonzero(md.mesh.lat <-80) 25 late=numpy.sum(md.mesh.lat[md.mesh.elements-1],axis=1)/3 26 longe=numpy.sum(md.mesh.long[md.mesh.elements-1],axis=1)/3 27 pos=numpy.nonzero(late <-80) 26 28 md.slr.deltathickness[pos]=-100 27 #greenl nd28 pos=numpy.nonzero(logical_and_n( md.mesh.lat > 70,md.mesh.lat < 80,md.mesh.long>-60,md.mesh.long<-30))29 #greenland 30 pos=numpy.nonzero(logical_and_n(late > 70,late < 80,longe>-60,longe<-30)) 29 31 md.slr.deltathickness[pos]=-100 30 32
Note:
See TracChangeset
for help on using the changeset viewer.