c parameter(im=360,jm=180) c CFS operational daily forecast SST c 1. 4 runs / day c 2. 270 days time serires forecast starts from c the 2nd day (24hr) of the initial date (00UTC) c 3. grid resolution 1.0deg longitude by 1.000deg latitude c starts from (89.5S,0.5E),eastward and northward c 4. undefined value=-99999.00 c 5. format: binary in big_endian c 6. file name: csstf.yyyyddmm.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 sst1(im,jm) real*4 sst2(im,jm) real*4 sst3(im,jm) real*4 sst4(im,jm) c open(10, & file='csstf.20101020.gr', & form='unformatted', & access='direct',recl=im*jm*4) c open(20, & file='csstf.20101020.txt.sample', & form='formatted',access='sequential') c write(20,*) &'fcst days: member1(K): member2(K): member3(K): member4:(K)' do ifcst=1,270 irec=(ifcst-1)*4+1 read(10,rec=irec) sst1 irec=(ifcst-1)*4+2 read(10,rec=irec) sst2 irec=(ifcst-1)*4+3 read(10,rec=irec) sst3 irec=(ifcst-1)*4+4 read(10,rec=irec) sst4 write(20,'(i10,4f12.2)') ifcst,sst1(180,90),sst2(180,90), & sst3(180,90),sst4(180,90) enddo c stop end cccccccccccccccccccccccccccccccc