Changeset 2814

Timestamp:

Mar 6, 2017, 5:38:36 PM (8 years ago)

Author:

jyg

Message:

Cosmetic changes to fisrtilp

File:

• LMDZ5/trunk/libf/phylmd/fisrtilp.F90 (modified) (6 diffs)

LMDZ5/trunk/libf/phylmd/fisrtilp.F90

@@ -53 +53 @@
   ! Principaux inputs:
   !
-  REAL dtime ! intervalle du temps (s)
-  REAL paprs(klon,klev+1) ! pression a inter-couche
-  REAL pplay(klon,klev) ! pression au milieu de couche
-  REAL t(klon,klev) ! temperature (K)
-  REAL q(klon,klev) ! humidite specifique (kg/kg)
+  REAL, INTENT(IN)                              :: dtime  ! intervalle du temps (s)
+  REAL, DIMENSION(klon,klev+1),    INTENT(IN)   :: paprs  ! pression a inter-couche
+  REAL, DIMENSION(klon,klev),      INTENT(IN)   :: pplay  ! pression au milieu de couche
+  REAL, DIMENSION(klon,klev),      INTENT(IN)   :: t      ! temperature (K)
+  REAL, DIMENSION(klon,klev),      INTENT(IN)   :: q      ! humidite specifique (kg/kg)
+  LOGICAL, DIMENSION(klon,klev),   INTENT(IN)   :: ptconv ! points ou le schema de conv. prof. est actif
+  INTEGER,                         INTENT(IN)   :: iflag_cld_th
+  INTEGER,                         INTENT(IN)   :: iflag_ice_thermo
+  !
+  ! Inputs lies aux thermiques
+  ! 
+  REAL, DIMENSION(klon,klev),      INTENT(IN)   :: ztv
+  REAL, DIMENSION(klon,klev),      INTENT(IN)   :: zqta, fraca 
+  REAL, DIMENSION(klon,klev),      INTENT(IN)   :: zpspsk, ztla
+  REAL, DIMENSION(klon,klev),      INTENT(IN)   :: zthl
+  !
+  !  Input/output
+  REAL, DIMENSION(klon,klev),      INTENT(INOUT):: ratqs  ! determine la largeur de distribution de vapeur
   !
   ! Principaux outputs:
   !
-  REAL d_t(klon,klev) ! incrementation de la temperature (K)
-  REAL d_q(klon,klev) ! incrementation de la vapeur d'eau
-  REAL d_ql(klon,klev) ! incrementation de l'eau liquide
-  REAL d_qi(klon,klev) ! incrementation de l'eau glace
-  REAL rneb(klon,klev) ! fraction nuageuse
-  REAL radliq(klon,klev) ! eau liquide utilisee dans rayonnements
-  REAL rhcl(klon,klev) ! humidite relative en ciel clair
-  REAL rain(klon) ! pluies (mm/s)
-  REAL snow(klon) ! neige (mm/s)
-  REAL prfl(klon,klev+1) ! flux d'eau precipitante aux interfaces (kg/m2/s)
-  REAL psfl(klon,klev+1) ! flux d'eau precipitante aux interfaces (kg/m2/s) 
-  !
-  ! Autres arguments
-  ! 
-  REAL ztv(klon,klev)
-  REAL zqta(klon,klev),fraca(klon,klev) 
-  REAL sigma1(klon,klev),sigma2(klon,klev)
-  REAL qltot(klon,klev),ctot(klon,klev)
-  REAL zpspsk(klon,klev),ztla(klon,klev)
-  REAL zthl(klon,klev)
-  REAL ztfondue, qsl, qsi
-
-  logical lognormale(klon)
-  logical ice_thermo
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: d_t  ! incrementation de la temperature (K)
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: d_q  ! incrementation de la vapeur d'eau
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: d_ql ! incrementation de l'eau liquide
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: d_qi ! incrementation de l'eau glace
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: rneb ! fraction nuageuse
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: radliq ! eau liquide utilisee dans rayonnements
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: rhcl ! humidite relative en ciel clair
+  REAL, DIMENSION(klon),           INTENT(OUT)  :: rain
+  REAL, DIMENSION(klon),           INTENT(OUT)  :: snow
+  REAL, DIMENSION(klon,klev+1),    INTENT(OUT)  :: prfl
+  REAL, DIMENSION(klon,klev+1),    INTENT(OUT)  :: psfl 
 
   !AA
   ! Coeffients de fraction lessivee : pour OFF-LINE
   !
