Changeset 96

Timestamp:

08/03/12 23:47:57 (12 years ago)

Author:

dubos

Message:

Computation and output of omega=Dp/Dt?

Location:

codes/icosagcm/trunk/src

Files:

• caldyn_gcm.f90 (modified) (2 diffs)
• omega.f90 (added)

codes/icosagcm/trunk/src/caldyn_gcm.f90

@@ -641 +641 @@
     USE theta2theta_rhodz_mod
     USE pression_mod
+    USE omega_mod
     USE write_field
     TYPE(t_field),POINTER :: f_ps(:), f_phis(:), f_u(:), f_theta_rhodz(:), f_dps(:), &
@@ -650 +651 @@
     CALL vorticity(f_u,f_buf_v)
     CALL writefield("vort",f_buf_v)
+
+    CALL w_omega(f_ps, f_u, f_buf_i)
+    CALL writefield('omega', f_buf_i)
     
     ! Temperature