Changeset 3121

Timestamp:

Dec 13, 2017, 11:35:15 AM (7 years ago)

Author:

musat

Message:

Add LMDZ outputs: icc3dcon, icc3dstra corresponding to
CMIP6' clic, clis
IM

Location:

LMDZ6/trunk

Files:

• DefLists/CMIP6_ping_atmos.xml (modified) (1 diff)
• DefLists/field_def_lmdz.xml (modified) (1 diff)
• DefLists/file_def_histLES_lmdz.xml (modified) (1 diff)
• DefLists/file_def_histday_lmdz.xml (modified) (1 diff)
• DefLists/file_def_histins_lmdz.xml (modified) (1 diff)
• DefLists/file_def_histmth_lmdz.xml (modified) (1 diff)
• DefLists/file_def_histstn_lmdz.xml (modified) (1 diff)
• libf/phylmd/newmicro.F90 (modified) (3 diffs)
• libf/phylmd/phys_local_var_mod.F90 (modified) (3 diffs)
• libf/phylmd/phys_output_ctrlout_mod.F90 (modified) (2 diffs)
• libf/phylmd/phys_output_write_mod.F90 (modified) (4 diffs)

LMDZ6/trunk/DefLists/CMIP6_ping_atmos.xml

@@ -65 +65 @@
    <field id="CMIP6_clhcalipso"    field_ref="pclhcalipso"       /> <!-- P1 (%) cloud_area_fraction_in_atmosphere_layer : Percentage cloud cover in layer centred on 220hPa -->
    <field id="CMIP6_cli"           field_ref="iwcon"            /> <!-- P1 (kg kg-1) mass_fraction_of_cloud_ice_in_air : Includes both large-scale and convective cloud. This is calculated as the mass of cloud ice in the grid cell divided by the mass of air (including the water in all phases) in the grid cell. It includes precipitating hydrometeors ONLY if the precipitating hydrometeors affect the calculation of radiative transfer in model. -->
-   <field id="CMIP6_clic"          field_ref="dummy_XYA"        /> <!-- P2 (1.0) mass_fraction_of_convective_cloud_ice_in_air : Calculated as the mass of convective cloud ice  in the grid cell divided by the mass of air (including the water in all phases) in the grid cell.  This includes precipitating hydrometeors ONLY if the precipitating hydrometeors affect the calculation of radiative transfer in model. -->
+   <field id="CMIP6_clic"          field_ref="icc3dcon"        /> <!-- P2 (1.0) mass_fraction_of_convective_cloud_ice_in_air : Calculated as the mass of convective cloud ice  in the grid cell divided by the mass of air (including the water in all phases) in the grid cell.  This includes precipitating hydrometeors ONLY if the precipitating hydrometeors affect the calculation of radiative transfer in model. -->
    <field id="CMIP6_climodis"      field_ref="climodis"         /> <!-- P1 (%) ice_cloud_area_fraction : MODIS Ice Cloud Area Percentage -->
-   <field id="CMIP6_clis"          field_ref="dummy_XYA"        /> <!-- P2 (1.0) mass_fraction_of_stratiform_cloud_ice_in_air : Calculated as the mass of stratiform cloud ice  in the grid cell divided by the mass of air (including the water in all phases) in the grid cell.  This includes precipitating hydrometeors ONLY if the precipitating hydrometeors affect the calculation of radiative transfer in model. -->
+   <field id="CMIP6_clis"          field_ref="icc3dstra"        /> <!-- P2 (1.0) mass_fraction_of_stratiform_cloud_ice_in_air : Calculated as the mass of stratiform cloud ice  in the grid cell divided by the mass of air (including the water in all phases) in the grid cell.  This includes precipitating hydrometeors ONLY if the precipitating hydrometeors affect the calculation of radiative transfer in model. -->
    <field id="CMIP6_clisccp"       field_ref="clisccp2"         /> <!-- P1 (%) cloud_area_fraction_in_atmosphere_layer : Percentage cloud cover in optical depth categories. -->
    <field id="CMIP6_clivi"         field_ref="iwp"              /> <!-- P1 (kg m-2) atmosphere_cloud_ice_content : mass of ice water in the column divided by the area of the column (not just the area of the cloudy portion of the column). Includes precipitating frozen hydrometeors ONLY if the precipitating hydrometeor affects the calculation of radiative transfer in model. -->

LMDZ6/trunk/DefLists/field_def_lmdz.xml

