Changeset 121

Timestamp:

02/01/13 00:13:46 (11 years ago)

Author:

dubos

Message:

Options for calculation of Exner pressure

Location:

codes/icosagcm/trunk/src

Files:

• exner.f90 (modified) (2 diffs)
• geopotential_mod.f90 (modified) (1 diff)

codes/icosagcm/trunk/src/exner.f90

@@ -1 +1 @@
 MODULE exner_mod
 
+  INTEGER :: caldyn_exner
+
 CONTAINS
-
+  
   SUBROUTINE exner(f_ps,f_p,f_pks,f_pk)
   USE icosa
@@ -43 +45 @@
     REAL(rstd) :: delta 
     
-    IF(.FALSE.) THEN          ! LMD-Z style calculation of Exner pressure
+    IF(caldyn_exner == 1) THEN          ! LMD-Z style calculation of Exner pressure
        !! Compute Alpha and Beta
        

codes/icosagcm/trunk/src/geopotential_mod.f90

@@ -44 +44 @@
 
 !!! Compute geopotential
-
+  ! LMD-Z : dh=pi.dtheta-vdp  => (1/rho) dp = pi.dtheta - d(h=pi.theta) = -dtheta.gradpi
+  ! g.dz =  -(1/rho)dp =  theta.gradpi
+  ! direct (variational)
+  ! dz = -(1/rho.g)dp, rho = RT/p = (R/C_p).theta.pi/p
+   
   ! for first layer
    DO j=jj_begin-offset,jj_end+offset