program rdblend c******************************************************************** c Reads ASCII formatted blended CZCS/In situ data. Data are c seasonal climatologies of chlorophyll, units mg/m3. c 13 = Winter (Jan-Mar) c 14 = Spring (Apr-Jun) c 15 = Summer (Jul-Sep) c 16 = Autumn (Oct-Dec) c Data begin at -180 longitude, -90 latitude c Land flag = 999.999 c Missing data flag = -1.0 c Results from this program should agree with the file chk.dat. c Reference: Gregg, W.W. and M.E. Conkright, 2000. Global seasonal c climatologies of ocean chlorophyll: Blending In situ and satellite c data for the CZCS era. Journal of Geophysical Reasearch, in press. c********************************************************************* c Parameters: c igrd: longitude grids c jgrd: latitude grids c Character: c cfle: input file name c Real: c cblend: blended chlorophyll data c Data: c cfle: names of input files c********************************************************************* parameter(igrd=360,jgrd=180) character*50 cfle(4),testfile real cblend(igrd,jgrd) data cfle /'blend13.txt','blend14.txt','blend15.txt', *'blend16.txt'/ c********************************************************************* c Variables: c isea: identifies season c iprt: longitude grid number to print c jprt: latitude grid number to print c********************************************************************* c Read data do n = 1,4 open(4,file=cfle(n),status='old',form='formatted') do j = 1,jgrd read(4,10)(cblend(i,j),i=1,igrd) enddo close(4) isea = n + 12 iprt = 140 jprt = 40 write(6,*)'Season,i,j,cblend = ',isea,iprt,jprt,cblend(iprt,jprt) enddo 10 format(360f7.3) end