Changeset 174 in lmdz_wrf for branches/LMDZ_WRFmeas

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

Location:

branches/LMDZ_WRFmeas/WRFV3

Files:

• lmdz/thermcell_plume.F90 (modified) (1 diff)
• phys/module_lmdz_phys.F (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

branches/LMDZ_WRFmeas/WRFV3/phys/module_lmdz_phys.F

@@ -1081 +1081 @@
     preff=100000.
     CALL eta_to_pressure(wrf_halfeta, preff, wrf_Ptop, dimz, lmdzoutP)
+! L. Fita, LMD October 2014. LMDZ geopotential values are relative to the surface 
+    PRINT *,'  Lluis UBOUND GeopotValues:',UBOUND(GeopotValues), ' sfcZ:',      &
+      UBOUND(sfcGeopotValues)
+    PRINT *,'  Lluis: lp',lp,' z GeopotValues sfcGeopotValues ______'
+    DO iz=1, dimz
+      PRINT *,GeopotValues(lp,iz),sfcGeopotValues(lp),GeopotValues(lp,iz) -     &
+        sfcGeopotValues(lp)
+      GeopotValues(:,iz) = GeopotValues(:,iz) - sfcGeopotValues(:)
+    END DO
 
 ! We need to provide values to these variables in order to make it works