/[MITgcm]/MITgcm/pkg/obcs/obcs_init_fixed.F
ViewVC logotype

Diff of /MITgcm/pkg/obcs/obcs_init_fixed.F

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

revision 1.3 by heimbach, Mon May 14 21:36:45 2001 UTC revision 1.4 by jmc, Wed Jan 30 04:22:31 2002 UTC
# Line 17  C     === Global variables === Line 17  C     === Global variables ===
17  #include "EEPARAMS.h"  #include "EEPARAMS.h"
18  #include "PARAMS.h"  #include "PARAMS.h"
19  #include "OBCS.h"  #include "OBCS.h"
20    #ifdef NONLIN_FRSURF
21    #include "GRID.h"
22    #include "SURFACE.h"
23    #endif
24    
25  C     == Routine arguments ==  C     == Routine arguments ==
26  C     myThid -  Number of this instance of INI_DEPTHS  C     myThid -  Number of this instance of INI_DEPTHS
# Line 70  C OB_Jnorth(iGm) allows for the northern Line 74  C OB_Jnorth(iGm) allows for the northern
74           ENDDO           ENDDO
75          ENDDO          ENDDO
76    
77    #ifdef NONLIN_FRSURF
78    
79    C- save the initial hFacS at the N & S boundaries :
80            DO i=1-Olx,sNx+Olx
81              OBNhFac0(i,bi,bj)=0.
82              OBShFac0(i,bi,bj)=0.
83    C  Northern boundary
84              IF ( OB_Jn(i,bi,bj).NE.0 ) THEN
85                j = OB_Jn(i,bi,bj)
86                k = ksurfS(i,j,bi,bj)
87                IF (k.LE.Nr) OBNhFac0(i,bi,bj)=hFacS(i,j,k,bi,bj)
88              ENDIF
89    C  Southern boundary
90              IF ( OB_Js(i,bi,bj).NE.0 ) THEN
91                j = OB_Js(i,bi,bj)+1
92                k = ksurfS(i,j,bi,bj)
93                IF (k.LE.Nr) OBShFac0(i,bi,bj)=hFacS(i,j,k,bi,bj)
94              ENDIF
95            ENDDO    
96    
97    C- save the initial hFacW at the E & W boundaries :
98            DO j=1-Oly,sNy+Oly
99              OBEhFac0(j,bi,bj)=0.
100              OBWhFac0(j,bi,bj)=0.
101    C  Eastern boundary
102              IF ( OB_Ie(j,bi,bj).NE.0 ) THEN
103                i = OB_Ie(j,bi,bj)
104                k = ksurfW(i,j,bi,bj)
105                IF (k.LE.Nr) OBEhFac0(j,bi,bj)=hFacW(i,j,k,bi,bj)
106              ENDIF
107    C  Western boundary
108              IF ( OB_Iw(j,bi,bj).NE.0 ) THEN
109                i = OB_Iw(j,bi,bj)+1
110                k = ksurfW(i,j,bi,bj)
111                IF (k.LE.Nr) OBWhFac0(j,bi,bj)=hFacW(i,j,k,bi,bj)
112              ENDIF
113            ENDDO    
114    
115    #endif /* NONLIN_FRSURF */    
116    
117         ENDDO         ENDDO
118        ENDDO        ENDDO
119    
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22