![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
Changeset 22864 for issm/trunk-jpl/test/NightlyRun/test236.py
- Timestamp:
- 06/22/18 01:02:47 (7 years ago)
- File:
-
- 1 edited
-
issm/trunk-jpl/test/NightlyRun/test236.py (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test236.py
r21408 r22864 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) .reshape(-1,1)#-10*ones(md.mesh.numberofvertices,1)43 md.thermal.spctemperature=np.tile(md.thermal.spctemperature,( 1,md.timestepping.final_time/md.timestepping.time_step))42 md.thermal.spctemperature=np.mean(md.smb.temperatures_lgm[0:md.mesh.numberofvertices,:],axis=1) #-10*ones(md.mesh.numberofvertices,1) 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 md.thermal.spctemperature=np.vstack((md.thermal.spctemperature,itemp .reshape(1,-1)))45 md.thermal.spctemperature=np.vstack((md.thermal.spctemperature,itemp)) 46 46 47 47 md.initialization.temperature=md.smb.temperatures_lgm[0:md.mesh.numberofvertices,0] #*ones(md.mesh.numberofvertices,1) … … 52 52 md.smb.precipitations_lgm=np.zeros((md.mesh.numberofvertices+1,12)) 53 53 for imonth in xrange(0,12): 54 55 56 57 54 md.smb.precipitations_presentday[0:md.mesh.numberofvertices,imonth]=-0.4*10**(-6)*md.mesh.y+0.5 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.5 57 md.smb.precipitations_lgm[md.mesh.numberofvertices,imonth]=((float(imonth)+1.)/12.) 58 58 59 59 # Interpolation factors 60 md.smb.Tdiff=0.5*np.ones((2, md.timestepping.final_time))61 md.smb.sealev=0.5*np.ones((2, md.timestepping.final_time))60 md.smb.Tdiff=0.5*np.ones((2,int(md.timestepping.final_time))) 61 md.smb.sealev=0.5*np.ones((2,int(md.timestepping.final_time))) 62 62 # Year of each data point 63 md.smb.Tdiff[1,1: md.timestepping.final_time]=np.arange(1.,md.timestepping.final_time)64 md.smb.sealev[1,1: md.timestepping.final_time]=np.arange(1.,md.timestepping.final_time)63 md.smb.Tdiff[1,1:int(md.timestepping.final_time)]=np.arange(1.,int(md.timestepping.final_time)) 64 md.smb.sealev[1,1:int(md.timestepping.final_time)]=np.arange(1.,int(md.timestepping.final_time)) 65 65 66 66 # time steps and resolution … … 69 69 70 70 71 # 71 # 72 72 md.transient.requested_outputs=['default','SmbMonthlytemperatures'] 73 73 md=setflowequation(md,'SSA','all') … … 76 76 77 77 #Fields and tolerances to track changes 78 field_names =['Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SmbMonthlytemperatures1','SmbMassBalance1',\ 79 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SmbMonthlytemperatures2','SmbMassBalance2',\ 80 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SmbMonthlytemperatures3','SmbMassBalance3'] 81 field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\ 82 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\ 83 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13] 84 field_values=[\ 85 md.results.TransientSolution[0].Vx,\ 86 md.results.TransientSolution[0].Vy,\ 87 md.results.TransientSolution[0].Vel,\ 88 md.results.TransientSolution[0].Pressure,\ 89 md.results.TransientSolution[0].Base,\ 90 md.results.TransientSolution[0].Surface,\ 91 md.results.TransientSolution[0].Thickness,\ 92 md.results.TransientSolution[0].SmbMonthlytemperatures,\ 93 md.results.TransientSolution[0].SmbMassBalance,\ 94 md.results.TransientSolution[1].Vx,\ 95 md.results.TransientSolution[1].Vy,\ 96 md.results.TransientSolution[1].Vel,\ 97 md.results.TransientSolution[1].Pressure,\ 98 md.results.TransientSolution[1].Base,\ 99 md.results.TransientSolution[1].Surface,\ 100 md.results.TransientSolution[1].Thickness,\ 101 md.results.TransientSolution[1].SmbMonthlytemperatures,\ 102 md.results.TransientSolution[1].SmbMassBalance,\ 103 md.results.TransientSolution[2].Vx,\ 104 md.results.TransientSolution[2].Vy,\ 105 md.results.TransientSolution[2].Vel,\ 106 md.results.TransientSolution[2].Pressure,\ 107 md.results.TransientSolution[2].Base,\ 108 md.results.TransientSolution[2].Surface,\ 109 md.results.TransientSolution[2].Thickness,\ 110 md.results.TransientSolution[2].SmbMonthlytemperatures,\ 111 md.results.TransientSolution[2].SmbMassBalance,\ 112 ] 78 field_names =['Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SmbMonthlytemperatures1','SmbMassBalance1', 79 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SmbMonthlytemperatures2','SmbMassBalance2', 80 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SmbMonthlytemperatures3','SmbMassBalance3'] 81 field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13, 82 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13, 83 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13] 84 field_values=[md.results.TransientSolution[0].Vx, 85 md.results.TransientSolution[0].Vy, 86 md.results.TransientSolution[0].Vel, 87 md.results.TransientSolution[0].Pressure, 88 md.results.TransientSolution[0].Base, 89 md.results.TransientSolution[0].Surface, 90 md.results.TransientSolution[0].Thickness, 91 md.results.TransientSolution[0].SmbMonthlytemperatures, 92 md.results.TransientSolution[0].SmbMassBalance, 93 md.results.TransientSolution[1].Vx, 94 md.results.TransientSolution[1].Vy, 95 md.results.TransientSolution[1].Vel, 96 md.results.TransientSolution[1].Pressure, 97 md.results.TransientSolution[1].Base, 98 md.results.TransientSolution[1].Surface, 99 md.results.TransientSolution[1].Thickness, 100 md.results.TransientSolution[1].SmbMonthlytemperatures, 101 md.results.TransientSolution[1].SmbMassBalance, 102 md.results.TransientSolution[2].Vx, 103 md.results.TransientSolution[2].Vy, 104 md.results.TransientSolution[2].Vel, 105 md.results.TransientSolution[2].Pressure, 106 md.results.TransientSolution[2].Base, 107 md.results.TransientSolution[2].Surface, 108 md.results.TransientSolution[2].Thickness, 109 md.results.TransientSolution[2].SmbMonthlytemperatures, 110 md.results.TransientSolution[2].SmbMassBalance]
Note:
See TracChangeset
for help on using the changeset viewer.
