Changeset 2993 for LMDZ5/trunk


Ignore:
Timestamp:
Sep 25, 2017, 11:41:28 AM (7 years ago)
Author:
oboucher
Message:

Adding tropospheric and stratospheric

Location:
LMDZ5/trunk
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/DefLists/field_def_lmdz.xml

    r2992 r2993  
    464464        <field id="z_tropopause"  long_name="Tropopause height" unit="m"      detect_missing_value=".true."  />
    465465        <field id="t_tropopause"  long_name="Tropopause temperature" unit="K" detect_missing_value=".true."  />
     466        <field id="colO3_strat"    long_name="Ozone stratospheric column"     unit="DU" />             <!-- Added ThL -->
     467        <field id="colO3_trop"     long_name="Ozone tropospheric column"      unit="DU" />             <!-- Added ThL -->
    466468    </field_group>
    467469
  • LMDZ5/trunk/DefLists/file_def_histday_lmdz.xml

    r2958 r2993  
    238238                <field field_ref="alp_bl_conv" level="1" />
    239239                <field field_ref="alp_bl_stat" level="1" />
     240                <field field_ref="colO3_strat" level="6" />     <!-- Added ThL -->
     241                <field field_ref="colO3_trop"  level="6" />     <!-- Added ThL -->
    240242
    241243                <field_group operation="average" freq_op="1ts" detect_missing_value=".true.">
  • LMDZ5/trunk/DefLists/file_def_histmth_lmdz.xml

    r2977 r2993  
    258258                <field field_ref="alp_bl_conv" level="1" />
    259259                <field field_ref="alp_bl_stat" level="1" />
     260                <field field_ref="colO3_strat" level="2" />     <!-- Added ThL -->
     261                <field field_ref="colO3_trop"  level="2" />     <!-- Added ThL -->
    260262
    261263                <field_group operation="average" freq_op="1ts" detect_missing_value=".true.">
  • LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

    r2992 r2993  
    12341234  TYPE(ctrl_out), SAVE :: o_t_tropopause = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
    12351235    't_tropopause', 'Tropopause temperature', 'K', (/ ('', i=1, 10) /))
     1236!--Added ThL
     1237  TYPE(ctrl_out), SAVE :: o_col_O3_strato = ctrl_out((/2, 6, 10, 10, 10, 10, 11, 11, 11, 11/), &
     1238    'colO3_strat','Ozone stratospheric column', 'DU', (/('', i=1, 10) /))
     1239  TYPE(ctrl_out), SAVE :: o_col_O3_tropo = ctrl_out((/2, 6, 10, 10, 10, 10, 11, 11, 11, 11/), &
     1240    'colO3_trop','Ozone tropospheric column', 'DU', (/('', i=1, 10) /))
     1241!--end add ThL
    12361242
    12371243#ifdef CPP_StratAer
  • LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

    r2992 r2993  
    186186         o_alt_tropo, &
    187187! Tropopause
    188          o_p_tropopause, o_z_tropopause, o_t_tropopause
     188         o_p_tropopause, o_z_tropopause, o_t_tropopause,  &
     189         o_col_O3_strato, o_col_O3_tropo               ! Added ThL
     190
    189191
    190192#ifdef CPP_StratAer
     
    12281230       CALL histwrite_phy(o_z_tropopause, z_tropopause)
    12291231
     1232! ThL -- In the following, we assume read_climoz == 1
     1233       zx_tmp_fi2d = 0.0    ! Computation for strato, added ThL
     1234       DO k=1, klev
     1235          zx_tmp_fi2d(:) = zx_tmp_fi2d(:) + wo(:,k,1) * stratomask(:,k) * 1.e3
     1236       END DO
     1237       CALL histwrite_phy(o_col_O3_strato, zx_tmp_fi2d) ! Added ThL
     1238       zx_tmp_fi2d = 0.0    ! Computation for tropo, added ThL
     1239       DO k=1, klev
     1240          zx_tmp_fi2d(:) = zx_tmp_fi2d(:) + wo(:,k,1) * (1.0-stratomask(:,k)) * 1.e3
     1241       END DO
     1242       CALL histwrite_phy(o_col_O3_tropo, zx_tmp_fi2d)   ! Added ThL
     1243! end add ThL
     1244
    12301245#ifdef CPP_StratAer
    12311246       IF (type_trac=='coag') THEN
Note: See TracChangeset for help on using the changeset viewer.