anchor.diff:      use mod_gtl, only: &
          GTL_max_k, &
          GTL_min_k
  +    use mod_runconf, only: &
  +       DYN_DIV_NUM
       implicit none
   
       character(len=*), intent(in) :: hdiff_type      ! type of horizontal diffusion
  
positive.code:     real(RP) :: kh_max(ADM_kall)
      real(RP) :: d2T_dx2, coef
  
      integer :: g, k, l, nq, p
      !---------------------------------------------------------------------------
  
  
negative.code: 
      if( alpha_v > 0.0_RP ) NUMFILTER_DOdivdamp_v = .true.
  
      small_step_dt = TIME_DTS / real(DYN_DIV_NUM,kind=RP)
  
      divdamp_coef_v = -alpha_v * SOUND * SOUND * small_step_dt
  
      return
    end subroutine numfilter_divdamp_setup
  