Changeset 1220 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/radlwsw.F90

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/phylmd/radlwsw.F90 (modified) (1 diff)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/radlwsw.F90

@@ -137 +137 @@
 
 ! Local variables
-  REAL*8 ZFSUP(KDLON,KFLEV+1)
-  REAL*8 ZFSDN(KDLON,KFLEV+1)
-  REAL*8 ZFSUP0(KDLON,KFLEV+1)
-  REAL*8 ZFSDN0(KDLON,KFLEV+1)
-  REAL*8 ZFLUP(KDLON,KFLEV+1)
-  REAL*8 ZFLDN(KDLON,KFLEV+1)
-  REAL*8 ZFLUP0(KDLON,KFLEV+1)
-  REAL*8 ZFLDN0(KDLON,KFLEV+1)
-  REAL*8 zx_alpha1, zx_alpha2
+  REAL(KIND=8) ZFSUP(KDLON,KFLEV+1)
+  REAL(KIND=8) ZFSDN(KDLON,KFLEV+1)
+  REAL(KIND=8) ZFSUP0(KDLON,KFLEV+1)
+  REAL(KIND=8) ZFSDN0(KDLON,KFLEV+1)
+  REAL(KIND=8) ZFLUP(KDLON,KFLEV+1)
+  REAL(KIND=8) ZFLDN(KDLON,KFLEV+1)
+  REAL(KIND=8) ZFLUP0(KDLON,KFLEV+1)
+  REAL(KIND=8) ZFLDN0(KDLON,KFLEV+1)
+  REAL(KIND=8) zx_alpha1, zx_alpha2
   INTEGER k, kk, i, j, iof, nb_gr
-  REAL*8 PSCT
-  REAL*8 PALBD(kdlon,2), PALBP(kdlon,2)
-  REAL*8 PEMIS(kdlon), PDT0(kdlon), PVIEW(kdlon)
-  REAL*8 PPSOL(kdlon), PDP(kdlon,KLEV)
-  REAL*8 PTL(kdlon,kflev+1), PPMB(kdlon,kflev+1)
-  REAL*8 PTAVE(kdlon,kflev)
-  REAL*8 PWV(kdlon,kflev), PQS(kdlon,kflev)
+  REAL(KIND=8) PSCT
+  REAL(KIND=8) PALBD(kdlon,2), PALBP(kdlon,2)
+  REAL(KIND=8) PEMIS(kdlon), PDT0(kdlon), PVIEW(kdlon)
+  REAL(KIND=8) PPSOL(kdlon), PDP(kdlon,KLEV)
+  REAL(KIND=8) PTL(kdlon,kflev+1), PPMB(kdlon,kflev+1)
+  REAL(KIND=8) PTAVE(kdlon,kflev)
+  REAL(KIND=8) PWV(kdlon,kflev), PQS(kdlon,kflev)
   real(kind=8) POZON(kdlon,kflev) ! mass fraction of ozone
