Changeset 806

Timestamp:

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

Author:

jisesh

Message:

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

Location:

codes/icosagcm/devel/Python/test

Files:

• py/Baroclinic_3D_ullrich.py (modified) (3 diffs)
• xml/field_def_simple.xml (modified) (1 diff)
• xml/filedef_simple.xml (modified) (1 diff)

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)

codes/icosagcm/devel/Python/test/xml/field_def_simple.xml

@@ -22 +22 @@
       <field id="Ai" />
       <field id="phis" />
-      <field id="phi" />
 
+      <field_group axis_ref="levp1">
+        <field id="Phi" />
+      </field_group>
+ 
       <field_group axis_ref="lev"> 
         <field id="temp" />

codes/icosagcm/devel/Python/test/xml/filedef_simple.xml

@@ -19 +19 @@
       <field field_ref="uz" />
       <field field_ref="z" />
+      <field field_ref="Phi" />
       <field field_ref="theta" />
       <field field_ref="curl" />