c parameter(im=720,jm=360) c CFSv2 operational daily time series SST from seasonal forecast c 1. 4 runs / day from (00z, 06z, 12z, and 18z) c 2. 270 days time serires forecast starts from c the 2nd day (24hr) of the initial date (00UTC) c 3. grid resolution 0.5deg longitude by 0.5deg latitude c starts from (89.75S,-359.75E),eastward and northward c 4. undefined value=-99999.00 c 5. format: binary in big_endian c 6. file name: CFSv2.ocnsst.yyyymmdd.gr c yyyy=the year of initial date c mm=the month of initial date c dd=the day of initial date c 7. unit=K c real*4 sst00(im,jm) real*4 sst06(im,jm) real*4 sst12(im,jm) real*4 sst18(im,jm) c open(10, & file='CFSv2.ocnsst.20120615.gr', & form='unformatted', & access='direct',recl=im*jm*4) c open(20, & file='CFSv2.ocnsst.20120615.txt.sample', & form='formatted',access='sequential') c write(20,*) 'grid point (-179.75E,0.25N)' write(20,*) &'fcst days: member1(K): member2(K): member3(K): member4:(K)' write(20,*) &'fcst days: 00z: 06z: 12z: 18z' do ifcst=1,270 irec=(ifcst-1)*4+1 read(10,rec=irec) sst00 irec=(ifcst-1)*4+2 read(10,rec=irec) sst06 irec=(ifcst-1)*4+3 read(10,rec=irec) sst12 irec=(ifcst-1)*4+4 read(10,rec=irec) sst18 write(20,'(i10,4f12.2)') ifcst,sst00(361,181),sst06(361,181), & sst12(361,181),sst18(361,181) enddo c stop end cccccccccccccccccccccccccccccccc