Changeset 1735

Timestamp:

Mar 29, 2013, 2:21:19 PM (11 years ago)

Author:

idelkadi

Message:

Corrections dans le cas des simulations de sensibilite ou l'on desactive (partout ou seulement dans la couche limite) les effets radiatifs des nuages. Les parametres de controle de ces simulations sont contenus dans le fichier beta_crf.data

File:

• LMDZ5/trunk/libf/phylmd/physiq.F (modified) (7 diffs)

LMDZ5/trunk/libf/phylmd/physiq.F

@@ -1187 +1187 @@
 c$OMP THREADPRIVATE(mskocean_beta)
       REAL, dimension(klon, klev) :: beta         ! facteur sur cldtaurad et cldemirad pour evaluer les retros liees aux CRF
-      REAL, dimension(klon, klev) :: cldtaurad    ! epaisseur optique pour radlwsw,COSP
-      REAL, dimension(klon, klev) :: cldtaupirad  ! epaisseur optique pour radlwsw,COSP cas pre-industrial
-      REAL, dimension(klon, klev) :: cldemirad    ! emissivite pour radlwsw,COSP 
+      REAL, dimension(klon, klev) :: cldtaurad    ! epaisseur optique pour radlwsw pour tester "CRF off"
+      REAL, dimension(klon, klev) :: cldtaupirad  ! epaisseur optique pour radlwsw pour tester "CRF off"
+      REAL, dimension(klon, klev) :: cldemirad    ! emissivite pour radlwsw pour tester "CRF off"
+      REAL, dimension(klon, klev) :: cldfrarad    ! fraction nuageuse
+
       INTEGER :: nbtr_tmp ! Number of tracer inside concvl
       REAL, dimension(klon,klev) :: sh_in ! Specific humidity entering in phytrac
@@ -3176 +3178 @@
       cldtaupirad = cldtaupi
       cldemirad   = cldemi
+     
 c
       if(lon1_beta.EQ.-180..AND.lon2_beta.EQ.180..AND.
@@ -3195 +3198 @@
         cldtaupirad(i,k) = cldtaupi(i,k) * beta(i,k)
         cldemirad(i,k)   = cldemi(i,k) * beta(i,k)
+        cldfrarad(i,k)   = cldfra(i,k) * beta(i,k)
        ENDDO
        ENDDO 
@@ -3218 +3222 @@
         cldtaupirad(i,k) = cldtaupi(i,k) * beta(i,k)
         cldemirad(i,k)   = cldemi(i,k) * beta(i,k)
+        cldfrarad(i,k)   = cldfra(i,k) * beta(i,k)
         endif
 c
@@ -3253 +3258 @@
      e        paprs, pplay,zxtsol,albsol1, albsol2, t_seri,q_seri,
      e        wo(:, :, 1),
-     e        cldfra, cldemirad, cldtaurad,
+     e        cldfrarad, cldemirad, cldtaurad,
      s        heat,heat0,cool,cool0,radsol,albpla,
      s        topsw,toplw,solsw,sollw,
@@ -3288 +3293 @@
      e        paprs, pplay,zxtsol,albsol1, albsol2, 
      e        t_seri,q_seri,wo,
-     e        cldfra, cldemirad, cldtaurad,
+     e        cldfrarad, cldemirad, cldtaurad,
      e        ok_ade, ok_aie, flag_aerosol,
      e        tau_aero, piz_aero, cg_aero,
@@ -3597 +3602 @@
      $                   prfl(:,1:klev),psfl(:,1:klev),
      $                   pmflxr(:,1:klev),pmflxs(:,1:klev),
-     $                   mr_ozone,cldtaurad, cldemirad)
+     $                   mr_ozone,cldtau, cldemi)
 
 !     L          calipso2D,calipso3D,cfadlidar,parasolrefl,atb,betamol,