Changeset 179 in lmdz_wrf for trunk/WRFV3/phys

Timestamp:

Oct 10, 2014, 12:40:39 PM (10 years ago)

Author:

lfita

Message:

No Ale, Alp, ... values because pphi is at the intermediate level!!!
Removing raw prints of checkig for Ale, Alp

File:

• trunk/WRFV3/phys/module_lmdz_phys.F (modified) (6 diffs)

trunk/WRFV3/phys/module_lmdz_phys.F

@@ -134 +134 @@
        & wrfUtenddestagg, wrfVtenddestagg
     REAL, DIMENSION((wime-wims-2*wbdyw)*(wjme-wjms-2*wbdyw), wkms:wkme-1) ::         &
-       & TValues, TtendValues, HalfPressValues, WFluxMassValues
+       & GeopotValues, TValues, TtendValues, HalfPressValues, WFluxMassValues
     REAL, DIMENSION((wime-wims-2*wbdyw)*(wjme-wjms-2*wbdyw), wkms:wkme-1) ::         &  !! DE-STAGGED !!
        & UValues, UtendValues
@@ -140 +140 @@
        & VValues, VtendValues
     REAL, DIMENSION((wime-wims-2*wbdyw)*(wjme-wjms-2*wbdyw), wkms:wkme) ::           &  !! Z-STAGG !!
-       & FullPressValues, GeopotValues
+       & FullPressValues
     REAL, DIMENSION((wime-wims-2*wbdyw)*(wjme-wjms-2*wbdyw)) ::                      &  
        & sfcGeopotValues
@@ -148 +148 @@
     REAL, DIMENSION(wims:wime,wkms:wkme,wjms:wjme)       :: wrfWValues,              &
        & wrfPERMASSValues, wrfZFluxMassValues
-    REAL, DIMENSION(wims:wime,wkms:wkme,wjms:wjme)       :: wrf_geopot, w3Dmat
+    REAL, DIMENSION(wims:wime,wkms:wkme,wjms:wjme)       :: wrf_geopot, w3Dmat,&
+       & wrf_geopot_half
     REAL, DIMENSION(wims:wime,wjms:wjme)                 :: wrf_sfcgeopot
 
@@ -264 +265 @@
 ! wrf_basegeopot: WRF base of geopotential
 ! wrf_geopot: WRF geopotential
+! wrf_geopot_half: WRF geopotential at the half eta level
 ! wrf_psfctend: WRF surface pressure tendency
 ! wrf_sfcgeopot: WRF surface geopotential
@@ -982 +984 @@
     wrf_P = wrf_baseP + wrf_perP
     wrf_geopot = wrf_basegeopot + wrf_pergeopot
+    wrf_geopot_half(:,wkms:wkme-1,:) = ( wrf_geopot(:,wkms:wkme-1,:) +       &
+      wrf_geopot(:,wkms+1:wkme,:))/2.
     wrf_sfcgeopot = wrf_grid%ht*g
 ! de-staggering winds
@@ -1070 +1074 @@
     CALL mat_vect(wrf_P, dimx, dimy, dimz, wbdyw, HalfPressValues)
     CALL mat_vect_zstagg(wrf_fullpres, dimx, dimy, dimz+1, wbdyw, FullPressValues)
-    CALL mat_vect_zstagg(wrf_geopot, dimx, dimy, dimz+1, wbdyw, GeopotValues)
+    CALL mat_vect(wrf_geopot_half, dimx, dimy, dimz, wbdyw, GeopotValues)
     CALL mat_vect(wrfZFluxMassValues, dimx, dimy, dimz, wbdyw, WFluxMassValues)