anchor.diff:         lall,         &
          test_case,    &
          nicamcore,    &
  +       prs_rebuild,  &
          DIAG_var      )
       use mod_adm, only: &
          ADM_proc_stop, &
  
positive.code:                                       q(k)       )  ! [OUT] water vapor mixing ratio (kg/kg)
         enddo
  
         call diag_pressure( kdim, z, rho, tmp, q, ps, prs )
         call conv_vxvyvz ( kdim, lat, lon, wix, wiy, vx_local, vy_local, vz_local )
  
         do k = 1, kdim
  
negative.code:     character(len=*), intent(in)  :: test_case
      logical,          intent(in)  :: nicamcore
      logical,          intent(in)  :: chemtracer
      logical,          intent(in)  :: prs_rebuild
      real(RP),         intent(out) :: DIAG_var(ijdim,kdim,lall,6+TRC_VMAX)
  
      real(DP), parameter :: Xfact = 1.0_DP ! Earth scaling parameter
  