69 |
|
|
70 |
c == local variables == |
c == local variables == |
71 |
|
|
72 |
integer bi,bj |
INTEGER bi,bj |
73 |
integer i,j,k |
INTEGER i,j,ks |
74 |
INTEGER imin, imax |
INTEGER imin, imax |
75 |
INTEGER jmin, jmax |
INTEGER jmin, jmax |
76 |
PARAMETER ( imin = 1-OLx , imax = sNx+OLx ) |
PARAMETER ( imin = 1-OLx , imax = sNx+OLx ) |
78 |
|
|
79 |
c == end of interface == |
c == end of interface == |
80 |
|
|
81 |
|
C-- set surface level index: |
82 |
|
ks = 1 |
83 |
|
|
84 |
DO bj = myByLo(myThid),myByHi(myThid) |
DO bj = myByLo(myThid),myByHi(myThid) |
85 |
DO bi = myBxLo(myThid),myBxHi(myThid) |
DO bi = myBxLo(myThid),myBxHi(myThid) |
86 |
|
|
163 |
c Shift wind stresses calculated at Grid-center to W/S points |
c Shift wind stresses calculated at Grid-center to W/S points |
164 |
fu(i,j,bi,bj) = exf_outscal_ustress* |
fu(i,j,bi,bj) = exf_outscal_ustress* |
165 |
& (ustress(i,j,bi,bj)+ustress(i-1,j,bi,bj)) |
& (ustress(i,j,bi,bj)+ustress(i-1,j,bi,bj)) |
166 |
& *exf_half*maskW(i,j,1,bi,bj) |
& *exf_half*maskW(i,j,ks,bi,bj) |
167 |
enddo |
enddo |
168 |
enddo |
enddo |
169 |
ENDIF |
ENDIF |
201 |
c Shift wind stresses calculated at C-points to W/S points |
c Shift wind stresses calculated at C-points to W/S points |
202 |
fv(i,j,bi,bj) = exf_outscal_vstress* |
fv(i,j,bi,bj) = exf_outscal_vstress* |
203 |
& (vstress(i,j,bi,bj)+vstress(i,j-1,bi,bj)) |
& (vstress(i,j,bi,bj)+vstress(i,j-1,bi,bj)) |
204 |
& *exf_half*maskS(i,j,1,bi,bj) |
& *exf_half*maskS(i,j,ks,bi,bj) |
205 |
enddo |
enddo |
206 |
enddo |
enddo |
207 |
ENDIF |
ENDIF |