/[MITgcm]/MITgcm/pkg/ecco/cost_ssh_mean.F
ViewVC logotype

Diff of /MITgcm/pkg/ecco/cost_ssh_mean.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.4 by heimbach, Thu May 5 23:47:57 2005 UTC revision 1.5 by heimbach, Fri May 27 22:10:26 2005 UTC
# Line 41  c     == global variables == Line 41  c     == global variables ==
41  #include "PARAMS.h"  #include "PARAMS.h"
42    
43  #include "ecco_cost.h"  #include "ecco_cost.h"
44    #include "optim.h"
45    
46  #ifdef ALLOW_EGM96_ERROR_COV  #ifdef ALLOW_EGM96_ERROR_COV
47  # if (defined (ALLOW_USE_MPI) || defined (ALWAYS_USE_MPI))  # if (defined (ALLOW_USE_MPI) || defined (ALWAYS_USE_MPI))
# Line 68  c     == local variables == Line 69  c     == local variables ==
69        integer jmin,jmax        integer jmin,jmax
70        integer imin,imax        integer imin,imax
71    
72    cph(
73          _RL diagnosfld(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
74    cph)
75    
76  #ifdef ALLOW_EGM96_ERROR_COV  #ifdef ALLOW_EGM96_ERROR_COV
77        _RL cphi        _RL cphi
78        _RL xmean        _RL xmean
# Line 232  c--         --> no weighting by cosphi. Line 237  c--         --> no weighting by cosphi.
237       &                    * wp(i,j,bi,bj)*tpmeanmask(i,j,bi,bj)       &                    * wp(i,j,bi,bj)*tpmeanmask(i,j,bi,bj)
238                if ( wp(i,j,bi,bj)*tpmeanmask(i,j,bi,bj) .ne. 0. )                if ( wp(i,j,bi,bj)*tpmeanmask(i,j,bi,bj) .ne. 0. )
239       &             num_hmean = num_hmean + 1. _d 0       &             num_hmean = num_hmean + 1. _d 0
240                  diagnosfld(i,j,bi,bj) = diff*diff
241         &                    * wp(i,j,bi,bj)*tpmeanmask(i,j,bi,bj)
242              enddo              enddo
243            enddo            enddo
244          enddo          enddo
245        enddo        enddo
246    
247    cph(
248          CALL WRITE_FLD_XY_RL( 'DiagnosSSHmean', ' ', diagnosfld,
249         &                           optimcycle, mythid )
250    cph)
251    
252  c--   Do the global summation.  c--   Do the global summation.
253        _GLOBAL_SUM_R8( sumc, mythid )        _GLOBAL_SUM_R8( sumc, mythid )
254        _GLOBAL_SUM_R8( num_hmean, mythid )        _GLOBAL_SUM_R8( num_hmean, mythid )
255        costmean = sumc        costmean = sumc
256    
257  #endif  #endif
258    
259        end        end
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22