anchor.diff:      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]
  +    real(RP) :: precip2 (ijdim)
       integer  :: test
   
       ! for toy-chemistory
  
positive.code:                             pdel  (:,:),       & ! [INOUT] but not changed
                              rpdel (:,:),       & ! [INOUT] but not changed
                              ps    (:),         & ! [INOUT] but not changed
                              precip(:),         & ! [INOUT]
                              test,              & ! [IN]
                              SM_LargeScaleCond, & ! [IN]
                              SM_PBL_Bryan       ) ! [IN]
  
negative.code:           enddo
         enddo
  
         precip(:) = precip(:) + precip2(:)
  
      endif
  
      if ( USE_ToyChemistry ) then
  