/[MITgcm]/MITgcm/pkg/exf/exf_set_obcs.F
ViewVC logotype

Diff of /MITgcm/pkg/exf/exf_set_obcs.F

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

revision 1.4 by dimitri, Tue Feb 18 05:33:54 2003 UTC revision 1.5 by heimbach, Tue Jun 24 16:07:32 2003 UTC
# Line 140  c     == local variables == Line 140  c     == local variables ==
140    
141  c     == end of interface ==  c     == end of interface ==
142    
143        if ( first ) then        if ( obcs_file .NE. ' ' ) then
         if ( obcs_file .NE. ' ' )  
      &        call mdsreadfieldyz( obcs_file, exf_iprec, exf_yftype, Nr  
      &                     , obcs_yz_1, count0, mythid  
      &                     )  
       endif  
144    
145        if (( first ) .or. ( changed )) then           if ( first ) then
146          call exf_swapffields_yz( obcs_yz_0, obcs_yz_1, mythid )              call mdsreadfieldyz( obcs_file, exf_iprec, exf_yftype, Nr
147         &           , obcs_yz_1, count0, mythid
148         &           )
149             endif
150    
151             if (( first ) .or. ( changed )) then
152                call exf_swapffields_yz( obcs_yz_0, obcs_yz_1, mythid )
153    
154                call mdsreadfieldyz( obcs_file, exf_iprec, exf_yftype, Nr
155         &           , obcs_yz_1, count1, mythid
156         &           )
157             endif
158    
159             do bj = mybylo(mythid),mybyhi(mythid)
160                do bi = mybxlo(mythid),mybxhi(mythid)
161                   do k = 1,Nr
162                      do j = 1,sny
163                         obcs_fld_yz(j,k,bi,bj) =
164         &                    fac             *obcs_yz_0(j,k,bi,bj) +
165         &                    (exf_one - fac) *obcs_yz_1(j,k,bi,bj)
166                      enddo
167                   enddo
168                enddo
169             enddo
170    
         if ( obcs_file .NE. ' ' )  
      &       call mdsreadfieldyz( obcs_file, exf_iprec, exf_yftype, Nr  
      &                     , obcs_yz_1, count1, mythid  
      &                     )  
171        endif        endif
172    
       do bj = mybylo(mythid),mybyhi(mythid)  
         do bi = mybxlo(mythid),mybxhi(mythid)  
           do k = 1,Nr  
             do j = 1,sny  
               obcs_fld_yz(j,k,bi,bj) =  
      &              fac             *obcs_yz_0(j,k,bi,bj) +  
      &              (exf_one - fac) *obcs_yz_1(j,k,bi,bj)  
             enddo  
           enddo  
         enddo  
       enddo  
   
173  #endif  #endif
174    
175        end        end
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22