Changeset 780 for codes/icosagcm
- Timestamp:
- 11/20/18 12:44:31 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/Python/test/py/Baroclinic_3D_ullrich.py
r774 r780 64 64 def ulon(x,y,eta): 65 65 u = -u0*(sin(pi*y/Ly)**2)*log(eta)*(eta**(-log(eta)/(b*b))) 66 #u = u + up*exp(-(((x-xc)**2+(y-yc)**2)/lp**2))66 u = u + up*exp(-(((x-xc)**2+(y-yc)**2)/lp**2)) 67 67 return u 68 68 … … 134 134 135 135 def diagnose(Phi,S,m,W): 136 s=S/m ; s=.5*(s+abs(s))136 s=S/m 137 137 for l in range(llm): 138 138 v[:,l]=(Phi[:,l+1]-Phi[:,l])/(g*m[:,l]) … … 201 201 202 202 T = args.T 203 dt = 360.203 dt = 100.#180.#360. 204 204 dz = flow0[3].max()/(params.g*llm) 205 205 nt = int(math.ceil(T/dt)) … … 265 265 context.send_field_primal('theta', gas.s) 266 266 context.send_field_primal('uz', w) 267 context.send_field_primal('z', z) 267 268 268 269 print( 'ptop, model top (m) :', unst.getvar('ptop'), Phi.max()/unst.getvar('g'))
Note: See TracChangeset
for help on using the changeset viewer.