parameter (nyr=109)! precip three month mean!cho
      dimension help(144,73)
      dimension precip3(144,73,55:nyr,12)! observed field
      dimension zmember(144,73,56:nyr,12,12)! single member forecast
      dimension zmemberave(144,73,56:nyr,12,12)! the ens mean forecast
      dimension zmemberavea(144,73,56:nyr,12,12)! the ens mean forecast-climate

      dimension caclim(144,73,12,12)
       character*28 pname
c        pname='/tc2/yangbo/lukehenew_Dec06/'
c      open(29,file='/ptmp/wx51hd/test/prate3all_total'! obs
c     *,form='unformatted')
      open(61,file='caglobprecip3_m1',
     *form='unformatted')!input
      open(62,file='caglobprecip3_m2',
     *form='unformatted')!input
      open(63,file='caglobprecip3_m3',
     *form='unformatted')!input
      open(64,file='caglobprecip3_m4',
     *form='unformatted')!input
      open(65,file='caglobprecip3_m5',
     *form='unformatted')!input
      open(66,file='caglobprecip3_m6',
     *form='unformatted')!input
      open(67,file='caglobprecip3_m7',
     *form='unformatted')!input
      open(68,file='caglobprecip3_m8',
     *form='unformatted')!input
      open(69,file='caglobprecip3_m9',
     *form='unformatted')!input
      open(70,file='caglobprecip3_m10',
     *form='unformatted')!input
      open(71,file='caglobprecip3_m11',
     *form='unformatted')!input
      open(72,file='caglobprecip3_m12',
     *form='unformatted')!input


c      mfinal=8! cho; 8=JJA, last month of data
      mfinal=8! cho; 6=AMJ last month of data
      zmember=999999.
      zmemberave=0.

      do memb=1,12
          write(*,*) 'memb=',memb
      do m=5,5! initial seasons djf=1, etc
          write(*,*)'m=',m
       nyri=nyr
       if (m.gt.mfinal-1)nyri=nyr-1! 
       do lead=1,12! twelve leads by CPC definition
       do jear=56,nyri! 1956 to 2005/6 (years for initial state=season)
         read(60+memb)help
         if (memb.eq.12)zmember(:,:,jear,m,lead)=help(:,:)
         zmemberave(:,:,jear,m,lead)=zmemberave(:,:,jear,m,lead)
     *                               +help(:,:)/12.
      enddo
      enddo
      enddo
      enddo  
c
c* this concludes the reading of the frcts.
c* For m=1(2), DJF(JFM) is the last observed season, 
c* and the lead 1 seasonal forecast is for AMJ(MJJ). 
c* The initial years are 1956 (all cases) to 2005 (m=1-10) or 2006 (m=11 and 12)
c*
      calim=0
      do lead=1,12
c      do m=7,7
      do m=6,6
      do jear=71,100
      caclim(:,:,m,lead)=caclim(:,:,m,lead)+
     *zmemberave(:,:,jear,m,lead)/30.! lead dependent ca climo
      enddo
      enddo
      enddo

     
c      do m=7,7! initial seasons djf=1, etc
      do m=5,5! initial seasons djf=1, etc
       nyri=nyr
       if (m.gt.mfinal-1)nyri=nyr-1! 
       do lead=1,12! twelve leads by CPC definition
       do jear=56,nyri! 1956 to 2005/6 (years for initial state=season)
         zmemberavea(:,:,jear,m,lead)=zmemberave(:,:,jear,m,lead)
     *   -caclim(:,:,m,lead)                            
      enddo
      enddo
      enddo
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
        open(333,file='caforecast.dat',form='unformatted') 
c       jear=108
       jear=109
c       m=7    ! m=7(JJA),lead1=OND,lead2=NDJ,lead3=DJF,lead4=JFM,lead5=FMA,lead6=MAM
       m=5    ! m=7(JJA),lead1=OND,lead2=NDJ,lead3=DJF,lead4=JFM,lead5=FMA,lead6=MAM
        do lead=1,6!m=10(SON),lead1=JFM,lead2=FMA,lead3=MAM,lead4=AMJ,lead5=MJJ,lead6=JJA
        write(333) ((zmemberavea(i,j,jear,m,lead),i=1,144),j=1,73)
        write(333) ((zmemberave(i,j,jear,m,lead),i=1,144),j=1,73)

          write(*,*) ((zmemberavea(i,j,jear,m,lead),i=100,140),j=33,38)
        enddo
        close(33)



      stop
      end