C $Header: /u/gcmpack/MITgcm/pkg/cost/cost_state_final.F,v 1.7 2012/08/10 19:36:02 jmc Exp $ C $Name: $ #include "COST_OPTIONS.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) DO j=1,sNy DO i=1,sNx DO k=1,nr objf_state_final(i,j,bi,bj,0*Nr+k) = theta(i,j,k,bi,bj) objf_state_final(i,j,bi,bj,1*Nr+k) = salt(i,j,k,bi,bj) objf_state_final(i,j,bi,bj,2*Nr+k) = uvel(i,j,k,bi,bj) objf_state_final(i,j,bi,bj,3*Nr+k) = vvel(i,j,k,bi,bj) END
DO objf_state_final(i,j,bi,bj,4*Nr+1) = etan(i,j,bi,bj) END
DO END
DO END
DO END
DO cph fc = objf_state_final(45,4,1,1,1) cph _GLOBAL_SUM_RL( fc , myThid ) cph print *, 'fc for admtlm test at 83,33 : ', fc #endif end