- Timestamp:
- 02/12/19 06:10:51 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py
r21303 r23716 70 70 if np.ndim(md.inversion.vy_obs)==1: 71 71 md.inversion.vy_obs=md.inversion.vy_obs.reshape(-1,) 72 print " boundary conditions for stressbalance model: spc set as observed velocities"72 print(" boundary conditions for stressbalance model: spc set as observed velocities") 73 73 md.stressbalance.spcvx[pos]=md.inversion.vx_obs[pos] 74 74 md.stressbalance.spcvy[pos]=md.inversion.vy_obs[pos] 75 75 else: 76 print " boundary conditions for stressbalance model: spc set as zero"76 print(" boundary conditions for stressbalance model: spc set as zero") 77 77 78 78 #Create zeros basalforcings and smb … … 83 83 if np.all(np.isnan(md.balancethickness.thickening_rate)): 84 84 md.balancethickness.thickening_rate=np.zeros((md.mesh.numberofvertices)) 85 print " no balancethickness.thickening_rate specified: values set as zero"85 print(" no balancethickness.thickening_rate specified: values set as zero") 86 86 md.masstransport.spcthickness=float('nan')*np.ones((md.mesh.numberofvertices)) 87 87 md.balancethickness.spcthickness=float('nan')*np.ones((md.mesh.numberofvertices)) … … 96 96 md.basalforcings.geothermalflux=np.zeros((md.mesh.numberofvertices)) 97 97 else: 98 print " no thermal boundary conditions created: no observed temperature found"98 print(" no thermal boundary conditions created: no observed temperature found") 99 99 100 100 return md
Note:
See TracChangeset
for help on using the changeset viewer.