- Timestamp:
- 10/18/19 16:10:28 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/boundaryconditions/SetIceSheetBC.py
r24213 r24260 31 31 print(" boundary conditions for stressbalance model: spc set as zero") 32 32 33 #No ice front - 33 #No ice front -> do nothing 34 34 35 35 #Create zeros basalforcings and smb … … 51 51 md.thermal.spctemperature[pos] = md.initialization.temperature[pos] #impose observed temperature on surface 52 52 if not isinstance(md.basalforcings.geothermalflux, np.ndarray) or not np.size(md.basalforcings.geothermalflux) == md.mesh.numberofvertices: 53 md.basalforcings.geothermalflux = 50. * 10**- 53 md.basalforcings.geothermalflux = 50. * 10**-3 * np.ones((md.mesh.numberofvertices)) #50 mW / m^2 54 54 else: 55 55 print(" no thermal boundary conditions created: no observed temperature found")
Note:
See TracChangeset
for help on using the changeset viewer.