Changeset 669 for codes/icosagcm/trunk/src/vertical

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

File:

• codes/icosagcm/trunk/src/vertical/disvert_std.f90 (modified) (2 diffs)

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(:)