anchor.diff:  
   
             ! _________ SLEVE simple Implementation  _______________________
  -          SLEVE = .False.
  +          SLEVE = .True.
             
             if (SLEVE) then  ! should become sth like if (options%parameters%sleve) then
   
  
positive.code:         type(options_t), intent(in)     :: options
  
          real, allocatable :: temp(:,:,:), terrain_u(:,:), terrain_v(:,:)
          integer :: i, max_level
          real :: smooth_height, H, s
          logical :: SLEVE
  
  
negative.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 + terrain_u * (SINH( (H - dz(i)/2)/s ) / SINH(H/s) )
              z_v(:,i,:)  = dz(i)/2 + terrain_v * (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 )
  
              zr_u(:,i,:)  = (z_u(:,i,:) - terrain_u) / ( dz(i)/2 )
              zr_v(:,i,:)  = (z_v(:,i,:) - terrain_v) / (dz(i)/2 )
              
              ! xx=200 
              ! yy=1  ! center of ideal domain
              ! print*, "sinh term z_u: ", (SINH( (H - 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)
             
  
              ! - - - - -  higher k levels  - - - - - 
              do i = this%grid%kms+1, this%grid%kme
  