anchor.diff:  
   
             ! _________ SLEVE simple Implementation  _______________________
  -          SLEVE = .False.
  +          SLEVE = .True.
             
             if (SLEVE) then  ! should become sth like if (options%parameters%sleve) then
   
  
positive.code:             terrain_v =  z_v(:,i,:)  ! save for later on 
  
              ! Offset analogous to: z_u(:,i,:) = z_u(:,i,:) + dz(i) / 2 * zr_u(:,i,:)
              z_u(:,i,:)  = dz(i)/2 + z_u(:,i,:) * (SINH( (H - dz(i)/2)/s ) / SINH(H/s) )
              z_v(:,i,:)  = dz(i)/2 + z_v(:,i,:) * (SINH( (H - dz(i)/2)/s ) / SINH(H/s) )
              
              ! z_r_u       = (dz(i)/2 + terrain_u * (SINH( (H - dz(i)/2)/s ) / SINH(H/s) ) )  &
              !                 / ( terrain_u + dz(i) / 2 )  ! same as below
              zr_u(:,i,:)  = z_u(:,i,:) / (terrain_u + dz(i)/2 )
              zr_v(:,i,:)  = z_v(:,i,:) / (terrain_v + dz(i)/2 )
              
              
  
              ! - - - - -  higher k levels  - - - - - 
              do i = this%grid%kms+1, this%grid%kme
  
negative.code: 
                  ! - - - - -   u/v grid calculations - - - - -
                  z_u(:,i,:)  = sum(dz(1:(i-1))) + dz(i)/2  &
                                + terrain_u * (SINH( (H - sum(dz(1:(i-1))) - dz(i)/2 )/ s ) / SINH(H/s) )
  
                  z_v(:,i,:)  = sum(dz(1:(i-1))) + dz(i)/2  &
                                + terrain_v * (SINH( (H - sum(dz(1:(i-1))) - dz(i)/2)/s ) / SINH(H/s) )
                  
                  zr_u(:,i,:)  = (z_u(:,i,:) - z_u(:,i-1,:)) / (dz(i)/2 + dz(i-1)/2 )
                  zr_v(:,i,:)  = (z_v(:,i,:) - z_v(:,i-1,:)) / (dz(i)/2 + dz(i-1)/2 )
  
                  ! print*, ""
                  ! print*, "i - sum(dz(1:(i-1))) : ",i, " - ", sum(dz(1:(i-1)))
                  ! print*, "sinh term z_u: ", (SINH( (H - sum(dz(1:(i-1))) + dz(i)/2 )/ s ) / SINH(H/s) )
                  ! print*, "z_u (",xx,i,yy,") : ", z_u(xx,i,yy)
                  ! print*, "zr_u(",xx,i,yy,") : ", zr_u(xx,i,yy)
              enddo
              ! ____ end SLEVE simple Implementation  _______
              ! ############################################
  
  
              ! ! # # # # # #   BK debug wind field:    # # # # # # # # #
              ! if (options%parameters%debug ) then
              !     xx=200 
              !     yy=1  ! center of ideal domain
              !     print*, ""
              !     print*, "zru shape: ", SHAPE(zr_u)
              !     print*, "z_level_ratio shape: ", SHAPE(z_level_ratio)
              !     print*, ""
              !     print*, "z_level_ratio(",xx,",i,",yy,") : ", z_level_ratio(xx,:,yy)
              !     print*, "zr_u(",xx,",i,",yy,") : ", zr_u(xx,:,yy)
              !     ! print*, "z_level_ratio(150,i,150) : ", z_level_ratio(150,:,150)
              !     print*, "dzdx at  i,j,k = 200,1,k: " !, k, " : "
              !     print*, dzdx(xx,:,yy)
              !     ! print*, "dzdy at i,j,k = 150,150,k: "!, k, " : "
              !     ! print*, dzdy(150,:,150)
              ! endif 
  
  
              call io_write("zr_u.nc", "zr_u", zr_u(:,:,:) )
              call io_write("zrl.nc", "zrl", z_level_ratio(:,:,:) )
              call io_write("dzdx.nc", "dzdx", dzdx(:,:,:) )
  
  
            else  
              i = this%grid%kms
  