Changeset 24214 for issm/trunk-jpl/test/NightlyRun/test245.py
- Timestamp:
- 10/11/19 00:27:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test245.py
r23793 r24214 19 19 # initalize pdd fields 20 20 md.smb.initialize(md) 21 md.smb.s0p = md.geometry.surface.reshape(- 1, 1)22 md.smb.s0t = md.geometry.surface.reshape(- 1, 1)21 md.smb.s0p = md.geometry.surface.reshape(- 1, 1) 22 md.smb.s0t = md.geometry.surface.reshape(- 1, 1) 23 23 24 24 25 25 md.smb.monthlytemperatures = np.empty((md.mesh.numberofvertices + 1, 12)) 26 26 md.smb.precipitation = np.empty((md.mesh.numberofvertices + 1, 12)) 27 temp_ma_present = - 10. * np.ones((md.mesh.numberofvertices,)) - md.smb.rlaps * md.geometry.surface / 1000.28 temp_mj_present = 10. * np.ones((md.mesh.numberofvertices, )) - md.smb.rlaps * md.geometry.surface / 1000.29 precipitation = 5. * np.ones((md.mesh.numberofvertices, ))27 temp_ma_present = - 10. * np.ones((md.mesh.numberofvertices, )) - md.smb.rlaps * md.geometry.surface / 1000. 28 temp_mj_present = 10. * np.ones((md.mesh.numberofvertices, )) - md.smb.rlaps * md.geometry.surface / 1000. 29 precipitation = 5. * np.ones((md.mesh.numberofvertices, )) 30 30 for imonth in range(12): 31 31 md.smb.monthlytemperatures[0:md.mesh.numberofvertices, imonth] = md.materials.meltingpoint + temp_ma_present + (temp_mj_present - temp_ma_present) * np.sin((imonth + 1. - 4.) * np.pi / 6.0)
Note:
See TracChangeset
for help on using the changeset viewer.