地震のマグニチュード、エネルギー、原爆相当エネルギー、マグニチュードの差、エネルギーの差

\textup{log}_{10}E=4.8+1.5 \times M
M:マグニチュード , E:エネルギー(J)

参考資料: https://ja.wikipedia.org/wiki/%E3%82%A8%E3%83%8D%E3%83%AB%E3%82%AE%E3%83%BC%E3%81%AE%E6%AF%94%E8%BC%83

地震のマグニチュード-相当エネルギー
マグニチュード エネルギー(ペタジュール) 広島型原爆x発分相当 ツァーリ・ボンバx発分相当
6.00 0.06 1.15 0.00
6.10 0.09 1.62 0.00
6.20 0.13 2.29 0.00
6.30 0.18 3.23 0.00
6.40 0.25 4.57 0.00
6.50 0.35 6.45 0.00
6.60 0.50 9.11 0.00
6.70 0.71 12.87 0.00
6.80 1.00 18.18 0.00
6.90 1.41 25.68 0.01
7.00 2.00 36.28 0.01
7.10 2.82 51.24 0.01
7.20 3.98 72.38 0.02
7.30 5.62 102.24 0.03
7.40 7.94 144.42 0.04
マグニチュード エネルギー(ペタジュール) 広島型原爆x発分相当 ツァーリ・ボンバx発分相当
7.50 11.22 204.00 0.05
7.60 15.85 288.16 0.07
7.70 22.39 407.04 0.11
7.80 31.62 574.96 0.15
7.90 44.67 812.15 0.21
8.00 63.10 1147.19 0.30
8.10 89.13 1620.46 0.42
8.20 125.89 2288.95 0.60
8.30 177.83 3233.24 0.85
8.40 251.19 4567.07 1.20
8.50 354.81 6451.15 1.69
8.60 501.19 9112.50 2.39
8.70 707.95 12871.74 3.37
8.80 1000.00 18181.82 4.76
8.90 1412.54 25682.50 6.73
9.00 1995.26 36277.50 9.50
マグニチュード エネルギー(ペタジュール) 広島型原爆x発分相当 ツァーリ・ボンバx発分相当

地震のマグニチュードの差-相当エネルギーの差
マグニチュードの差(単位:マグニチュード) エネルギーの差(単位:倍)
0.10 1.41
0.20 2.00
0.30 2.82
0.40 3.98
0.50 5.62
0.60 7.94
0.70 11.22
0.80 15.85
0.90 22.39
1.00 31.62
1.10 44.67
1.20 63.10
1.30 89.13
1.40 125.89
マグニチュードの差(単位:マグニチュード) エネルギーの差(単位:倍)
1.50 177.83
1.60 251.19
1.70 354.81
1.80 501.19
1.90 707.95
2.00 1000.00
2.10 1412.54
2.20 1995.26
2.30 2818.38
2.40 3981.07
2.50 5623.41
2.60 7943.28
2.70 11220.18
2.80 15848.93
2.90 22387.21
3.00 31622.78
マグニチュードの差(単位:マグニチュード) エネルギーの差(単位:倍)

Script
library(xtable)
username<-Sys.info()['user']
pathOutput<-paste("C:/Users/",username,"/Desktop/R_Data_Write/",sep="")
setwd(pathOutput)
M<-seq(6,9.0,by=0.1)
E<-10^(4.8+1.5*M) # MJ:10^6,GJ:10^9,TJ:10^12,PJ:10^15,EJ(Exa):10^18
atomicbomb.HIROSHIMA<-55*10^12
atomicbomb.NAGASAKI<-84*10^12
atomicbomb.TsarBomba<-210*10^15
hiroshima<-round(E/atomicbomb.HIROSHIMA,3)
nagasaki<-round(E/atomicbomb.NAGASAKI,1)
tsar<-round(E/atomicbomb.TsarBomba,3)
dataset01<-data.frame(M,round(E/(10^15),2),hiroshima,tsar)
colnames(dataset01)<-c("マグニチュード","エネルギー(ペタジュール)","広島型原爆x発分相当","ツァーリ・ボンバx発分相当")
diffMagnitude<-seq(0.1,3,by=0.1)
diffEnergy<-10^(1.5*diffMagnitude)
dataset02<-data.frame(diffMagnitude,round(diffEnergy,2))
colnames(dataset02)<-c("マグニチュードの差(単位:マグニチュード)","エネルギーの差(単位:倍)")
htmlFile<-"magnitudeTable.html"
cat("",file=htmlFile,append=F)
cat("<b>地震のマグニチュード-相当エネルギー</b>",file=htmlFile,append=T)
print(xtable(dataset01),file=htmlFile,append=T,type="html",include.rownames=F)
cat("<hr>",file=htmlFile,append=T)
cat("<b>地震のマグニチュードの差-相当エネルギーの差</b>",file=htmlFile,append=T)
print(xtable(dataset02),file=htmlFile,append=T,type="html",include.rownames=F)
#library(DT)
#saveWidget(datatable(dataset01,options=list(autoWidth=T, paging=F,searching=F,info=F),rownames=F),"dataset01.html")
#saveWidget(datatable(dataset02,options=list(autoWidth=T, paging=F,searching=F,info=F),rownames=F),"dataset02.html")