Changeset 1220 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/sw_aeroAR4.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/sw_aeroAR4.F90 (modified) (3 diffs)

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

@@ -56 +56 @@
   !* ARGUMENTS:
   !
-  REAL*8 PSCT  ! constante solaire (valeur conseillee: 1370)
-
-  REAL*8 PPSOL(KDLON)        ! SURFACE PRESSURE (PA)
-  REAL*8 PDP(KDLON,KFLEV)    ! LAYER THICKNESS (PA)
-  REAL*8 PPMB(KDLON,KFLEV+1) ! HALF-LEVEL PRESSURE (MB)
-
-  REAL*8 PRMU0(KDLON)  ! COSINE OF ZENITHAL ANGLE
-  REAL*8 PFRAC(KDLON)  ! fraction de la journee
-
-  REAL*8 PTAVE(KDLON,KFLEV)  ! LAYER TEMPERATURE (K)
-  REAL*8 PWV(KDLON,KFLEV)    ! SPECIFI! HUMIDITY (KG/KG)
-  REAL*8 PQS(KDLON,KFLEV)    ! SATURATED WATER VAPOUR (KG/KG)
-  REAL*8 POZON(KDLON,KFLEV)  ! OZONE CONCENTRATION (KG/KG)
-  REAL*8 PAER(KDLON,KFLEV,5) ! AEROSOLS' OPTICAL THICKNESS
-
-  REAL*8 PALBD(KDLON,2)  ! albedo du sol (lumiere diffuse)
-  REAL*8 PALBP(KDLON,2)  ! albedo du sol (lumiere parallele)
-
-  REAL*8 PCLDSW(KDLON,KFLEV)    ! CLOUD FRACTION
-  REAL*8 PTAU(KDLON,2,KFLEV)    ! CLOUD OPTICAL THICKNESS
-  REAL*8 PCG(KDLON,2,KFLEV)     ! ASYMETRY FACTOR
-  REAL*8 POMEGA(KDLON,2,KFLEV)  ! SINGLE SCATTERING ALBEDO
-
-  REAL*8 PHEAT(KDLON,KFLEV) ! SHORTWAVE HEATING (K/DAY)
-  REAL*8 PHEAT0(KDLON,KFLEV)! SHORTWAVE HEATING (K/DAY) clear-sky
-  REAL*8 PALBPLA(KDLON)     ! PLANETARY ALBEDO
-  REAL*8 PTOPSW(KDLON)      ! SHORTWAVE FLUX AT T.O.A.
-  REAL*8 PSOLSW(KDLON)      ! SHORTWAVE FLUX AT SURFACE
-  REAL*8 PTOPSW0(KDLON)     ! SHORTWAVE FLUX AT T.O.A. (CLEAR-SKY)
-  REAL*8 PSOLSW0(KDLON)     ! SHORTWAVE FLUX AT SURFACE (CLEAR-SKY)
+  REAL(KIND=8) PSCT  ! constante solaire (valeur conseillee: 1370)
+
+  REAL(KIND=8) PPSOL(KDLON)        ! SURFACE PRESSURE (PA)
+  REAL(KIND=8) PDP(KDLON,KFLEV)    ! LAYER THICKNESS (PA)
+  REAL(KIND=8) PPMB(KDLON,KFLEV+1) ! HALF-LEVEL PRESSURE (MB)
+
+  REAL(KIND=8) PRMU0(KDLON)  ! COSINE OF ZENITHAL ANGLE
+  REAL(KIND=8) PFRAC(KDLON)  ! fraction de la journee
+
+  REAL(KIND=8) PTAVE(KDLON,KFLEV)  ! LAYER TEMPERATURE (K)
+  REAL(KIND=8) PWV(KDLON,KFLEV)    ! SPECIFI! HUMIDITY (KG/KG)
+  REAL(KIND=8) PQS(KDLON,KFLEV)    ! SATURATED WATER VAPOUR (KG/KG)
+  REAL(KIND=8) POZON(KDLON,KFLEV)  ! OZONE CONCENTRATION (KG/KG)
+  REAL(KIND=8) PAER(KDLON,KFLEV,5) ! AEROSOLS' OPTICAL THICKNESS
+
+  REAL(KIND=8) PALBD(KDLON,2)  ! albedo du sol (lumiere diffuse)
+  REAL(KIND=8) PALBP(KDLON,2)  ! albedo du sol (lumiere parallele)
+
+  REAL(KIND=8) PCLDSW(KDLON,KFLEV)    ! CLOUD FRACTION
+  REAL(KIND=8) PTAU(KDLON,2,KFLEV)    ! CLOUD OPTICAL THICKNESS
+  REAL(KIND=8) PCG(KDLON,2,KFLEV)     ! ASYMETRY FACTOR
+  REAL(KIND=8) POMEGA(KDLON,2,KFLEV)  ! SINGLE SCATTERING ALBEDO
+
+  REAL(KIND=8) PHEAT(KDLON,KFLEV) ! SHORTWAVE HEATING (K/DAY)
+  REAL(KIND=8) PHEAT0(KDLON,KFLEV)! SHORTWAVE HEATING (K/DAY) clear-sky
+  REAL(KIND=8) PALBPLA(KDLON)     ! PLANETARY ALBEDO
+  REAL(KIND=8) PTOPSW(KDLON)      ! SHORTWAVE FLUX AT T.O.A.
+  REAL(KIND=8) PSOLSW(KDLON)      ! SHORTWAVE FLUX AT SURFACE
+  REAL(KIND=8) PTOPSW0(KDLON)     ! SHORTWAVE FLUX AT T.O.A. (CLEAR-SKY)
+  REAL(KIND=8) PSOLSW0(KDLON)     ! SHORTWAVE FLUX AT SURFACE (CLEAR-SKY)
   !
   !* LOCAL VARIABLES:
