| 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 |