Changeset 3356 for LMDZ6/branches/DYNAMICO-conv/libf/phylmd/phys_output_mod.F90

Timestamp:

Jun 29, 2018, 12:31:11 PM (6 years ago)

Author:

Laurent Fairhead

Message:

First attempt at merging with trunk

Location:

LMDZ6/branches/DYNAMICO-conv

Files:

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

LMDZ6/branches/DYNAMICO-conv/libf/phylmd/phys_output_mod.F90

@@ -35 +35 @@
     USE iophy 
     USE dimphy
-    USE infotrac_phy, ONLY: nqtot, nqo, niadv, tname, ttext
+    USE infotrac_phy, ONLY: nqtot, nqo, niadv, tname, ttext, type_trac
     USE ioipsl
     USE phys_cal_mod, only : hour, calend
@@ -46 +46 @@
     USE vertical_layers_mod, ONLY: ap,bp,preff,presnivs, aps, bps, pseudoalt
     USE time_phylmdz_mod, ONLY: day_ini, itau_phy, start_time, annee_ref, day_ref
+#ifdef REPROBUS
+    USE chem_rep, ONLY: nbnas, tnamenas, ttextnas
+#endif
 #ifdef CPP_XIOS
     ! ug Pour les sorties XIOS
@@ -137 +140 @@
     REAL :: ini_heure
 #endif
+    INTEGER                         :: ISW
+    REAL, DIMENSION(NSW)            :: wl1_sun, wl2_sun !wavelength bounds (in um) for SW
+    REAL, DIMENSION(NSW)            :: wn1_sun, wn2_sun !wavenumber bounds (in m-1) for SW 
+    REAL, DIMENSION(NSW)            :: spectband  !mean wavenumb. of each sp.band
+    REAL, DIMENSION(NSW,2)          :: spbnds_sun !bounds of spectband
 
     WRITE(lunout,*) 'Debut phys_output_mod.F90'
@@ -151 +159 @@
     IF (.NOT. ALLOCATED(o_trac)) ALLOCATE(o_trac(nqtot))
     IF (.NOT. ALLOCATED(o_trac_cum)) ALLOCATE(o_trac_cum(nqtot))
+#ifdef REPROBUS
+    IF (.NOT. ALLOCATED(o_nas)) ALLOCATE(o_nas(nbnas))
+#endif
     ALLOCATE(o_dtr_the(nqtot),o_dtr_con(nqtot),o_dtr_lessi_impa(nqtot))
     ALLOCATE(o_dtr_lessi_nucl(nqtot),o_dtr_insc(nqtot),o_dtr_bcscav(nqtot))
@@ -354 +365 @@
      CALL wxios_add_vaxis("Alt", &
             levmax(iff) - levmin(iff) + 1, pseudoalt)
+
+    IF (NSW.EQ.6) THEN
+!
+!wl1_sun: minimum bound of wavelength (in um)
+!
+      wl1_sun(1)=0.180
+      wl1_sun(2)=0.250
+      wl1_sun(3)=0.440
+      wl1_sun(4)=0.690
+      wl1_sun(5)=1.190
+      wl1_sun(6)=2.380
+!
+!wl2_sun: maximum bound of wavelength (in um)
+!
+      wl2_sun(1)=0.250
+      wl2_sun(2)=0.440
+      wl2_sun(3)=0.690
+      wl2_sun(4)=1.190
+      wl2_sun(5)=2.380
+      wl2_sun(6)=4.000
+!
+    ELSE IF(NSW.EQ.2) THEN
+!
+!wl1_sun: minimum bound of wavelength (in um)
+!
+      wl1_sun(1)=0.250
+      wl1_sun(2)=0.690
+!
+!wl2_sun: maximum bound of wavelength (in um)
+!
+      wl2_sun(1)=0.690
+      wl2_sun(2)=4.000
+    ENDIF
+
+    DO ISW=1, NSW
+     wn1_sun(ISW)=1.e+6/wl1_sun(ISW) 
+     wn2_sun(ISW)=1.e+6/wl2_sun(ISW) 
+     spbnds_sun(ISW,1)=wn2_sun(ISW)
+     spbnds_sun(ISW,2)=wn1_sun(ISW)
+     spectband(ISW)=(wn1_sun(ISW)+wn2_sun(ISW))/2
+    ENDDO
+!
+!!! ajout axe vertical spectband : solar band number
+    CALL wxios_add_vaxis("spectband", NSW, spectband, positif='down')
    ELSE
     ! NMC files
@@ -523 +578 @@
             ENDDO
       ENDIF
+      IF (type_trac=='repr') THEN
+#ifdef REPROBUS
+         DO iiq=1,nbnas
+            o_nas(iiq) = ctrl_out((/ 4, 5, 5, 5, 10, 10, 11, 11, 11, 11 /), &
+                 tnamenas(iiq),ttextnas(iiq), "-", &
+                 (/ '', '', '', '', '', '', '', '', '', '' /))
+         ENDDO
+#endif
+      ENDIF
 
    ENDDO !  iff
@@ -536 +600 @@
 
     IF (prt_level >= 10) THEN
+      WRITE(lunout,*)'swaerofree_diag=',swaerofree_diag
       WRITE(lunout,*)'swaero_diag=',swaero_diag
       WRITE(lunout,*)'dryaod_diag=',dryaod_diag