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

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

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

revision 1.4 by heimbach, Mon Mar 28 19:40:59 2005 UTC revision 1.5 by heimbach, Mon Mar 28 23:49:49 2005 UTC
# Line 159  cph Line 159  cph
159  cph               if ( rC(K) .GT. -1000. ) then  cph               if ( rC(K) .GT. -1000. ) then
160  cph)  cph)
161                do j = jmin,jmax                do j = jmin,jmax
162                  do i = imin,imax                                 do i = imin,imax                
163                     if (argosobs(i,j,k,bi,bj) .ge. spval .and.                  if ( (argosobs(i,j,k,bi,bj) .ge. spval) .and.
164       $                 argosobs(i,j,k,bi,bj) .lt. spmax .and.       $              (argosobs(i,j,k,bi,bj) .lt. spmax) .and.
165       $                  _hFacC(i,j,13,bi,bj) .ne. 0. .and.       $              (_hFacC(i,j,13,bi,bj) .ne. 0.) .and.
166       $                  _hFacC(i,j,k,bi,bj) .ne. 0.) then       $              (_hFacC(i,j,k,bi,bj) .ne. 0.) ) then
167                        fctile_argos = fctile_argos +                   fctile_argos = fctile_argos +
168       &           wsalt2(i,j,k,bi,bj)*cosphi(i,j,bi,bj)*       &             wsalt2(i,j,k,bi,bj)*cosphi(i,j,bi,bj)*
169       &                 (sbar(i,j,k,bi,bj)-argosobs(i,j,k,bi,bj))*       &             (sbar(i,j,k,bi,bj)-argosobs(i,j,k,bi,bj))*
170       &                 (sbar(i,j,k,bi,bj)-argosobs(i,j,k,bi,bj))       &             (sbar(i,j,k,bi,bj)-argosobs(i,j,k,bi,bj))
171                    endif                                     if ( wsalt2(i,j,k,bi,bj)*cosphi(i,j,bi,bj) .ne. 0. )
172                  enddo       &             num_argos(bi,bj) = num_argos(bi,bj) + 1. _d 0
173                    endif                  
174                   enddo
175                enddo                enddo
176              enddo              enddo
177  c--         End of loop over layers.  c--         End of loop over layers.
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22