Changeset 114 in lmdz_wrf

Timestamp:

Jul 25, 2014, 6:15:29 PM (10 years ago)

Author:

lfita

Message:

Adding ZLEV < 0. for alim_star computations

File:

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

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_plume.F90

@@ -666 +666 @@
          do ig=1,ngrid
             if (ztv(ig,l)> ztv(ig,l+1) .and. ztv(ig,1)>=ztv(ig,l) ) then
-               alim_star(ig,l)=MAX((ztv(ig,l)-ztv(ig,l+1)),0.)  &
+! L. Fita, LMD July 2014. Adding IF for zlev == 0.
+               IF ( zlev(ig,l+1) < 0.) THEN
+                 alim_star(ig,l) == 0. 
+               ELSE
+                 alim_star(ig,l)=MAX((ztv(ig,l)-ztv(ig,l+1)),0.)  &
      &                       *sqrt(zlev(ig,l+1)) 
-               lalim(ig)=l+1
+                 lalim(ig)=l+1
+               END IF
                alim_star_tot(ig)=alim_star_tot(ig)+alim_star(ig,l)
             IF (zlev(ig,l+1) < 0. .OR. alim_star_tot(ig) == 0.) THEN