Changeset 32 for codes/icosagcm/trunk/src/earth_const.f90

Timestamp:

07/27/12 16:04:22 (12 years ago)

Author:

ymipsl

Message:

externalize earth constant in run.def and add scalling factor functionnality

YM

File:

• codes/icosagcm/trunk/src/earth_const.f90 (modified) (1 diff)

codes/icosagcm/trunk/src/earth_const.f90

@@ -3 +3 @@
   USE math_const
   
-  REAL(rstd),PARAMETER :: radius=6.37122E6
-!  REAL(rstd),PARAMETER :: radius=1.
-  REAL(rstd),PARAMETER :: g=9.80616
+  REAL(rstd),SAVE :: radius=6.37122E6
+  REAL(rstd),SAVE :: g=9.80616
   REAL(rstd),PARAMETER :: daysec=86400
-!  REAL(rstd),PARAMETER :: omega=2*Pi/daysec
-  REAL(rstd),PARAMETER :: omega=7.292E-5
-!  REAL(rstd),PARAMETER :: omega=0
-  REAL(rstd),PARAMETER :: kappa=0.2857143
-  REAL(rstd),PARAMETER :: cpp=1004.70885
-  REAL(rstd),PARAMETER :: preff=101325.
-  REAL(rstd),PARAMETER :: pa=50000.  
+  REAL(rstd),SAVE :: omega=7.292E-5
+  REAL(rstd),SAVE :: kappa=0.2857143
+  REAL(rstd),SAVE :: cpp=1004.70885
+  REAL(rstd),SAVE :: preff=101325.
+  REAL(rstd),SAVE :: pa=50000.
+  REAL(rstd),SAVE :: scale_factor=1.
 
+CONTAINS
+  
+  SUBROUTINE init_earth_const
+  USE ioipsl
+  IMPLICIT NONE
+  REAL(rstd) :: X=1
+  
+    CALL getin("radius",radius)
+    CALL getin("g",g)
+    CALL getin("X",scale_factor)
+    CALL getin("omega",omega)  
+    CALL getin("kappa",kappa)  
+    CALL getin("cpp",cpp)  
+    CALL getin("preff",preff)  
+    
+    radius=radius/scale_factor
+    omega=omega/scale_factor
+    PRINT *,"radius = ",radius
+    
+  END SUBROUTINE init_earth_const
+  
+  
 END MODULE earth_const