Changeset 174 in lmdz_wrf for branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_plume.F90

Timestamp:

Oct 2, 2014, 11:32:36 AM (10 years ago)

Author:

lfita

Message:

NaN values due to negative topography (topography below sea level) due to:

• pphi (Geopotential in LMDZ) is relative to the surface. Thus, one should pass

pphi = P + PB - g*HGT

File:

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

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_plume.F90

@@ -603 +603 @@
          do ig=1,ngrid
             if (ztv(ig,l)> ztv(ig,l+1) .and. ztv(ig,1)>=ztv(ig,l) ) then
-! 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
-               END IF
+!!! 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
+!!               END IF
                alim_star_tot(ig)=alim_star_tot(ig)+alim_star(ig,l)
             endif