Changeset 1561

Timestamp:

Aug 23, 2011, 9:10:27 AM (13 years ago)

Author:

Ehouarn Millour

Message:

Bug correction on the initialiaation of surface pressure in "academic" and aqua/terra-planet cases.
EM
PS: Don't forget that the aqua-terra planet case is so far only implemented in the seqiuential case.

Location:

LMDZ5/trunk/libf

Files:

• dyn3d/iniacademic.F90 (modified) (1 diff)
• dyn3dpar/iniacademic.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/dyn3d/iniacademic.F90

@@ -209 +209 @@
         ! surface pressure
         if (iflag_phys>2) then
+           ! specific value for CMIP5 aqua/terra planets
+           ! "Specify the initial dry mass to be equivalent to
+           !  a global mean surface pressure (101325 minus 245) Pa."
+           ps(:)=101080.  
+        else
+           ! use reference surface pressure
            ps(:)=preff
-        else
-           ps(:)=101080.
         endif
+        
         ! ground geopotential
         phis(:)=0.

LMDZ5/trunk/libf/dyn3dpar/iniacademic.F90

@@ -115 +115 @@
   endif
 
-  academic_case: if (iflag_phys == 2) then
+  academic_case: if (iflag_phys >= 2) then
      ! initializations
 
@@ -208 +208 @@
      IF (.NOT. read_start) THEN
         ! surface pressure
-        ps(:)=preff
+        if (iflag_phys>2) then
+           ! specific value for CMIP5 aqua/terra planets
+           ! "Specify the initial dry mass to be equivalent to
+           !  a global mean surface pressure (101325 minus 245) Pa."
+           ps(:)=101080.  
+        else
+           ! use reference surface pressure
+           ps(:)=preff
+        endif
+        
         ! ground geopotential
         phis(:)=0.