地震シュミレーションプログラム

ブラウザで表示させると何故かプログラムが欠落しますので、ソースを表示させで
保存、そしてＭＳ−ＤＯＳ　ＢＡＳＩＣ（Ｎ８８ー日本語ＢＡＳＩＣ（８６）６．
０）　にて実行させてください。（１００から７００までの６０行がプログラム本体。）
なおネットスケープではソースの編集ができないためプログラムの移植ができません。
ＩＥを使って下さい。またはコピーしてエディターで移植して下さい。

ＩＥまたはマイクロソフトのエディターではＢＡＳＩＣのソースを扱えないようで、
しばしば問題が発生しています。このプログラムでは行番号５７０のソースが原因
とみられる問題が発生しましたので５７０を削除してあります。前ページの５７０
を挿入して走らせて下さい。

プログラムの解説はプログラコードの最後にあります。/A>



100 REM EARTH   G426C  5/13       11/13 1996
110 DEFINT I,J,K,L,M,N  :SCREEN 3:CONSOLE ,,,1 :RP=1.1 :RESTORE 5710:COLOR 7
120 ND=49:HB=.4 :SP1=128!:SP2=28!:GP1=139!:GP2=35: MIRAI=100:GP=1/2:TP=1/1
130 P1=3.3E+11:P2=1.41      :DIP= 0 :STR=45:JX=(GP1-SP1)/HB:JY=(GP2-SP2)/HB
140 DIM AVE(ND,ND),BASE(ND,ND),BLOCK(ND,ND),GDATA(500,4),FLD(ND,ND)
150 DIM BLOCKA(ND,ND),SDMAXA(ND,ND),SDMAX(ND,ND)
160 DEF FNFS(X)=10^(1.02*X-4)*P2:DEF FNGME(X)=10^(11.8+1.5*X):CD=1E+10
170 DEF FNGTR(X,R)=2.32*2.72^(-X*X/(2*R*R)):P3=1/1
180 FOR I=1 TO 1000:FOR J=0 TO 4:READ GDATA(I,J):GDATA(I,J)=GDATA(I,J)/10
190 IF GDATA(I,J)=0 THEN GOTO *OW ELSE NEXT J,I
200 *OW:KOSU=I-1:KIKAN=(GDATA(KOSU,1)-GDATA(1,1))*10:IFY=GDATA(1,1)*10
210 ILY=GDATA(KOSU,1)*10:DP=.0175*DIP :AA=COS(DP):SR=.0175*STR
220 FOR I=1 TO KOSU:TPGE=0:GM=GDATA(I,4):EM=10^(1.5*GM+16.1):SM=FNFS(GM)
230 RG=SQR(SM/3.14)/100:SD=EM/(P1*SM*CD):PRINT GDATA(I,0)*10;
240 FOR K=0 TO ND:FOR L=0 TO  ND: GOSUB *DAS:IF JDS<0 THEN GOTO *FAR1
250 IF RSD<HB THEN RRG=(RG+RG*AA)/2  ELSE RRG=RSA*RG/RSD
260 RRG=RRG*P3:SE=SD*FNGTR(RSA,RRG)
270 FLD(K,L)=FLD(K,L)+1:BLOCK(K,L)=BLOCK(K,L)+SE
280 *FAR1:NEXT L,K,I
290 FOR K=0 TO ND:FOR L=0 TO ND:AVE(K,L)=BLOCK(K,L)/KIKAN:DLF=FLD(K,L)
300 BLOCK(K,L)=0:NEXT L,K     :I=1:PRINT:IW=1:IY=IFY:CLS 2
310 *RT:FOR K=0 TO ND:FOR L=0 TO ND:IF IW=1 THEN BLOCK(K,L)=BLOCK(K,L)+AVE(K,L)
320 NEXT L,K:JG=GDATA(I,1)*10 :IF JG <>IY THEN GOTO *GN
330 IF GDATA(I,1)=GDATA(I+1,1) THEN IY=IY-1:IW=10   ELSE IW=1
340 GM=GDATA(I,4):EM=10^(1.5*GM+16.1):PRINT JG;:COLOR 7
350 SM=FNFS(GM):RG=SQR(SM/3.14)/100:SD=EM/(P1*SM*CD)
360 FOR K=0 TO ND:FOR L=0 TO ND:GOSUB *DAS:IF JDS<0 THEN GOTO *FAR2
370 IF RSD<HB THEN RRG=(RG+RG*AA)/2  ELSE RRG=RSA*RG/RSD
380 RRG=RRG*P3:SE=SD*FNGTR(RSA,RRG)
390 BLOCK(K,L)=BLOCK(K,L)-SE:BASE(K,L)=BLOCK(K,L):COLOR 2
400 IF SDMAX(K,L)<SE THEN SDMAX(K,L)=SE
410 *FAR2:NEXT L,K      :I=I+1
420 *GN:YI=BLOCK(JX,JY)*GP:XI=(IY-IFY)*TP :LINE (XI,200)-(XI,YI+200),2
430 IY=IY+1:IF IY>ILY THEN GOTO *KO ELSE GOTO *RT
440 *KO:FOR K=0 TO ND:FOR L=0 TO ND:BLOCK(K,L)=BLOCK(K,L)-BASE(K,L):SDMAXA(K,L)=SDMAX(K,L):BLOCKA(K,L)=BLOCK(K,L):NEXT L,K:PK=SP1+JX*HB:COLOR 7
450 PI=SP2+JY*HB:PRINT"IDO and KEIDO of the point";USING "######.#"; PI;PK
460 PRINT" S OF NOW    MAX        AVE/YEAR   N"
470 PRINT USING"#######.#";BLOCK(JX,JY);SDMAX(JX,JY);AVE(JX,JY) ;FLD(JX,JY):
480 *NEX:GOSUB *YOST: PRINT "FROM";ILY;"TO";ILY+MIRAI; "IDO=";SP2;"-";SP2+ND*HB;"  KEIDO=";SP1;"-";SP1+ND*HB;"   RP=";RP
490 PRINT " YEAR       IDO    KEIDO     MG1     MG2     LENGTH of FAULT (Km) "
500 FOR IY =ILY+1 TO ILY+MIRAI
510 FOR K=0 TO ND:FOR L=0 TO ND:BLOCK(K,L)=BLOCK(K,L)+AVE(K,L):NEXT L,K
520 FOR KK=0 TO ND:FOR LL=0 TO ND:   IF BLOCK(KK,LL)>SDMAX(KK,LL) AND FLD(KK,LL)<1.5 THEN  BLOCK(KK,LL)=SDMAX(KK,LL)
530 IF BLOCK(KK,LL)>SDMAX(KK,LL) AND FLD(KK,LL)>1.5 THEN  XX=KK*HB:YY=LL*HB:SD=BLOCK(KK,LL) ELSE GOTO *BUGI
540 GM=(.434*LOG(SD/2.32*P1*CD*P2)-20.1)/.48:RG=(SQR(FNFS(GM)/3.14))/100:GMT=GM
550 TGE=0:FOR K=0 TO ND:FOR L=0 TO ND:DX=K*HB-XX:DY=L*HB-YY:GOSUB *DAS2
560 IF JDS<0 THEN GOTO *FAR3
570 IF RSD<HB THEN RRG=(RG+RG*AA)/2  ELSE RRG=RSA*RG/RSD
580 RRG=RRG*P3: POE=FNGTR(RSA,RRG)/2.32   :SE=SD        *POE
590 PGE=BLOCK(K,L)*P1*HB*HB*1E+14*POE/AA:TGE=TGE+PGE:BLOCK(K,L)=BLOCK(K,L)-SE
600  IF BLOCK(K,L)<0 THEN BLOCK(K,L)=0
610 *FAR3:NEXT L,K
620 GM=(.434*LOG(TGE)-16.1)/1.5:PKEIDO=SP1+XX: PIDO=SP2+YY
630 PRINT USING "######.#";IY,PIDO;PKEIDO;GM;GMT;:PRINT USING"########";RG*200
640 *BUGI:NEXT LL,KK:YI=(BLOCK(JX,JY)+BASE(JX,JY))*GP:XI=(IY-IFY)*TP
650 LINE(XI,200)-(XI,YI+200),KC:NEXT IY:GOTO *NEX                         :END
660 *YOST:INPUT "RP=";RP:FOR K=0 TO ND:FOR L=0 TO ND:SDMAX(K,L)=SDMAXA(K,L)*RP
670 BLOCK(K,L)=BLOCKA(K,L):NEXT L,K:KC=KC+1:RETURN
680 *DAS:PX=SP1+K*HB:PY=SP2+L*HB:X=GDATA(I,2):Y=GDATA(I,3):DX=PX-X:DY=PY-Y
690 *DAS2:RSA=SQR(DX*DX+DY*DY):RR=((DX*COS(SR)-DY*SIN(SR))*(DX*COS(SR)-DY*SIN(SR))/(AA*AA)+(DX*SIN(SR)+DY*COS(SR))*(DX*SIN(SR)+DY*COS(SR))):RSD=SQR(RR):IF RSD<=RG THEN JDS=1:RETURN ELSE IF RSD>RG THEN JDS=-1:RETURN
700 *ZERO:FOR K=0 TO ND:FOR L=0 TO ND:AVE(K,L)=0:BASE(K,L)=0:BLOCK(K,L)=0:FLD(K,L)=0:BLOCKA(K,L)=0:SDMAX(K,L)=0:NEXT L,K:RETURN