@@ -568 +568 @@
         <field id="lcc3dcon"    long_name="Convective cloud liquid fraction"    unit="1" />
         <field id="lcc3dstra"    long_name="Stratiform cloud liquid fraction"    unit="1" />
+        <field id="icc3dcon"    long_name="Mass Fraction of Convective Cloud Ice"    unit="kg/kg" />
+        <field id="icc3dstra"    long_name="Mass Fraction of Stratiform Cloud Ice "    unit="kg/kg" />
         <field id="clwcon"    long_name="Convective Cloud Liquid water content"    unit="kg/kg" />
         <field id="Ma"    long_name="undilute adiab updraft"    unit="kg/m2/s" />

LMDZ6/trunk/DefLists/file_def_histLES_lmdz.xml

@@ -499 +499 @@
                 <field field_ref="lcc3dcon" level="10" />
                 <field field_ref="lcc3dstra" level="10" />
+                <field field_ref="icc3dcon" level="10" />
+                <field field_ref="icc3dstra" level="10" />
                 <field field_ref="clwcon" level="10" />
                 <field field_ref="Ma" level="10" />

LMDZ6/trunk/DefLists/file_def_histday_lmdz.xml

@@ -512 +512 @@
                 <field field_ref="lcc3dcon" level="6" />
                 <field field_ref="lcc3dstra" level="6" />
+                <field field_ref="icc3dcon" level="6" />
+                <field field_ref="icc3dstra" level="6" />
                 <field field_ref="clwcon" level="10" />
                 <field field_ref="Ma" level="10" />

LMDZ6/trunk/DefLists/file_def_histins_lmdz.xml

@@ -499 +499 @@
                 <field field_ref="lcc3dcon" level="10" />
                 <field field_ref="lcc3dstra" level="10" />
+                <field field_ref="icc3dcon" level="10" />
+                <field field_ref="icc3dstra" level="10" />
                 <field field_ref="clwcon" level="10" />
                 <field field_ref="Ma" level="10" />

LMDZ6/trunk/DefLists/file_def_histmth_lmdz.xml

@@ -554 +554 @@
                 <field field_ref="lcc3dcon" level="2" />
                 <field field_ref="lcc3dstra" level="2" />
+                <field field_ref="icc3dcon" level="2" />
+                <field field_ref="icc3dstra" level="2" />
                 <field field_ref="clwcon" level="4" />
                 <field field_ref="Ma" level="4" />

LMDZ6/trunk/DefLists/file_def_histstn_lmdz.xml

@@ -499 +499 @@
                 <field field_ref="lcc3dcon" level="10" />
                 <field field_ref="lcc3dstra" level="10" />
+                <field field_ref="icc3dcon" level="10" />
+                <field field_ref="icc3dstra" level="10" />
                 <field field_ref="clwcon" level="10" />
                 <field field_ref="Ma" level="10" />

LMDZ6/trunk/libf/phylmd/newmicro.F90

@@ -10 +10 @@
   USE dimphy
   USE phys_local_var_mod, ONLY: scdnc, cldncl, reffclwtop, lcc, reffclws, &
-    reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra
+    reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra, icc3dcon, icc3dstra
   USE phys_state_var_mod, ONLY: rnebcon, clwcon
   USE icefrac_lsc_mod ! computes ice fraction (JBM 3/14)
@@ -586 +586 @@
         lcc3dstra(i, k) = lcc3dstra(i, k) - lcc3dcon(i, k) ! eau liquide stratiforme
         lcc3dstra(i, k) = max(lcc3dstra(i,k), 0.0)
