mixiユーザー(id:2145604)

2012年12月24日10:04

500 view

近隣恒星相互の距離と明るさ

大学卒業の年、卒研でコンピューターを触れたのでお遊びで入力
プリントアウトの1部を柴野さんへ送ったところ野田宏一郎氏が喜びそうと廻されたとのこと…
--------------------------------------------------------------------
C____ MUTUAL DISTANCE OF NEARIST STARS
_____ DIMENSION STAR(15),H(15),DIG(15),YE(15),ZL(15)
_____ DIMENSION DIS(15,15),AKA(15,15),ITH(15,15),ITM(15,15)
_____ RADF(A,P,Q)=Q*P*(INTF(A)/180.+(A-INTF(A))/108.)
_____ P=3.1416
_____ READ 501,(STAR(I),H(I),DIG(I),YE(I),ZL(I),I=1,15)
__501 FORMAT(A8,F7.3,F7.2,F6.3,F6.2)
_____ DO 100 I=1,15
_____ PSC=1./YE(I)
_____ DIS(1,I)=PSC*3.26
_____ DIS(I,1)=DIS(1,I)
_____ FI=RADF(H(I),P,15.)
_____ THI=RADF(DIG(I),P,1.)
_____ IF(I-1)100,100,60
___60 K=I+1
_____ DO 300 J=K,15
_____ FJ=RADF(H(J),9,15.)
_____ THJ=RADF(DIG(J),P,1.)
_____ PAS=1./YE(J)
_____ XSUB=(SINF(THI)*PSC-SINF(THJ)*PAS)**2
_____ YSUB=(PSC*COSF(THI))**2+PAS*COSF(THJ))**2
_____ ZSUB=2.PSC*COSF(THI)*PAS*COSF(THJ)*COSF(FI-FJ)
_____ DIS(I,J)=SQRF(XSUB+YSUB-ZSUB)*3.26
_____ FAI=(FJ-FI)*12./P+24.
_____ TA=SQRTF(XSUB/(YSUB-ZSUB))
_____ THT=SIGNF(ATANF(TA),THJ-THI)*180./P
_____ DIC(I,J)=INTF(THT)+MODF(THT,1.)*0.6
_____ ITH(I,J)=XMODF(XFIXF(FAI),24)
_____ IHM(I,J)=MODF(FAI,1.)*60.+100.5
_____ DIC(J,I)=-DIC(I,J)
_____ IHT(J,I)=XMODF(ITH(I,J)+12,24)
_____ ITM(J,I)=ITM(I,J)
__300 CONTINE
__100 CONTINE
_____ DO 200 I=1,15
_____ DO 400 J=1,15
_____ IF(I-J)30,400,30
___30 AKA(I,J)=ZL(J)-5.+5.*LOGF(DIS(I,J)/3.26)/LOGF(10.)
__400 CONTINUE
_____ AKA(I,I)=1000.0003
_____ DIS(I,I)=1000.003
_____ DIC(I,1)=-DIG(I)
_____ DIC(1,I)=DIG(I)
_____ DIC(I,I)=1000.003
_____ ITH(1,I)=H(I)
_____ ITM(1,I)=MODF(H(I),1.)*100.+100.5
_____ ITM(I,1)=ITM(1,I)
_____ ITH(I,1)=XMODF(ITH(1,I)+36,24)
_____ ITH(I,I)=900
_____ ITM(I,I)=900
__200 CONTINUE
_____ PLINT 601 ,(I,I=1,15)
__601 FORMAT(1H1,50X,35H MUTUAL DISTANCE OF NEARIST STARS //1H ,11X,15(I5,3X))
_____ DO 102 I=1,15
_____ PLINT 602,I,STAR(I),(DIS(I,J),J=1,1)
__602 FORMAT(1H0,I2,1X,A8,,15(F5.1,3HLY ))
_____ PLINT 603,(AKA(I,J),J=1,15)
__603 FORMAT(1H ,11X,15(F6.2,2X))
__102 CONTINUE
_____ PLINT 601,(I,I=1,15)
_____ PLINT 604,(I,STAR(I),DIC(I,J),J=1,15),(ITH(I,J),ITM(I,J),J=1,15),I=1,15)
__604 FORMAT(1H0,I2,1X,A8,15(F6.2,2X)/1H 11X15(I3,1H:I2,2X))
_____ STOP
_____ END(1,0,0,1,0)
--------------------------------------------------------------------
僕が中学の頃に買った百科事典の記載から採取したデータが元でかなり怪しく今wikiなどに記載されているデータで再計算しなくては(笑)
計算結果でいま話題のタウケチからのデータを見ると
SOL_____ 10.9LY 2.5等星
A-CEN.__ 12.5LY 2.6等星
BARNARD_ 14.7LY 11.4等星
LUY726-8 3.0LY 10.4等星
WOLF-359 17.8LY 15.2等星
LA-21185 18.1LY 9.1等星
SIRIUS__ 11.6LY -0.8等星
ROSS-154 15.0LY 11.5等星
LUY789-6 8.0LY 11.9等星
ROSS-248 11.7LY 12.5等星
E-ERI.__ 5.1LY 2.2等星
ROSS-128 21.1LY 12.6等星
61-CYG._ 14.8LY 6.2等星
PROCYON_ 15.7LY 1.2等星
となってます。半角大文字英数字しか使えなかったので星名は略してます...
というか言語も死語のフォートランあせあせ(飛び散る汗)しかもJIS非準拠の大学独自方言(笑)


しかし、学生にコンピューターを(自由に)触らせると【今も昔も】遊びにしか使わないな(笑)
3 5

コメント

mixiユーザー

ログインしてコメントを確認・投稿する

<2012年12月>
      1
2345678
9101112131415
16171819202122
23242526272829
3031     

最近の日記