-  REAL pfrac_nucl(klon,klev)
-  REAL pfrac_1nucl(klon,klev)
-  REAL pfrac_impa(klon,klev)
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: pfrac_nucl
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: pfrac_1nucl
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: pfrac_impa
   !
   ! Fraction d'aerosols lessivee par impaction et par nucleation
   ! POur ON-LINE
   !
-  REAL frac_impa(klon,klev)
-  REAL frac_nucl(klon,klev)
-  real zct      ,zcl
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: frac_impa
+  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: frac_nucl
   !AA
+  ! --------------------------------------------------------------------------------
   !
   ! Options du programme:
@@ -107 +107 @@
 
   INTEGER ninter ! sous-intervals pour la precipitation
-  INTEGER ncoreczq  
-  INTEGER iflag_cld_th
-  INTEGER iflag_ice_thermo
   PARAMETER (ninter=5)
   LOGICAL evap_prec ! evaporation de la pluie
   PARAMETER (evap_prec=.TRUE.)
-  REAL ratqs(klon,klev) ! determine la largeur de distribution de vapeur
-  logical ptconv(klon,klev) ! determine la largeur de distribution de vapeur
+  !
+  LOGICAL cpartiel ! condensation partielle
+  PARAMETER (cpartiel=.TRUE.)
+  REAL t_coup
+  PARAMETER (t_coup=234.0)
+  REAL DDT0
+  PARAMETER (DDT0=.01)
+  REAL ztfondue
+  PARAMETER (ztfondue=278.15)
+  ! --------------------------------------------------------------------------------
+  !
+  ! Variables locales:
+  !
+  INTEGER i, k, n, kk
+  REAL qsl, qsi
+  real zct      ,zcl
+  INTEGER ncoreczq  
+  REAL ctot(klon,klev)
+  REAL zqs(klon), zdqs(klon), zdelta, zcor, zcvm5  
+  REAL zdqsdT_raw(klon)
+  REAL Tbef(klon),qlbef(klon),DT(klon),num,denom
+
+  logical lognormale(klon)
+  logical ice_thermo
+  LOGICAL convergence(klon)
+  INTEGER n_i(klon), iter
+  REAL cste
 
   real zpdf_sig(klon),zpdf_k(klon),zpdf_delta(klon)
@@ -120 +142 @@
   real erf   
   REAL qcloud(klon)
-  !
-  LOGICAL cpartiel ! condensation partielle
-  PARAMETER (cpartiel=.TRUE.)
-  REAL t_coup
-  PARAMETER (t_coup=234.0)
-  !
-  ! Variables locales:
-  !
-  INTEGER i, k, n, kk
-  REAL zqs(klon), zdqs(klon), zdelta, zcor, zcvm5  
-  REAL zdqsdT_raw(klon)
-  REAL Tbef(klon),qlbef(klon),DT(klon),num,denom
-  LOGICAL convergence(klon)
-  REAL DDT0
-  PARAMETER (DDT0=.01)
-  INTEGER n_i(klon), iter
-  REAL cste
  
   REAL zrfl(klon), zrfln(klon), zqev, zqevt 
@@ -150 +155 @@
   REAL zdz(klon),zrho(klon),ztot      , zrhol(klon)
   REAL zchau      ,zfroi      ,zfice(klon),zneb(klon)
-  REAL zmelt, zpluie, zice, zcondold
-  PARAMETER (ztfondue=278.15)
+  REAL zmelt, zpluie, zice
   REAL dzfice(klon)
   REAL zsolid
@@ -434 +438 @@
                 zt(i) = zt(i) - zqev &
                      * RLVTT/RCPD/(1.0+RVTMP2*(zq(i)+zmqc(i)))
-!!JLD debut de partie a supprimer a therme
+!!JLD debut de partie a supprimer a terme
             else ! if  (fl_cor_ebil .GT. 0)
                 ! Calcul de l'evaporation du flux de precip herite
@@ -457 +461 @@
                      * RLVTT/RCPD/(1.0+RVTMP2*zq(i))
               end if ! if  (fl_cor_ebil .GT. 0)
-!!JLD fin de partie a supprimer a therme
+!!JLD fin de partie a supprimer a terme
                 zrfl(i) = zrfln(i)
                 zifl(i) = 0.