Changeset 4800 for LMDZ6/branches/SETHET_DECOUPLE/libf/phylmd/phys_output_write_mod.F90

Timestamp:

Feb 2, 2024, 5:11:05 PM (4 months ago)

Author:

nfevrier

Message:

First save of N. Février's modifications

File:

• LMDZ6/branches/SETHET_DECOUPLE/libf/phylmd/phys_output_write_mod.F90 (modified) (4 diffs)

LMDZ6/branches/SETHET_DECOUPLE/libf/phylmd/phys_output_write_mod.F90

@@ -413 +413 @@
     ! Local
     INTEGER :: itau_w
-    INTEGER :: i, iinit, iinitend=1, iff, iq, iiq, nsrf, k, ll, naero
+    INTEGER :: i, iinit, iinitend=1, iff, iq, iiq, nsrf, k, ll, naero, j
     REAL, DIMENSION (klon) :: zx_tmp_fi2d, zpt_conv2d, wind100m
     REAL, DIMENSION (klon,klev) :: zx_tmp_fi3d, zpt_conv
@@ -446 +446 @@
     REAL,DIMENSION(klon,klev) :: z, dz
     REAL,DIMENSION(klon)      :: zrho, zt
+
+    INTEGER,DIMENSION(2) :: index_bc=[106,116]
 
     ! On calcul le nouveau tau:
@@ -2368 +2370 @@
              CALL histwrite_phy(o_dtr_uscav(iq-nqo),d_tr_uscav(:,:,iq-nqo))
             !--2D fields
-             CALL histwrite_phy(o_dtr_dry(iq-nqo), flux_tr_dry(:,iq-nqo))
+             !CALL histwrite_phy(o_dtr_dry(iq-nqo), flux_tr_dry(:,iq-nqo))
              zx_tmp_fi2d=0.
              IF (vars_defined) THEN
@@ -2378 +2380 @@
            ENDDO !--iq
          ENDIF   !--type_trac
+
+         ! Pour bilan
+         DO j=1,size(index_bc)
+             iq=index_bc(j)
+            ! WRITE(lunout,*) "WRITE_TEST_1"," / iq=",iq
+             CALL histwrite_phy(o_dtr_vdf(iq-nqo),d_tr_cl(:,:,iq-nqo))
+             !WRITE(lunout,*) "WRITE_TEST_2"
+         END DO
+         
 !
          IF (type_trac == 'co2i') THEN