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" |
56 |
|
#include "exf_clim_param.h" |
57 |
#include "exf_clim_fields.h" |
#include "exf_clim_fields.h" |
58 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
59 |
# include "tamc.h" |
# include "tamc.h" |
68 |
c == local variables == |
c == local variables == |
69 |
|
|
70 |
integer bi,bj |
integer bi,bj |
71 |
integer i,j |
integer i,j,k |
72 |
integer jtlo |
integer jtlo |
73 |
integer jthi |
integer jthi |
74 |
integer itlo |
integer itlo |
126 |
do j = jmin,jmax |
do j = jmin,jmax |
127 |
do i = imin,imax |
do i = imin,imax |
128 |
c Zonal wind stress. |
c Zonal wind stress. |
129 |
if (ustress(i,j,bi,bj).gt.2.0D0) then |
if (ustress(i,j,bi,bj).gt.windstressmax) then |
130 |
ustress(i,j,bi,bj)=2.0D0 |
ustress(i,j,bi,bj)=windstressmax |
131 |
endif |
endif |
132 |
enddo |
enddo |
133 |
enddo |
enddo |
136 |
#endif |
#endif |
137 |
do j = jmin,jmax |
do j = jmin,jmax |
138 |
do i = imin,imax |
do i = imin,imax |
139 |
if (ustress(i,j,bi,bj).lt.-2.0D0) then |
if (ustress(i,j,bi,bj).lt.-windstressmax) then |
140 |
ustress(i,j,bi,bj)=-2.0D0 |
ustress(i,j,bi,bj)=-windstressmax |
141 |
endif |
endif |
142 |
enddo |
enddo |
143 |
enddo |
enddo |
160 |
do j = jmin,jmax |
do j = jmin,jmax |
161 |
do i = imin,imax |
do i = imin,imax |
162 |
c Meridional wind stress. |
c Meridional wind stress. |
163 |
if (vstress(i,j,bi,bj).gt.2.0D0) then |
if (vstress(i,j,bi,bj).gt.windstressmax) then |
164 |
vstress(i,j,bi,bj)=2.0D0 |
vstress(i,j,bi,bj)=windstressmax |
165 |
endif |
endif |
166 |
enddo |
enddo |
167 |
enddo |
enddo |
170 |
#endif |
#endif |
171 |
do j = jmin,jmax |
do j = jmin,jmax |
172 |
do i = imin,imax |
do i = imin,imax |
173 |
if (vstress(i,j,bi,bj).lt.-2.0D0) then |
if (vstress(i,j,bi,bj).lt.-windstressmax) then |
174 |
vstress(i,j,bi,bj)=-2.0D0 |
vstress(i,j,bi,bj)=-windstressmax |
175 |
endif |
endif |
176 |
enddo |
enddo |
177 |
enddo |
enddo |
213 |
enddo |
enddo |
214 |
#endif |
#endif |
215 |
|
|
216 |
|
#ifdef ALLOW_CLIMTEMP_RELAXATION |
217 |
|
if ( climtempfile .NE. ' ' ) then |
218 |
|
do k = 1, Nr |
219 |
|
do j = jmin,jmax |
220 |
|
do i = imin,imax |
221 |
|
thetaStar(i,j,k,bi,bj) = climtemp(i,j,k,bi,bj) |
222 |
|
enddo |
223 |
|
enddo |
224 |
|
enddo |
225 |
|
endif |
226 |
|
#endif |
227 |
|
|
228 |
|
#ifdef ALLOW_CLIMSALT_RELAXATION |
229 |
|
if ( climsaltfile .NE. ' ' ) then |
230 |
|
do k = 1, Nr |
231 |
|
do j = jmin,jmax |
232 |
|
do i = imin,imax |
233 |
|
saltStar(i,j,k,bi,bj) = climsalt(i,j,k,bi,bj) |
234 |
|
enddo |
235 |
|
enddo |
236 |
|
enddo |
237 |
|
endif |
238 |
|
#endif |
239 |
|
|
240 |
#ifdef ATMOSPHERIC_LOADING |
#ifdef ATMOSPHERIC_LOADING |
241 |
do j = jmin,jmax |
do j = jmin,jmax |
242 |
do i = imin,imax |
do i = imin,imax |
264 |
#ifdef ALLOW_CLIMSSS_RELAXATION |
#ifdef ALLOW_CLIMSSS_RELAXATION |
265 |
_EXCH_XY_R4( sss, mythid ) |
_EXCH_XY_R4( sss, mythid ) |
266 |
#endif |
#endif |
267 |
|
#ifdef ALLOW_CLIMTEMP_RELAXATION |
268 |
|
_EXCH_XYZ_R4( thetaStar, mythid ) |
269 |
|
#endif |
270 |
|
#ifdef ALLOW_CLIMSALT_RELAXATION |
271 |
|
_EXCH_XYZ_R4( saltStar, mythid ) |
272 |
|
#endif |
273 |
#ifdef ATMOSPHERIC_LOADING |
#ifdef ATMOSPHERIC_LOADING |
274 |
_EXCH_XY_R4( pload, mythid ) |
_EXCH_XY_R4( pload, mythid ) |
275 |
#endif |
#endif |