source: trunk/LMDZ.COMMON/libf/misc/ran1.F @ 3553

Last change on this file since 3553 was 1441, checked in by emillour, 10 years ago

Updates in common dynamics (seq and ) to keep up with updates
in LMDZ5 (up to LMDZ5 trunk, rev 2250):

  • compilation:
  • added test in grid/dimension/makdim to check that # of longitudes is a multiple of 8
  • dyn3d_common:

Bug correction concerning zoom (cf LMDZ5 rev 2218)

  • coefpoly.F becomes coefpoly_m.F90 (in misc)
  • fxhyp.F => fxhyp_m.F90 , fyhyp.F => fyhyp_m.F90
  • new routines for zoom: invert_zoom_x_m.F90 and principal_cshift_m.F90
  • inigeom.F adapted to new zoom definition routines
  • fluxstokenc.F : got rid of calls to initial0()
  • dyn3d:
  • advtrac.F90 : got rid of calls to initial0()
  • conf_gcm.F90 : cosmetic changes and change in default dzoomx,dzoomy values
  • guide_mod.F90 : followed updates from Earth Model
  • gcm.F is now gcm.F90
  • dyn3dpar:
  • advtrac_p.F90, covcont_p.F90, mod_hallo.F90 : cosmetic changes
  • conf_gcm.F90 : cosmetic and changed in default dzoomx,dzoomy values
  • parallel_lmdz.F90 : updates to keep up with Earth model
  • misc:
  • arth.F90 becomes arth_m.F90
  • wxios.F90 updated wrt Earth model changes
  • nrtype.F90 and coefpoly_m.F90 added
  • ran1.F, sort.F, minmax.F, minmax2.F, juldate.F moved over from dyn3d_common

EM
File size: 1.0 KB
Line 
1!
2! $Id: ran1.F 2197 2015-02-09 07:13:05Z emillour $
3!
4      FUNCTION RAN1(IDUM)
5      IMPLICIT NONE
6      REAL RAN1
7      REAL,SAVE :: R(97)
8      REAL,PARAMETER :: RM1=3.8580247E-6,RM2=7.4373773E-6
9      INTEGER,SAVE :: IFF=0
10      integer,save :: ix1,ix2,ix3
11      INTEGER,PARAMETER :: M1=259200,IA1=7141,IC1=54773
12      INTEGER,PARAMETER :: M2=134456,IA2=8121,IC2=28411
13      INTEGER,PARAMETER :: M3=243000,IA3=4561,IC3=51349
14      INTEGER :: IDUM,J
15
16      IF (IDUM.LT.0.OR.IFF.EQ.0) THEN
17        IFF=1
18        IX1=MOD(IC1-IDUM,M1)
19        IX1=MOD(IA1*IX1+IC1,M1)
20        IX2=MOD(IX1,M2)
21        IX1=MOD(IA1*IX1+IC1,M1)
22        IX3=MOD(IX1,M3)
23        DO 11 J=1,97
24          IX1=MOD(IA1*IX1+IC1,M1)
25          IX2=MOD(IA2*IX2+IC2,M2)
26          R(J)=(REAL(IX1)+REAL(IX2)*RM2)*RM1
2711      CONTINUE
28        IDUM=1
29      ENDIF
30      IX1=MOD(IA1*IX1+IC1,M1)
31      IX2=MOD(IA2*IX2+IC2,M2)
32      IX3=MOD(IA3*IX3+IC3,M3)
33      J=1+(97*IX3)/M3
34      IF(J.GT.97.OR.J.LT.1) stop 1
35      RAN1=R(J)
36      R(J)=(REAL(IX1)+REAL(IX2)*RM2)*RM1
37      RETURN
38      END
Note: See TracBrowser for help on using the repository browser.