142 |
TY = TAN(ASIN(AL_P(3))) * SIN(AL_P(4)) |
TY = TAN(ASIN(AL_P(3))) * SIN(AL_P(4)) |
143 |
DO I = 1, NPLA |
DO I = 1, NPLA |
144 |
NN = 0 |
NN = 0 |
145 |
C????? IF (M.EQ.2) NN = 1 |
IF (M.EQ.2) NN = 1 |
146 |
IF (MOD(I,2).EQ.NN) THEN |
IF (MOD(I,2).EQ.NN) THEN |
147 |
SHIFT = +0.5 |
SHIFT = +0.5 |
148 |
ELSE |
ELSE |
253 |
C |
C |
254 |
DO M = 1,2 |
DO M = 1,2 |
255 |
C |
C |
256 |
|
c print *,' ax ',ax,' ay ',ay |
257 |
|
c print *,' bx ',bx,' by ',by |
258 |
IF (NPCFIT(M).GE.2) THEN |
IF (NPCFIT(M).GE.2) THEN |
259 |
IF (M.EQ.1) THEN |
IF (M.EQ.1) THEN |
260 |
VARCFIT(1) = CHI2X |
VARCFIT(1) = CHI2X |
261 |
IMPX = 10. * ( AX + 12.1 ) |
IMPX = AX ! PAMELA REF |
262 |
TANX = BX |
TANX = BX |
263 |
ELSE |
ELSE |
264 |
VARCFIT(2) = CHI2Y |
VARCFIT(2) = CHI2Y |
265 |
IMPY = 10. * ( AY + 12.2 ) |
IMPY = AY ! PAMELA REF |
266 |
TANY = BY |
TANY = BY |
267 |
ENDIF |
ENDIF |
268 |
C |
C |
269 |
DO I = 1,NPLA |
DO I = 1,NPLA |
270 |
NN = 0 |
NN = 0 |
271 |
C????? IF (M.EQ.2) NN = 1 |
IF (M.EQ.2) NN = 1 |
272 |
IF (MOD(I,2).EQ.NN) THEN |
IF (MOD(I,2).EQ.NN) THEN |
273 |
SHIFT = +0.5 |
SHIFT = +0.5 |
274 |
ELSE |
ELSE |
277 |
C |
C |
278 |
IF (M.EQ.1) THEN |
IF (M.EQ.1) THEN |
279 |
DISTX = PIANO(I) - 5.81 |
DISTX = PIANO(I) - 5.81 |
280 |
Y(I) = DISTX * TANX + CX |
Y(I) = (DISTX * TANX) + AX - XALIG |
|
BAR(M,I) = Y(I) |
|
281 |
c CBAR(M,I) = Y(I) |
c CBAR(M,I) = Y(I) |
282 |
CBAR(M,I) = (Y(I) - XALIG)/10. |
BAR(M,I) = Y(I) |
283 |
|
CBAR(M,I) = (Y(I) + XALIG)/10. |
284 |
IF (I.EQ.22) MX=ABS(Y(1)-Y(22))/ABS(ZIN(1)-ZIN(22)) |
IF (I.EQ.22) MX=ABS(Y(1)-Y(22))/ABS(ZIN(1)-ZIN(22)) |
285 |
C |
C |
286 |
ELSE |
ELSE |
287 |
DISTY = PIANO(I) |
DISTY = PIANO(I) |
288 |
YY(I) = DISTY * TANY + CY |
YY(I) = (DISTY * TANY) + AY - YALIG |
|
BAR(M,I) = YY(I) |
|
289 |
c CBAR(M,I) = YY(I) |
c CBAR(M,I) = YY(I) |
290 |
CBAR(M,I) = (-YALIG + YY(I))/10. |
BAR(M,I) = YY(I) |
291 |
|
CBAR(M,I) = (YY(I) + YALIG)/10. |
292 |
IF (I.EQ.22) MY=ABS(Y(1)-Y(22))/ABS(ZIN(1)-ZIN(22)) |
IF (I.EQ.22) MY=ABS(Y(1)-Y(22))/ABS(ZIN(1)-ZIN(22)) |
293 |
C |
C |
294 |
ENDIF |
ENDIF |
616 |
ENDIF |
ENDIF |
617 |
ENDDO |
ENDDO |
618 |
C |
C |
|
EINF = EMIN |
|
|
ESUP = 150. |
|
|
C |
|
|
IF (TRIGTY.GE.2.AND.HZN.NE.0) THEN |
|
|
EINF = 50. |
|
|
ESUP = 1500. |
|
|
ENDIF |
|
619 |
C |
C |
620 |
C CALCULATE PLANETOT AND QMEAN |
C CALCULATE PLANETOT AND QMEAN |
621 |
C |
C |
626 |
NPIANI = 5 |
NPIANI = 5 |
627 |
QMEAN = 0. |
QMEAN = 0. |
628 |
INDEX = 0 |
INDEX = 0 |
629 |
CALL ELIO(RPIANO,NPIANI,QMEAN,NTOT,INDEX) |
C |
630 |
PLANETOT = RPIANO(1) + RPIANO(2) |
IF (TRIGTY.GE.2.AND.HZN.NE.0) THEN |
631 |
|
EINF = 50. |
632 |
|
ESUP = 15000. |
633 |
|
CALL NUCLEI(RPIANO,NPIANI,QMEAN,NTOT,INDEX) |
634 |
|
PLANETOT = RPIANO(1) + RPIANO(2) |
635 |
|
ELSE |
636 |
|
EINF = EMIN |
637 |
|
ESUP = 15000. |
638 |
|
CALL ELIO(RPIANO,NPIANI,QMEAN,NTOT,INDEX) |
639 |
|
PLANETOT = RPIANO(1) + RPIANO(2) |
640 |
|
ENDIF |
641 |
C |
C |
642 |
50 CONTINUE |
50 CONTINUE |
643 |
C |
C |