Changeset 2224 for LMDZ5/trunk/libf

Timestamp:

Mar 11, 2015, 12:23:40 PM (9 years ago)

Author:

crio

Message:

cv3p1_closure: Bug fix in the toothpaste closure computation (dimension of klfc).
physiq: check number of H2O tracers if thermodynamical effect of ice is included. Stop simulation if ice water is missing.

Location:

LMDZ5/trunk/libf/phylmd

Files:

• cv3p1_closure.F90 (modified) (3 diffs)
• physiq.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/cv3p1_closure.F90

@@ -53 +53 @@
 
   ! local variables:
-  INTEGER il, i, j, k, icbmax, i0(nloc), klfc
+  INTEGER il, i, j, k, icbmax, i0(nloc), klfc(nloc)
   REAL deltap, fac, w, amu
   REAL rhodp
@@ -525 +525 @@
 
 !CR:Compute k at plfc
+  DO il=1,ncum
+           klfc(il)=nl
+  ENDDO
   DO k=1,nl
      DO il=1,ncum
         if ((plfc(il).lt.ph(il,k)).and.(plfc(il).ge.ph(il,k+1))) then
-           klfc=k
+           klfc(il)=k
         endif
      ENDDO
@@ -540 +543 @@
 !CR: Add large-scale component to the mass-flux
 !encore connu sous le nom "Experience du tube de dentifrice"
-    if (coef_clos_ls.gt.0.) then 
-       cbmf1(il) = cbmf1(il) - coef_clos_ls*min(0.,1./RG*omega(il,klfc))
+    if ((coef_clos_ls.gt.0.).and.(plfc(il).gt.0.)) then 
+       cbmf1(il) = cbmf1(il) - coef_clos_ls*min(0.,1./RG*omega(il,klfc(il)))
     endif
 !RC

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -924 +924 @@
      print*, '================================================='
      !
+!CR: check sur le nb de traceurs de l eau
+     if ((iflag_ice_thermo.gt.0).and.(nqo==2)) then
+          WRITE (lunout, *) ' iflag_ice_thermo==1 requires 3 H2O tracers (H2Ov, H2Ol, H2Oi)', ' but nqo=', nqo, &
+          '. Might as well stop here.'
+          STOP
+     endif
+
      dnwd0=0.0
      ftd=0.0