anchor.diff:          allocate(day_frac(ims:ime))
           allocate(solar_elevation(ims:ime))
           allocate(albedo(ims:ime,jms:jme))
  -        allocate(cldfra3d(ims:ime,kms:kme,jms:jme))
   
           ! Note, need to link NoahMP to update albedo
  -        
  +
           qc = 0
           qi = 0
           qs = 0
   
  +        cldfra=0
           albedo=0.17
           F_QI=.false.
           F_QC=.false.
  
positive.code:         integer :: its, ite, jts, jte, kts, kte
          integer :: ids, ide, jds, jde, kds, kde
  
          real, dimension(:,:,:,:), pointer :: tauaer_sw=>null(), ssaaer_sw=>null(), asyaer_sw=>null() 
          real, allocatable :: cldfra3d(:,:,:)
          real, allocatable :: day_frac(:), solar_elevation(:)
          real, allocatable:: albedo(:,:)
          integer :: j
  
negative.code:         endif
          update_interval=options%rad_options%update_interval_rrtmg ! 30 min, 1800 s   600 ! 10 min (600 s)
          last_model_time=-999
  
      end subroutine radiation_init
  
  
  