Changeset 132 in lmdz_wrf

Timestamp:

Jul 28, 2014, 6:19:53 PM (10 years ago)

Author:

lfita

Message:

Adding checks for rhobarz0

File:

• branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_main.F90 (modified) (1 diff)

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_main.F90

@@ -953 +953 @@
     do ig =1,ngrid
      zmax(ig)=pphi(ig,lmax(ig))/rg
-     if (ok_lcl(ig)) then 
-      rhobarz0(ig)=rhobarz(ig,klcl(ig))+(rhobarz(ig,klcl(ig)+1) &
+!     if (ok_lcl(ig)) then 
+! L. Fita, LMD July 2014. Adding avoiding divisons by zero...
+     if ( ok_lcl(ig) .and. (klcl(ig).gt.0) .and. (klcl(ig)+1.le.nlay-1) .and.        &
+       (klcl(ig)+1.gt.0) ) then
+       rhobarz0(ig)=rhobarz(ig,klcl(ig))+(rhobarz(ig,klcl(ig)+1) &
  &               -rhobarz(ig,klcl(ig)))*interp(ig)
-      zlcl(ig)=(pplev(ig,1)-pcon(ig))/(rhobarz0(ig)*RG)
-      zlcl(ig)=min(zlcl(ig),zmax(ig))   ! Si zlcl > zmax alors on pose zlcl = zmax
+       zlcl(ig)=(pplev(ig,1)-pcon(ig))/(rhobarz0(ig)*RG)
+       zlcl(ig)=min(zlcl(ig),zmax(ig))   ! Si zlcl > zmax alors on pose zlcl = zmax
      else
       rhobarz0(ig)=0.