anchor.diff:      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
  
positive.code: 
    !-----------------------------------------------------------------------------
    subroutine diag_pressure( &
        kdim, &  !--- IN : # of z dimension
        z,    &  !--- IN : height information
        rho,  &  !--- IN : density
        t,    &  !--- IN : temperature
        q,    &  !--- IN : water vapor content
        ps,   &  !--- IN : surface pressure
        prs   )  !--- OUT: pressure
      !
      implicit none
      integer, intent(in)   :: kdim
      real(RP), intent(in)  :: z(:)
      real(RP), intent(in)  :: rho(:)
      real(RP), intent(in)  :: t(:)
      real(RP), intent(in)  :: q(:)
      real(RP), intent(in)  :: ps
      real(RP), intent(out) :: prs(:)
      !
      integer :: k
      !-----
      do k=1, kdim
         if ( z(k) < 1.0D-7 ) then
            prs(k) = ps
         else
            prs(k) = rho(k) * t(k) * ( (1.0_RP-q(k))*Rd + q(k)*Rv )
         endif
      enddo
  
      return
    end subroutine diag_pressure
  
negative.code:     integer,          intent(in)  :: lall
      character(len=*), intent(in)  :: test_case
      logical,          intent(in)  :: nicamcore
      logical,          intent(in)  :: prs_rebuild
      real(RP),         intent(out) :: DIAG_var(ijdim,kdim,lall,6+TRC_VMAX)
  
      real(DP) :: lat, lon               ! latitude, longitude on Icosahedral grid
  