--- MITgcm/pkg/cost/cost_final.F 2003/01/21 19:20:52 1.5 +++ MITgcm/pkg/cost/cost_final.F 2003/10/02 21:34:45 1.6 @@ -1,4 +1,4 @@ -C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/cost/cost_final.F,v 1.5 2003/01/21 19:20:52 heimbach Exp $ +C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/cost/cost_final.F,v 1.6 2003/10/02 21:34:45 heimbach Exp $ #include "COST_CPPOPTIONS.h" @@ -52,18 +52,13 @@ itlo = mybxlo(mythid) ithi = mybxhi(mythid) -#ifdef ALLOW_COST_VECTOR +#if (defined (ALLOW_COST_VECTOR)) CALL COST_VECTOR (myThid) - do bj = jtlo,jthi - do bi = itlo,ithi - do i = 1,sNx - print*,' --> objf_vector(i,bi,bj) = ', - & objf_vector(i,bi,bj) - end do - end do - end do +#elif (defined (ALLOW_COST_STATE_FINAL)) + + CALL COST_STATE_FINAL (myThid) #else /* ALLOW_COST_VECTOR undef */ @@ -90,13 +85,15 @@ enddo enddo - print*,' fc = ', fc + print*,' local fc = ', fc c-- Do global summation. _GLOBAL_SUM_R8( fc , myThid ) #endif /* ALLOW_COST_VECTOR */ + print*,' global fc = ', fc + c-- set averaging freq. to zero to avoid re-write of c-- averaged fields in reverse checkpointing loops taveFreq = 0.