| 859 |
if( nsatstrips(icx).gt.0 )then |
if( nsatstrips(icx).gt.0 )then |
| 860 |
stripx = stripx + cog(4,icx) |
stripx = stripx + cog(4,icx) |
| 861 |
resxPAM = pitchX*1e-4/sqrt(12.) |
resxPAM = pitchX*1e-4/sqrt(12.) |
| 862 |
cc=cog(4,icx) |
cc cc=cog(4,icx) |
| 863 |
c$$$ print*,icx,' *** ',cc |
c$$$ print*,icx,' *** ',cc |
| 864 |
c$$$ print*,icx,' *** ',resxPAM |
c$$$ print*,icx,' *** ',resxPAM |
| 865 |
endif |
endif |
| 999 |
if( nsatstrips(icy).gt.0 )then |
if( nsatstrips(icy).gt.0 )then |
| 1000 |
stripy = stripy + cog(4,icy) |
stripy = stripy + cog(4,icy) |
| 1001 |
resyPAM = pitchY*1e-4/sqrt(12.) |
resyPAM = pitchY*1e-4/sqrt(12.) |
| 1002 |
cc=cog(4,icy) |
cc cc=cog(4,icy) |
| 1003 |
c$$$ print*,icy,' *** ',cc |
c$$$ print*,icy,' *** ',cc |
| 1004 |
c$$$ print*,icy,' *** ',resyPAM |
c$$$ print*,icy,' *** ',resyPAM |
| 1005 |
endif |
endif |
| 1253 |
c$$$ PFAx = 'COG4'!PFA |
c$$$ PFAx = 'COG4'!PFA |
| 1254 |
c$$$ PFAy = 'COG4'!PFA |
c$$$ PFAy = 'COG4'!PFA |
| 1255 |
|
|
| 1256 |
|
|
| 1257 |
if(icx.gt.nclstr1.or.icy.gt.nclstr1)then |
if(icx.gt.nclstr1.or.icy.gt.nclstr1)then |
| 1258 |
print*,'xyzpam: ***WARNING*** clusters ',icx,icy |
print*,'xyzpam: ***WARNING*** clusters ',icx,icy |
| 1259 |
$ ,' does not exists (nclstr1=',nclstr1,')' |
$ ,' does not exists (nclstr1=',nclstr1,')' |
| 1266 |
call idtoc(pfaid,PFAx) |
call idtoc(pfaid,PFAx) |
| 1267 |
call idtoc(pfaid,PFAy) |
call idtoc(pfaid,PFAy) |
| 1268 |
|
|
| 1269 |
|
cc print*,PFAx,PFAy |
| 1270 |
|
|
| 1271 |
c$$$ call xyz_PAM(icx,icy,sensor,PFAx,PFAy,ax,ay,bfx,bfy) |
c$$$ call xyz_PAM(icx,icy,sensor,PFAx,PFAy,ax,ay,bfx,bfy) |
| 1272 |
|
|
| 1273 |
c$$$ print*,icx,icy,sensor,PFAx,PFAy,ax,ay,bfx,bfy |
c$$$ print*,icx,icy,sensor,PFAx,PFAy,ax,ay,bfx,bfy |