183 |
if(mod(int(VIEW(ic)),2).eq.1)then !Y-view |
if(mod(int(VIEW(ic)),2).eq.1)then !Y-view |
184 |
|
|
185 |
if( abs(angle).ge.e2fay.and.abs(angle).le.e2tay )then |
if( abs(angle).ge.e2fay.and.abs(angle).le.e2tay )then |
186 |
pfaeta = pfaeta2(ic,angle)+pfcorr(ic,angle) |
pfaeta = pfaeta2(ic,angle)+pfacorr(ic,angle) |
187 |
elseif( abs(angle).ge.e3fay.and.abs(angle).le.e3tay )then |
elseif( abs(angle).ge.e3fay.and.abs(angle).le.e3tay )then |
188 |
pfaeta = pfaeta3(ic,angle)+pfcorr(ic,angle) |
pfaeta = pfaeta3(ic,angle)+pfacorr(ic,angle) |
189 |
elseif( abs(angle).ge.e4fay.and.abs(angle).le.e4tay )then |
elseif( abs(angle).ge.e4fay.and.abs(angle).le.e4tay )then |
190 |
pfaeta = pfaeta4(ic,angle)+pfcorr(ic,angle) |
pfaeta = pfaeta4(ic,angle)+pfacorr(ic,angle) |
191 |
else |
else |
192 |
pfaeta = cog(4,ic) |
pfaeta = cog(4,ic) |
193 |
endif |
endif |
195 |
else !X-view |
else !X-view |
196 |
|
|
197 |
if( abs(angle).ge.e2fax.and.abs(angle).le.e2tax )then |
if( abs(angle).ge.e2fax.and.abs(angle).le.e2tax )then |
198 |
pfaeta = pfaeta2(ic,angle)+pfcorr(ic,angle) |
pfaeta = pfaeta2(ic,angle)+pfacorr(ic,angle) |
199 |
elseif( abs(angle).ge.e3fax.and.abs(angle).le.e3tax )then |
elseif( abs(angle).ge.e3fax.and.abs(angle).le.e3tax )then |
200 |
pfaeta = pfaeta3(ic,angle)+pfcorr(ic,angle) |
pfaeta = pfaeta3(ic,angle)+pfacorr(ic,angle) |
201 |
elseif( abs(angle).ge.e4fax.and.abs(angle).le.e4tax )then |
elseif( abs(angle).ge.e4fax.and.abs(angle).le.e4tax )then |
202 |
pfaeta = pfaeta4(ic,angle)+pfcorr(ic,angle) |
pfaeta = pfaeta4(ic,angle)+pfacorr(ic,angle) |
203 |
else |
else |
204 |
pfaeta = cog(4,ic) |
pfaeta = cog(4,ic) |
205 |
endif |
endif |