/[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.3 by heimbach, Tue Nov 12 20:34:41 2002 UTC revision 1.4 by dimitri, Tue Feb 18 05:33:54 2003 UTC
# Line 14  c Line 14  c
14  c     o set open boundary conditions  c     o set open boundary conditions
15  c  c
16  c     started: heimbach@mit.edu 01-May-2001  c     started: heimbach@mit.edu 01-May-2001
17    c     mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002
18    
19  c     ==================================================================  c     ==================================================================
20  c     SUBROUTINE exf_set_obcs_xz  c     SUBROUTINE exf_set_obcs_xz
# Line 53  c     == local variables == Line 54  c     == local variables ==
54    
55  c     == end of interface ==  c     == end of interface ==
56    
57        if ( first ) then        if ( obcs_file .NE. ' ' ) then
         if ( obcs_file .NE. ' ' )  
      &        call mdsreadfieldxz( obcs_file, exf_iprec, exf_yftype, Nr  
      &                     , obcs_xz_1, count0, mythid  
      &                     )  
       endif  
58    
59        if (( first ) .or. ( changed )) then           if ( first ) then
60          call exf_swapffields_xz( obcs_xz_0, obcs_xz_1, mythid )              call mdsreadfieldxz( obcs_file, exf_iprec, exf_yftype, Nr
61         &           , obcs_xz_1, count0, mythid
62         &           )
63             endif
64    
65             if (( first ) .or. ( changed )) then
66                call exf_swapffields_xz( obcs_xz_0, obcs_xz_1, mythid )
67    
68                call mdsreadfieldxz( obcs_file, exf_iprec, exf_yftype, Nr
69         &           , obcs_xz_1, count1, mythid
70         &           )
71             endif
72    
73             do bj = mybylo(mythid),mybyhi(mythid)
74                do bi = mybxlo(mythid),mybxhi(mythid)
75                   do k = 1,Nr
76                      do i = 1,snx
77                         obcs_fld_xz(i,k,bi,bj) =
78         &                    fac * obcs_xz_0(i,k,bi,bj) +
79         &                    (exf_one - fac) * obcs_xz_1(i,k,bi,bj)
80                      enddo
81                   enddo
82                enddo
83             enddo
84    
         if ( obcs_file .NE. ' ' )  
      &       call mdsreadfieldxz( obcs_file, exf_iprec, exf_yftype, Nr  
      &                     , obcs_xz_1, count1, mythid  
      &                     )  
85        endif        endif
86    
       do bj = mybylo(mythid),mybyhi(mythid)  
         do bi = mybxlo(mythid),mybxhi(mythid)  
           do k = 1,Nr  
             do i = 1,snx  
               obcs_fld_xz(i,k,bi,bj) =  
      &              fac             *obcs_xz_0(i,k,bi,bj) +  
      &              (exf_one - fac) *obcs_xz_1(i,k,bi,bj)  
             enddo  
           enddo  
         enddo  
       enddo  
   
87  #endif  #endif
88    
89        end        end
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22