Changeset 1924

Timestamp:

Jan 8, 2014, 3:41:30 PM (11 years ago)

Author:

idelkadi

Message:

• Rajout de nouveaux diagnostiques pour la comparaison entre les fractions nuageuses sorties derictement dans LMDZ et celles sorties via le simulateur ISCCP. Routines concernees : physiq.F90, phys_output_ctrlout_mod.F90, phys_local_var_mod.F90 et phys_output_write_mod.F90

• Rajout de la fraction d'ensoleillement par jour aux arguments de l'interface avec Cosp (entree pour le simulateur ISCCP). Routine concernée : physiq.F90

Location:

LMDZ5/trunk/libf/phylmd

Files:

• phys_local_var_mod.F90 (modified) (5 diffs)
• phys_output_ctrlout_mod.F90 (modified) (2 diffs)
• phys_output_write_mod.F90 (modified) (6 diffs)
• physiq.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -200 +200 @@
       REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: cldh, cldl, cldm, cldq, cldt, qsat2m, qsol
 !$OMP THREADPRIVATE(cldh, cldl, cldm, cldq, cldt, qsat2m, qsol)
+      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: cldhjn, cldljn, cldmjn,cldtjn
+!$OMP THREADPRIVATE(cldhjn, cldljn, cldmjn, cldtjn)
+      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: JrNt
+!$OMP THREADPRIVATE(JrNt)
       REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: dthmin, evap, fder, plcl, plfc, prw
 !$OMP THREADPRIVATE(dthmin, evap, fder, plcl, plfc, prw)
@@ -263 +267 @@
       REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: beta_prec
 !$OMP THREADPRIVATE(beta_prec)
-      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: rneb
-!$OMP THREADPRIVATE(rneb)
+      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: rneb,rnebjn
+!$OMP THREADPRIVATE(rneb,rnebjn)
 
 ! variables de sorties MM
@@ -373 +377 @@
       ALLOCATE(cdragm(klon), cdragh(klon), cldh(klon), cldl(klon))
       ALLOCATE(cldm(klon), cldq(klon), cldt(klon), qsat2m(klon), qsol(klon))
+      ALLOCATE(cldhjn(klon), cldljn(klon), cldmjn(klon), cldtjn(klon))
+      ALLOCATE(JrNt(klon))
       ALLOCATE(dthmin(klon), evap(klon), fder(klon), plcl(klon), plfc(klon))
       ALLOCATE(prw(klon), zustar(klon), zu10m(klon), zv10m(klon), rh2m(klon))
@@ -420 +426 @@
       ALLOCATE(wdtrainA(klon,klev),wdtrainM(klon,klev))
       ALLOCATE(beta_prec(klon,klev))
-      ALLOCATE(rneb(klon,klev))
+      ALLOCATE(rneb(klon,klev),rnebjn(klon,klev))
 
 
@@ -524 +530 @@
       DEALLOCATE(cdragm, cdragh, cldh, cldl)
       DEALLOCATE(cldm, cldq, cldt, qsat2m, qsol)
+      DEALLOCATE(cldljn, cldmjn, cldhjn, cldtjn, JrNt)
       DEALLOCATE(dthmin, evap, fder, plcl, plfc)
       DEALLOCATE(prw, zustar, zu10m, zv10m, rh2m, s_lcl)

LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -391 +391 @@
   TYPE(ctrl_out), SAVE :: o_cldt = ctrl_out((/ 1, 1, 2, 10, 5, 10, 11, 11, 11 /), &
     'cldt', 'Total cloudiness', '-', (/ ('', i=1, 9) /))
+  TYPE(ctrl_out), SAVE :: o_JrNt = ctrl_out((/ 1, 1, 10, 7, 10, 10, 11, 11, 11 /), &
+    'JrNt', '1 if Day 0 if Night', '-', (/ ('', i=1, 9) /))                                                                                  
+  TYPE(ctrl_out), SAVE :: o_cldhjn = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11 /), &
+    'cldhjn', 'High-level cloudiness Day', '-', (/ ('', i=1, 9) /))                                                                                   
+  TYPE(ctrl_out), SAVE :: o_cldmjn = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11/), &      
+    'cldmjn', 'Mid-level cloudiness day', '-', (/ ('', i=1, 9) /))                                                           
+  TYPE(ctrl_out), SAVE :: o_cldljn = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11/), &       
+    'cldljn', 'Low-level cloudiness day', '-', (/ ('', i=1, 9) /))
+  TYPE(ctrl_out), SAVE :: o_cldtjn = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11/), &     
+    'cldtjn', 'Total cloudiness day', '-', (/ ('', i=1, 9) /))
+                                                       
   TYPE(ctrl_out), SAVE :: o_cldq = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
     'cldq', 'Cloud liquid water path', 'kg/m2', (/ ('', i=1, 9) /))
@@ -852 +863 @@
   TYPE(ctrl_out), SAVE :: o_rneb = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11 /), &
     'rneb', 'Cloud fraction', '-', (/ ('', i=1, 9) /))