5000 DATA 1, 416,1000,300,70
5010 DATA 2, 599,1000,300,70
5020 DATA 3, 679,1000,300,70
5030 DATA 4, 684,1000,300,82
5040 DATA 5, 701,1000,300,70
5050 DATA 6, 715,1378,351,70
5060 DATA 7, 715,1374,348,68
5070 DATA 8, 734,1000,300,70
5080 DATA 9, 745,1365,354,79
5090 DATA 10,762,1000,300,70
5100 DATA 11,818,1000,300,76
5110 DATA 12,827,1358,350,68
5120 DATA 13,830,1401,398,73
5130 DATA 14,841,1380,362,66
5140 DATA 15,841,1389,351,70
5150 DATA 16,850,1397,390,70
5160 DATA 17,856,1000,300,63
5170 DATA 18,863,1000,300,70
5180 DATA 19,868,1348,348,71
5190 DATA 20,869,1000,300,83
5200 DATA 21,878,1393,355,74
5210 DATA 22,880,1332,354,70
5220 DATA 23,881,1000,300,64
5230 DATA 24,887,1350,330,83
5240 DATA 25,890,1000,300,60
5250 DATA 26,934,1000,300,60
5260 DATA 27,938,1358,350,70
5270 DATA 28,976,1358,349,68
5280 DATA 29,1038,1356,343,70
5290 DATA 30,1041,1000,300,70
5300 DATA 31,1070,1358,348,63
5310 DATA 32,1091,1358,347,64
5320 DATA 33,1093,1000,300,62
5330 DATA 34,1096,1000,300,83
5340 DATA 35,1099,1000,300,82
5350 DATA 36,1177,1358,347,63
5360 DATA 37,1185,1358,350,74
5370 DATA 38,1213,1000,300,70
5380 DATA 39,1227,1000,300,70
5390 DATA 40,1230,1000,300,70
5400 DATA 41,1240,1000,300,70
5410 DATA 42,1241,1000,300,70
5420 DATA 43,1245,1000,300,70
5430 DATA 44,1257,1395,352,73
5440 DATA 45,1293,1000,300,70
5450 DATA 46,1317,1358,350,68
5460 DATA 47,1325,1361,356,65
5470 DATA 48,1331,1352,337,71
5480 DATA 49,1350,1358,350,60
5490 DATA 50,1360,1362,334,78
5500 DATA 51,1361,1000,300,70
5510 DATA 52,1361,1350,330,84
5520 DATA 53,1408,1360,330,75
5530 DATA 54,1425,1358,350,60
5540 DATA 55,1433,1395,349,71
5550 DATA 56,1449,1358,350,62
5560 DATA 57,1456,1000,300,70
5570 DATA 58,1466,1000,300,70
5580 DATA 59,1494,1357,346,60
5590 DATA 60,1498,1323,330,73
5600 DATA 61,1498,1380,340,83
5610 DATA 62,1502,1382,372,68
5620 DATA 63,1510,1356,346,68
5630 DATA 64,1517,1000,300,70
5640 DATA 65,1520,1360,330,74
5650 DATA 66,1525,1000,300,70
5660 DATA 67,1579,1355,347,60
5670 DATA 68,1586,1369,360,78
5680 DATA 69,1589,1382,348,67
5690 DATA 70,1596,1316,333,70
5700 DATA 71,1596,1356,347,75
5710 DATA 72,1605,1385,335,79
5720 DATA 72,1605,1349,330,79
5730 DATA 73,1611,1398,376,69
5740 DATA 74,1611,1444,390,81
5750 DATA 75,1614,1000,300,70
5760 DATA 76,1615,1397,357,65
5770 DATA 77,1616,1420,381,70
5780 DATA 78,1619,1306,325,60
5790 DATA 79,1625,1000,300,70
5800 DATA 80,1625,1306,328,55
5810 DATA 81,1627,1382,366,60
5820 DATA 82,1628,1000,300,60
5830 DATA 83,1630,1398,358,63
5840 DATA 84,1633,1392,352,70
5850 DATA 85,1635,1398,358,60
5860 DATA 86,1640,1407,421,70
5870 DATA 87,1640,1362,363,65
5880 DATA 88,1644,1400,394,65
5890 DATA 89,1646,1407,381,66
5900 DATA 90,1646,1000,300,70
5910 DATA 91,1647,1000,300,65
5920 DATA 92,1648,1392,352,70
5930 DATA 93,1649,1325,337,70
5940 DATA 94,1649,1395,358,70
5950 DATA 95,1649,1397,355,64
5960 DATA 96,1650,1000,300,63
5970 DATA 97,1659,1398,371,69
5980 DATA 98,1662,1360,352,74
5990 DATA 99,1662,1320,317,77
6000 DATA 100,1664,1000,300,59
6010 DATA 101,1664,1000,300,70
6020 DATA 102,1664,1000,300,70
6030 DATA 103,1665,1000,300,60
6040 DATA 104,1666,1382,371,67
6050 DATA 105,1667,1000,300,70
6060 DATA 106,1667,1416,406,62
6070 DATA 107,1668,1000,300,70
6080 DATA 108,1668,1000,300,59
6090 DATA 109,1670,1394,377,67
6100 DATA 110,1671,1000,300,70
6110 DATA 111,1674,1416,406,60
6120 DATA 112,1676,1318,345,65
6130 DATA 113,1677,1423,410,74
6140 DATA 114,1677,1420,355,80
6150 DATA 115,1678,1425,390,75
6160 DATA 116,1683,1396,367,63
6170 DATA 117,1683,1397,368,68
6180 DATA 118,1683,1397,369,70
6190 DATA 119,1685,1000,300,65
6200 DATA 120,1686,1326,340,72
6210 DATA 121,1686,1374,347,68
6220 DATA 122,1694,1401,402,70
6230 DATA 123,1694,1000,300,70
6240 DATA 124,1696,1000,300,70
6250 DATA 125,1697,1396,354,65
6260 DATA 126,1698,1315,331,60
6270 DATA 127,1700,1296,339,70
6280 DATA 128,1703,1314,333,65
6290 DATA 129,1703,1398,347,81
6300 DATA 130,1704,1400,404,70
6310 DATA 131,1705,1312,330,70
6320 DATA 132,1706,1399,386,58
6330 DATA 133,1706,1398,356,58
6340 DATA 134,1707,1359,332,84
6350 DATA 135,1708,1000,300,70
6360 DATA 136,1710,1415,370,65
6370 DATA 137,1710,1337,355,65
6380 DATA 138,1711,1338,352,63
6390 DATA 139,1714,1379,367,63
6400 DATA 140,1715,1366,354,68
6410 DATA 141,1717,1425,385,75
6420 DATA 142,1718,1376,353,70
6430 DATA 143,1723,1306,329,65
6440 DATA 144,1725,1397,363,60
6450 DATA 145,1725,1364,364,60
6460 DATA 146,1725,1381,360,63
6470 DATA 147,1729,1000,300,70
6480 DATA 148,1729,1371,374,68
6490 DATA 149,1731,1406,380,65
6500 DATA 150,1731,1000,300,70
6510 DATA 151,1733,1000,300,66
6520 DATA 152,1735,1000,300,70
6530 DATA 153,1736,1000,300,60
6540 DATA 154,1738,1387,370,55
6550 DATA 155,1739,1000,300,70
6560 DATA 156,1741,1394,416,10
6570 DATA 157,1746,1000,300,70
6580 DATA 158,1749,1326,332,68
6590 DATA 159,1751,1358,350,58
6600 DATA 160,1751,1382,371,72
6610 DATA 161,1755,1000,300,70
6620 DATA 162,1755,1396,368,70
6630 DATA 163,1756,1409,357,58
6640 DATA 164,1760,1000,300,70
6650 DATA 165,1762,1000,300,70
6660 DATA 166,1762,1387,381,70
6670 DATA 167,1763,1423,410,74
6680 DATA 168,1763,1420,410,73
6690 DATA 169,1763,1420,410,70
6700 DATA 170,1766,1405,407,73
6710 DATA 171,1767,1000,300,70
6720 DATA 172,1767,1398,357,60
6730 DATA 173,1768,1275,262,70
6740 DATA 174,1768,1000,300,70
6750 DATA 175,1769,1416,406,70
6760 DATA 176,1769,1321,330,78
6770 DATA 177,1771,1243,240,74
6780 DATA 178,1772,1435,395,74
6790 DATA 179,1778,1320,346,65
6800 DATA 180,1780,1399,389,50
6810 DATA 181,1782,1391,354,70
6820 DATA 182,1786,1391,352,53
6830 DATA 183,1789,1343,337,70
6840 DATA 184,1791,1396,358,63
6850 DATA 185,1791,1380,362,68
6860 DATA 186,1792,1303,328,64
6870 DATA 187,1792,1400,438,71
6880 DATA 188,1793,1315,341,64
6890 DATA 189,1793,1400,409,70
6900 DATA 190,1793,1445,385,82
6910 DATA 191,1794,1000,300,70
6920 DATA 192,1796,1343,357,55
6930 DATA 193,1799,1366,366,60
6940 DATA 194,1801,1401,353,70
6950 DATA 195,1802,1365,352,68
6960 DATA 196,1802,1384,378,68
6970 DATA 197,1804,1400,391,70
6980 DATA 198,1810,1399,399,65
6990 DADA 199,1811,1000,300,10
7000 DATA 200,1812,1335,335,70
7010 DATA 201,1812,1397,355,63
7020 DATA 202,1815,1365,364,60
7030 DATA 203,1817,1391,352,60
7040 DATA 204,1819,1363,352,73
7050 DATA 205,1821,1000,300,70
7060 DATA 206,1821,1396,375,58
7070 DATA 207,1823,1411,400,59
7080 DATA 208,1826,1373,362,60
7090 DATA 209,1828,1299,326,60
7100 DATA 210,1828,1389,376,69
7110 DATA 211,1830,1356,351,65
7120 DATA 212,1831,1303,332,61
7130 DATA 213,1832,1416,407,65
7140 DATA 214,1833,1366,355,63
7150 DATA 215,1833,1393,389,75
7160 DATA 216,1834,1414,433,64
7170 DATA 217,1835,1326,351,55
7180 DATA 218,1835,1425,385,70
7190 DATA 219,1836,1392,344,55
7200 DATA 220,1839,1000,300,70
7210 DATA 221,1841,1385,350,63
7220 DATA 222,1841,1324,332,60
7230 DATA 223,1842,1000,300,70
7240 DATA 224,1843,1391,354,65
7250 DATA 225,1843,1460,420,75
7260 DATA 226,1844,1313,330,50
7270 DATA 227,1847,1000,300,70
7280 DATA 228,1847,1382,367,74
7290 DATA 229,1847,1383,372,65
7300 DATA 230,1848,1304,332,59
7310 DATA 231,1848,1406,407,60
7320 DATA 232,1848,1307,329,50
7330 DATA 233,1853,1381,366,65
7340 DATA 234,1853,1392,353,67
7350 DATA 235,1854,1360,348,73
7360 DATA 236,1854,1416,406,65
7370 DATA 237,1854,1378,340,84
7380 DATA 238,1854,1350,330,84
7390 DATA 239,1854,1320,333,74
7400 DATA 240,1855,1000,300,70
7410 DATA 241,1855,1369,363,67
7420 DATA 242,1855,1000,300,70
7430 DATA 243,1855,1420,381,73
7440 DATA 244,1855,1378,345,73
7450 DATA 245,1855,1398,357,69
7460 DATA 246,1856,1423,410,75
7470 DATA 247,1856,1395,357,63
7480 DATA 248,1857,1314,344,60
7490 DATA 249,1857,1382,348,63
7500 DATA 250,1857,1325,340,73
7510 DATA 251,1858,1372,364,71
7520 DATA 252,1858,1000,300,50
7530 DATA 253,1858,1379,366,57
7540 DATA 254,1858,1420,408,73
7550 DATA 255,1858,1000,300,50
7560 DATA 256,1858,1408,409,60
7570 DATA 257,1859,1319,348,62
7580 DATA 258,1859,1397,359,60
7590 DATA 259,1859,1320,345,63
7600 DATA 260,1861,1412,386,64
7610 DATA 261,1864,1000,300,63
7620 DATA 262,1866,1000,300,70
7630 DATA 263,1870,1391,353,63
7640 DATA 264,1872,1321,352,71
7650 DATA 265,1880,1398,354,58
7660 DATA 266,1881,1473,433,70
7670 DATA 267,1882,1000,300,50
7680 DATA 268,1884,1398,357,50
7690 DATA 269,1886,1385,371,53
7700 DATA 270,1889,1307,328,63
7710 DATA 271,1890,1380,365,62
7720 DATA 272,1890,1393,342,68
7730 DATA 273,1891,1366,356,80
7740 DATA 274,1892,1399,357,62
7750 DATA 275,1892,1367,371,64
7760 DATA 276,1893,1480,435,78
7770 DATA 277,1893,1305,314,53
7780 DATA 278,1894,1460,425,79
7790 DATA 279,1894,1398,357,70
7800 DATA 280,1894,1310,328,63
7810 DATA 281,1894,1399,389,70
7820 DATA 282,1895,1404,361,72
7830 DATA 283,1895,1310,328,63
7840 DATA 284,1896,1410,365,70
7850 DATA 285,1896,1373,375,57
7860 DATA 286,1896,1440,395,85
7870 DATA 287,1896,1407,395,72
7880 DATA 288,1897,1383,367,52
7890 DATA 289,1897,1419,381,74
7900 DATA 290,1897,1433,383,77
7910 DATA 291,1898,1312,346,62
7920 DATA 292,1898,1420,386,72
7930 DATA 293,1898,1389,370,61
7940 DATA 294,1898,1302,336,60
7950 DATA 295,1898,1247,245,70
7960 DATA 296,1899,1361,341,70
7970 DATA 297,1899,1320,319,71
7980 DATA 297,1899,1323,327,69
7990 DATA 298,1900,1362,358,58
8000 DATA 299,1900,1411,387,70
8010 DATA 300,1900,1394,339,66
8020 'DATA 0,0,0,0,0
8030 DATA 301,1901,1300,280,75
8040 DATA 302,1901,1425,405,72
8050 DATA 303,1902,1413,405,70
8060 DATA 304,1903,1376,362,55
8070 DATA 305,1904,1389,371,61
8080 DATA 306,1905,1325,341,73
8090 DATA 307,1905,1393,348,58
8100 DATA 308,1909,1415,345,67
8110 DATA 308,1909,1415,345,75
8120 DATA 309,1909,1363,354,68
8130 DATA 310,1909,1280,260,62
8140 DATA 311,1909,1311,323,76
8150 DATA 312,1910,1409,425,51
8160 DATA 313,1910,1416,442,53
8170 DATA 314,1911,1300,280,80
8180 DATA 315,1912,1383,364,51
8190 DATA 316,1913,1303,316,57
8200 DATA 317,1914,1306,316,71
8210 DATA 318,1914,1404,395,71
8220 DATA 319,1914,1404,392,61
8230 DATA 320,1915,1436,421,70
8240 DATA 321,1915,1403,354,60
8250 DATA 322,1916,1385,365,62
8260 DATA 323,1916,1350,346,61
8270 DATA 324,1917,1381,350,63
8280 DATA 325,1918,1520,455,80
8290 DATA 326,1918,1505,445,77
8300 DATA 327,1918,1379,365,61
8310 DATA 327,1918,1379,365,65
8320 DATA 328,1921,1402,360,70
8330 DATA 329,1922,1398,352,68
8340 DATA 330,1922,1301,327,69
8350 DATA 330,1922,1301,327,65
8360 DATA 331,1923,1312,306,71
8370 DATA 332,1923,1395,351,79
8380 DATA 333,1923,1389,355,68
8390 DATA 334,1924,1392,355,73
8400 DATA 335,1925,1348,356,68
8410 DATA 336,1927,1352,355,73
8420 DATA 337,1927,1387,376,52
8430 DATA 338,1930,1406,366,65
8440 DATA 339,1930,1391,351,73
8450 DATA 340,1931,1392,362,69
8460 DATA 341,1931,1326,323,71
8470 DATA 342,1933,1445,392,81
8480 DATA 343,1933,1368,371,60
8490 DATA 344,1934,1370,356,63
8500 DATA 345,1935,1384,350,64
8510 DATA 346,1936,1357,346,64
8520 DATA 347,1936,1421,382,75
8530 DATA 348,1936,1390,344,63
8540 DATA 349,1938,1351,336,68
8550 DATA 350,1938,1416,367,70
8560 DATA 351,1938,1445,436,61
8570 DATA 352,1938,1252,253,67
8580 DATA 353,1938,1422,373,75
8590 DATA 354,1939,1320,323,65
8600 DATA 355,1939,1395,401,68
8610 DATA 356,1940,1395,443,75
8620 DATA 357,1941,1382,367,61
8630 DATA 358,1941,1321,320,72
8640 DATA 359,1943,1342,354,62
8650 DATA 360,1943,1399,373,62
8660 DATA 361,1943,1341,355,72
8670 DATA 362,1943,1381,368,59
8680 DATA 363,1944,1366,338,79
8690 DATA 364,1945,1371,347,68
8700 DATA 365,1946,1356,330,80
8710 DATA 366,1947,1232,247,74
8720 DATA 367,1947,1410,438,67
8730 DATA 368,1948,1354,338,67
8740 DATA 369,1948,1362,362,71
8750 DATA 370,1949,1327,341,62
8760 DATA 371,1949,1396,367,62
8770 DATA 371,1949,1398,366,64
8780 DATA 372,1952,1441,418,82
8790 DATA 373,1952,1362,365,65
8800 DATA 374,1952,1358,345,68
8810 DATA 375,1953,1417,340,74
8820 DATA 376,1955,1343,337,64
8830 DATA 377,1955,1402,403,59
8840 DATA 378,1956,1441,443,63
8850 DATA 379,1956,1406,380,60
8860 DATA 380,1956,1401,357,63
8870 DATA 381,1957,1393,342,60
8880 DATA 382,1958,1485,443,81
8890 DATA 383,1959,1446,433,63
8900 DATA 383,1959,1445,434,61
8910 DATA 384,1960,1434,398,72
8920 DATA 385,1961,1388,375,52
8930 DATA 386,1961,1319,316,70
8940 DATA 387,1961,1456,429,72
8950 DATA 388,1961,1368,360,70
8960 DATA 389,1962,1439,422,70
8970 DATA 390,1962,1411,387,65
8980 DATA 391,1962,1395,341,59
8990 DATA 392,1963,1358,358,69
9000 DATA 393,1963,1500,438,81
9010 DATA 394,1964,1390,403,69
9020 DATA 395,1964,1392,384,75
9030 DATA 396,1965,1383,349,61
9040 DATA 397,1965,1000,300,64
9050 DATA 398,1967,1392,342,53
9060 DATA 399,1968,1307,320,61
9070 DATA 400,1968,1325,323,75
9080 DATA 401,1968,1436,407,79
9090 DATA 402,1968,1394,360,61
9100 DATA 403,1968,1324,333,66
9110 DATA 404,1968,1383,368,53
9120 DATA 405,1969,1371,358,66
9130 DATA 406,1970,1292,284,61
9140 DATA 407,1970,1431,424,67
9150 DATA 408,1970,1320,321,67
9160 DATA 409,1970,1408,392,62
9170 DATA 410,1971,1384,371,55
9180 'DATA 0,0,0,0,0
9190 DATA 411,1972,1413,332,71
9200 DATA 412,1972,1411,332,72
9210 DATA 413,1973,1460,430,74
9220 DATA 414,1974,1388,346,69
9230 DATA 415,1975,1311,330,61
9240 DATA 416,1975,1313,331,64
9250 DATA 417,1976,1390,355,55
9260 DATA 418,1978,1393,348,70
9270 DATA 419,1978,1327,351,61
9280 DATA 420,1978,1422,382,74
9290 DATA 421,1980,1392,349,67
9300 DATA 422,1982,1426,421,71
9310 DATA 423,1983,1391,404,77
9320 DATA 424,1983,1390,355,60
9330 DATA 425,1984,1322,324,71
9340 DATA 426,1984,1376,358,68
9350 DATA 427,1987,1321,320,66
9360 DATA 428,1987,1405,354,67
9370 DATA 429,1993,1444,429,78
9380 DATA 430,1993,1392,428,78
9390 DATA
9400 DATA
9410 DATA
9420 DATA
9430 DATA
9440 DATA 0,0,0,0
9450 FOR L=ND TO 0 STEP -1:PRINT:FOR K=0 TO ND
9460 PRINT USING"###";FLD(K,L);:NEXT:NEXT:STOP
9470 DATA      0,0,0,0,0
　　

