program read_sst c c Sample program provides read statements to read AVHRR MCSST data. c c The MCSST data are on a 1024 X 512 element grid with the origin c at 89.824 degrees North latitude and 0.176 degrees East longitude c (i.e., the grid ranges from 89.824 degrees North to 89.824 degrees c South in latitude and 0.176 degrees East to 359.824 degrees East c in longitude). There are 12 monthly data files for 1994. Note c that the 12 months do not entirely correspond to calendar months. c The days over which the data were collected are stored in the files. c c Declarations: integer*2 sst(1024,512) ! SST data in 1/10s degree C ! -9999 indicates land ! 9999 indicates no data integer*2 yr,mon ! year and month of data integer*2 smon,sday ! month and day of start of defined 'month' integer*2 emon,eday ! month and day of end of defined 'month' c Open & read a data file open(8,file='sst9401.dat',form='unformatted') read(8) yr,mon write(*,*) yr,mon ! Should be 94, 1 for first data file read(8) smon, sday, emon, eday write(*,*) smon,sday,emon,eday ! should have values of 1,6,2,2 ! for first data file. read(8) sst close(8) c Print values and convert SST to degrees Celsius. do 100 i=1,5 write(*,*) sst(i,256)/10. 100 continue c c Example Output from file: sst9401.dat c c yr,mon 94 1 c smon,sday,emon,eday 1 6 2 2 c c sst(1,256) = 28.30000 c sst(2,256) = 29.30000 c sst(3,256) = 28.30000 c sst(4,256) = 27.70000 c sst(5,256) = 27.40000 c stop end