anchor.diff:          if (options%physics%radiation==kRA_RRTMG) then
               do j = jms,jme
                   solar_elevation  = calc_solar_elevation(date=domain%model_time, lon=domain%longitude%data_2d, &
  -                                j=j, ims=ims,ime=ime,jms=jms,jme=jme,its=its,ite=ite,day_frac=day_frac)                
  +                                j=j, ims=ims,ime=ime,jms=jms,jme=jme,its=its,ite=ite,day_frac=day_frac)
                   domain%cosine_zenith_angle%data_2d(its:ite,j)=sin(solar_elevation(its:ite))
  -            enddo  
  +            enddo
   
               if (last_model_time==-999) then
                   last_model_time = domain%model_time%seconds()-update_interval
  
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:         allocate(day_frac(ims:ime))
          allocate(solar_elevation(ims:ime))
          allocate(albedo(ims:ime,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.
  