+  TYPE(ctrl_out), SAVE :: o_rnebjn = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11,11 /), &      
+    'rnebjn', 'Cloud fraction in day', '-', (/ ('', i=1, 9) /))
   TYPE(ctrl_out), SAVE :: o_rnebcon = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11 /), &
     'rnebcon', 'Convective Cloud Fraction', '-', (/ ('', i=1, 9) /))

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -58 +58 @@
          o_tke_srf, o_tke_max_srf, o_wstar, &
          o_cdrm, o_cdrh, o_cldl, o_cldm, o_cldh, &
-         o_cldt, o_cldq, o_lwp, o_iwp, o_ue, &
-         o_ve, o_uq, o_vq, o_cape, o_pbase, &
+         o_cldt, o_JrNt, o_cldljn, o_cldmjn, &
+         o_cldhjn, o_cldtjn, o_cldq, o_lwp, o_iwp, &
+         o_ue, o_ve, o_uq, o_vq, o_cape, o_pbase, &
          o_ptop, o_fbase, o_plcl, o_plfc, &
          o_wbeff, o_cape_max, o_upwd, o_Ma, &
@@ -106 +107 @@
          o_ovapinit, o_ovap, o_oliq, o_geop, &
          o_vitu, o_vitv, o_vitw, o_pres, o_paprs, &
-         o_zfull, o_zhalf, o_rneb, o_rnebcon, &
+         o_zfull, o_zhalf, o_rneb, o_rnebjn, o_rnebcon, &
          o_rnebls, o_rhum, o_ozone, o_ozone_light, &
          o_dtphy, o_dqphy, o_albe_srf, o_rugs_srf, &
@@ -176 +177 @@
          t2m, fevap, fluxt, fluxlat, fsollw, fsolsw, &
          wfbils, wfbilo, cdragm, cdragh, cldl, cldm, &
-         cldh, cldt, cldq, flwp, fiwp, ue, ve, uq, vq, &
+         cldh, cldt, JrNt, cldljn, cldmjn, cldhjn, &
+         cldtjn, cldq, flwp, fiwp, ue, ve, uq, vq, &
          plcl, plfc, wbeff, upwd, dnwd, dnwd0, prw, &
          s_pblh, s_pblt, s_lcl, s_therm, uwriteSTD, &
@@ -201 +203 @@
          ec550aer, flwc, fiwc, t_seri, theta, q_seri, &
          ql_seri, zphi, u_seri, v_seri, omega, cldfra, &
-         rneb, zx_rh, frugs, agesno, d_t_dyn, d_q_dyn, &
+         rneb, rnebjn, zx_rh, frugs, agesno, d_t_dyn, d_q_dyn, &
          d_u_dyn, d_v_dyn, d_t_con, d_t_ajsb, d_t_ajs, &
          d_u_ajsb, d_u_ajs, d_v_ajsb, d_v_ajs, &
@@ -551 +553 @@
        CALL histwrite_phy(o_cldh, cldh)
        CALL histwrite_phy(o_cldt, cldt)
+       CALL histwrite_phy(o_JrNt, JrNt)
+       CALL histwrite_phy(o_cldljn, cldl*JrNt)
+       CALL histwrite_phy(o_cldmjn, cldm*JrNt)
+       CALL histwrite_phy(o_cldhjn, cldh*JrNt)
+       CALL histwrite_phy(o_cldtjn, cldt*JrNt)
        CALL histwrite_phy(o_cldq, cldq)
        IF (vars_defined)       zx_tmp_fi2d(1:klon) = flwp(1:klon)
@@ -885 +892 @@
        CALL histwrite_phy(o_rnebcon, rnebcon)
        CALL histwrite_phy(o_rnebls, rneb)
+IF (vars_defined)  THEN
+      DO k=1, klev
+         DO i=1, klon
+            zx_tmp_fi3d(i,k)=cldfra(i,k)*JrNt(i)
+         ENDDO
+      ENDDO
+ENDIF
+       CALL histwrite_phy(o_rnebjn, zx_tmp_fi3d)
        CALL histwrite_phy(o_rhum, zx_rh)
        CALL histwrite_phy(o_ozone, &

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -1692 +1692 @@
   ENDIF
 
+! AI Janv 2014 
+        do i = 1, klon
+         if (fract(i).le.0.) then
+            JrNt(i)=0.
+         else
+            JrNt(i)=1.
+         endif
+        enddo
+
   if (mydebug) then
      call writefield_phy('u_seri',u_seri,llm)
@@ -3311 +3320 @@
              ecrit_mth,ecrit_day,ecrit_hf, &
              klon,klev,rlon,rlat,presnivs,overlap, &
-             ref_liq,ref_ice, &
+             fract,ref_liq,ref_ice, &
              pctsrf(:,is_ter)+pctsrf(:,is_lic), &
              zu10m,zv10m,pphis, &