52 |
itlo = mybxlo(mythid) |
itlo = mybxlo(mythid) |
53 |
ithi = mybxhi(mythid) |
ithi = mybxhi(mythid) |
54 |
|
|
55 |
#ifdef ALLOW_COST_VECTOR |
#if (defined (ALLOW_COST_VECTOR)) |
56 |
|
|
57 |
CALL COST_VECTOR (myThid) |
CALL COST_VECTOR (myThid) |
58 |
|
|
59 |
do bj = jtlo,jthi |
#elif (defined (ALLOW_COST_STATE_FINAL)) |
60 |
do bi = itlo,ithi |
|
61 |
do i = 1,sNx |
CALL COST_STATE_FINAL (myThid) |
|
print*,' --> objf_vector(i,bi,bj) = ', |
|
|
& objf_vector(i,bi,bj) |
|
|
end do |
|
|
end do |
|
|
end do |
|
62 |
|
|
63 |
#else /* ALLOW_COST_VECTOR undef */ |
#else /* ALLOW_COST_VECTOR undef */ |
64 |
|
|
85 |
enddo |
enddo |
86 |
enddo |
enddo |
87 |
|
|
88 |
print*,' fc = ', fc |
print*,' local fc = ', fc |
89 |
|
|
90 |
c-- Do global summation. |
c-- Do global summation. |
91 |
_GLOBAL_SUM_R8( fc , myThid ) |
_GLOBAL_SUM_R8( fc , myThid ) |
92 |
|
|
93 |
#endif /* ALLOW_COST_VECTOR */ |
#endif /* ALLOW_COST_VECTOR */ |
94 |
|
|
95 |
|
print*,' global fc = ', fc |
96 |
|
|
97 |
c-- set averaging freq. to zero to avoid re-write of |
c-- set averaging freq. to zero to avoid re-write of |
98 |
c-- averaged fields in reverse checkpointing loops |
c-- averaged fields in reverse checkpointing loops |
99 |
taveFreq = 0. |
taveFreq = 0. |