Changeset 7086

Timestamp:

2021-02-16T16:25:05+01:00 (4 years ago)

Author:

nicolas.vuichard

Message:

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

File:

• branches/ORCHIDEE_Quest/ORCHIDEE/src_sechiba/enerbil.f90 (modified) (2 diffs)

branches/ORCHIDEE_Quest/ORCHIDEE/src_sechiba/enerbil.f90

@@ -1209 +1209 @@
     REAL(r_std), DIMENSION (kjpindex)                            :: fluxsens_tmp
     REAL(r_std), DIMENSION (kjpindex)                            :: fluxlat_tmp
-
+    REAL(r_std), DIMENSION (kjpindex)                            :: betatot           !! beta used in the computation of qsurf  (-)
 !_ ================================================================================================================================
     
@@ -1265 +1265 @@
       !!     \input{enerbilflux6.tex}
       !! \endlatexonly
-      qsurf(ji) = (vbeta1(ji) * (un - vbeta5(ji)) + vbeta5(ji)) * qsol_sat_new(ji) + &
-           & (un - vbeta1(ji))*(un - vbeta5(ji)) * vbeta(ji) * qsol_sat_new(ji)
-      
-      qsurf(ji)=qsurf(ji) + (1 - (vbeta1(ji) * (un - vbeta5(ji)) + vbeta5(ji)))*qair(ji) + &                          
-           & (1- (un - vbeta1(ji))*(un - vbeta5(ji)) * vbeta(ji) )*qair(ji)
-      
+      
+       betatot(ji) = vbeta1(ji) * (un - vbeta5(ji)) + vbeta5(ji) + &
+            & (un - vbeta1(ji))*(un - vbeta5(ji)) * vbeta(ji)
+ 
+       qsurf(ji) = betatot(ji) * qsol_sat_new(ji) + (un - betatot(ji)) * qair(ji)    
+
       ! When no turbulence and no wind (qc small) the near surface air moisture is equal to qair.
       ! This is related to problem with interpolation in stdlevvar in LMDZ.