2009.01.13 Each climatology data file: 4 variables: fclm: 26-yr monthly mean climatology of the ensemble foreacst 15 members: f06-f20 ensemble: e15b oclm: 26-yr monthly mean climatology of observations fstd: standard deviation of hindcasts of all 15 members (26 yrs) ostd: standard deviation of observation Some CFS hindacst statistics are available at http://www.cpc.ncep.noaa.gov/products/people/wwang/cfs_metrics/ ---------------------------------------- Each data file contains 20 CFS individual runs, 3 different ensembles and the observation (total 24 variables). Each variable: 12 initial months (t=1-12) 9 montn foreacsts (z=1-9) 12x9 monthly mean data Initial conditions for 20 members: last 2 days of previous month (2) first 3 days of the initial month (3) 9-13 of the initial month (5) 19-23 of the initial month (5) last 2 days of the initial month (2) first 3 days of the following month (3) Totol: 20 For t=2, z=3: It is a 3-month foreacst with Feb initial conditions (20): Jan: 30, 31 Feb: 1, 2, 3, 9, 10, 11, 12, 13, 19, 20, 21, 22, 23, 27, 28, Mar: 1, 2, 3 for the same target month: May (z=3) Note that for Jan 30 and 31 initial conditions, it is actually a 4-month foreacst for May. -------------------------------------------- C Read SST PARAMETER ( NI = 180, NJ = 139, nv=9) real a(ni,nj,nv,24,12) open(11,file='SST2x11997.gr',form='unformatted', 1 access='direct', recl=ni*nj*4) c iv=1,24 c iv: 1-20 20 members c iv: 21-23 3 ensembles ensm20 ensm15a ensm15b c iv: 24 observations irec=0 do im=1,12 do iv=1,24 do k=1,nv irec=irec+1 read(11,rec=irec) ((a(i,j,k,iv,im),i=1,ni),j=1,nj) end do end do end do close (11) C------------------------------------ C Read Precip. PARAMETER ( NI = 192, NJ = 94, nv=9) real a(ni,nj,nv,24,12) open(11,file='PrateT621997.gr',form='unformatted', 1 access='direct', recl=ni*nj*4) irec=0 do im=1,12 do iv=1,24 do k=1,nv irec=irec+1 read(11,rec=irec) ((a(i,j,k,iv,im),i=1,ni),j=1,nj) end do end do end do close (11) C--------------------------------------- C Read Z500 or U200 C For height and wind, the data format is different from C SST and precip, which is indicated in the .ctl file C "options yrev sequential" C 1) It is 'sequential' C 2) 'yrev', data in meridional direction start from C the North Pole and end at the South Pole PARAMETER ( NI = 144, NJ = 73, nv=9) real a(ni,nj,nv,24,12) open(11,file='U2001997.gr', 1 form='unformatted',access='sequential') do im=1,12 do iv=1,24 do k=1,nv read(11) ((a(i,j,k,iv,im),i=1,ni),j=1,nj) end do end do end do close (11)