anchor.diff:      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)
   
       integer :: n, k, l
       !---------------------------------------------------------------------------
   
       DIAG_var(:,:,:,:) = 0.0_RP
       p = 0.0_RP
  -
  -!    moist = 0
  -!    pertt = 0
   
       write(ADM_LOG_FID,*) "### DO NOT INPUT ANY TOPOGRAPHY ###"
       if ( NQW_MAX < 3 ) then
  
positive.code:        I_QV,      &
         NQW_MAX,   &
         TRC_vmax
  !       NCHEM_MAX, &
  !       NCHEM_STR, &
  !       NCHEM_END
      implicit none
  
      integer,          intent(in)  :: ijdim
  
negative.code:                                       z(k),      &  ! [INOUT]
                                        zcoords,   &  ! [IN ]
                                        wix(k),    &  ! [OUT] Zonal wind (m s^-1)
                                        wiy(k),    &  ! [OUT] Meridional wind (m s^-1)
                                        tmp(k),    &  ! [OUT] Temperature (K)
                                        thetav(k), &  ! [OUT] Virtual potential Temperature (K)
                                        phis,      &  ! [OUT] Surface Geopotential (m^2 s^-2)
  