Changeset 1625 for LMDZ5/trunk/libf/dyn3dpar/etat0_netcdf.F90

Timestamp:

May 9, 2012, 3:14:48 PM (13 years ago)

Author:

lguez

Message:

Created logical variable "pressure_exner". "pressure_exner" replaces
"disvert_type" to choose between calls to "exner_hyb" and
"exner_milieu". If "pressure_exner" is true, pressure inside layers is
computed from Exner function ("exner_hyb"), else it is the mean of
pressure values at interfaces ("exner_milieu"). "disvert_type" now
only chooses between "disvert" and "disvert_noterre". Default value
for "pressure_exner" is true if "disvert_type" equals 1.

File:

• LMDZ5/trunk/libf/dyn3dpar/etat0_netcdf.F90 (modified) (1 diff)

LMDZ5/trunk/libf/dyn3dpar/etat0_netcdf.F90

@@ -251 +251 @@
 !*******************************************************************************
   CALL pression(ip1jmp1, ap, bp, psol, p3d)
-  if (disvert_type.eq.1) then
+  if (pressure_exner) then
     CALL exner_hyb(ip1jmp1, psol, p3d, alpha, beta, pks, pk, y)
-  else ! we assume that we are in the disvert_type==2 case
+  else
     CALL exner_milieu(ip1jmp1,psol,p3d,beta,pks,pk,y)
   endif