anchor.diff:  enddo
   enddo
   
  -write(dayfile,'(a,a,a)') basedir,'/PROCESSING/overpass_corr/TRAD1_',tile_num
  -open(10,file=dayfile,form='unformatted',access='direct',recl=dx*dy*4)
  +!write(dayfile,'(a,a,a)') basedir,'/PROCESSING/overpass_corr/TRAD1_',tile_num
  +open(10,file=outdayfile,form='unformatted',access='direct',recl=dx*dy*4)
   write(10,rec=1) day_corr
   close(10)
   
  
positive.code: rday=int(((offset+nweek1)*7)+1)
  write(6,*) yyyyddd, rday
  
  write(dztimefile,'(a,a,a,a)') basedir,'/PROCESSING/overpass_corr/CURRENT_DAY_ZTIME_',tile_num,'.dat'
  open(10,file=dztimefile,form='unformatted',access='direct',recl=dx*dy*4)
  read(10,rec=1) dztime
  close(10)
  
  write(dayfile,'(a,a,a,a)') basedir,'/PROCESSING/overpass_corr/RAW_TRAD1_',tile_num,'.dat'
  write(avgfile,'(a,a,a,a)') basedir,'/PROCESSING/overpass_corr/CURRENT_DTRAD_AVG_',tile_num,'.dat'
  
  open(10,file=dayfile,form='unformatted',access='direct',recl=dx*dy*4)
  open(14,file=avgfile,form='unformatted',access='direct',recl=dx*dy*4)
  
negative.code: !character (len=*), allocatable :: basedir
  character(len=4) :: tile_num
  integer, parameter :: dx=3750, dy=3750
  character(len=400) :: dayfile, outdayfile, dtimefile, ntimefile
  character(len=400) :: avgfile, dztimefile, nztimefile
  character(len=255) :: arg4, basedir
  character(len=4) :: arg1, arg3
  character(len=3) :: arg2
  