Changeset 24214 for issm/trunk-jpl/test/Par/RoundSheetEISMINT.py

Timestamp:

10/11/19 00:27:00 (5 years ago)

Author:

bdef

Message:

CHG: syntax cahnge to meet most of Pep8 requirement

File:

• issm/trunk-jpl/test/Par/RoundSheetEISMINT.py (modified) (1 diff)

issm/trunk-jpl/test/Par/RoundSheetEISMINT.py

@@ -4 +4 @@
 #Ok, start defining model parameters here
 print("      creating thickness")
-md.geometry.thickness=10.*numpy.ones((md.mesh.numberofvertices))
-md.geometry.base=numpy.zeros((md.mesh.numberofvertices))
-md.geometry.surface=md.geometry.base+md.geometry.thickness
+md.geometry.thickness = 10. * numpy.ones((md.mesh.numberofvertices))
+md.geometry.base = numpy.zeros((md.mesh.numberofvertices))
+md.geometry.surface = md.geometry.base + md.geometry.thickness
 
 print("      creating drag")
-md.friction.coefficient=20.*numpy.ones((md.mesh.numberofvertices)) 
-md.friction.p=numpy.ones((md.mesh.numberofelements))
-md.friction.q=numpy.ones((md.mesh.numberofelements))
+md.friction.coefficient = 20. * numpy.ones((md.mesh.numberofvertices))
+md.friction.p = numpy.ones((md.mesh.numberofelements))
+md.friction.q = numpy.ones((md.mesh.numberofelements))
 
 print("      creating temperatures")
-tmin=238.15    #K
-st=1.67*10**-2/1000.    #k/m
-radius=numpy.sqrt((md.mesh.x)**2+(md.mesh.y)**2)
-md.initialization.temperature=tmin+st*radius
-md.basalforcings.geothermalflux=4.2*10**-2*numpy.ones((md.mesh.numberofvertices))
+tmin = 238.15  #K
+st = 1.67 * 10**- 2 / 1000.  #k / m
+radius = numpy.sqrt((md.mesh.x)**2 + (md.mesh.y)**2)
+md.initialization.temperature = tmin + st * radius
+md.basalforcings.geothermalflux = 4.2 * 10**- 2 * numpy.ones((md.mesh.numberofvertices))
 
 print("      creating flow law parameter")
-md.materials.rheology_B=6.81*10**7*numpy.ones((md.mesh.numberofvertices))    #to have the same B as the analytical solution 
-md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements))
+md.materials.rheology_B = 6.81 * 10**7 * numpy.ones((md.mesh.numberofvertices))  #to have the same B as the analytical solution
+md.materials.rheology_n = 3. * numpy.ones((md.mesh.numberofelements))
 
 print("      creating surface mass balance")
-smb_max=0.5    #m/yr
-sb=10**-2/1000.    #m/yr/m
-rel=450.*1000.    #m
-md.smb.mass_balance=numpy.minimum(smb_max*numpy.ones_like(radius),sb*(rel-radius))
+smb_max = 0.5  #m / yr
+sb = 10**- 2 / 1000.  #m / yr / m
+rel = 450. * 1000.  #m
+md.smb.mass_balance = numpy.minimum(smb_max * numpy.ones_like(radius), sb * (rel - radius))
 
 print("      creating velocities")
-constant=0.3
-md.inversion.vx_obs=constant/2.*md.mesh.x*(md.geometry.thickness)**-1
-md.inversion.vy_obs=constant/2.*md.mesh.y*(md.geometry.thickness)**-1
-md.inversion.vel_obs=numpy.sqrt((md.inversion.vx_obs)**2+(md.inversion.vy_obs)**2)
-md.initialization.vx=numpy.zeros((md.mesh.numberofvertices))
-md.initialization.vy=numpy.zeros((md.mesh.numberofvertices))
-md.initialization.vz=numpy.zeros((md.mesh.numberofvertices))
-md.initialization.pressure=numpy.zeros((md.mesh.numberofvertices))
+constant = 0.3
+md.inversion.vx_obs = constant / 2. * md.mesh.x * (md.geometry.thickness)**- 1
+md.inversion.vy_obs = constant / 2. * md.mesh.y * (md.geometry.thickness)**- 1
+md.inversion.vel_obs = numpy.sqrt((md.inversion.vx_obs)**2 + (md.inversion.vy_obs)**2)
+md.initialization.vx = numpy.zeros((md.mesh.numberofvertices))
+md.initialization.vy = numpy.zeros((md.mesh.numberofvertices))
+md.initialization.vz = numpy.zeros((md.mesh.numberofvertices))
+md.initialization.pressure = numpy.zeros((md.mesh.numberofvertices))
 
 #Deal with boundary conditions:
 print("      boundary conditions for stressbalance model:")
-md=SetMarineIceSheetBC(md,'../Exp/RoundFrontEISMINT.exp')
+md = SetMarineIceSheetBC(md, '../Exp/RoundFrontEISMINT.exp')
 
-radius=numpy.sqrt((md.mesh.x)**2+(md.mesh.y)**2)
-pos=numpy.nonzero(radius==numpy.min(radius))[0]
-md.mesh.x[pos]=0.
-md.mesh.y[pos]=0.    #the closest node to the center is changed to be exactly at the center
+radius = numpy.sqrt((md.mesh.x)**2 + (md.mesh.y)**2)
+pos = numpy.nonzero(radius == numpy.min(radius))[0]
+md.mesh.x[pos] = 0.
+md.mesh.y[pos] = 0.  #the closest node to the center is changed to be exactly at the center
 
-md.stressbalance.spcvx[pos]=0.
-md.stressbalance.spcvy[pos]=0.
-md.stressbalance.spcvz[pos]=0.
+md.stressbalance.spcvx[pos] = 0.
+md.stressbalance.spcvy[pos] = 0.
+md.stressbalance.spcvz[pos] = 0.
 
 #parallel options
-md.timestepping.final_time=50000.
+md.timestepping.final_time = 50000.
 
 #Constants
-md.materials.rho_ice=910.
-md.materials.thermalconductivity=2.1
-md.materials.latentheat=3.35*10**5
-md.materials.beta=8.66*10**-4/(md.materials.rho_ice*md.constants.g)    #conversion from K/m to K/Pa
-md.constants.yts=31556926.
+md.materials.rho_ice = 910.
+md.materials.thermalconductivity = 2.1
+md.materials.latentheat = 3.35 * 10**5
+md.materials.beta = 8.66 * 10**- 4 / (md.materials.rho_ice * md.constants.g)  #conversion from K / m to K / Pa
+md.constants.yts = 31556926.