Changeset 23634
- Timestamp:
- 01/16/19 08:19:51 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/materials/cuffey.py
r23633 r23634 27 27 rigidity=np.zeros_like(T) 28 28 pos=np.nonzero(T<=-45) 29 print pos30 29 if len(pos): 31 30 rigidity[pos]=10**8*(-0.000396645116301*(T[pos]+50)**3+ 0.013345579471334*(T[pos]+50)**2 -0.356868703259105*(T[pos]+50)+7.272363035371383)
Note:
See TracChangeset
for help on using the changeset viewer.