anchor.diff:                      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,     &
  
positive.code:         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
  
          albedo=0.17
          F_QI=.false.
          F_QC=.false.
  
negative.code:         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)
                  domain%cosine_zenith_angle%data_2d(its:ite,j)=sin(solar_elevation(its:ite))
              enddo
  
              if (last_model_time==-999) then
                  last_model_time = domain%model_time%seconds()-update_interval
  