anchor.diff:  
               enddo  ! ____ end SLEVE simple Implementation  _______
               
  -
  -            ! if (options%parameters%debug) then 
  -            !   ! Write the z_level_ratio for analysis:
  -            !   write(file_idn, *) options%parameters%sleve_n
  -            !   write(file_idS, *) options%parameters%sleve_decay_factor
  -            !   filename = 'zlr_n_' // trim(adjustl(file_idn)) // '_S_'//trim(adjustl(file_idS)) //'.nc'
  -            !   call io_write(filename, "z_level_ratio", z_level_ratio(:,:,:) )
  -            !   ! call io_write("zlr.nc", "z_level_ratio", z_level_ratio(:,:,:) )
  -            !   ! call io_write("dzdx.nc", "dzdx", dzdx(:,:,:) )
  -            ! endif
  -
   
             else  !. i.e. no sleve coordinates
               i = this%grid%kms
  
positive.code: 
  
  
      !> -------------------------------
      !! Initialize various domain variables, mostly z, dz, etc.
      !!
      !! -------------------------------
      ! subroutine initialize_core_variables(this, options)
  ! =======
  ! >>>>>>> feature/coarray
  !         implicit none
  !         type(options_t), intent(in) :: options
  !         integer,         intent(in) :: nz
  !         real,            intent(in) :: dz(:)
  !         integer :: max_level
  
  !         integer :: j
  !         real :: height
  
  !         if (options%parameters%flat_z_height > nz) then
  !             if (this_image()==1) write(*,*) "Treating flat_z_height as specified in meters above mean terrain height: ", options%parameters%flat_z_height," meters"
  !             height = 0
  !             do j = 1, nz
  !                 if (height <= options%parameters%flat_z_height) then
  !                     height = height + dz(j)
  !                     max_level = j
  !                 endif
  !             enddo
  
  ! <<<<<<< HEAD
  !         real, allocatable :: temp(:,:,:), terrain_u(:,:), terrain_v(:,:)
  !         real, allocatable :: dz_scl(:)
  !         integer :: i, max_level, xx, yy
  !         real :: smooth_height, H, s, h_avg
  !         logical :: SLEVE
  
  ! =======
  !         elseif (options%parameters%flat_z_height <= 0) then
  !             if (this_image()==1) write(*,*) "Treating flat_z_height as counting levels down from the model top: ", options%parameters%flat_z_height," levels"
  !             max_level = nz + options%parameters%flat_z_height
  ! >>>>>>> feature/coarray
  
  !         else
  !             if (this_image()==1) write(*,*) "Treating flat_z_height as counting levels up from the ground: ", options%parameters%flat_z_height," levels"
  !             max_level = options%parameters%flat_z_height
  !         endif
  
  !     end function find_flat_model_level
  
  
      subroutine allocate_z_arrays(this)
          implicit none
  
negative.code:               write(*,*) "sum(dz) ", sum(dz(1:max_level))
                print*, ""
              endif
  
  
              i=kms 
  