anchor.diff:          real, allocatable :: dz_scl(:)
           integer :: i, max_level 
           real :: smooth_height, H, s, n
  -        logical :: SLEVE
  +        logical :: SLEVE, use_delta_terrain
           character :: filename, file_idS, file_idn
   
   
  
positive.code:                                           this% kds : this% kde,   &
                                            this% jds : this% jde)   )
  
          allocate(this%zfr_u(this% ims : this% ime, &
                                      this% kms : this% kme, &
                                      this% jms : this% jme) )
  
negative.code: 
          call io_write("zfr_u.nc", "zfr_u", zfr_u(:,:,:) ) ! check in plot
          
          print *, "zfr_u min/max: ", minval(zfr_u), " / ", maxval(zfr_u)
  
          ! =>  But how to scale the winds with this??
  
  