33 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
34 |
#include "SIZE.h" |
#include "SIZE.h" |
35 |
#include "DYNVARS.h" |
#include "DYNVARS.h" |
36 |
|
#include "TR1.h" |
37 |
|
|
38 |
#include "ctrl.h" |
#include "ctrl.h" |
39 |
#include "ctrl_dummy.h" |
#include "ctrl_dummy.h" |
40 |
|
#include "optim.h" |
41 |
|
|
42 |
c == routine arguments == |
c == routine arguments == |
43 |
|
|
60 |
|
|
61 |
character*( 80) fnametheta |
character*( 80) fnametheta |
62 |
character*( 80) fnamesalt |
character*( 80) fnamesalt |
63 |
|
character*( 80) fnametr1 |
64 |
|
|
65 |
c == external == |
c == external == |
66 |
|
|
139 |
enddo |
enddo |
140 |
#endif |
#endif |
141 |
|
|
142 |
|
#ifdef ALLOW_TR10_CONTROL |
143 |
|
c-- Temperature field. |
144 |
|
il=ilnblnk( xx_tr1_file ) |
145 |
|
write(fnametr1(1:80),'(2a,i10.10)') |
146 |
|
& xx_tr1_file(1:il),'.',optimcycle |
147 |
|
call active_read_xyz( fnametr1, tmpfld3d, 1, |
148 |
|
& doglobalread, ladinit, optimcycle, |
149 |
|
& mythid, xx_tr1_dummy ) |
150 |
|
|
151 |
|
do bj = jtlo,jthi |
152 |
|
do bi = itlo,ithi |
153 |
|
do k = 1,nr |
154 |
|
do j = jmin,jmax |
155 |
|
do i = imin,imax |
156 |
|
tr1(i,j,k,bi,bj) = tr1(i,j,k,bi,bj) + |
157 |
|
& fac*tmpfld3d(i,j,k,bi,bj) |
158 |
|
gsNm1(i,j,k,bi,bj) = gsNm1(i,j,k,bi,bj) + |
159 |
|
& fac*tmpfld3d(i,j,k,bi,bj) |
160 |
|
enddo |
161 |
|
enddo |
162 |
|
enddo |
163 |
|
enddo |
164 |
|
enddo |
165 |
|
#endif |
166 |
|
|
167 |
|
|
168 |
c-- Update the tile edges. |
c-- Update the tile edges. |
169 |
|
|
175 |
_EXCH_XYZ_R8( salt, mythid ) |
_EXCH_XYZ_R8( salt, mythid ) |
176 |
_EXCH_XYZ_R8( gsNm1, mythid ) |
_EXCH_XYZ_R8( gsNm1, mythid ) |
177 |
#endif |
#endif |
178 |
|
#ifdef ALLOW_TR10_CONTROL |
179 |
|
_EXCH_XYZ_R8( tr1, mythid ) |
180 |
|
_EXCH_XYZ_R8( gTr1Nm1, mythid ) |
181 |
|
#endif |
182 |
|
|
183 |
return |
return |
184 |
end |
end |