Changeset 3125 for LMDZ6/trunk/libf/phylmd/phys_output_mod.F90

Timestamp:

Dec 14, 2017, 8:55:21 AM (6 years ago)

Author:

acozic

Message:

Update some routines to coupled LMDZ6 with REPROBUS

File:

• LMDZ6/trunk/libf/phylmd/phys_output_mod.F90 (modified) (4 diffs)

LMDZ6/trunk/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
@@ -156 +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))
@@ -572 +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