anchor.diff:                  call RRTMG_SWRAD(rthratensw=domain%tend%th_swrad,         &
   !                swupt, swuptc, swuptcln, swdnt, swdntc, swdntcln, &
   !                swupb, swupbc, swupbcln, swdnb, swdnbc, swdnbcln, &
  -!                      swupflx, swupflxc, swdnflx, swdnflxc,      &            
  +!                      swupflx, swupflxc, swdnflx, swdnflxc,      &
                       swcf = domain%shortwave_cloud_forcing%data_2d,        &
                       gsw = domain%shortwave%data_2d,                       &
  -                    xtime = 0., gmt = 0.,                                 &  ! not used  
  +                    xtime = 0., gmt = 0.,                                 &  ! not used
                       xlat = domain%latitude%data_2d,                       &  ! not used
                       xlong = domain%longitude%data_2d,                     &  ! not used
  -                    radt = 0., degrad = 0., declin = 0.,                  &  ! not used                        
  -                    coszr = domain%cosine_zenith_angle%data_2d,           & 
  +                    radt = 0., degrad = 0., declin = 0.,                  &  ! not used
  +                    coszr = domain%cosine_zenith_angle%data_2d,           &
                       julday = 0,                                           &  ! not used
                       solcon = solar_constant,                              &
                       albedo = albedo,                                      &
                       t3d = domain%temperature%data_3d,                     &
  -                    t8w = domain%temperature_interface%data_3d,           & 
  +                    t8w = domain%temperature_interface%data_3d,           &
                       tsk = domain%skin_temperature%data_2d,                &
                       p3d = domain%pressure%data_3d,                        &
                       p8w = domain%pressure_interface%data_3d,              &
  
positive.code:                     yr=domain%model_time%year,                            &
                      julian=domain%model_time%day_of_year()                &
                                                     )
                        
                  call RRTMG_LWRAD(rthratenlw=domain%tend%th_lwrad,                 &
  !                           lwupt, lwuptc, lwuptcln, lwdnt, lwdntc, lwdntcln,     &        !if lwupt defined, all MUST be defined
  !                           lwupb, lwupbc, lwupbcln, lwdnb, lwdnbc, lwdnbcln,     &
  
negative.code:                                     Dz_1d(k) = domain%dz_interface%data_3d(i,k,j)
                                      cf_1d(k) = cldfra(i,k,j)
                                  ENDDO
                                  gridkm = domain%dx/1000
                                  XLAND = domain%land_mask
                                  CALL cal_cldfra3(cf_1d, qv_1d, qc_1d, qi_1d, qs_1d, Dz_1d, &
                   &                              p_1d, t_1d, XLAND(i,j), gridkm,        &
                   &                              .false., 1.5, kms, kme)
  
                                  DO k = kts,kte
                                      ! qc, qi and qs are locally recalculated in cal_cldfra3 base on RH to account for subgrid clouds                                     qc(i,k,j) = qc_1d(k)
  