PROGRAM convert
! -------------------------------------------
! read in Dudley Chelton's ascii data file 
! convert to 1x1 deg global gridded data will -1 for missing C and L_r
! output format is obvious from WRITE statement below 
!
! Coded by Rob Scott
! do with it whatever you want.
! -------------------------------------------
REAL, DIMENSION(360,180):: lr, c, lat, lon

REAL :: raw_lat, raw_lon, raw_c, raw_lr

OPEN(UNIT=10,FILE='rossby_radius.dat',STATUS='OLD')
OPEN(UNIT=20,FILE='rossby_radius_map.dat',STATUS='NEW')

lr = -1
c = -1
lon = 999
lat = 999
DO n = 1,31927
   READ(10,*) raw_lat, raw_lon, raw_c, raw_lr
   grid_i = CEILING(raw_lon) 
   grid_j = CEILING(raw_lat) + 90
   lr(grid_i,grid_j) = raw_lr
   c(grid_i,grid_j) = raw_c
   lat(grid_i,grid_j) = raw_lat
   lon(grid_i,grid_j) = raw_lon
END DO

DO j=1,180
   DO i=1,360
      WRITE(UNIT=20,'(2(1x,I4),4(1x,E12.5))') j,i,lat(i,j), lon(i,j), c(i,j), lr(i,j)
   END DO
END DO

STOP
END PROGRAM convert