Ignore:
Timestamp:
Apr 22, 2024, 9:06:52 PM (4 weeks ago)
Author:
evignon
Message:

correction formulation de l'erosion de la neige soufflee suite aux investigations de Nicolas Chiabrando,
prise en compte d'un temps d'erosion de la neige fraiche nouvellement accumulee

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/phylmd/lmdz_blowing_snow_sublim_sedim.F90

    r4835 r4916  
    7474
    7575
    76 Ka=2.4e-2      ! thermal conductivity of the air, SI
    7776p0=101325.0    ! ref pressure
    7877
    … …  
    229228              IF (zqv(i) .LT. qsi(i)) THEN
    230229                 rhoair=zpres(i)/ztemp(i)/RD
    231                  !diffusivity of water vapor
    232230                 Dv=0.0001*0.211*(p0/zpres(i))*((ztemp(i)/RTT)**1.94) ! water vapor diffusivity in air, SI
    233231                 Ka=(5.69+0.017*(ztemp(i)-RTT))*1.e-5*100.*4.184                ! thermal conductivity of the air, SI
    … …  
    241239                 dqbsub=(-b_p+sqrt(delta_p))/(2.*a_p) - zqb(i)       
    242240                 dqbsub = MIN(0.0,MAX(dqbsub,-zqb(i)))
    243 
    244241                 ! Sublimation limit: we ensure that the whole mesh does not exceed saturation wrt ice
    245242                 maxdqbsub = MAX(0.0, qsi(i)-zqv(i))
Note: See TracChangeset for help on using the changeset viewer.