Changeset 2958 for LMDZ5/trunk


Ignore:
Timestamp:
Jul 18, 2017, 2:21:01 PM (7 years ago)
Author:
jbmadeleine
Message:

Added new outputs to cloudth_vert (level 10):
s_th, Thermal plume saturation deficit, kg/kg
s_env, Environment saturation deficit, kg/kg
sigma_th, Thermal plume gauss variance, kg/kg
sigma_env, Environment gauss variance, kg/kg

Location:
LMDZ5/trunk
Files:
11 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/DefLists/field_def_lmdz.xml

    r2955 r2958  
    596596        <field id="q_th"    long_name="Thermal plume total humidity"    unit="kg/kg" />
    597597        <field id="a_th"    long_name="Thermal plume fraction"    unit="-" />
     598        <field id="s_th"    long_name="Thermal plume saturation deficit"    unit="kg/kg" />
     599        <field id="s_env"    long_name="Environment saturation deficit"    unit="kg/kg" />
     600        <field id="sigma_th"    long_name="Thermal plume gauss variance"    unit="kg/kg" />
     601        <field id="sigma_env"    long_name="Environment gauss variance"    unit="kg/kg" />
    598602        <field id="d_th"    long_name="Thermal plume detrainment"    unit="K/s" />
    599603        <field id="dqthe"    long_name="Thermal dQ"    unit="(kg/kg)/s" />

  • LMDZ5/trunk/DefLists/file_def_histLES_lmdz.xml

    r2945 r2958  
    548548                <field field_ref="q_th" level="10" />
    549549                <field field_ref="a_th" level="10" />
     550                <field field_ref="s_th" level="10" />
     551                <field field_ref="s_env" level="10" />
     552                <field field_ref="sigma_th" level="10" />
     553                <field field_ref="sigma_env" level="10" />
    550554                <field field_ref="d_th" level="10" />
    551555                <field field_ref="dqthe" level="10" />

  • LMDZ5/trunk/DefLists/file_def_histday_lmdz.xml

    r2952 r2958  
    559559                <field field_ref="q_th" level="10" />
    560560                <field field_ref="a_th" level="10" />
     561                <field field_ref="s_th" level="10" />
     562                <field field_ref="s_env" level="10" />
     563                <field field_ref="sigma_th" level="10" />
     564                <field field_ref="sigma_env" level="10" />
    561565                <field field_ref="d_th" level="10" />
    562566                <field field_ref="dqthe" level="10" />

  • LMDZ5/trunk/DefLists/file_def_histhf_lmdz.xml

    r2952 r2958  
    582582                <field field_ref="q_th" level="10" />
    583583                <field field_ref="a_th" level="10" />
     584                <field field_ref="s_th" level="10" />
     585                <field field_ref="s_env" level="10" />
     586                <field field_ref="sigma_th" level="10" />
     587                <field field_ref="sigma_env" level="10" />
    584588                <field field_ref="d_th" level="10" />
    585589                <field field_ref="dqthe" level="10" />

  • LMDZ5/trunk/DefLists/file_def_histins_lmdz.xml

    r2945 r2958  
    548548                <field field_ref="q_th" level="10" />
    549549                <field field_ref="a_th" level="10" />
     550                <field field_ref="s_th" level="10" />
     551                <field field_ref="s_env" level="10" />
     552                <field field_ref="sigma_th" level="10" />
     553                <field field_ref="sigma_env" level="10" />
    550554                <field field_ref="d_th" level="10" />
    551555                <field field_ref="dqthe" level="10" />

  • LMDZ5/trunk/DefLists/file_def_histmth_lmdz.xml

    r2952 r2958  
    592592                <field field_ref="q_th" level="4" />
    593593                <field field_ref="a_th" level="4" />
     594                <field field_ref="s_th" level="4" />
     595                <field field_ref="s_env" level="4" />
     596                <field field_ref="sigma_th" level="4" />
     597                <field field_ref="sigma_env" level="4" />
    594598                <field field_ref="d_th" level="4" />
    595599                <field field_ref="dqthe" level="4" />

  • LMDZ5/trunk/DefLists/file_def_histstn_lmdz.xml

    r2945 r2958  
    548548                <field field_ref="q_th" level="10" />
    549549                <field field_ref="a_th" level="10" />
     550                <field field_ref="s_th" level="10" />
     551                <field field_ref="s_env" level="10" />
     552                <field field_ref="sigma_th" level="10" />
     553                <field field_ref="sigma_env" level="10" />
    550554                <field field_ref="d_th" level="10" />
    551555                <field field_ref="dqthe" level="10" />

  • LMDZ5/trunk/libf/phylmd/cloudth_mod.F90

    r2957 r2958  
    819819
    820820      USE ioipsl_getin_p_mod, ONLY : getin_p
     821      USE phys_output_var_mod, ONLY : cloudth_sth,cloudth_senv, &
     822     &                                cloudth_sigmath,cloudth_sigmaenv
    821823
    822824      IMPLICIT NONE
    … …  
    11451147      alenv=(0.622*Lv*zqsatenv(ind1,ind2))/(rdd*zthl(ind1,ind2)**2) 
    11461148      aenv=1./(1.+(alenv*Lv/cppd))
    1147       senv=aenv*(po(ind1)-zqsatenv(ind1,ind2)) 
    1148      
     1149      senv=aenv*(po(ind1)-zqsatenv(ind1,ind2))
     1150      sth=0.
    11491151
    11501152      sigma1s=ratqs(ind1,ind2)*zqenv(ind1)
     1153      sigma2s=0.
    11511154
    11521155      xenv=senv/(sqrt2*sigma1s)
    … …  
    11661169 
    11671170      endif       ! From the separation (thermal/envrionnement) et (environnement) only, l.335 et l.492
     1171      ! Outputs used to check the PDFs
     1172      cloudth_senv(ind1,ind2) = senv
     1173      cloudth_sth(ind1,ind2) = sth
     1174      cloudth_sigmaenv(ind1,ind2) = sigma1s
     1175      cloudth_sigmath(ind1,ind2) = sigma2s
     1176
    11681177      enddo       ! from the loop on ngrid l.333
    11691178     

  • LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

    r2952 r2958  
    15711571  TYPE(ctrl_out), SAVE :: o_a_th = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
    15721572    'a_th', "Thermal plume fraction", "", (/ ('', i=1, 10) /))
     1573
     1574  TYPE(ctrl_out), SAVE :: o_cloudth_sth = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     1575    's_th', "Thermal plume saturation deficit", "kg/kg", (/ ('', i=1, 10) /))
     1576  TYPE(ctrl_out), SAVE :: o_cloudth_senv = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     1577    's_env', "Environment saturation deficit", "kg/kg", (/ ('', i=1, 10) /))
     1578  TYPE(ctrl_out), SAVE :: o_cloudth_sigmath = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     1579    'sigma_th', "Thermal plume gauss variance", "kg/kg", (/ ('', i=1, 10) /))
     1580  TYPE(ctrl_out), SAVE :: o_cloudth_sigmaenv = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     1581    'sigma_env', "Environment gauss variance", "kg/kg", (/ ('', i=1, 10) /))
     1582
    15731583  TYPE(ctrl_out), SAVE :: o_d_th = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
    15741584    'd_th', 'Thermal plume detrainment', 'K/s', (/ ('', i=1, 10) /))

  • LMDZ5/trunk/libf/phylmd/phys_output_var_mod.F90

    r2854 r2958  
    3838  !$OMP THREADPRIVATE(d_h_qw_col, d_h_ql_col, d_h_qs_col, d_h_col)
    3939
     40  ! Outputs used in cloudth_vert to extract the moments of the horizontal and
     41  ! vertical PDFs
     42  REAL, SAVE, ALLOCATABLE :: cloudth_sth(:,:),cloudth_senv(:,:)
     43  !$OMP THREADPRIVATE(cloudth_sth,cloudth_senv)
     44  REAL, SAVE, ALLOCATABLE :: cloudth_sigmath(:,:),cloudth_sigmaenv(:,:)
     45  !$OMP THREADPRIVATE(cloudth_sigmath,cloudth_sigmaenv)
     46
    4047! Marine
    4148! Variables de sortie du simulateur AIRS
    … …  
    141148    d_h_qw_col=0. ; d_h_ql_col=0. ; d_h_qs_col=0. ; d_h_col=0.
    142149
     150    ! Outputs used in cloudth_vert
     151    allocate(cloudth_sth(klon,klev))
     152    allocate(cloudth_senv(klon,klev))
     153    cloudth_sth = 0. ; cloudth_senv = 0.
     154    allocate(cloudth_sigmath(klon,klev))
     155    allocate(cloudth_sigmaenv(klon,klev))
     156    cloudth_sigmath = 0. ; cloudth_sigmaenv = 0.
     157
    143158! Marine
    144159! Variables de sortie simulateur AIRS
    … …  
    177192  &           , d_h_qw_col, d_h_ql_col, d_h_qs_col, d_h_col)
    178193
     194    ! Outputs used in cloudth_vert
     195    deallocate(cloudth_sth)
     196    deallocate(cloudth_senv)
     197    deallocate(cloudth_sigmath)
     198    deallocate(cloudth_sigmaenv)
     199
    179200! Marine
    180201! Variables de sortie simulateur AIRS

  • LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

    r2952 r2958  
    141141         o_duthe, o_dvthe, o_ftime_th, &
    142142         o_f_th, o_e_th, o_w_th, o_q_th, &
    143          o_a_th, o_d_th, o_f0_th, o_zmax_th, &
     143         o_a_th, o_cloudth_sth, o_cloudth_senv, &
     144         o_cloudth_sigmath, o_cloudth_sigmaenv, &
     145         o_d_th, o_f0_th, o_zmax_th, &
    144146         o_dqthe, o_dqthe2d, o_dtajs, o_dqajs, o_dqajs2d, o_dtswr, &
    145147         o_dtsw0, o_dtlwr, o_dtlw0, o_dtec, &
    … …  
    306308         zustr_gwd_front, zvstr_gwd_front,     &
    307309         sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o, &
     310         cloudth_sth,cloudth_senv,cloudth_sigmath,cloudth_sigmaenv, &
    308311! Marine
    309312         map_prop_hc, map_prop_hist, &
    … …  
    15701573          CALL histwrite_phy(o_q_th, zqasc)
    15711574          CALL histwrite_phy(o_a_th, fraca)
     1575          CALL histwrite_phy(o_cloudth_sth, cloudth_sth)
     1576          CALL histwrite_phy(o_cloudth_senv, cloudth_senv)
     1577          CALL histwrite_phy(o_cloudth_sigmath, cloudth_sigmath)
     1578          CALL histwrite_phy(o_cloudth_sigmaenv, cloudth_sigmaenv)
    15721579          CALL histwrite_phy(o_d_th, detr_therm)
    15731580          CALL histwrite_phy(o_f0_th, f0)
Note: See TracChangeset for help on using the changeset viewer.