Changeset 1220 for LMDZ4/branches/LMDZ4-dev/libf/dyn3d/fyhyp.F

Timestamp:

Aug 5, 2009, 4:38:34 PM (15 years ago)

Author:

lguez

Message:

-- Replaced "integer*4" declarations by "integer", "real*8" by

"real(kind=8)" and "real*4" by "real". Note that these are the only
modifications in the files "radiation_AR4.F" and "sw_aeroAR4.F90".

-- Corrected the kind of arguments to "max" and "min".

-- Replaced "nH" edit descriptors, which is a deleted feature in

Fortran 95, by character strings.

-- "regr_lat_time_climoz" now allows the pressure coordinate in the

input file to be in descending order.

-- Replaced call to not standard function "float" by call to intrinsic

function "real".

-- Included file "radepsi.h" in "physiq" was not used. Removed it.

The following set of modifications is related to the management of time.

-- In "gcm", "leapfrog" and "sortvarc0", "day_ini" was defined as 1

plus number of days between the reference date "(annee_ref,
day_ref)" and the first day of the current simulation. Changed
definition: "(annee_ref, day_ini)" is the first day of the current
simulation. There is an accompanying modification for "day_end".

-- Corrected bug in call to "ioconf_startdate" in "gcm".

-- Added call to "ioconf_calendar" in "create_etat0_limit".

File:

• LMDZ4/branches/LMDZ4-dev/libf/dyn3d/fyhyp.F (modified) (2 diffs)

LMDZ4/branches/LMDZ4-dev/libf/dyn3d/fyhyp.F

@@ -50 +50 @@
      
        REAL   dzoom
-       REAL*8 ylat(jjp1), yprim(jjp1)
-       REAL*8 yuv
-       REAL*8 yt(0:nmax2)
-       REAL*8 fhyp(0:nmax2),beta,Ytprim(0:nmax2),fxm(0:nmax2)
+       REAL(KIND=8) ylat(jjp1), yprim(jjp1)
+       REAL(KIND=8) yuv
+       REAL(KIND=8) yt(0:nmax2)
+       REAL(KIND=8) fhyp(0:nmax2),beta,Ytprim(0:nmax2),fxm(0:nmax2)
        SAVE Ytprim, yt,Yf
-       REAL*8 Yf(0:nmax2),yypr(0:nmax2)
-       REAL*8 yvrai(jjp1), yprimm(jjp1),ylatt(jjp1)
-       REAL*8 pi,depi,pis2,epsilon,y0,pisjm
-       REAL*8 yo1,yi,ylon2,ymoy,Yprimin,champmin,champmax
-       REAL*8 yfi,Yf1,ffdy
-       REAL*8 ypn,deply,y00
+       REAL(KIND=8) Yf(0:nmax2),yypr(0:nmax2)
+       REAL(KIND=8) yvrai(jjp1), yprimm(jjp1),ylatt(jjp1)
+       REAL(KIND=8) pi,depi,pis2,epsilon,y0,pisjm
+       REAL(KIND=8) yo1,yi,ylon2,ymoy,Yprimin,champmin,champmax
+       REAL(KIND=8) yfi,Yf1,ffdy
+       REAL(KIND=8) ypn,deply,y00
        SAVE y00, deply
 
@@ -66 +66 @@
        INTEGER jpn,jjpn
        SAVE jpn
-       REAL*8 a0,a1,a2,a3,yi2,heavyy0,heavyy0m
-       REAL*8 fa(0:nmax2),fb(0:nmax2)
+       REAL(KIND=8) a0,a1,a2,a3,yi2,heavyy0,heavyy0m
+       REAL(KIND=8) fa(0:nmax2),fb(0:nmax2)
        REAL y0min,y0max
 
-       REAL*8     heavyside
+       REAL(KIND=8)     heavyside
 
        pi       = 2. * ASIN(1.)