25 |
#include "PARAMS.h" |
#include "PARAMS.h" |
26 |
#include "DYNVARS.h" |
#include "DYNVARS.h" |
27 |
#include "GRID.h" |
#include "GRID.h" |
|
#include "TR1.h" |
|
28 |
#include "ctrl.h" |
#include "ctrl.h" |
29 |
#include "ctrl_dummy.h" |
#include "ctrl_dummy.h" |
30 |
#include "optim.h" |
#include "optim.h" |
31 |
|
#ifdef ALLOW_PTRACERS |
32 |
|
# include "PTRACERS_SIZE.h" |
33 |
|
# include "PTRACERS.h" |
34 |
|
#endif |
35 |
#ifdef ALLOW_ECCO |
#ifdef ALLOW_ECCO |
36 |
# include "ecco_cost.h" |
# include "ecco_cost.h" |
37 |
#endif |
#endif |
157 |
#endif |
#endif |
158 |
|
|
159 |
#ifdef ALLOW_TR10_CONTROL |
#ifdef ALLOW_TR10_CONTROL |
160 |
|
#ifdef ALLOW_PTRACERS |
161 |
c-- Temperature field. |
c-- Temperature field. |
162 |
il=ilnblnk( xx_tr1_file ) |
il=ilnblnk( xx_tr1_file ) |
163 |
write(fnametr1(1:80),'(2a,i10.10)') |
write(fnametr1(1:80),'(2a,i10.10)') |
171 |
do k = 1,nr |
do k = 1,nr |
172 |
do j = jmin,jmax |
do j = jmin,jmax |
173 |
do i = imin,imax |
do i = imin,imax |
174 |
tr1(i,j,k,bi,bj) = tr1(i,j,k,bi,bj) + |
ptracer(i,j,k,bi,bj,1) = ptracer(i,j,k,bi,bj,1) + |
175 |
& fac*tmpfld3d(i,j,k,bi,bj) |
& fac*tmpfld3d(i,j,k,bi,bj) |
176 |
enddo |
enddo |
177 |
enddo |
enddo |
179 |
enddo |
enddo |
180 |
enddo |
enddo |
181 |
#endif |
#endif |
182 |
|
#endif |
183 |
|
|
184 |
#ifdef ALLOW_DIFFKR_CONTROL |
#ifdef ALLOW_DIFFKR_CONTROL |
185 |
c-- diffkr. |
c-- diffkr. |
306 |
|
|
307 |
#ifdef ALLOW_THETA0_CONTROL |
#ifdef ALLOW_THETA0_CONTROL |
308 |
_EXCH_XYZ_R8( theta, mythid ) |
_EXCH_XYZ_R8( theta, mythid ) |
|
cph _EXCH_XYZ_R8( gtNm1, mythid ) |
|
309 |
#endif |
#endif |
310 |
#ifdef ALLOW_SALT0_CONTROL |
#ifdef ALLOW_SALT0_CONTROL |
311 |
_EXCH_XYZ_R8( salt, mythid ) |
_EXCH_XYZ_R8( salt, mythid ) |
|
cph _EXCH_XYZ_R8( gsNm1, mythid ) |
|
312 |
#endif |
#endif |
313 |
#ifdef ALLOW_TR10_CONTROL |
#ifdef ALLOW_TR10_CONTROL |
314 |
_EXCH_XYZ_R8( tr1, mythid ) |
#ifdef ALLOW_PTRACERS |
315 |
cph _EXCH_XYZ_R8( gTr1Nm1, mythid ) |
_EXCH_XYZ_R8(pTracer(1-Olx,1-Oly,1,1,1,1),myThid) |
316 |
|
#endif |
317 |
#endif |
#endif |
318 |
#ifdef ALLOW_DIFFKR_CONTROL |
#ifdef ALLOW_DIFFKR_CONTROL |
319 |
_EXCH_XYZ_R8( diffkr, mythid) |
_EXCH_XYZ_R8( diffkr, mythid) |