Changeset 13857 for issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py

Timestamp:

10/30/12 15:03:19 (12 years ago)

Author:

jschierm

Message:

NEW: Added python model.extract and contourenvelope (plus matlab cosmetic changes).

File:

• issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py (modified) (4 diffs)

issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py

@@ -38 +38 @@
         pos=numpy.nonzero(numpy.logical_and(md.mesh.vertexonboundary,numpy.logical_not(vertexonicefront)))[0]
         if not numpy.size(pos):
-                print("SetMarineIceSheetBC warning: ice front all around the glacier, no dirichlet found. Dirichlet must be added manually.")
+                print "SetMarineIceSheetBC warning: ice front all around the glacier, no dirichlet found. Dirichlet must be added manually."
 
         md.diagnostic.spcvx=float('nan')*numpy.ones((md.mesh.numberofvertices,1))
@@ -50 +50 @@
         #Dirichlet Values
         if isinstance(md.inversion.vx_obs,numpy.ndarray) and numpy.size(md.inversion.vx_obs,axis=0)==md.mesh.numberofvertices and isinstance(md.inversion.vy_obs,numpy.ndarray) and numpy.size(md.inversion.vy_obs,axis=0)==md.mesh.numberofvertices:
-                print("      boundary conditions for diagnostic model: spc set as observed velocities")
+                print "      boundary conditions for diagnostic model: spc set as observed velocities"
                 md.diagnostic.spcvx[pos]=md.inversion.vx_obs[pos]
                 md.diagnostic.spcvy[pos]=md.inversion.vy_obs[pos]
         else:
-                print("      boundary conditions for diagnostic model: spc set as zero")
+                print "      boundary conditions for diagnostic model: spc set as zero"
 
         md.hydrology.spcwatercolumn=numpy.zeros((md.mesh.numberofvertices,2))
@@ -83 +83 @@
         if numpy.all(numpy.isnan(md.surfaceforcings.precipitation)) and (md.surfaceforcings.ispdd==1):
                 md.surfaceforcings.precipitation=numpy.zeros((md.mesh.numberofvertices,1))
-                print("      no surfaceforcings.precipitation specified: values set as zero")
+                print "      no surfaceforcings.precipitation specified: values set as zero"
         if numpy.all(numpy.isnan(md.surfaceforcings.mass_balance)) and (md.surfaceforcings.ispdd==0):
                 md.surfaceforcings.mass_balance=numpy.zeros((md.mesh.numberofvertices,1))
-                print("      no surfaceforcings.mass_balance specified: values set as zero")
+                print "      no surfaceforcings.mass_balance specified: values set as zero"
         if numpy.all(numpy.isnan(md.basalforcings.melting_rate)):
                 md.basalforcings.melting_rate=numpy.zeros((md.mesh.numberofvertices,1))
-                print("      no basalforcings.melting_rate specified: values set as zero")
+                print "      no basalforcings.melting_rate specified: values set as zero"
         if numpy.all(numpy.isnan(md.balancethickness.thickening_rate)):
                 md.balancethickness.thickening_rate=numpy.zeros((md.mesh.numberofvertices,1))
-                print("      no balancethickness.thickening_rate specified: values set as zero")
+                print "      no balancethickness.thickening_rate specified: values set as zero"
 
         md.prognostic.spcthickness=float('nan')*numpy.ones((md.mesh.numberofvertices,1))
@@ -104 +104 @@
                 if not isinstance(md.basalforcings.geothermalflux,numpy.ndarray) or not numpy.size(md.basalforcings.geothermalflux,axis=0)==md.mesh.numberofvertices:
                         md.basalforcings.geothermalflux=numpy.zeros((md.mesh.numberofvertices,1))
-                        md.basalforcings.geothermalflux[numpy.nonzero(md.mask.vertexongroundedice)]=50.*10.**-3; #50mW/m2
+                        md.basalforcings.geothermalflux[numpy.nonzero(md.mask.vertexongroundedice)]=50.*10.**-3    #50mW/m2
         else:
-                print("      no thermal boundary conditions created: no observed temperature found");
+                print "      no thermal boundary conditions created: no observed temperature found"
 
         return md