　　　　　　　　　　　　　　　　　　　　　　　　プログラム　終了

<A NAME="解説"></A>
<FONT SIZE=5>プログラムの解説</FONT>
北緯２８ー４８度、東経１２８ー１４８度、の範囲を５０Ｘ５０のメッシュ（ブロ
ック）に分け、１９９６年からの１００年について、西暦１６０５年からの地震記録
（理科年表）をデータとして予想（シュミレーション）するように設定されています。
また（３５、１３９）の地点のブロックの歪みをグラフ表示します。
２時間ほどで過去のデータ処理が終わりますので、１　と入力してください。以後未
来のシュミレーションにおいて、過去の破壊限界平均値の１００パーセントまで歪み
が溜ると地震が起きます。０．９と入力すると９０パーセントで地震が起きます。
　


始めにデーターを読み込み、配列ＧＤＡＴＡに入れます。次に処理します。処理済
みデータは７２　７３　７４．．．．と表示していきます。この部分、メモリー及
び計算時間の節約のためにブラックボックス化してしまっています。不規則に起こ
っている地震を処理して配列に格納するのは３次元配列を使用すれば容易ですが、
パソコンではメモリー不足です。地震記録データにダミーを入れて整えるという方
法もあるでしょう。
次に過去の地震データの地震をシュミレーションすることにより、各ブロックの毎
年溜るであろう歪み量（配列ＡＶＥ）から、破壊限界の歪み量（配列ＳＤＭＡＸ）
などの推定値を計算します。地下３０キロメートルの地震断層で破壊が起こるとし
て、こう性率を（３．３Ｅ＋１１）として歪み量はセンチメートルで表します。
処理済みデータの西暦を表示し計算範囲内の場合は赤で表示します。また始めに指
定した地点の歪み量をグラフィック表示していきます。
次に得られたＳＤＭＡＸやＡＶＥを用いて未来の地震をシュミレートしていきます。
ＲＰ＝　？　と聞いてきますので、１　と入力します。これは破壊限界の歪み量（
過去の歪み量の最大値）の１００パーセントに達したら地震を起こすという意味で、
ＲＰの値を変えて何回でも未来の地震を起こしてみることができるようになってい
ます。０．８と入力すれば破壊限界の８０パーセントで地震が起こりますので、よ
り早期に規模が小さい地震が多数起こります。
ＮＥＣインタープレタ・ベーシック、１００ＭＨＺのＣＰＵを用いた場合３０分程
かかります。

得られる地震予報の精度は例えば次のようにして推定します。西暦１５００年から
１９００年までの４００年分のデータを用いて、１９０１年から１９９６年までを
シュミレート予報してみる。そして現実に１９０１年から１９９６年までに起こっ
た地震と比較して、予報誤差を推定・計算する。
見事に予報が当たる地域も有りますし、全くダメな地域もあります。このプログラ
ムでは地震断層の走行・傾斜も指定できるようになっていますが、日本全国一律の
理論、パラメーター、走行・傾斜、地震データではなく、その地域の特性を考慮に
いれた地震モデル・プログラム・及びデータを構築する必要があります。

なおデータとして始めはサンプルデータを用いてプログラムのチェックを行うわ
けですが、自然地震データでなけれな自然な地震は起こりません。このプログラム
では前震や余震も起こりますが本震だけをデータに入れた場合は、始めは本震だけ
が起こります。