anchor.diff:          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%terrain) )                    call setup(this%forcing_terrain,          this%grid2d,    forcing_var=opt%parameters%hgtvar, list=this%variables_to_force)
  -        ! if (0<opt%vars_to_allocate( kVARS%terrain) )                    call setup(this%forcing_terrain_u,        this%u_grid2d) !,  forcing_var=opt%parameters%hgtvar, list=this%variables_to_force)
  -        ! if (0<opt%vars_to_allocate( kVARS%terrain) )                    call setup(this%forcing_terrain_v,        this%v_grid,  forcing_var=opt%parameters%hgtvar, list=this%variables_to_force)
  -        ! if (0<opt%vars_to_allocate( kVARS%terrain) )                    call setup(this%delta_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)
  
positive.code: 
              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
  
negative.code: 
  
  
  
      subroutine allocate_z_arrays(this)
          implicit none
  