#include "COST_CPPOPTIONS.h"

      subroutine COST_STATE_FINAL( myThid )
C     /==========================================================\
C     | subroutine cost_state_final                              |
C     | o This routine assigns final T,S to cost function        |
C     \==========================================================/
       implicit none

C     == Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "DYNVARS.h"

#include "cost.h"

C     ======== Routine arguments ======================
C     myThid - Thread number for this instance of the routine.
      integer myThid

#ifdef ALLOW_COST_STATE_FINAL
C     ========= Local variables =========================
      integer i, j, k
      integer bi, bj

      DO bj=myByLo(myThid),myByHi(myThid)
       DO bi=myBxLo(myThid),myBxHi(myThid)
cph        DO k=1,nr
         DO j=1,sNy
          DO i=1,sNx
           objf_state_final(i,j,bi,bj,1) = theta(i,j,1,bi,bj)
           objf_state_final(i,j,bi,bj,2) = salt(i,j,1,bi,bj)
cph           objf_state_final(i,j,k,bi,bj,1) = theta(i,j,k,bi,bj)
cph           objf_state_final(i,j,k,bi,bj,2) = salt(i,j,k,bi,bj)
          END
 DO
         END
 DO
cph        END DO
       END
 DO
      END
 DO

cph      fc = objf_state_final(45,4,1,1,1)
cph      _GLOBAL_SUM_R8( fc , myThid )
cph      print *, 'fc for admtlm test at 83,33 : ', fc

#endif
      
      end