program readmon dimension sst(360,180),ice(360,180),err(360,180) dimension iTAGLS(360,180) character fname*64 c---------------------------------------------------------------------------- c this program reads the version 2 weekly ASCII oi SST data c c VARIABLE LIST c ------------- c sst - sea surface temperature array (deg C) c err - normalized error variance c ice - ice concentration array (%) (0-100, >100 = land or coast) c iyrst - year of start date of analysis c imst - month of start date of analysis c idst - day of start date of analysis c iyrnd - year of end date of analysis c imnd - month of end date of analysis c idnd - day of end date of analysis c ndays - number of days in analysis (start date thru enddate) c indx - analysis version for reference c itagls - land/sea tag array (0=land, 1=water) c c c GRID ORIENTATION c ---------------- c The first gridbox of each array is centered on 0.5E, 89.5S. The points c move eastward to 359.5E, then northward to 89.5N. c c So, the orientation of the sst and cice data arrays is : c c sst(1,1) = 0.5E, 89.5S c sst(2,1) = 1.5E, 89.5S c sst(360,1) = 359.5E, 89.5S c sst(1,2) = 0.5E, 88.5S c sst(1,180) = 0.5E, 89.5N c sst(360,180) = 359.5E, 89.5N c c c---------------------------------------------------------------------------- c uncomment out the lines below if you want to read the land-sea mask c c c Read in land sea tags (1 for ocean; 0 for land) c open (11,file='land.sea.mask.v2.asc',form='formatted') c READ (11,'(360i1)') iTAGLS c close(11) write(*,'("Enter Input File Name : ",$)') read(*,'(a)') fname open(10,file=fname) 10 read(10,6,end=99) iyrst,imst,idst,iyrnd,imnd,idnd,ndays,indx 6 format(7I5,i10) | read(10,8) sst read(10,7) err read(10,9) ice write(*,11) iyrst,imst,idst,iyrnd,imnd,idnd,ndays,indx 11 format(7i6,i11) 8 format(20f4.1) 9 format(26i3) 7 FORMAT(16f5.3) goto 10 99 continue end