-  REAL*8 PAER(kdlon,kflev,5)
-  REAL*8 PCLDLD(kdlon,kflev)
-  REAL*8 PCLDLU(kdlon,kflev)
-  REAL*8 PCLDSW(kdlon,kflev)
-  REAL*8 PTAU(kdlon,2,kflev)
-  REAL*8 POMEGA(kdlon,2,kflev)
-  REAL*8 PCG(kdlon,2,kflev)
-  REAL*8 zfract(kdlon), zrmu0(kdlon), zdist
-  REAL*8 zheat(kdlon,kflev), zcool(kdlon,kflev)
-  REAL*8 zheat0(kdlon,kflev), zcool0(kdlon,kflev)
-  REAL*8 ztopsw(kdlon), ztoplw(kdlon)
-  REAL*8 zsolsw(kdlon), zsollw(kdlon), zalbpla(kdlon)
-  REAL*8 zsollwdown(kdlon)
-  REAL*8 ztopsw0(kdlon), ztoplw0(kdlon)
-  REAL*8 zsolsw0(kdlon), zsollw0(kdlon)
-  REAL*8 zznormcp
-  REAL*8 tauaero(kdlon,kflev,9,2)                     ! aer opt properties
-  REAL*8 pizaero(kdlon,kflev,9,2)
-  REAL*8 cgaero(kdlon,kflev,9,2)
-  REAL*8 PTAUA(kdlon,2,kflev)                         ! present-day value of cloud opt thickness (PTAU is pre-industrial value), local use
-  REAL*8 POMEGAA(kdlon,2,kflev)                       ! dito for single scatt albedo
-  REAL*8 ztopswadaero(kdlon), zsolswadaero(kdlon)     ! Aerosol direct forcing at TOAand surface
-  REAL*8 ztopswad0aero(kdlon), zsolswad0aero(kdlon)   ! Aerosol direct forcing at TOAand surface
-  REAL*8 ztopswaiaero(kdlon), zsolswaiaero(kdlon)     ! dito, indirect
-  REAL*8 ztopsw_aero(kdlon,9), ztopsw0_aero(kdlon,9)
-  REAL*8 zsolsw_aero(kdlon,9), zsolsw0_aero(kdlon,9)
+  REAL(KIND=8) PAER(kdlon,kflev,5)
+  REAL(KIND=8) PCLDLD(kdlon,kflev)
+  REAL(KIND=8) PCLDLU(kdlon,kflev)
+  REAL(KIND=8) PCLDSW(kdlon,kflev)
+  REAL(KIND=8) PTAU(kdlon,2,kflev)
+  REAL(KIND=8) POMEGA(kdlon,2,kflev)
+  REAL(KIND=8) PCG(kdlon,2,kflev)
+  REAL(KIND=8) zfract(kdlon), zrmu0(kdlon), zdist
+  REAL(KIND=8) zheat(kdlon,kflev), zcool(kdlon,kflev)
+  REAL(KIND=8) zheat0(kdlon,kflev), zcool0(kdlon,kflev)
+  REAL(KIND=8) ztopsw(kdlon), ztoplw(kdlon)
+  REAL(KIND=8) zsolsw(kdlon), zsollw(kdlon), zalbpla(kdlon)
+  REAL(KIND=8) zsollwdown(kdlon)
+  REAL(KIND=8) ztopsw0(kdlon), ztoplw0(kdlon)
+  REAL(KIND=8) zsolsw0(kdlon), zsollw0(kdlon)
+  REAL(KIND=8) zznormcp
+  REAL(KIND=8) tauaero(kdlon,kflev,9,2)                     ! aer opt properties
+  REAL(KIND=8) pizaero(kdlon,kflev,9,2)
+  REAL(KIND=8) cgaero(kdlon,kflev,9,2)
+  REAL(KIND=8) PTAUA(kdlon,2,kflev)                         ! present-day value of cloud opt thickness (PTAU is pre-industrial value), local use
+  REAL(KIND=8) POMEGAA(kdlon,2,kflev)                       ! dito for single scatt albedo
+  REAL(KIND=8) ztopswadaero(kdlon), zsolswadaero(kdlon)     ! Aerosol direct forcing at TOAand surface
+  REAL(KIND=8) ztopswad0aero(kdlon), zsolswad0aero(kdlon)   ! Aerosol direct forcing at TOAand surface
+  REAL(KIND=8) ztopswaiaero(kdlon), zsolswaiaero(kdlon)     ! dito, indirect
+  REAL(KIND=8) ztopsw_aero(kdlon,9), ztopsw0_aero(kdlon,9)
+  REAL(KIND=8) zsolsw_aero(kdlon,9), zsolsw0_aero(kdlon,9)
   real, parameter:: dobson_u = 2.1415e-05 ! Dobson unit, in kg m-2
-
 
   ! initialisation