148 |
if( check ) then |
if( check ) then |
149 |
undef = getcon('UNDEF') |
undef = getcon('UNDEF') |
150 |
do i= 1,Lena |
do i= 1,Lena |
151 |
jindx = 1 + (indx(i+offset-1)-1/sNx) |
jindx = 1.0001 + (indx(i+offset-1)-1/sNx) |
152 |
newindx = indx(i+offset-1) + (jindx-1)*Olx+(sNx+Olx)*Oly+Olx |
newindx = indx(i+offset-1)+(jindx-1)*2*Olx+(sNx+2*Olx)*Oly+Olx |
153 |
if(qdiag(newindx,1,kd,bi,bj).eq.undef |
if(qdiag(newindx,1,kd,bi,bj).eq.undef |
154 |
. .or.field(i,k).eq.undef)then |
. .or.field(i,k).eq.undef)then |
155 |
qdiag(newindx,1,kd,bi,bj) = undef |
qdiag(newindx,1,kd,bi,bj) = undef |
160 |
enddo |
enddo |
161 |
else |
else |
162 |
do i= 1,Lena |
do i= 1,Lena |
163 |
jindx = 1 + (indx(i+offset-1)-1/sNx) |
jindx = 1.0001 + (indx(i+offset-1)-1/sNx) |
164 |
newindx = indx(i+offset-1) + (jindx-1)*Olx+(sNx+Olx)*Oly+Olx |
newindx = indx(i+offset-1)+(jindx-1)*2*Olx+(sNx+2*Olx)*Oly+Olx |
165 |
qdiag(newindx,1,kd,bi,bj)=qdiag(newindx,1,kd,bi,bj)+ |
qdiag(newindx,1,kd,bi,bj)=qdiag(newindx,1,kd,bi,bj)+ |
166 |
. field(i,k)*chfr(i) |
. field(i,k)*chfr(i) |
167 |
enddo |
enddo |