Changeset 2784 for LMDZ5/trunk/libf/phylmd

Timestamp:

Jan 27, 2017, 9:41:14 AM (8 years ago)

Author:

acozic

Message:

Add new variable in inca (Strat Chemistry)

Location:

LMDZ5/trunk/libf/phylmd

Files:

• physiq_mod.F90 (modified) (3 diffs)
• phytrac_mod.F90 (modified) (4 diffs)
• tracinca_mod.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -1080 +1080 @@
     INTEGER :: nbtr_tmp ! Number of tracer inside concvl
     REAL, dimension(klon,klev) :: sh_in ! Specific humidity entering in phytrac
+    REAL, dimension(klon,klev) :: ch_in ! Condensed humidity entering in phytrac (eau liquide)
     integer iostat
 
@@ -4060 +4061 @@
     ELSE
        sh_in(:,:) = qx(:,:,ivap)
+       ch_in(:,:) = qx(:,:,iliq)
     ENDIF
 
@@ -4092 +4094 @@
          frac_impa,frac_nucl, beta_prec_fisrt,beta_prec, &
          presnivs, pphis,     pphi,     albsol1, &
-         sh_in,    rhcl,      cldfra,   rneb, &
+         sh_in,   ch_in,    rhcl,      cldfra,   rneb, &
          diafra,   cldliq,    itop_con, ibas_con, &
          pmflxr,   pmflxs,    prfl,     psfl, &

LMDZ5/trunk/libf/phylmd/phytrac_mod.F90

@@ -66 +66 @@
        frac_impa,frac_nucl,beta_fisrt,beta_v1,        &
        presnivs,  pphis,    pphi,     albsol,         &
-       sh,        rh,       cldfra,   rneb,           &
+       sh,        ch, rh,   cldfra,   rneb,           &
        diafra,    cldliq,   itop_con, ibas_con,       &
        pmflxr,    pmflxs,   prfl,     psfl,           &
@@ -139 +139 @@
     REAL,DIMENSION(klon,klev),INTENT(IN)   :: sh      ! humidite specifique
     REAL,DIMENSION(klon,klev),INTENT(IN)   :: rh      ! humidite relative
+    REAL,DIMENSION(klon,klev),INTENT(IN)   :: ch      ! eau liquide
     REAL,DIMENSION(klon,klev+1),INTENT(IN) :: paprs   ! pression pour chaque inter-couche (en Pa)
     REAL,DIMENSION(klon,klev),INTENT(IN)   :: pplay   ! pression pour le mileu de chaque couche (en Pa)
@@ -573 +574 @@
        write(lunout,*)  'flag_cvltr    = ', flag_cvltr
 
-       IF (lessivage.AND.config_inca.EQ.'inca') THEN
+       IF (lessivage .AND. type_trac .EQ. 'inca') THEN
           CALL abort_physic('phytrac', 'lessivage=T config_inca=inca impossible',1)
           STOP
@@ -1022 +1023 @@
             pdtphys,  t_seri,   paprs,          pplay,     &
             pmfu,     upwd,     ftsol,  pctsrf, pphis,     &
-            pphi,     albsol,   sh,             rh,        &
+            pphi,     albsol,   sh,    ch,     rh,        &
             cldfra,   rneb,     diafra,         cldliq,    &
             itop_con, ibas_con, pmflxr,         pmflxs,    &

LMDZ5/trunk/libf/phylmd/tracinca_mod.F90

@@ -36 +36 @@
        pdtphys,  t_seri,   paprs,          pplay,     &
        pmfu,     upwd,     ftsol,  pctsrf, pphis,     &
-       pphi,     albsol,   sh,             rh,        &
+       pphi,     albsol,   sh,             ch, rh,    &
        cldfra,   rneb,     diafra,         cldliq,    &
        itop_con, ibas_con, pmflxr,         pmflxs,    &
@@ -76 +76 @@
     REAL,DIMENSION(klon,klev),INTENT(IN)   :: t_seri  ! Temperature
     REAL,DIMENSION(klon,klev),INTENT(IN)   :: sh      ! humidite specifique
+    REAL,DIMENSION(klon,klev),INTENT(IN)   :: ch      ! eau liquide
     REAL,DIMENSION(klon,klev),INTENT(IN)   :: rh      ! humidite relative
     REAL,DIMENSION(klon,klev+1),INTENT(IN) :: paprs   ! pression pour chaque inter-couche (en Pa)
@@ -185 +186 @@
          t_seri,     & !tfld
          sh,         & !sh
+         ch,         & !ql
          rh,         & !rh
          nbp_lon,    & !nx