anchor.diff:      logical,          intent(in)  :: nicamcore
       real(RP),         intent(out) :: DIAG_var(ijdim,kdim,lall,6+TRC_VMAX)
   
  -!    real(DP), parameter :: Xfact = 1.0_DP ! Earth scaling parameter
       real(DP) :: lat, lon               ! latitude, longitude on Icosahedral grid
       real(DP) :: prs(kdim), tmp(kdim)   ! presssure and temperature in ICO-grid field
       real(DP) :: wix(kdim),   wiy(kdim) ! zonal/meridional wind components in ICO-grid field
  
positive.code:     real(RP) :: vy_local(kdim)
      real(RP) :: vz_local(kdim)
      real(DP) :: ps
  !, phis
  
      real(RP) :: cl, cl2
  !    real(RP) :: uave, dz, f_cf
  
negative.code: 
      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            :: pert       ! type of perturbation (0 = no perturbation, 1 = perturbation)
  
      integer :: n, k, l
  