Changeset 34

Timestamp:

07/27/12 16:09:21 (12 years ago)

Author:

ymipsl

Message:

Some correction due to externalisation of earth const parameters

YM

File:

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

codes/icosagcm/trunk/src/etat0_ncar.f90

@@ -3 +3 @@
   PRIVATE
 
-  REAL(rstd), PARAMETER :: h0=1.
-  REAL(rstd), PARAMETER :: lon0=3*pi/2
-  REAL(rstd), PARAMETER :: lat0=0.0
-  REAL(rstd), PARAMETER :: alpha=0.0
-  REAL(rstd), PARAMETER :: R0 = radius*0.5
-  REAL(rstd), PARAMETER :: lat1=0.
-  REAL(rstd), PARAMETER :: lat2=0.
-  REAL(rstd), PARAMETER :: lon1=pi/6
-  REAL(rstd), PARAMETER :: lon2=-pi/6
-  REAL(rstd), PARAMETER :: latc1=0.
-  REAL(rstd), PARAMETER :: latc2=0.
-  REAL(rstd), PARAMETER :: lonc1=5*pi/6
-  REAL(rstd), PARAMETER :: lonc2=7*pi/6
-  REAL(rstd), PARAMETER :: zt=1000.0
-  REAL(rstd), PARAMETER :: rt=radius*0.5
-  REAL(rstd), PARAMETER :: zc=5000.0
+  REAL(rstd), SAVE  :: h0=1.
+  REAL(rstd), SAVE  :: lon0=3*pi/2
+  REAL(rstd), SAVE  :: lat0=0.0
+  REAL(rstd), SAVE  :: alpha=0.0
+  REAL(rstd), SAVE  :: R0 
+  REAL(rstd), SAVE  :: lat1=0.
+  REAL(rstd), SAVE  :: lat2=0.
+  REAL(rstd), SAVE  :: lon1=pi/6
+  REAL(rstd), SAVE  :: lon2=-pi/6
+  REAL(rstd), SAVE  :: latc1=0.
+  REAL(rstd), SAVE  :: latc2=0.
+  REAL(rstd), SAVE  :: lonc1=5*pi/6
+  REAL(rstd), SAVE  :: lonc2=7*pi/6
+  REAL(rstd), SAVE  :: zt=1000.0
+  REAL(rstd), SAVE  :: rt
+  REAL(rstd), SAVE  :: zc=5000.0
 
   PUBLIC etat0
@@ -39 +39 @@
     REAL(rstd),POINTER :: q(:,:,:)
     INTEGER :: ind
+    
+    R0=radius*0.5
+    rt=radius*0.5
     
     DO ind=1,ndomain