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

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

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

revision 1.7 by heimbach, Mon May 14 22:06:55 2007 UTC revision 1.8 by gforget, Wed Jun 20 19:37:35 2007 UTC
# Line 116  c--         Determine the weights to be Line 116  c--         Determine the weights to be
116                do i = imin,imax                do i = imin,imax
117                  if (_hFacC(i,j,k,bi,bj) .ne. 0.) then                  if (_hFacC(i,j,k,bi,bj) .ne. 0.) then
118                    tmpx = tmpfld3d(i,j,k,bi,bj)                    tmpx = tmpfld3d(i,j,k,bi,bj)
119                    if ( ABS(R_low(i,j,bi,bj)) .LT. 100. )  #ifndef ALLOW_SMOOTH_CORREL3D
120                      if ( ABS(R_low(i,j,bi,bj)) .LT. 100. )
121       &              tmpx = tmpx*ABS(R_low(i,j,bi,bj))/100.       &              tmpx = tmpx*ABS(R_low(i,j,bi,bj))/100.
122                    fctile = fctile                    fctile = fctile
123       &                 + wthetaLev(i,j,k,bi,bj)*cosphi(i,j,bi,bj)       &                 + wthetaLev(i,j,k,bi,bj)*cosphi(i,j,bi,bj)
124       &                 *tmpx*tmpx       &                 *tmpx*tmpx
125    #else
126                      fctile = fctile + tmpx*tmpx
127    #endif
128                    if ( wthetaLev(i,j,k,bi,bj)*cosphi(i,j,bi,bj).ne.0. )                    if ( wthetaLev(i,j,k,bi,bj)*cosphi(i,j,bi,bj).ne.0. )
129       &                 num_temp0(bi,bj) = num_temp0(bi,bj) + 1. _d 0       &                 num_temp0(bi,bj) = num_temp0(bi,bj) + 1. _d 0
130                  endif                  endif
# Line 133  c--         Determine the weights to be Line 137  c--         Determine the weights to be
137            enddo            enddo
138          enddo          enddo
139    
140    #ifndef ALLOW_SMOOTH_CORREL3D
141  #ifdef ALLOW_SMOOTH_IC_COST_CONTRIBUTION  #ifdef ALLOW_SMOOTH_IC_COST_CONTRIBUTION
142    
143        call active_read_xyz(        call active_read_xyz(
# Line 159  c--     Loop over this thread's tiles. Line 164  c--     Loop over this thread's tiles.
164       &                *maskS(i,j+1,k,bi,bj)*maskS(i,j+2,k,bi,bj)       &                *maskS(i,j+1,k,bi,bj)*maskS(i,j+2,k,bi,bj)
165       &            + ( tmpfld3d(i,j+1,k,bi,bj)-tmpfld3d(i,j,k,bi,bj) )       &            + ( tmpfld3d(i,j+1,k,bi,bj)-tmpfld3d(i,j,k,bi,bj) )
166       &                *maskS(i,j+1,k,bi,bj)       &                *maskS(i,j+1,k,bi,bj)
167                    if ( ABS(R_low(i,j,bi,bj)) .LT. 100. )                    if ( ABS(R_low(i,j,bi,bj)) .LT. 100. )
168       &              tmpx = tmpx*ABS(R_low(i,j,bi,bj))/100.       &              tmpx = tmpx*ABS(R_low(i,j,bi,bj))/100.
169                    fctile = fctile                    fctile = fctile
170       &               + wthetaLev(i,j,k,bi,bj)*cosphi(i,j,bi,bj)       &               + wthetaLev(i,j,k,bi,bj)*cosphi(i,j,bi,bj)
# Line 175  c--     Loop over this thread's tiles. Line 180  c--     Loop over this thread's tiles.
180            enddo            enddo
181          enddo          enddo
182  #endif  #endif
183    #endif
184    
185  #endif  #endif
186    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22