comment out "extending to surface" to get zero-diff

GEOS-ESM · Jul 26, 2023 · be8b1e3 · be8b1e3
1 parent 0a5ab38
commit be8b1e3
Showing 1 changed file with 8 additions and 8 deletions.
diff --git a/base/MAPL_VerticalMethods.F90 b/base/MAPL_VerticalMethods.F90
@@ -309,10 +309,10 @@ subroutine init_indices(rc)
             do lev =1, levo
                pp = flip_sign*this%interp_levels(lev)
                pb = flip_sign*ple3d(:,:,km_e)
-               where (pp > pb)
-                  this%ks_e(:,:,lev) = km_e
-                  this%weight_e(:,:,lev) = 1.
-               endwhere
+               !where (pp > pb)
+               !   this%ks_e(:,:,lev) = km_e
+               !   this%weight_e(:,:,lev) = 1.
+               !endwhere
                do k = km_e-1, 1,-1 ! levels of input
                  if(all(pb<pp)) exit
                  pt = flip_sign*ple3d(:,:,k)
@@ -360,10 +360,10 @@ subroutine regrid_eta_to_pressure(this,ptrin,ptrout,rc)
             do i = 1, D1
                k = ks_(i,j,lev)
                if (k == -1) cycle
-               if (k == km) then
-                  ptrout(i,j,lev) = ptrin(i,j,k)
-                  cycle
-               endif
+               !if (k == km) then
+               !   ptrout(i,j,lev) = ptrin(i,j,k)
+               !   cycle
+               !endif
                weight = weights_(i,j,lev)
                if (ptrin(i,j,k)   == MAPL_UNDEF) then
                   ptrout(i,j,lev) = ptrin(i,j,k+1)