Changeset 24214 for issm/trunk-jpl/test/NightlyRun/test236.py
- Timestamp:
- 10/11/19 00:27:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test236.py
r23793 r24214 35 35 md.smb.temperatures_presentday[0:md.mesh.numberofvertices, imonth] = tmonth[imonth] 36 36 md.smb.temperatures_lgm[0:md.mesh.numberofvertices, imonth] = tmonth[imonth] - 20. 37 37 # Time for the last line: 38 38 md.smb.temperatures_presentday[md.mesh.numberofvertices, imonth] = ((float(imonth) + 1.) / 12.) 39 39 md.smb.temperatures_lgm[md.mesh.numberofvertices, imonth] = ((float(imonth) + 1.) / 12.) 40 40 41 41 # creating initialization and spc temperatures initialization and spc 42 md.thermal.spctemperature = np.mean(md.smb.temperatures_lgm[0:md.mesh.numberofvertices, :], axis=1) #-10*ones(md.mesh.numberofvertices, 1)42 md.thermal.spctemperature = np.mean(md.smb.temperatures_lgm[0:md.mesh.numberofvertices, :], axis=1) # - 10 * ones(md.mesh.numberofvertices, 1) 43 43 md.thermal.spctemperature = np.tile(md.thermal.spctemperature, (int(md.timestepping.final_time / md.timestepping.time_step), 1)).T 44 44 itemp = np.arange(0, md.timestepping.final_time, md.timestepping.time_step) 45 45 md.thermal.spctemperature = np.vstack((md.thermal.spctemperature, itemp)) 46 46 47 md.initialization.temperature = md.smb.temperatures_lgm[0:md.mesh.numberofvertices, 0] #*ones(md.mesh.numberofvertices, 1)47 md.initialization.temperature = md.smb.temperatures_lgm[0:md.mesh.numberofvertices, 0] # * ones(md.mesh.numberofvertices, 1) 48 48 md.smb.initialize(md) 49 49 … … 52 52 md.smb.precipitations_lgm = np.zeros((md.mesh.numberofvertices + 1, 12)) 53 53 for imonth in range(0, 12): 54 md.smb.precipitations_presentday[0:md.mesh.numberofvertices, imonth] = - 0.4 * 10**(-6) * md.mesh.y + 0.554 md.smb.precipitations_presentday[0:md.mesh.numberofvertices, imonth] = - 0.4 * 10**(- 6) * md.mesh.y + 0.5 55 55 md.smb.precipitations_presentday[md.mesh.numberofvertices, imonth] = ((float(imonth) + 1.) / 12.) 56 md.smb.precipitations_lgm[0:md.mesh.numberofvertices, imonth] = - 0.4 * 10**(-6) * md.mesh.y + 0.556 md.smb.precipitations_lgm[0:md.mesh.numberofvertices, imonth] = - 0.4 * 10**(- 6) * md.mesh.y + 0.5 57 57 md.smb.precipitations_lgm[md.mesh.numberofvertices, imonth] = ((float(imonth) + 1.) / 12.) 58 58
Note:
See TracChangeset
for help on using the changeset viewer.