Changeset 2993 for LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

Timestamp:

Sep 25, 2017, 11:41:28 AM (7 years ago)

Author:

oboucher

Message:

Adding tropospheric and stratospheric

File:

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

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -186 +186 @@
          o_alt_tropo, & 
 ! Tropopause
-         o_p_tropopause, o_z_tropopause, o_t_tropopause
+         o_p_tropopause, o_z_tropopause, o_t_tropopause,  &
+         o_col_O3_strato, o_col_O3_tropo               ! Added ThL
+
 
 #ifdef CPP_StratAer
@@ -1228 +1230 @@
        CALL histwrite_phy(o_z_tropopause, z_tropopause)
 
+! ThL -- In the following, we assume read_climoz == 1
+       zx_tmp_fi2d = 0.0    ! Computation for strato, added ThL
+       DO k=1, klev
+          zx_tmp_fi2d(:) = zx_tmp_fi2d(:) + wo(:,k,1) * stratomask(:,k) * 1.e3
+       END DO
+       CALL histwrite_phy(o_col_O3_strato, zx_tmp_fi2d) ! Added ThL
+       zx_tmp_fi2d = 0.0    ! Computation for tropo, added ThL
+       DO k=1, klev
+          zx_tmp_fi2d(:) = zx_tmp_fi2d(:) + wo(:,k,1) * (1.0-stratomask(:,k)) * 1.e3
+       END DO
+       CALL histwrite_phy(o_col_O3_tropo, zx_tmp_fi2d)   ! Added ThL
+! end add ThL
+
 #ifdef CPP_StratAer
        IF (type_trac=='coag') THEN