--- MITgcm/pkg/cost/cost_test.F 2001/11/20 21:06:48 1.3 +++ MITgcm/pkg/cost/cost_test.F 2002/03/24 02:14:21 1.4 @@ -1,4 +1,4 @@ -C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/cost/cost_test.F,v 1.3 2001/11/20 21:06:48 heimbach Exp $ +C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/cost/cost_test.F,v 1.4 2002/03/24 02:14:21 heimbach Exp $ #include "CPP_OPTIONS.h" @@ -42,8 +42,8 @@ itlo = mybxlo(mythid) ithi = mybxhi(mythid) - iLocOut = 6 - jLocOut = 35 + iLocOut = 80 + jLocOut = 30 kLocOut = 1 ce some reference temperature @@ -64,6 +64,10 @@ write(*,'(a,3(x,i4),a,4(x,i4))') & 'COST ',ig,jg,kLocOut,' TILE ',i,j,bi,bj objf_test(bi,bj) = theta(i,j,kLocOut,bi,bj) +cph( + print *, 'ph-cost ', ig, jg, kLocOut, + & theta(i,j,kLocOut,bi,bj) +cph) endif end do