program z_cntrl
!
!**********************************************************************
! *              .      .    .
! * PROGRAM:    z_cntrl
! * PRGMMR:     M. ZUPANSKI     ORG: CIRA/CSU   DATE: 2003-09-03
! *
! * ABSTRACT:   INITIALIZE CNTRL VARIABLE IN ENSEMBLE SPACE 
! *             Assign zero value to initial z_cntrl file
! *
! *
! * PROGRAM LOG:
! *
! * 09/03/2003    ..... M. ZUPANSKI:
! * 10/14/2003    ..... M. ZUPANSKI:
! *
! **********************************************************************
! OUTPUT:    z(i)     = control variable (ensemble subspace)
!--------------------------------------------------------------------
!-----
  integer,parameter::icntrl=21       ! cntrl vrbl (ensemble space)
!-----

 integer :: NENS
 integer :: NENS_START
 real,allocatable::z(:)

!
!      DECLARE NAMELIST
!
      NAMELIST /ENSEMBLE_SIZE/ NENS_START,NENS

!=========================================================
!==============start calculation===================
        write(*,*) "Start z_cntrl"

!-- read ensemble size

       write(*,*) "z_cntrl: start reading input 15"

       REWIND 15
       READ(15,ENSEMBLE_SIZE)

       write(*,*) "z_cntrl: NENS_START, NENS=",NENS_START,NENS

!    write the cntrl vrbl in ensemble space (z)

        allocate(z(NENS_START:NENS))
        z(:)=0.0

      CLOSE(icntrl)
      OPEN(UNIT=icntrl,FILE='zcntrl_file',FORM='UNFORMATTED',IOSTAT=IER)
      IF(IER.NE.0) WRITE(*,*)' icntrl OPEN UNIT ERROR IER=',IER
      WRITE(icntrl) z
      CLOSE(icntrl,status='KEEP')

        deallocate(z)

!====================================================================

 stop
 end