Changeset 366 for codes/icosagcm/trunk/src/caldyn_kernels_base.f90

Timestamp:

10/30/15 15:41:06 (9 years ago)

Author:

dubos

Message:

Progress towards NH

File:

• codes/icosagcm/trunk/src/caldyn_kernels_base.f90 (modified) (3 diffs)

codes/icosagcm/trunk/src/caldyn_kernels_base.f90

@@ -18 +18 @@
   !$OMP THREADPRIVATE(caldyn_conserv) 
 
-  TYPE(t_message),PUBLIC :: req_ps, req_mass, req_theta_rhodz, req_u, req_qu
+  TYPE(t_message),PUBLIC :: req_ps, req_mass, req_theta_rhodz, req_u, req_qu, req_geopot, req_w
 
   PUBLIC :: compute_geopot, compute_caldyn_vert
@@ -51 +51 @@
 
     IF(caldyn_eta==eta_mass .AND. .NOT. DEC) THEN
-
 !!! Compute exner function and geopotential
        DO l = 1,llm
@@ -64 +63 @@
        ENDDO
     ELSE 
-       ! We are using a Lagrangian vertical coordinate
-       ! Pressure must be computed first top-down (temporarily stored in pk)
+       ! We are using DEC or a Lagrangian vertical coordinate
+       ! Pressure is computed first top-down (temporarily stored in pk)
        ! Then Exner pressure and geopotential are computed bottom-up
-       ! Notice that the computation below should work also when caldyn_eta=eta_mass          
+       ! Works also when caldyn_eta=eta_mass          
 
        IF(boussinesq) THEN ! compute geopotential and pk=Lagrange multiplier