+        !FC pour la glace (CAUSES)
+        icc3dcon(i, k) = rnebcon(i, k)*(1-phase3d(i, k))*clwcon(i, k) !  glace convective
+        icc3dstra(i, k)= pclc(i, k)*pqlwp(i, k)*(1-phase3d(i, k))
+        icc3dstra(i, k) = icc3dstra(i, k) - icc3dcon(i, k) ! glace stratiforme
+        icc3dstra(i, k) = max( icc3dstra(i, k), 0.0)
+        !FC (CAUSES)
+
         ! Compute cloud droplet radius as above in meter
         radius = 1.1*((pqlwp(i,k)*pplay(i,k)/(rd*t(i,k)))/(4./3*rpi*1000.* &
@@ -626 +633 @@
         IF (lcc3dcon(i,k)<=0.0) lcc3dcon(i, k) = 0.0
         IF (lcc3dstra(i,k)<=0.0) lcc3dstra(i, k) = 0.0
+!FC (CAUSES)
+        IF (icc3dcon(i,k)<=0.0) icc3dcon(i, k) = 0.0
+        IF (icc3dstra(i,k)<=0.0) icc3dstra(i, k) = 0.0
+!FC (CAUSES)
       END DO
       IF (reffclwtop(i)<=0.0) reffclwtop(i) = 0.0

LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -169 +169 @@
       REAL, SAVE, ALLOCATABLE :: lcc3dstra(:,:)
       !$OMP THREADPRIVATE(lcc3dstra)
+      REAL, SAVE, ALLOCATABLE :: icc3dcon(:,:)
+      !$OMP THREADPRIVATE(icc3dcon)
+      REAL, SAVE, ALLOCATABLE :: icc3dstra(:,:)
+      !$OMP THREADPRIVATE(icc3dstra)
       REAL, SAVE, ALLOCATABLE :: od443aer(:) 
       !$OMP THREADPRIVATE(od443aer) 
@@ -593 +597 @@
       ALLOCATE(lcc3dcon(klon, klev))
       ALLOCATE(lcc3dstra(klon, klev))
+      ALLOCATE(icc3dcon(klon, klev))
+      ALLOCATE(icc3dstra(klon, klev))
       ALLOCATE(od443aer(klon))
       ALLOCATE(od550aer(klon))
@@ -868 +874 @@
       DEALLOCATE(lcc3dcon)
       DEALLOCATE(lcc3dstra)
+      DEALLOCATE(icc3dcon)
+      DEALLOCATE(icc3dstra)
       DEALLOCATE(od443aer)
       DEALLOCATE(od550aer)

LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -809 +809 @@
     's2', 'Surface moyenne des panaches de type 2', 'm2', (/ ('', i=1, 10) /))
               
-!-------Déclenchement stochastiqu
+!-------Déclenchement stochastique
   TYPE(ctrl_out), SAVE :: o_proba_notrig = ctrl_out((/ 1, 6, 6, 6, 10, 10, 11, 11, 11, 11/), &
     'proba_notrig', 'Probabilite de non-declenchement', ' ', (/ ('', i=1, 10) /))
@@ -1447 +1447 @@
   TYPE(ctrl_out), SAVE :: o_lcc3dstra = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lcc3dstra', 'Stratiform cloud liquid fraction', '1', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_icc3dcon = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'icc3dcon', 'Convective cloud ice fraction', '1', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_icc3dstra = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'icc3dstra', 'Stratiform cloud ice fraction', '1', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_solbnd = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'solbnd', 'Top-of-Atmosphere Solar Insolation for each band', 'W m-2', (/ ('', i=1, 10) /))

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -14 +14 @@
 CONTAINS 
 
-  ! ug Routine pour définir (los du premier passageà) ET sortir les variables
+  ! ug Routine pour définir (lors du premier passage) ET sortir les variables
   SUBROUTINE phys_output_write(itap, pdtphys, paprs, pphis, &
        pplay, lmax_th, aerosol_couple,         &
@@ -122 +122 @@
          o_cldncl, o_reffclws, o_reffclwc, o_solbnd, o_stratomask,&
          o_cldnvi, o_lcc, o_lcc3d, o_lcc3dcon, &
-         o_lcc3dstra, o_reffclwtop, o_ec550aer, &
+         o_lcc3dstra, o_icc3dcon, o_icc3dstra, &
+         o_reffclwtop, o_ec550aer, &
          o_lwcon, o_iwcon, o_temp, o_theta, &
          o_ovapinit, o_ovap, o_oliq, o_ocond, o_geop, &
@@ -278 +279 @@
          sollwad0_aero, toplwai_aero, sollwai_aero, &
          scdnc, cldncl, reffclws, reffclwc, cldnvi, stratomask,&
-         lcc, lcc3d, lcc3dcon, lcc3dstra, reffclwtop, &
+         lcc, lcc3d, lcc3dcon, lcc3dstra, &
+         icc3dcon, icc3dstra, reffclwtop, &
          ec550aer, flwc, fiwc, t_seri, theta, q_seri, &
          ql_seri, qs_seri, tr_seri, &
@@ -1406 +1408 @@
           CALL histwrite_phy(o_lcc3dcon, lcc3dcon)
           CALL histwrite_phy(o_lcc3dstra, lcc3dstra)
+          CALL histwrite_phy(o_icc3dcon, icc3dcon)
+          CALL histwrite_phy(o_icc3dstra, icc3dstra)
           CALL histwrite_phy(o_reffclwtop, reffclwtop)
        ENDIF