anchor.diff:      real(DP) :: ps
   
       real(RP) :: cl, cl2
  -!    real(RP) :: uave, dz, f_cf
  -    real(RP) :: tv(kdim)
   
       integer, parameter :: deep    = 0 ! deep atmosphere (1 = yes or 0 = no)
       integer, parameter :: zcoords = 1 ! 1 if z is specified, 0 if p is specified
  
positive.code:     real(DP) :: ps, phis
  
      real(RP) :: cl, cl2
  !    real(RP) :: uave, dz, f_cf
      real(RP) :: tv(kdim)
  
      integer, parameter :: deep    = 0 ! deep atmosphere (1 = yes or 0 = no)
      integer, parameter :: zcoords = 1 ! 1 if z is specified, 0 if p is specified
  
negative.code:     real(DP) :: ps, phis
  
      real(RP) :: cl, cl2
  
      integer, parameter :: deep    = 0 ! deep atmosphere (1 = yes or 0 = no)
      integer, parameter :: zcoords = 1 ! 1 if z is specified, 0 if p is specified
      integer            :: moist       ! include moisture (1 = yes or 0 = no)
      integer            :: pertt       ! type of perturbation (0 = exponential, 1 = stream function)
      logical, parameter :: prs_dry = .false.
  
      integer :: n, k, l
      !---------------------------------------------------------------------------
  