2 |
|
|
3 |
#include "EXF_OPTIONS.h" |
#include "EXF_OPTIONS.h" |
4 |
|
|
|
|
|
5 |
subroutine exf_mapfields( mythid ) |
subroutine exf_mapfields( mythid ) |
6 |
|
|
7 |
c ================================================================== |
c ================================================================== |
48 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
49 |
#include "SIZE.h" |
#include "SIZE.h" |
50 |
#include "FFIELDS.h" |
#include "FFIELDS.h" |
51 |
|
#include "GRID.h" |
52 |
|
|
53 |
#include "exf_param.h" |
#include "exf_param.h" |
54 |
#include "exf_constants.h" |
#include "exf_constants.h" |
55 |
#include "exf_fields.h" |
#include "exf_fields.h" |
141 |
enddo |
enddo |
142 |
enddo |
enddo |
143 |
do j = jmin,jmax |
do j = jmin,jmax |
144 |
do i = imin,imax |
do i = imin+1,imax |
145 |
fu(i,j,bi,bj) = exf_outscal_ustress*ustress(i,j,bi,bj) |
#if (defined (ALLOW_ATM_TEMP) || defined (ALLOW_ATM_WIND)) |
146 |
|
c Shift wind stresses calculated at C-points to W/S points |
147 |
|
fu(i,j,bi,bj) = exf_outscal_ustress* |
148 |
|
& (ustress(i,j,bi,bj)+ustress(i-1,j,bi,bj))/2.* |
149 |
|
& maskW(i,j,1,bi,bj) |
150 |
|
#else |
151 |
|
fu(i,j,bi,bj) = exf_outscal_ustress*ustress(i,j,bi,bj) |
152 |
|
#endif |
153 |
enddo |
enddo |
154 |
enddo |
enddo |
155 |
|
|
174 |
endif |
endif |
175 |
enddo |
enddo |
176 |
enddo |
enddo |
177 |
do j = jmin,jmax |
do j = jmin+1,jmax |
178 |
do i = imin,imax |
do i = imin,imax |
179 |
fv(i,j,bi,bj) = exf_outscal_vstress*vstress(i,j,bi,bj) |
#if (defined (ALLOW_ATM_TEMP) || defined (ALLOW_ATM_WIND)) |
180 |
|
c Shift wind stresses calculated at C-points to W/S points |
181 |
|
fv(i,j,bi,bj) = exf_outscal_vstress* |
182 |
|
& (vstress(i,j,bi,bj)+vstress(i,j-1,bi,bj))/2.* |
183 |
|
& maskS(i,j,1,bi,bj) |
184 |
|
#else |
185 |
|
fv(i,j,bi,bj) = exf_outscal_vstress*vstress(i,j,bi,bj) |
186 |
|
#endif |
187 |
enddo |
enddo |
188 |
enddo |
enddo |
189 |
|
|