Changeset 3411 for LMDZ6/branches/DYNAMICO-conv/libf/phylmd/cosp/cosp_output_mod.F90

Timestamp:

Nov 5, 2018, 3:24:59 PM (6 years ago)

Author:

Laurent Fairhead

Message:

Undoing merge with trunk (r3356) to properly register Yann's latest modifications

Location:

LMDZ6/branches/DYNAMICO-conv

Files:

• . (modified) (1 prop)
• libf/phylmd/cosp/cosp_output_mod.F90 (modified) (5 diffs)

LMDZ6/branches/DYNAMICO-conv/libf/phylmd/cosp/cosp_output_mod.F90

@@ -12 +12 @@
                             numMODISReffLiqBins, reffLIQ_binCenters
 
-     IMPLICIT NONE
 ! cosp_output_mod
       INTEGER, PRIVATE             :: i
@@ -255 +254 @@
   real                      :: zjulian,zjulian_start
   real,dimension(Ncolumns)  :: column_ax
+  real,dimension(2,SR_BINS) :: sratio_bounds
+  real,dimension(SR_BINS)   ::  sratio_ax
   real,dimension(DBZE_BINS) ::  dbze_ax
   CHARACTER(LEN=20), DIMENSION(3)  :: chfreq = (/ '1day', '1d  ', '3h  ' /)            
-  real,parameter,dimension(SR_BINS) :: sratio_ax = (/0.005, &
-                                                  0.605,2.09,4.,6., & 
-                                          8.5,12.5,17.5,22.5,27.5,35.,45.,55.,70.,50040./)
 
 !!! Variables d'entree
@@ -284 +282 @@
     enddo
  
+
+    sratio_bounds(2,:)=stlidar%srbval(:) ! srbval contains the upper
+!                                         limits from lmd_ipsl_stats.f90
+    sratio_bounds(1,2:SR_BINS) = stlidar%srbval(1:SR_BINS-1)
+    sratio_bounds(1,1)         = 0.0
+    sratio_bounds(2,SR_BINS)   = 1.e5 ! This matches with Chepfer et al., JGR,
+!                                    ! 2009. However, it is not consistent 
+                                     ! with the upper limit in
+                                     ! lmd_ipsl_stats.f90, which is
+                                     ! LIDAR_UNDEF-1=998.999
+     sratio_ax(:) = (sratio_bounds(1,:)+sratio_bounds(2,:))/2.0 
+
     cosp_outfilenames(1) = 'histmthCOSP'
     cosp_outfilenames(2) = 'histdayCOSP'
@@ -342 +352 @@
 ! AI nov 2015
    CALL wxios_add_vaxis("temp", LIDAR_NTEMP, LIDAR_PHASE_TEMP)
-   CALL wxios_add_vaxis("cth", MISR_N_CTH, MISR_CTH)
+   CALL wxios_add_vaxis("cth16", MISR_N_CTH, MISR_CTH)
    CALL wxios_add_vaxis("dbze", DBZE_BINS, dbze_ax) 
    CALL wxios_add_vaxis("scatratio", SR_BINS, sratio_ax)
@@ -390 +400 @@
       CALL histvert(cosp_nidfiles(iff),"temp","temperature","C",LIDAR_NTEMP,LIDAR_PHASE_TEMP,nverttemp(iff))
 
-      CALL histvert(cosp_nidfiles(iff),"cth","altitude","m",MISR_N_CTH,MISR_CTH,nvertmisr(iff))
+      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, &