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

Timestamp:

Jan 28, 2019, 7:31:11 PM (5 years ago)

Author:

oboucher

Message:

Adding some diagnostics for type_trac=co2i

File:

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

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -196 +196 @@
 ! Tropopause
          o_p_tropopause, o_z_tropopause, o_t_tropopause,  &
-         o_col_O3_strato, o_col_O3_tropo               ! Added ThL
+         o_col_O3_strato, o_col_O3_tropo,                 & 
+!--interactive CO2
+         o_flx_co2_ocean, o_flx_co2_land, o_flx_co2_ff, o_flx_co2_bb
 
 
@@ -327 +329 @@
     USE CHEM_REP, ONLY : nas, nbnas, tnamenas, ttextnas
 #endif
+
+    USE carbon_cycle_mod, ONLY: fco2_ff, fco2_bb, fco2_land, fco2_ocean
 
     USE phys_output_var_mod, ONLY: vars_defined, snow_o, zfra_o, bils_diss, &
@@ -1253 +1257 @@
           IF (slab_gm) THEN
              CALL histwrite_phy(o_slab_gm, dt_gm(:,1:nslay))
-          END IF
+          ENDIF
           IF (slab_hdiff) THEN
             IF (nslay.EQ.1) THEN
@@ -1663 +1667 @@
        ENDIF
        CALL histwrite_phy(o_dtcon, zx_tmp_fi3d)
-       if(iflag_thermals.eq.0)then
+       IF (iflag_thermals.EQ.0) THEN
           IF (vars_defined) THEN
              zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys + &
@@ -1669 +1673 @@
           ENDIF
           CALL histwrite_phy(o_tntc, zx_tmp_fi3d)
-       else if(iflag_thermals.ge.1.and.iflag_wake.EQ.1)then
+       ELSE IF(iflag_thermals.GE.1.AND.iflag_wake.EQ.1) THEN
           IF (vars_defined) THEN
              zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys + &
@@ -1676 +1680 @@
           ENDIF
           CALL histwrite_phy(o_tntc, zx_tmp_fi3d)
-       endif
+       ENDIF
        IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys
        CALL histwrite_phy(o_ducon, zx_tmp_fi3d)
@@ -1727 +1731 @@
           CALL histwrite_phy(o_plulst, plul_st)
           IF (vars_defined) THEN
-             do i=1,klon
+             DO i=1,klon
                 zx_tmp_fi2d(1:klon)=lmax_th(:)
-             enddo
+             ENDDO
           ENDIF
           CALL histwrite_phy(o_lmaxth, zx_tmp_fi2d)
@@ -2192 +2196 @@
             IF (type_trac == 'lmdz' .OR. type_trac == 'repr' .OR. type_trac == 'coag') THEN
              !--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))
-!             CALL histwrite_phy(o_dtr_lessi_impa(iq-nqo),d_tr_lessi_impa(:,:,iq-nqo))
-!             CALL histwrite_phy(o_dtr_lessi_nucl(iq-nqo),d_tr_lessi_nucl(:,:,iq-nqo))
-!             CALL histwrite_phy(o_dtr_insc(iq-nqo),d_tr_insc(:,:,iq-nqo))
-!             CALL histwrite_phy(o_dtr_bcscav(iq-nqo),d_tr_bcscav(:,:,iq-nqo))
-!             CALL histwrite_phy(o_dtr_evapls(iq-nqo),d_tr_evapls(:,:,iq-nqo))
-!             CALL histwrite_phy(o_dtr_ls(iq-nqo),d_tr_ls(:,:,iq-nqo))
-!             CALL histwrite_phy(o_dtr_trsp(iq-nqo),d_tr_trsp(:,:,iq-nqo))
-!             CALL histwrite_phy(o_dtr_sscav(iq-nqo),d_tr_sscav(:,:,iq-nqo))
-!             CALL histwrite_phy(o_dtr_sat(iq-nqo),d_tr_sat(:,:,iq-nqo))
-!             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_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))
+             CALL histwrite_phy(o_dtr_lessi_impa(iq-nqo),d_tr_lessi_impa(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_lessi_nucl(iq-nqo),d_tr_lessi_nucl(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_insc(iq-nqo),d_tr_insc(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_bcscav(iq-nqo),d_tr_bcscav(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_evapls(iq-nqo),d_tr_evapls(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_ls(iq-nqo),d_tr_ls(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_trsp(iq-nqo),d_tr_trsp(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_sscav(iq-nqo),d_tr_sscav(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_sat(iq-nqo),d_tr_sat(:,:,iq-nqo))
+             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))
              zx_tmp_fi2d=0.
              IF (vars_defined) THEN
@@ -2217 +2221 @@
              CALL histwrite_phy(o_trac_cum(iq-nqo), zx_tmp_fi2d)
 #endif
-            ENDIF
-          ENDDO
-       ENDIF
+            ENDIF   !--type_trac
+!
+            IF (type_trac == 'co2i') THEN
+             !--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)
+             !--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_ff,    fco2_ff)
+             CALL histwrite_phy(o_flx_co2_bb,    fco2_bb)
+            ENDIF !--type_trac co2i
+
+          ENDDO !--iq
+       ENDIF !--nqtot > nqo+1
 
        IF (type_trac == 'repr') THEN
@@ -2230 +2257 @@
 
        ENDIF   !(iflag_phytrac==1)
-
 
        IF (.NOT.vars_defined) THEN