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

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

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

revision 1.2 by heimbach, Mon Oct 11 16:38:53 2004 UTC revision 1.3 by heimbach, Mon Mar 28 23:49:49 2005 UTC
# Line 186  c--       Loop over the model layers Line 186  c--       Loop over the model layers
186  c--         Compute model misfit and cost function term for  c--         Compute model misfit and cost function term for
187  c           the vertical velovity field. The error is 1e-4 m/s.  c           the vertical velovity field. The error is 1e-4 m/s.
188              do j = jmin,jmax              do j = jmin,jmax
189                do i = imin,imax               do i = imin,imax
190                    if (_hFacC(i,j,k,bi,bj) .ne. 0.) then                if (_hFacC(i,j,k,bi,bj) .ne. 0.) then
191                       fctilew = fctilew +                       fctilew = fctilew +
192       &                  (2.5e11*cosphi(i,j,bi,bj)*       &                  (2.5e11*cosphi(i,j,bi,bj)*
193       &          (wlmean(i,j,k,bi,bj)/nl - wfmean(i,j,k,bi,bj)/nf)*       &          (wlmean(i,j,k,bi,bj)/nl - wfmean(i,j,k,bi,bj)/nf)*
194       &          (wlmean(i,j,k,bi,bj)/nl - wfmean(i,j,k,bi,bj)/nf))       &          (wlmean(i,j,k,bi,bj)/nl - wfmean(i,j,k,bi,bj)/nf))
195                    endif                       if ( cosphi(i,j,bi,bj) .ne. 0. )
196                enddo       &                 num_wdrift(bi,bj) = num_wdrift(bi,bj) + 1. _d 0
197                  endif
198                 enddo
199              enddo              enddo
200    
201            enddo            enddo
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22