@@ -91 +91 @@
   real, parameter:: dobson_u = 2.1415e-05 ! Dobson unit, in kg m-2
 
-  REAL(kind=8) ZOZ(KDLON,KFLEV)
+  REAL(KIND=8) ZOZ(KDLON,KFLEV)
   ! column-density of ozone in layer, in kilo-Dobsons
 
-  REAL*8 ZAKI(KDLON,2)     
-  REAL*8 ZCLD(KDLON,KFLEV)
-  REAL*8 ZCLEAR(KDLON) 
-  REAL*8 ZDSIG(KDLON,KFLEV)
-  REAL*8 ZFACT(KDLON)
-  REAL*8 ZFD(KDLON,KFLEV+1)
-  REAL*8 ZFDOWN(KDLON,KFLEV+1)
-  REAL*8 ZFU(KDLON,KFLEV+1)
-  REAL*8 ZFUP(KDLON,KFLEV+1)
-  REAL*8 ZRMU(KDLON)
-  REAL*8 ZSEC(KDLON)
-  REAL*8 ZUD(KDLON,5,KFLEV+1)
-  REAL*8 ZCLDSW0(KDLON,KFLEV)
-
-  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(KIND=8) ZAKI(KDLON,2)     
+  REAL(KIND=8) ZCLD(KDLON,KFLEV)
+  REAL(KIND=8) ZCLEAR(KDLON) 
+  REAL(KIND=8) ZDSIG(KDLON,KFLEV)
+  REAL(KIND=8) ZFACT(KDLON)
+  REAL(KIND=8) ZFD(KDLON,KFLEV+1)
+  REAL(KIND=8) ZFDOWN(KDLON,KFLEV+1)
+  REAL(KIND=8) ZFU(KDLON,KFLEV+1)
+  REAL(KIND=8) ZFUP(KDLON,KFLEV+1)
+  REAL(KIND=8) ZRMU(KDLON)
+  REAL(KIND=8) ZSEC(KDLON)
+  REAL(KIND=8) ZUD(KDLON,5,KFLEV+1)
+  REAL(KIND=8) ZCLDSW0(KDLON,KFLEV)
+
+  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)
 
   INTEGER inu, jl, jk, i, k, kpl1
@@ -126 +126 @@
 
   !jq-Introduced for aerosol forcings
-  REAL*8, SAVE :: flag_aer
+  REAL(KIND=8), SAVE :: flag_aer
   !$OMP THREADPRIVATE(flag_aer)
 
   LOGICAL ok_ade, ok_aie    ! use aerosol forcings or not?
