Changeset 1391 for trunk/LMDZ.COMMON/libf/dyn3d_common/ran1.F

Timestamp:

Mar 6, 2015, 3:12:12 PM (10 years ago)

Author:

emillour

Message:

Common dynamical core:
Updates in the dynamics to keeup up with updates in LMDZ5
(up to LMDZ5 trunk rev 2200):

• compilation:
• create_make_gcm : added processing of .f & .f90 files (not just .F and .F90)
• makelmdz: add "mix" option for -io (ouptut with both IOIPSL and XIOS)
• makelmdz_fcm: add "mix" option for -io

• filtrez:
• acc.F and eigen.F : add "implicit none" and variable declarations

• bibio:
• handle_err_m.F90: replace "stop" with call to abort_gcm()
• i1mach.F, j4save.F: add "implicit none" and variable declarations
• xercnt.F, xermsg.F, xerprn.F, xersve.F, xgetua.F: add "implicit none" and variable declarations

• dyn3d_common:
• disvert.F90 : added comments on meaning of "pa" variable
• grid_atob.F : better control on level of default ouputs
• infotrac.F90: update Earth-specific stuff (nqo water tracers)
• interpre.F: correction on the size of input array w
• juldate.F, massbar.F, ppm3d.F, ran1.F: add "implicit none" and variable declarations
• sortvarc.F: code cleanup
• iniacademic.F90: cleanup and extra sanity check.

• dyn3d:
• abort_gcm.F: additions for XIOS
• conf_gcm.F90: transformed to free form from conf_gcm.F
• gcm.F: added test to check that iphysiq is a multiple of iperiod
• getparam.F90, guidz_mod.F: update from LMDZ5
• integrd.F: replace stop with call_abort()

• dyn3dpar:
• abort_gcm.F: minor cleanup
• gcm.F: added test to check that iphysiq is a multiple of iperiod
• getparam.F90, guide_p_mod.F90: update from LMDZ5
• integrd_p.F: abort with call_abort when there is negative surface pressure
• leapfrog_p.F: add INCA specific stuff to keep up with current LMDZ5
• conf_gcm.F90: transformed to free form from conf_gcm.F

EM

File:

• trunk/LMDZ.COMMON/libf/dyn3d_common/ran1.F (modified) (1 diff)

trunk/LMDZ.COMMON/libf/dyn3d_common/ran1.F

@@ -1 +1 @@
 !
-! $Id: ran1.F 1403 2010-07-01 09:02:53Z fairhead $
+! $Id: ran1.F 2197 2015-02-09 07:13:05Z emillour $
 !
       FUNCTION RAN1(IDUM)
-      DIMENSION R(97)
-      save r
-      save iff,ix1,ix2,ix3
-      PARAMETER (M1=259200,IA1=7141,IC1=54773,RM1=3.8580247E-6)
-      PARAMETER (M2=134456,IA2=8121,IC2=28411,RM2=7.4373773E-6)
-      PARAMETER (M3=243000,IA3=4561,IC3=51349)
-      DATA IFF /0/
+      IMPLICIT NONE
+      REAL RAN1
+      REAL,SAVE :: R(97)
+      REAL,PARAMETER :: RM1=3.8580247E-6,RM2=7.4373773E-6
+      INTEGER,SAVE :: IFF=0
+      integer,save :: ix1,ix2,ix3
+      INTEGER,PARAMETER :: M1=259200,IA1=7141,IC1=54773
+      INTEGER,PARAMETER :: M2=134456,IA2=8121,IC2=28411
+      INTEGER,PARAMETER :: M3=243000,IA3=4561,IC3=51349
+      INTEGER :: IDUM,J
+
       IF (IDUM.LT.0.OR.IFF.EQ.0) THEN
         IFF=1