1 |
C $Header$ |
C $Header$ |
2 |
|
|
3 |
#include "CTRL_CPPOPTIONS.h" |
#include "CTRL_CPPOPTIONS.h" |
4 |
|
#ifdef ALLOW_PTRACERS |
5 |
|
# include "PTRACERS_OPTIONS.h" |
6 |
|
#endif |
7 |
|
|
8 |
CBOP |
CBOP |
9 |
C !ROUTINE: ctrl_map_ini |
C !ROUTINE: ctrl_map_ini |
28 |
#include "PARAMS.h" |
#include "PARAMS.h" |
29 |
#include "DYNVARS.h" |
#include "DYNVARS.h" |
30 |
#include "GRID.h" |
#include "GRID.h" |
31 |
#include "TR1.h" |
#ifdef ALLOW_PASSIVE_TRACER |
32 |
|
# include "TR1.h" |
33 |
|
#endif |
34 |
|
#ifdef ALLOW_PTRACERS |
35 |
|
# include "PTRACERS.h" |
36 |
|
#endif |
37 |
|
|
38 |
#include "ctrl.h" |
#include "ctrl.h" |
39 |
#include "ctrl_dummy.h" |
#include "ctrl_dummy.h" |
40 |
#include "optim.h" |
#include "optim.h" |
158 |
do k = 1,nr |
do k = 1,nr |
159 |
do j = jmin,jmax |
do j = jmin,jmax |
160 |
do i = imin,imax |
do i = imin,imax |
161 |
|
#if (defined (ALLOW_PASSIVE_TRACER)) |
162 |
tr1(i,j,k,bi,bj) = tr1(i,j,k,bi,bj) + |
tr1(i,j,k,bi,bj) = tr1(i,j,k,bi,bj) + |
163 |
& fac*tmpfld3d(i,j,k,bi,bj) |
& fac*tmpfld3d(i,j,k,bi,bj) |
164 |
|
#elif (defined (ALLOW_PTRACERS)) |
165 |
|
IF ( NUMBER_OF_PTRACERS .GT. 1 ) STOP |
166 |
|
& 'ALLOW_TR10_CONTROL with ALLOW_PTRACERS implemented for 1 tracer' |
167 |
|
ptracer(i,j,k,bi,bj,1) = ptracer(i,j,k,bi,bj,1) + |
168 |
|
& fac*tmpfld3d(i,j,k,bi,bj) |
169 |
|
#endif |
170 |
enddo |
enddo |
171 |
enddo |
enddo |
172 |
enddo |
enddo |
299 |
|
|
300 |
#ifdef ALLOW_THETA0_CONTROL |
#ifdef ALLOW_THETA0_CONTROL |
301 |
_EXCH_XYZ_R8( theta, mythid ) |
_EXCH_XYZ_R8( theta, mythid ) |
|
cph _EXCH_XYZ_R8( gtNm1, mythid ) |
|
302 |
#endif |
#endif |
303 |
#ifdef ALLOW_SALT0_CONTROL |
#ifdef ALLOW_SALT0_CONTROL |
304 |
_EXCH_XYZ_R8( salt, mythid ) |
_EXCH_XYZ_R8( salt, mythid ) |
|
cph _EXCH_XYZ_R8( gsNm1, mythid ) |
|
305 |
#endif |
#endif |
306 |
#ifdef ALLOW_TR10_CONTROL |
#ifdef ALLOW_TR10_CONTROL |
307 |
|
# if (defined (ALLOW_PASSIVE_TRACER)) |
308 |
_EXCH_XYZ_R8( tr1, mythid ) |
_EXCH_XYZ_R8( tr1, mythid ) |
309 |
cph _EXCH_XYZ_R8( gTr1Nm1, mythid ) |
# elif (defined (ALLOW_PTRACERS)) |
310 |
|
_EXCH_XYZ_R8(pTracer(1-Olx,1-Oly,1,1,1,1),myThid) |
311 |
|
# endif |
312 |
#endif |
#endif |
313 |
#ifdef ALLOW_DIFFKR_CONTROL |
#ifdef ALLOW_DIFFKR_CONTROL |
314 |
_EXCH_XYZ_R8( diffkr, mythid) |
_EXCH_XYZ_R8( diffkr, mythid) |