Changeset 5631 for LMDZ6


Ignore:
Timestamp:
Apr 26, 2025, 1:32:12 AM (6 weeks ago)
Author:
aborella
Message:

Corrections to the deep convection-stratiform clouds coupling

Location:
LMDZ6/branches/contrails/libf/phylmd
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp.f90

    r5626 r5631  
    12311231        !--the sink of condensed water from precipitation
    12321232        IF ( ptconv(i,k) ) THEN
    1233           qvcon_old(i,k) = qvcon(i,k)
    1234           qccon_old(i,k) = qccon(i,k) * zcond(i) / zoliq(i)
     1233          IF ( zoliq(i) .GT. 0. ) THEN
     1234            qvcon_old(i,k) = qvcon(i,k)
     1235            qccon_old(i,k) = qccon(i,k) * zcond(i) / zoliq(i)
     1236          ELSE
     1237            qvcon_old(i,k) = 0.
     1238            qccon_old(i,k) = 0.
     1239          ENDIF
    12351240        ELSE
    12361241          qvcon_old(i,k) = 0.

  • LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_condensation.f90

    r5626 r5631  
    515515
    516516        qvapincld = qvc(i) / cldfra(i)
     517        IF ( qvapincld .GT. gamma_cond(i) * qsat(i) ) THEN
     518          qvapincld = gamma_cond(i) * qsat(i)
     519          qvc(i) = qvapincld * cldfra(i)
     520        ENDIF
    517521        qiceincld = ( qcld(i) / cldfra(i) - qvapincld )
    518522       
    … …  
    576580          !------------------------------------
    577581          !--Additionally to a minimum in cloud water vapor, we impose a minimum
    578           !--in in-cloud ice water content. It is calculated following
     582          !--on the in-cloud ice water content. It is calculated following
    579583          !--Marti and Mauersberger (1993), see also Schiller et al. (2008)
    580584          qiceincld_min = qsati_diss(i) - qsat(i)

  • LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_ini.f90

    r5625 r5631  
    487487    CALL getin_p('depo_coef_cirrus',depo_coef_cirrus)
    488488    CALL getin_p('capa_cond_cirrus',capa_cond_cirrus)
     489    CALL getin_p('N_ice_volume',N_ice_volume)
    489490    CALL getin_p('nu_iwc_pdf_lscp',nu_iwc_pdf_lscp)
    490491    CALL getin_p('beta_pdf_lscp',beta_pdf_lscp)
    … …  
    586587    WRITE(lunout,*) 'lscp_ini, depo_coef_cirrus:', depo_coef_cirrus
    587588    WRITE(lunout,*) 'lscp_ini, capa_cond_cirrus:', capa_cond_cirrus
     589    WRITE(lunout,*) 'lscp_ini, N_ice_volume:', N_ice_volume
    588590    WRITE(lunout,*) 'lscp_ini, nu_iwc_pdf_lscp:', nu_iwc_pdf_lscp
    589591    WRITE(lunout,*) 'lscp_ini, beta_pdf_lscp:', beta_pdf_lscp
Note: See TracChangeset for help on using the changeset viewer.