Changeset 669 for codes/icosagcm/trunk/src

Timestamp:

01/30/18 15:28:08 (6 years ago)

Author:

dubos

Message:

trunk/disvert : set pa=preff/2 rather than hard-coded pa=50000

Location:

codes/icosagcm/trunk/src

Files:

• base/earth_const.f90 (modified) (1 diff)
• vertical/disvert_std.f90 (modified) (2 diffs)

codes/icosagcm/trunk/src/base/earth_const.f90

@@ -13 +13 @@
   REAL(rstd),SAVE :: Treff=273.
   REAL(rstd),SAVE :: preff=101325.
-  REAL(rstd),SAVE :: pa=50000.
+  REAL(rstd),SAVE :: pa=50000. ! default value set to preff/2 by disvert_std
   REAL(rstd),SAVE :: scale_height=8000. ! atmospheric scale height (m)
   REAL(rstd),SAVE :: scale_factor=1.

codes/icosagcm/trunk/src/vertical/disvert_std.f90

@@ -16 +16 @@
 
   SUBROUTINE init_disvert
-  USE icosa
-  IMPLICIT NONE
-  
+    pa = .5*preff ! differs from 50000 if preff differs from 1e5
+    CALL getin('pa',pa)
+
     ALLOCATE(ap(llm+1))
     ALLOCATE(bp(llm+1))
@@ -29 +29 @@
 
   SUBROUTINE disvert(ap,bp,presnivs)
-  USE icosa
   USE mpipara
-  USE earth_const
-  IMPLICIT NONE
   REAL(rstd),INTENT(OUT) :: ap(:)
   REAL(rstd),INTENT(OUT) :: bp(:)