write out a file containing the summed footprint

prep_regions/prep_regions.f90

@@ -130,6 +130,15 @@ program prep_regions
   varname = 'regions'
   call write_ncdf(files, filename, varname, nbox_xy)
 
+  ! write the average fooptrint for all receptors
+  open(50,file=trim(files%path_output)//'all_footprint_average.txt')
+  do ix=1,nxregrid
+      do jy=1,nyregrid
+          write(50,*)  reg_lon(ix), reg_lat(jy), surfinf(ix,jy)
+       end do
+  end do
+  close(50)
+
   ! write lsm 
   allocate( mask(nxregrid,nyregrid), stat=ierr)
   mask=lsm