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:               write(*,*) "sum(dz) ", sum(dz(1:max_level))
                print*, ""
              endif
  
              !     call geo_LUT(this%geo_u, forcing%geo_u)
              ! call geo_LUT(this%geo_v, forcing%geo_v)
              ! call geo_LUT(this%geo,   forcing%geo)
              ! print*, "this%geo size: ", size(this%geo)  
              ! print*," forcing%geo size: ",  size(forcing%geo)
  
  
  
              i=kms 
  
negative.code:         if (0<opt%vars_to_allocate( kVARS%v_latitude) )                 call setup(this%v_latitude,               this%v_grid2d)
          if (0<opt%vars_to_allocate( kVARS%v_longitude) )                call setup(this%v_longitude,              this%v_grid2d)
          if (0<opt%vars_to_allocate( kVARS%terrain) )                    call setup(this%terrain,                  this%grid2d)
          if (0<opt%vars_to_allocate( kVARS%sensible_heat) )              call setup(this%sensible_heat,            this%grid2d)
          if (0<opt%vars_to_allocate( kVARS%latent_heat) )                call setup(this%latent_heat,              this%grid2d)
          if (0<opt%vars_to_allocate( kVARS%u_10m) )                      call setup(this%u_10m,                    this%grid2d)
  