C
C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/ctrl/ctrl_init_ctrlvar.F,v 1.3 2004/11/11 06:28:33 heimbach Exp $
C $Name:  $

#include "CTRL_CPPOPTIONS.h"

      subroutine ctrl_init_ctrlvar( 
     &       xx_fname
     &     , ivarindex
     &     , ncvarindex_loc
     &     , ncvarrecs_loc
     &     , ncvarrecstart_loc
     &     , ncvarrecsend_loc
     &     , ncvarxmax_loc
     &     , ncvarymax_loc
     &     , ncvarnrmax_loc
     &     , ncvargrd_loc
     &     , whichxyz
     &     , mythid
     &     )

c     ==================================================================
c     SUBROUTINE ctrl_init_ctrlvar
c     ==================================================================
c
c     started: heimbach@mit.edu 28-Oct-2003
c
c     ==================================================================
c     SUBROUTINE ctrl_init_ctrlvar
c     ==================================================================

      implicit none

c     == global variables ==

#include "EEPARAMS.h"
#include "SIZE.h"
#include "PARAMS.h"
#include "GRID.h"
#include "ctrl.h"

#ifdef ALLOW_ECCO_OPTIMIZATION
# include "optim.h"
#endif

c     == routine arguments ==
      character*(MAX_LEN_FNAM) xx_fname
      integer ivarindex
      integer ncvarindex_loc
      integer ncvarrecs_loc
      integer ncvarrecstart_loc
      integer ncvarrecsend_loc
      integer ncvarxmax_loc
      integer ncvarymax_loc
      integer ncvarnrmax_loc
      character ncvargrd_loc*1
      character whichxyz*2
      integer mythid

c     == local variables ==
#ifndef ALLOW_ECCO_OPTIMIZATION
      integer optimcycle
#endif
      character*( 80) fname(2)
      logical  exst

#ifndef ALLOW_ECCO_OPTIMIZATION
      optimcycle = 0
#endif

      _BEGIN_MASTER( mythid )

      ncvarindex(ivarindex)    = ncvarindex_loc
      ncvarrecs (ivarindex)    = ncvarrecs_loc
      ncvarrecstart(ivarindex) = ncvarrecstart_loc
      ncvarrecsend(ivarindex)  = ncvarrecsend_loc
      ncvarxmax (ivarindex)    = ncvarxmax_loc
      ncvarymax (ivarindex)    = ncvarymax_loc
      ncvarnrmax(ivarindex)    = ncvarnrmax_loc
      ncvargrd  (ivarindex)    = ncvargrd_loc

      call ctrl_set_fname( xx_fname, fname, mythid )

      inquire( file='costfinal', exist=exst )
      IF ( .NOT. exst) THEN

         if ( whichxyz .EQ. '3d') then
            call ctrl_set_globfld_xyz( fname(2), ivarindex, mythid)
            if ( doInitXX .AND. optimcycle.eq.0 )
     &           call ctrl_set_globfld_xyz( fname(1), ivarindex, mythid)
         else if ( whichxyz .EQ. 'xy') then
            call ctrl_set_globfld_xy( fname(2), ivarindex, mythid)
            if ( doInitXX .AND. optimcycle.eq.0 )
     &           call ctrl_set_globfld_xy( fname(1), ivarindex, mythid )
         else if ( whichxyz .EQ. 'xz') then
            call ctrl_set_globfld_xz( fname(2), ivarindex, mythid)
            if ( doInitXX .AND. optimcycle.eq.0 )
     &           call ctrl_set_globfld_xz( fname(1), ivarindex, mythid )
         else if ( whichxyz .EQ. 'yz') then
            call ctrl_set_globfld_yz( fname(2), ivarindex, mythid)
            if ( doInitXX .AND. optimcycle.eq.0 )
     &           call ctrl_set_globfld_yz( fname(1), ivarindex, mythid )
         else
            STOP 'whichxyz option not implemented'
         end if

      ENDIF

      _END_MASTER( mythid )

      return
      end