Changeset 5620 for LMDZ6


Ignore:
Timestamp:
Apr 15, 2025, 1:46:48 PM (2 months ago)
Author:
evignon
Message:

correction du calcul de la temperature virtuelle dans thermcell_plume_6A.
Attention! Cette commission conduit une perte de convergence numérique

File:
1 edited

Legend:

Unmodified
Added
Removed

  • TabularUnified LMDZ6/trunk/libf/phylmd/lmdz_thermcell_plume_6A.f90

    r5512 r5620  
    224224        zta_est(ig,l)=ztva_est(ig,l)
    225225        ztva_est(ig,l) = ztva_est(ig,l)/zpspsk(ig,l)
    226         ztva_est(ig,l) = ztva_est(ig,l)*(1.+RETV*(zqta(ig,l-1)  &
    227      &      -zqla_est(ig,l))-zqla_est(ig,l))
     226        ztva_est(ig,l) = ztva_est(ig,l)*(1.+RETV*(zqta(ig,l-1)-zqla_est(ig,l)))
    228227 
    229228
    … …  
    566565!on rajoute le calcul de zha pour diagnostiques (temp potentielle)
    567566           zha(ig,l) = ztva(ig,l)
    568            ztva(ig,l) = ztva(ig,l)*(1.+RETV*(zqta(ig,l)  &
    569      &              -zqla(ig,l))-zqla(ig,l))
     567           ztva(ig,l) = ztva(ig,l)*(1.+RETV*(zqta(ig,l)-zqla(ig,l)))
    570568           zbuoy(ig,l)=RG*(ztva(ig,l)-ztv(ig,l))/ztv(ig,l)
    571569           zdz=zlev(ig,l+1)-zlev(ig,l)
Note: See TracChangeset for help on using the changeset viewer.