218 |
|
|
219 |
if ( (obcswfirst) .or. (obcswchanged)) then |
if ( (obcswfirst) .or. (obcswchanged)) then |
220 |
|
|
|
cgg( This is a terribly long way to do it. However, the dimensions do not exactly |
|
|
cgg match up. I will blame Fortran for the ugliness. |
|
|
|
|
|
do bj = jtlo,jthi |
|
|
do bi = itlo,ithi |
|
|
do k = 1,nr |
|
|
do j = jmin,jmax |
|
|
tmpfldyz(j,k,bi,bj) = xx_obcsw1(j,k,bi,bj,iobcs) |
|
|
enddo |
|
|
enddo |
|
|
enddo |
|
|
enddo |
|
|
|
|
|
call exf_swapffields_yz( tmpfldyz2, tmpfldyz, mythid) |
|
|
|
|
221 |
do bj = jtlo,jthi |
do bj = jtlo,jthi |
222 |
do bi = itlo,ithi |
do bi = itlo,ithi |
223 |
do k = 1,nr |
do k = 1,nr |
224 |
do j = jmin,jmax |
do j = jmin,jmax |
225 |
xx_obcsw0(j,k,bi,bj,iobcs) = tmpfldyz2(j,k,bi,bj) |
xx_obcsw0(j,k,bi,bj,iobcs) = xx_obcsw1(j,k,bi,bj,iobcs) |
226 |
enddo |
tmpfldyz (j,k,bi,bj) = 0. _d 0 |
227 |
enddo |
enddo |
228 |
enddo |
enddo |
229 |
|
enddo |
230 |
enddo |
enddo |
231 |
|
|
232 |
call active_read_yz( fnameobcsw, tmpfldyz, |
call active_read_yz( fnameobcsw, tmpfldyz, |