Changeset 5791 for LMDZ6/branches/contrails/libf/phylmd/phys_output_mod.F90

Timestamp:

Jul 28, 2025, 7:23:15 PM (6 days ago)

Author:

aborella

Message:

Merge with trunk r5789

Location:

LMDZ6/branches/contrails

Files:

• . (modified) (1 prop)
• libf/phylmd/phys_output_mod.F90 (modified) (13 diffs)

LMDZ6/branches/contrails/libf/phylmd/phys_output_mod.F90

@@ -6 +6 @@
   USE phys_output_var_mod
   USE phys_output_write_mod, ONLY : phys_output_write
+
   REAL, DIMENSION(nfiles),SAVE :: ecrit_files
+  LOGICAL, DIMENSION(nfiles)   :: phys_out_filekeys
 
 
@@ -51 +53 @@
     USE vertical_layers_mod, ONLY: ap,bp,preff,presnivs, aps, bps, pseudoalt, presinter
     USE time_phylmdz_mod, ONLY: day_ini, itau_phy, start_time, annee_ref, day_ref
+
+    USE AERO_MOD, ONLY : nbands_lw_rrtm !FC
+
     ! ug Pour les sorties XIOS
     use wxios_mod
@@ -123 +128 @@
     CHARACTER(LEN=20), DIMENSION(nfiles)  :: chtimestep = (/ 'Default', 'Default', 'Default', 'Default', 'Default', &
                                                              'Default', 'Default', 'Default', 'Default', 'Default' /)
-    LOGICAL, DIMENSION(nfiles)            :: phys_out_filekeys
     LOGICAL, DIMENSION(nfiles)            :: phys_out_filestations
 
@@ -155 +159 @@
     REAL, DIMENSION(NSW)            :: spectband  !mean wavenumb. of each sp.band
     REAL, DIMENSION(NSW,2)          :: spbnds_sun !bounds of spectband
+!FC
+    INTEGER                         :: ILW
+    REAL,  DIMENSION(nbands_lw_rrtm)  :: wl1_lw, wl2_lw
 
     WRITE(lunout,*) 'Debut phys_output_mod.F90'
