Changeset 3247 for LMDZ6/trunk/libf/phylmd/cosp/cosp_output_mod.F90

Timestamp:

Mar 12, 2018, 9:52:27 AM (7 years ago)

Author:

idelkadi

Message:

Recriture de la routine de gestion des cles de sorties Cosp :

ecriture d'un module cosp_read_otputkeys.F90 incluant 3 routines :

• une subroutine pour initialiser au 1er appel a Cosp les Cles
• une subroutine appelee au 2e passage dans Cosp, permettant de piloter les cles de sorties en fonction de ce qui est demande dans les fichiers xml (dans le cas ou les sorties sont geres par Xios)
• une routine appelee au 2e passage dans Cosp, permettant de lire les cles dans le ficher cosp_output_nl.txt dans le cas ou les sorties ne sont pas geres par Xios

Nettoyage dans le module cosp_output_write.F90

File:

• LMDZ6/trunk/libf/phylmd/cosp/cosp_output_mod.F90 (modified) (2 diffs)

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

@@ -255 +255 @@
   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
@@ -283 +284 @@
     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'