/[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.3 by heimbach, Mon Mar 28 23:49:49 2005 UTC revision 1.4 by heimbach, Sat Oct 15 18:27:32 2005 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. 1000. )                    if ( ABS(R_low(i,j,bi,bj)) .LT. 100. )
120       &              tmpx = tmpx*ABS(R_low(i,j,bi,bj))/1000.       &              tmpx = tmpx*ABS(R_low(i,j,bi,bj))/100.
121                    fctile = fctile                    fctile = fctile
122       &                 + wtheta(k,bi,bj)*cosphi(i,j,bi,bj)       &                 + wthetaLev(i,j,k,bi,bj)*cosphi(i,j,bi,bj)
123       &                 *tmpx*tmpx       &                 *tmpx*tmpx
124                    if ( wtheta(k,bi,bj)*cosphi(i,j,bi,bj) .ne. 0. )                    if ( wthetaLev(i,j,k,bi,bj)*cosphi(i,j,bi,bj).ne.0. )
125       &                 num_temp0(bi,bj) = num_temp0(bi,bj) + 1. _d 0       &                 num_temp0(bi,bj) = num_temp0(bi,bj) + 1. _d 0
126                  endif                  endif
127                enddo                enddo
# Line 169  c--     Print cost function for all tile Line 169  c--     Print cost function for all tile
169       &                      SQUEEZE_RIGHT , mythid)       &                      SQUEEZE_RIGHT , mythid)
170  #endif  #endif
171    
 #else  
 c--   Do not enter the calculation of the salinity increment  
 c--   contribution to the final cost function.  
   
       fctile   = 0. _d 0  
       fcthread = 0. _d 0  
   
 #ifdef ECCO_VERBOSE  
       _BEGIN_MASTER( mythid )  
         write(msgbuf,'(a)') ' '  
         call print_message( msgbuf, standardmessageunit,  
      &                      SQUEEZE_RIGHT , mythid)  
         write(msgbuf,'(a,a)')  
      &    ' cost_theta0 : no contribution of the I.C. in salin. ',  
      &                    ' to cost function.'  
         call print_message( msgbuf, standardmessageunit,  
      &                      SQUEEZE_RIGHT , mythid)  
         write(msgbuf,'(a)') ' '  
         call print_message( msgbuf, standardmessageunit,  
      &                      SQUEEZE_RIGHT , mythid)  
       _END_MASTER( mythid )  
 #endif  
   
172  #endif  #endif
173    
174        return        return
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22