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:                                z(k),      &  ! [INOUT]
                                 zcoords,   &  ! [IN ]
                                 wix(k),    &  ! [OUT] Zonal wind (m s^-1)
                                 wix(k),    &  ! [OUT] Meridional wind (m s^-1)
                                 tmp(k),    &  ! [OUT] Temperature (K)
                                 thetav(k), &  ! [OUT] Virtual potential Temperature (K)
                                 ps,        &  ! [OUT] Surface Pressure (Pa)
  
negative.code:        I_QV,      &
         NQW_MAX,   &
         TRC_vmax
      implicit none
  
      integer,          intent(in)  :: ijdim
  