-    ! Initialisations (Valeurs par defaut
+! Initialisations (Valeurs par defaut
+!FC
+       wl1_lw = [ 10., 250., 500., 630., 700., 820., 980.,1080. &
+     & ,1180.,1390.,1480.,1800.,2080.,2250.,2380.,2600.]
+       wl2_lw = [250., 500., 630., 700., 820., 980.,1080.,1180. &
+     & ,1390.,1480.,1800.,2080.,2250.,2380.,2600.,3000.]
+!         print*, 'avant boucle', nbands_lw_rrtm
+!        DO ILW=1,nbands_lw_rrtm
+!        spectbandLW(ILW)= (wl1_lw(ilw) + wl2_lw(ILW) )/2
+!        print*, 'on a les canaux ? ',ILW,spectbandLW(ILW),wl1_lw(ILW),wl2_lw(ilw)
+!        enddo
+!        print*, 'spectbandLW',spectbandLW
+
+!FC
+
 
     DO ilev=1,klev
@@ -229 +250 @@
     type_ecri(10)= 'ave(X)'
 
-    clef_files(1) = ok_mensuel
-    clef_files(2) = ok_journe
-    clef_files(3) = ok_hf
-    clef_files(4) = ok_instan
-    clef_files(5) = ok_LES
-    clef_files(6) = ok_instan
-    clef_files(7) = ok_histNMC(1)
-    clef_files(8) = ok_histNMC(2)
-    clef_files(9) = ok_histNMC(3)
+    clef_files(1:3) = .TRUE.
+    clef_files(4:10) = .FALSE.
     IF (CPPKEY_STRATAER) THEN
       clef_files(10)= .TRUE.
@@ -245 +259 @@
 
     !sortir des fichiers "stations" si clef_stations(:)=.TRUE.
-    clef_stations(1) = .FALSE.
-    clef_stations(2) = .FALSE.
-    clef_stations(3) = .FALSE.
-    clef_stations(4) = .FALSE.
-    clef_stations(5) = .FALSE.
-    clef_stations(6) = .FALSE.
-    clef_stations(7) = .FALSE.
-    clef_stations(8) = .FALSE.
-    clef_stations(9) = .FALSE.
-    clef_stations(10)= .FALSE.
-
-    lev_files(1) = lev_histmth
-    lev_files(2) = lev_histday
-    lev_files(3) = lev_histhf
-    lev_files(4) = lev_histins
-    lev_files(5) = lev_histLES
-    lev_files(6) = lev_histins
-    lev_files(7) = levout_histNMC(1)
-    lev_files(8) = levout_histNMC(2)
-    lev_files(9) = levout_histNMC(3)
-    lev_files(10)= 5
-
+    clef_stations(1:10) = .FALSE.
+    lev_files(1:10) = 5
+
+    print*,'A ecrit_mth=',ecrit_mth
+    ! Frequencies of the history files;
+    ! Defaut
+    ! overwritten by chtimestep given to convers_timesteps
     ecrit_files(1) = ecrit_mth
     ecrit_files(2) = ecrit_day
@@ -272 +271 @@
     ecrit_files(4) = ecrit_ins
     ecrit_files(5) = ecrit_LES
-    ecrit_files(6) = ecrit_ins
-    ecrit_files(7) = freq_outNMC(1)
-    ecrit_files(8) = freq_outNMC(2)
-    ecrit_files(9) = freq_outNMC(3)
-    ecrit_files(10)= ecrit_mth
+    ecrit_files(6:10) = ecrit_ins
 
     !! Lectures des parametres de sorties dans physiq.def
@@ -367 +362 @@
       ELSE
        CALL convers_timesteps(chtimestep(iff),dtime,ecrit_files(iff)) 
+      print*,'Dans phys_output_open, iff=',iff,' ecrit_files=',ecrit_files(iff)
       ENDIF
+      ! ecrit_files contains frequency of file iif in seconds 
+
 
        WRITE(lunout,*)'ecrit_files(',iff,')= ',ecrit_files(iff)
@@ -502 +500 @@
           CALL histvert(nid_files(iff), "Alt","Height approx for scale heigh of 8km at levels", "Km", &                       
                levmax(iff) - levmin(iff) + 1,pseudoalt,nvertAlt(iff))
+!FC
+!          CALL histvert(nid_files(iff), "spectbandLW"," LW bands ", "cm-1", &
+!                  nbands_lw_rrtm ,wl1_lw, ncanaux(iff))
+!          print *, ' apres ncanaux = ' , ncanaux(iff),iff
+!FC
+
 
           ELSE
@@ -624 +628 @@
    ! Updated write frequencies due to phys_out_filetimesteps. 
     ! Write frequencies are now in seconds.  
-    ecrit_mth = ecrit_files(1)
-    ecrit_day = ecrit_files(2)
-    ecrit_hf  = ecrit_files(3)
-    ecrit_ins = ecrit_files(4)
-    ecrit_LES = ecrit_files(5)
-    ecrit_ins = ecrit_files(6)
+! WHY CHANGING ecrit_mth ?
+! For Cosp ?
+!    ecrit_mth = ecrit_files(1)
+!    print*,'B ecrit_mth=',ecrit_mth
+!    ecrit_day = ecrit_files(2)
+!    ecrit_hf  = ecrit_files(3)
+!    ecrit_ins = ecrit_files(4)
+!    ecrit_LES = ecrit_files(5)
+!    ecrit_ins = ecrit_files(6)
 
     IF (prt_level >= 10) THEN
@@ -638 +645 @@
       WRITE(lunout,*)'phys_output_open: ends here'
     ENDIF
+    print*,'Dans phys_output_open,ecrit_files B',ecrit_files(1:6)
 
 !  DO iq=1,nqtot
@@ -679 +687 @@
     IF ( type == 'day'.or.type == 'days'.or.type == 'd'.or.type == 'jours'.or.type == 'jour'.or.type == 'j' )&
            &  timestep = ttt * dayseconde
-    IF ( type == 'months'.or.type == 'month'.or.type == 'mth'.or.type == 'mois'.or.type == 'm' ) THEN
+    IF ( type == 'months'.or.type == 'month'.or.type == 'mth'.or.type == 'mois'.or.type == 'm' .or.type == 'mo' ) THEN
        WRITE(lunout,*)'annee_ref,day_ref mon_len',annee_ref,day_ref,mth_len
        timestep = ttt * dayseconde * mth_len
@@ -687 +695 @@
     IF ( type == 'mn'.or.type == 'minutes'.or.type == 'minute'.or.type == 'm' ) timestep = ttt * 60.
     IF ( type == 's'.or.type == 'sec'.or.type == 'secondes'.or.type =='seconde'   ) timestep = ttt
-    IF ( type == 'TS' ) timestep = ttt * dtime
+    IF ( type == 'TS' .or. type == 'ts' ) timestep = ttt * dtime
 
     WRITE(lunout,*)'type =      ',type