anchor.diff:      fe (:,:)   = 0.0_RP
       fq (:,:,:) = 0.0_RP
   
  -    precip(:) = 0.0_RP
  +    precip (:) = 0.0_RP
  +    precip2(:) = 0.0_RP
   
       if ( USE_Kessler ) then
          do ij = 1, ijdim
  
positive.code:     logical :: SET_DCMIP2016_13 = .false.
      logical :: SET_DCMIP2016_DRY = .false.
      logical :: SET_DCMIP2016_LSC = .false. !large scale condensation
  
      namelist /FORCING_DCMIP_PARAM/ &
         SET_RJ2012,          &
  
negative.code:        USE_SimpleMicrophys = .true.
         SM_Latdepend_SST    = .true.
         SM_LargeScaleCond   = .true.
         !SM_PBL_Bryan        = .false.
         USE_Kessler         = .false.
         USE_ToyChemistry    = .false.
      elseif( SET_DCMIP2016_11 ) then
         write(ADM_LOG_FID,*) '*** Force setting of DCMIP2016 Case 1-1 (Moist baroclinic wave with terminator chemistry)'
         USE_Kessler         = .true.
         SM_LargeScaleCond   = .false.
         !SM_PBL_Bryan        = .false.
         USE_ToyChemistry    = .true.
         if ( SET_DCMIP2016_NOSST ) then
            write(ADM_LOG_FID,*) '*** Force setting of DCMIP2016: Only Precipitation'
            USE_SimpleMicrophys = .false.
            SM_Latdepend_SST    = .false.
         else
            USE_SimpleMicrophys = .true.
            SM_Latdepend_SST    = .true.
         endif
         if ( SET_DCMIP2016_LSC ) then
            write(ADM_LOG_FID,*) '*** Force setting of DCMIP2016: USE SM_LargeScaleCond'
            USE_Kessler      = .false.
  