anchor.diff:  !                   aer_ra_feedback,                                      &
   !jdfcz              progn,prescribe,                                      &
                       calc_clean_atm_diag=0,                                &
  -                    qndrop3d=domain%cloud_number%data_3d,                 &
  +!                    qndrop3d=domain%cloud_number%data_3d,                 &
                       f_qndrop=f_qndrop,                                    & !czhao
                       mp_physics=0,                                         & !wang 2014/12
                       ids=ids, ide=ide, jds=jds, jde=jde, kds=kds, kde=kde, &
  
positive.code:         F_QR=associated(domain%rain_mass%data_3d )
          F_QS=associated(domain%snow_mass%data_3d )
          F_QV=associated(domain%water_vapor%data_3d )
          F_QG=associated(domain%graupel_mass%data_3d )
          F_QNDROP=associated(domain%cloud_number%data_3d)
  
          if (F_QC) qc(:,:,:) = domain%cloud_water_mass%data_3d
  
negative.code:         where(domain%cloud_water_mass%data_3d+domain%snow_mass%data_3d+domain%cloud_ice_mass%data_3d>0) cldfra3d=1
  
          ! Note, need to link NoahMP to update albedo
          
          qc = 0
          qi = 0
          qs = 0
  
          albedo=0.17
          F_QI=.false.
          F_QC=.false.
  