anchor.diff:            enddo
          enddo
   
  +       precip(:) = precip(:) + precip2(:)
  +
       endif
   
       if ( USE_ToyChemistry ) then
  
positive.code:     real(RP) :: pdel (ijdim,vlayer)   ! Layer thickness (Pa)
      real(RP) :: rpdel(ijdim,vlayer)   ! Reciprocal of layer thickness (1/Pa)
      real(RP) :: ps   (ijdim)          ! surface pressure output [dummy]
      integer  :: test
  
      ! for toy-chemistory
  
negative.code:                             pdel  (:,:),       & ! [INOUT] but not changed
                              rpdel (:,:),       & ! [INOUT] but not changed
                              ps    (:),         & ! [INOUT] but not changed
                              precip2(:),        & ! [OUT]
                              test,              & ! [IN]
                              SM_LargeScaleCond, & ! [IN]
                              SM_PBL_Bryan       ) ! [IN]
  