Changeset 3865 for LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

Timestamp:

Mar 23, 2021, 4:14:07 PM (3 years ago)

Author:

lmdz-users

Message:

Modifications from Thibaut to create an ESM with interactive CO2 + INCA aerosols

File:

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

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -453 +453 @@
     REAL,DIMENSION(klon)      :: zrho, zt
 
+    INTEGER :: nqup
+
     ! On calcul le nouveau tau:
     itau_w = itau_phy + itap 
@@ -2436 +2438 @@
          ENDIF !--type_trac co2i
 
+         IF (type_trac == 'inco') THEN
+           nqup = nqo+1
+           DO iq=nqo+1, nqup
+             !--3D fields
+             CALL histwrite_phy(o_trac(iq-nqo), tr_seri(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_vdf(iq-nqo),d_tr_cl(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_the(iq-nqo),d_tr_th(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_con(iq-nqo),d_tr_cv(:,:,iq-nqo))
+             !--2D fields
+             !--CO2 burden
+             zx_tmp_fi2d=0.
+             IF (vars_defined) THEN
+                DO k=1,klev
+                   zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+zmasse(:,k)*tr_seri(:,k,iq-nqo)
+                ENDDO
+             ENDIF
+             CALL histwrite_phy(o_trac_cum(iq-nqo), zx_tmp_fi2d)
+           ENDDO !--iq
+           !--CO2 net fluxes
+           CALL histwrite_phy(o_flx_co2_land,  fco2_land)
+           CALL histwrite_phy(o_flx_co2_ocean, fco2_ocean)
+           CALL histwrite_phy(o_flx_co2_ocean_cor, fco2_ocean_cor)
+           CALL histwrite_phy(o_flx_co2_land_cor, fco2_land_cor)
+           CALL histwrite_phy(o_flx_co2_ff,    fco2_ff)
+           CALL histwrite_phy(o_flx_co2_bb,    fco2_bb)
+         ENDIF !--type_trac inco
+
        ENDIF   !(iflag_phytrac==1)