Changeset 2713 for LMDZ5/trunk/libf/phylmd/cosp/cosp_output_mod.F90

Timestamp:

Nov 24, 2016, 4:02:04 PM (8 years ago)

Author:

idelkadi

Message:

Mise a jour du simulateur COSP.
Passage de la version 1.4.0 a la version 1.4.1
(Version suggeree pour CFMIP3/CMIP6)

File:

• LMDZ5/trunk/libf/phylmd/cosp/cosp_output_mod.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/cosp/cosp_output_mod.F90

@@ -9 +9 @@
   USE MOD_COSP_TYPES
   use MOD_COSP_Modis_Simulator, only : cosp_modis
+  use mod_modis_sim, only : numMODISReffIceBins, reffICE_binCenters, &
+                            numMODISReffLiqBins, reffLIQ_binCenters
 
 ! cosp_output_mod
@@ -17 +19 @@
 !$OMP THREADPRIVATE(cosp_outfilekeys, cosp_nidfiles)
       INTEGER, DIMENSION(3), SAVE  :: nhoricosp,nvert,nvertmcosp,nvertcol,nvertbze, &
-                                      nvertsratio,nvertisccp,nvertp,nverttemp,nvertmisr
+                                      nvertsratio,nvertisccp,nvertp,nverttemp,nvertmisr, &
+                                      nvertReffIce,nvertReffLiq
       REAL, DIMENSION(3), SAVE                :: zoutm_cosp
 !$OMP THREADPRIVATE(nhoricosp, nvert,nvertmcosp,nvertcol,nvertsratio,nvertbze,nvertisccp,nvertp,zoutm_cosp,nverttemp,nvertmisr)
+!$OMP THREADPRIVATE(nvertReffIce,nvertReffLiq)
       REAL, SAVE                   :: zdtimemoy_cosp
 !$OMP THREADPRIVATE(zdtimemoy_cosp) 
@@ -176 +180 @@
   TYPE(ctrl_outcosp), SAVE :: o_clmodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
          "clmodis", "MODIS Cloud Area Fraction", "%", (/ ('', i=1, 3) /))
+  TYPE(ctrl_outcosp), SAVE :: o_crimodis = ctrl_outcosp((/ .TRUE., .TRUE., .TRUE. /), &
+         "crimodis", "Optical_Thickness_vs_ReffIce from Modis", "%", (/ ('',i=1, 3) /))          
+  TYPE(ctrl_outcosp), SAVE :: o_crlmodis = ctrl_outcosp((/ .TRUE., .TRUE.,.TRUE. /), &
+         "crlmodis", "Optical_Thickness_vs_ReffLiq from Modis", "%", (/ ('',i=1, 3) /))         
 
 ! Rttovs simulator
@@ -325 +333 @@
 !   CALL wxios_add_vaxis("dbze", DBZE_BINS, dbze_ax) 
 !   CALL wxios_add_vaxis("scatratio", SR_BINS, sratio_ax)
+   CALL wxios_add_vaxis("ReffIce", numMODISReffIceBins, reffICE_binCenters)
+   CALL wxios_add_vaxis("ReffLiq", numMODISReffLiqBins, reffLIQ_binCenters)
+
 #endif
    
@@ -366 +377 @@
 
       CALL histvert(cosp_nidfiles(iff),"cth16","altitude","m",MISR_N_CTH,MISR_CTH,nvertmisr(iff))
+  
+      CALL histvert(cosp_nidfiles(iff),"ReffIce","Effective_particle_size_Ice","microns",numMODISReffIceBins, reffICE_binCenters, &
+                    nvertReffIce(iff))                                         
+     
+      CALL histvert(cosp_nidfiles(iff),"ReffLiq","Effective_particle_size_Liq","microns",numMODISReffLiqBins, reffLIQ_binCenters, &                                  
+                    nvertReffLiq(iff))
 
 !      CALL histvert(cosp_nidfiles(iff),"dbze","equivalent_reflectivity_factor","dBZ",DBZE_BINS,dbze_ax,nvertbze(iff))