anchor.diff:                            - MAXVAL(h2)/s2 * COSH(smooth_height/s2)/SINH(smooth_height/s2)
   
               ! with the new (leuenberger et al 2010) Sleve formulation, the inveribiltiy criterion is as follows:
  +            ! ( Although an argument could be made to calculate this on the offset (u/v) grid b/c that is most
  +            !   relevant for advection? In reality this is probably a sufficient approximation, as long as we
  +            !   aren't pushing the gamma factor too close to zero )
               allocate(gamma_n(this%kds : this%kde+1))
               i=kms                          
               gamma_n(i) =  1                                                     &    
  
positive.code:                 print*, "    Using a sleve_n of ", options%parameters%sleve_n
                  write(*,*) "    Smooth height is ", smooth_height, "m.a.s.l     (model top ", sum(dz(1:nz)), "m.a.s.l.)"
                  write(*,*) "    invertibility parameter gamma is: ", gamma_min
                  if(gamma_min <= 0) print*, " CAUTION: coordinate transformation is not invertible (gamma <= 0 ) !!! reduce decay rate(s)!"
                  if(options%parameters%debug)  write(*,*) "   (btw: for (debugging) reference, 'gamma(n=1)'= ", gamma,")"
                  print*, ""
              endif
  
  
negative.code:                     zr_u(:,i,:)  = (z_u(:,i,:) - z_u(:,i-1,:)) / (dz_scl(i)/2 + dz_scl(i-1)/2 )  ! if dz_scl(i-1) = 0 (and no error)  k=1 can be included
                      zr_v(:,i,:)  = (z_v(:,i,:) - z_v(:,i-1,:)) / (dz_scl(i)/2 + dz_scl(i-1)/2 )
  
                      
                  else ! above the flat_z_height
  
                      zr_u(:,i,:) = 1
  