142 |
enddo |
enddo |
143 |
do j = jmin,jmax |
do j = jmin,jmax |
144 |
do i = imin+1,imax |
do i = imin+1,imax |
145 |
#if (defined (ALLOW_ATM_TEMP) || defined (ALLOW_ATM_WIND)) |
#if (defined (ALLOW_BULKFORMULAE) || defined (USE_EXF_INTERPOLATION)) |
146 |
c Shift wind stresses calculated at C-points to W/S points |
c Shift wind stresses calculated at C-points to W/S points |
147 |
fu(i,j,bi,bj) = exf_outscal_ustress* |
fu(i,j,bi,bj) = exf_outscal_ustress* |
148 |
& (ustress(i,j,bi,bj)+ustress(i-1,j,bi,bj))/2.* |
& (ustress(i,j,bi,bj)+ustress(i-1,j,bi,bj))/2.* |
176 |
enddo |
enddo |
177 |
do j = jmin+1,jmax |
do j = jmin+1,jmax |
178 |
do i = imin,imax |
do i = imin,imax |
179 |
#if (defined (ALLOW_ATM_TEMP) || defined (ALLOW_ATM_WIND)) |
#if (defined (ALLOW_BULKFORMULAE) || defined (USE_EXF_INTERPOLATION)) |
180 |
c Shift wind stresses calculated at C-points to W/S points |
c Shift wind stresses calculated at C-points to W/S points |
181 |
fv(i,j,bi,bj) = exf_outscal_vstress* |
fv(i,j,bi,bj) = exf_outscal_vstress* |
182 |
& (vstress(i,j,bi,bj)+vstress(i,j-1,bi,bj))/2.* |
& (vstress(i,j,bi,bj)+vstress(i,j-1,bi,bj))/2.* |