Changeset 7086


Ignore:
Timestamp:
02/16/21 16:25:05 (5 months ago)
Author:
nicolas.vuichard
Message:

bug correction on qsurf introduced in rev [6019]. See ticket #762

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/ORCHIDEE_Quest/ORCHIDEE/src_sechiba/enerbil.f90

    r6019 r7086  
    12091209    REAL(r_std), DIMENSION (kjpindex)                            :: fluxsens_tmp 
    12101210    REAL(r_std), DIMENSION (kjpindex)                            :: fluxlat_tmp 
    1211  
     1211    REAL(r_std), DIMENSION (kjpindex)                            :: betatot           !! beta used in the computation of qsurf  (-) 
    12121212!_ ================================================================================================================================ 
    12131213     
     
    12651265      !!     \input{enerbilflux6.tex} 
    12661266      !! \endlatexonly 
    1267       qsurf(ji) = (vbeta1(ji) * (un - vbeta5(ji)) + vbeta5(ji)) * qsol_sat_new(ji) + & 
    1268            & (un - vbeta1(ji))*(un - vbeta5(ji)) * vbeta(ji) * qsol_sat_new(ji) 
    1269        
    1270       qsurf(ji)=qsurf(ji) + (1 - (vbeta1(ji) * (un - vbeta5(ji)) + vbeta5(ji)))*qair(ji) + &                           
    1271            & (1- (un - vbeta1(ji))*(un - vbeta5(ji)) * vbeta(ji) )*qair(ji) 
    1272        
     1267       
     1268       betatot(ji) = vbeta1(ji) * (un - vbeta5(ji)) + vbeta5(ji) + & 
     1269            & (un - vbeta1(ji))*(un - vbeta5(ji)) * vbeta(ji) 
     1270  
     1271       qsurf(ji) = betatot(ji) * qsol_sat_new(ji) + (un - betatot(ji)) * qair(ji)     
     1272 
    12731273      ! When no turbulence and no wind (qc small) the near surface air moisture is equal to qair. 
    12741274      ! This is related to problem with interpolation in stdlevvar in LMDZ.  
Note: See TracChangeset for help on using the changeset viewer.