Changeset 736 for trunk/LMDZ.GENERIC/libf

Timestamp:

Jul 24, 2012, 6:50:33 PM (12 years ago)

Author:

jleconte

Message:

24/07/2012 == JL

• Correction of a bug in turbulent diffusion (turbdiff.F90)

=> This solves a water conservation problem arising when the code tries to
evaporate over dry land.

File:

• trunk/LMDZ.GENERIC/libf/phystd/turbdiff.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phystd/turbdiff.F90

@@ -580 +580 @@
                       !recompute q1 with new water flux from surface  
                       zq(ig,1,iq) = (zmass(ig,1)*(pq(ig,1,iq)+ptimestep*pdqfi(ig,1,iq))  &
-                                            +zfluxq(ig,2)*zcq(ig,1)-dqsdif_total(ig))     &
+                                            +zfluxq(ig,2)*zcq(ig,2)-dqsdif_total(ig))     &
                                  / (zmass(ig,1)+(1.-zdq(ig,2))*zfluxq(ig,2))                 
                   end if