Changeset 806 for codes/icosagcm/devel/Python/test/py

Timestamp:

02/08/19 18:51:47 (5 years ago)

Author:

jisesh

Message:

devel/Python/Baroclinic_3D_ullrich : output Phi at half-levels

File:

• codes/icosagcm/devel/Python/test/py/Baroclinic_3D_ullrich.py (modified) (3 diffs)

codes/icosagcm/devel/Python/test/py/Baroclinic_3D_ullrich.py

@@ -38 +38 @@
 from dynamico import xios
 from dynamico import precision as prec
-from dynamico.meshes import Cartesian_mesh as Mesh
 from dynamico.kernels import grad, curl, div, KE
 from dynamico.LAM import Davies
@@ -158 +157 @@
     for l in range(llm):
         v[:,l]=(Phi[:,l+1]-Phi[:,l])/(g*m[:,l])
-        w[:,l]=.5*params.g*(W[:,l+1]+W[:,l])/m[:,l]
-        z[:,l]=.5*(Phi[:,l+1]+Phi[:,l])/params.g
+        w[:,l]=.5*g*(W[:,l+1]+W[:,l])/m[:,l]
+        z[:,l]=.5*(Phi[:,l+1]+Phi[:,l])/g
         un[:,l]=u[:,l]/mesh.de
         abs_vort_vk[:,l]=curl_vk[:,l] + mesh.fv
@@ -282 +281 @@
 
             context.send_field_primal('uz', w)
+            context.send_field_primal('Phi', Phi)
             context.send_field_primal('z', z)
             context.send_field_primal('div_fast', div_fast)