Changeset 22267 for issm/trunk-jpl/test/Par/SquareThermal.py
- Timestamp:
- 11/18/17 15:19:15 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/Par/SquareThermal.py
r21409 r22267 26 26 print " creating temperatures" 27 27 md.initialization.temperature=(273.-20.)*numpy.ones((md.mesh.numberofvertices)) 28 md.initialization.pressure=numpy.zeros((md.mesh.numberofvertices,)) 29 md.initialization.waterfraction=numpy.zeros((md.mesh.numberofvertices,)) 30 md.initialization.watercolumn=numpy.zeros((md.mesh.numberofvertices,)) 28 31 29 32 print " creating flow law parameter" … … 33 36 print " creating surface mass balance" 34 37 md.smb.mass_balance=numpy.ones((md.mesh.numberofvertices))/md.constants.yts #1m/a 35 md.basalforcings.melting_rate=0.*numpy.ones((md.mesh.numberofvertices))/md.constants.yts #1m/a 38 #md.basalforcings.melting_rate=0.*numpy.ones((md.mesh.numberofvertices))/md.constants.yts #1m/a 39 md.basalforcings.groundedice_melting_rate=0.*numpy.ones((md.mesh.numberofvertices))/md.constants.yts #1m/a 40 md.basalforcings.floatingice_melting_rate=0.*numpy.ones((md.mesh.numberofvertices))/md.constants.yts #1m/a 36 41 37 42 #Deal with boundary conditions: … … 41 46 42 47 print " boundary conditions for thermal model" 43 md.thermal.spctemperature =md.initialization.temperature48 md.thermal.spctemperature[:]=md.initialization.temperature 44 49 md.basalforcings.geothermalflux=numpy.zeros((md.mesh.numberofvertices)) 45 50 md.basalforcings.geothermalflux[numpy.nonzero(md.mask.groundedice_levelset>0.)[0]]=1.*10**-3 #1 mW/m^2
Note:
See TracChangeset
for help on using the changeset viewer.