-  REAL*8 tauaero(kdlon,kflev,9,2)  ! aerosol optical properties
-  REAL*8 pizaero(kdlon,kflev,9,2)  ! (see aeropt.F)
-  REAL*8 cgaero(kdlon,kflev,9,2)   ! -"-
-  REAL*8 PTAUA(KDLON,2,KFLEV)    ! CLOUD OPTICAL THICKNESS (pre-industrial value)
-  REAL*8 POMEGAA(KDLON,2,KFLEV)  ! SINGLE SCATTERING ALBEDO
-  REAL*8 PTOPSWADAERO(KDLON)     ! SHORTWAVE FLUX AT T.O.A.(+AEROSOL DIR)
-  REAL*8 PSOLSWADAERO(KDLON)     ! SHORTWAVE FLUX AT SURFACE(+AEROSOL DIR)
-  REAL*8 PTOPSWAD0AERO(KDLON)     ! SHORTWAVE FLUX AT T.O.A.(+AEROSOL DIR)
-  REAL*8 PSOLSWAD0AERO(KDLON)     ! SHORTWAVE FLUX AT SURFACE(+AEROSOL DIR)
-  REAL*8 PTOPSWAIAERO(KDLON)     ! SHORTWAVE FLUX AT T.O.A.(+AEROSOL IND)
-  REAL*8 PSOLSWAIAERO(KDLON)     ! SHORTWAVE FLUX AT SURFACE(+AEROSOL IND)
-  REAL*8 PTOPSWAERO(KDLON,9)
-  REAL*8 PTOPSW0AERO(KDLON,9)
-  REAL*8 PSOLSWAERO(KDLON,9)
-  REAL*8 PSOLSW0AERO(KDLON,9)
+  REAL(KIND=8) tauaero(kdlon,kflev,9,2)  ! aerosol optical properties
+  REAL(KIND=8) pizaero(kdlon,kflev,9,2)  ! (see aeropt.F)
+  REAL(KIND=8) cgaero(kdlon,kflev,9,2)   ! -"-
+  REAL(KIND=8) PTAUA(KDLON,2,KFLEV)    ! CLOUD OPTICAL THICKNESS (pre-industrial value)
+  REAL(KIND=8) POMEGAA(KDLON,2,KFLEV)  ! SINGLE SCATTERING ALBEDO
+  REAL(KIND=8) PTOPSWADAERO(KDLON)     ! SHORTWAVE FLUX AT T.O.A.(+AEROSOL DIR)
+  REAL(KIND=8) PSOLSWADAERO(KDLON)     ! SHORTWAVE FLUX AT SURFACE(+AEROSOL DIR)
+  REAL(KIND=8) PTOPSWAD0AERO(KDLON)     ! SHORTWAVE FLUX AT T.O.A.(+AEROSOL DIR)
+  REAL(KIND=8) PSOLSWAD0AERO(KDLON)     ! SHORTWAVE FLUX AT SURFACE(+AEROSOL DIR)
+  REAL(KIND=8) PTOPSWAIAERO(KDLON)     ! SHORTWAVE FLUX AT T.O.A.(+AEROSOL IND)
+  REAL(KIND=8) PSOLSWAIAERO(KDLON)     ! SHORTWAVE FLUX AT SURFACE(+AEROSOL IND)
+  REAL(KIND=8) PTOPSWAERO(KDLON,9)
+  REAL(KIND=8) PTOPSW0AERO(KDLON,9)
+  REAL(KIND=8) PSOLSWAERO(KDLON,9)
+  REAL(KIND=8) PSOLSW0AERO(KDLON,9)
 
   !jq - Fluxes including aerosol effects
-  REAL*8,ALLOCATABLE,SAVE :: ZFSUPAD_AERO(:,:)
+  REAL(KIND=8),ALLOCATABLE,SAVE :: ZFSUPAD_AERO(:,:)
   !$OMP THREADPRIVATE(ZFSUPAD_AERO)
-  REAL*8,ALLOCATABLE,SAVE :: ZFSDNAD_AERO(:,:)
+  REAL(KIND=8),ALLOCATABLE,SAVE :: ZFSDNAD_AERO(:,:)
   !$OMP THREADPRIVATE(ZFSDNAD_AERO)
   !jq - Fluxes including aerosol effects
-  REAL*8,ALLOCATABLE,SAVE :: ZFSUPAD0_AERO(:,:)
+  REAL(KIND=8),ALLOCATABLE,SAVE :: ZFSUPAD0_AERO(:,:)
   !$OMP THREADPRIVATE(ZFSUPAD0_AERO)
-  REAL*8,ALLOCATABLE,SAVE :: ZFSDNAD0_AERO(:,:)
+  REAL(KIND=8),ALLOCATABLE,SAVE :: ZFSDNAD0_AERO(:,:)
   !$OMP THREADPRIVATE(ZFSDNAD0_AERO)
-  REAL*8,ALLOCATABLE,SAVE :: ZFSUPAI_AERO(:,:)
+  REAL(KIND=8),ALLOCATABLE,SAVE :: ZFSUPAI_AERO(:,:)
   !$OMP THREADPRIVATE(ZFSUPAI_AERO)
-  REAL*8,ALLOCATABLE,SAVE :: ZFSDNAI_AERO(:,:)
+  REAL(KIND=8),ALLOCATABLE,SAVE :: ZFSDNAI_AERO(:,:)
   !$OMP THREADPRIVATE(ZFSDNAI_AERO)
-  REAL*8,ALLOCATABLE,SAVE ::  ZFSUP_AERO(:,:,:)
+  REAL(KIND=8),ALLOCATABLE,SAVE ::  ZFSUP_AERO(:,:,:)
   !$OMP THREADPRIVATE(ZFSUP_AERO)
-  REAL*8,ALLOCATABLE,SAVE ::  ZFSDN_AERO(:,:,:)
+  REAL(KIND=8),ALLOCATABLE,SAVE ::  ZFSDN_AERO(:,:,:)
   !$OMP THREADPRIVATE(ZFSDN_AERO)
-  REAL*8,ALLOCATABLE,SAVE ::  ZFSUP0_AERO(:,:,:)
+  REAL(KIND=8),ALLOCATABLE,SAVE ::  ZFSUP0_AERO(:,:,:)
   !$OMP THREADPRIVATE(ZFSUP0_AERO)
-  REAL*8,ALLOCATABLE,SAVE ::  ZFSDN0_AERO(:,:,:)
+  REAL(KIND=8),ALLOCATABLE,SAVE ::  ZFSDN0_AERO(:,:,:)
   !$OMP THREADPRIVATE(ZFSDN0_AERO)