안녕하세요 C언어 쌩초보입니다... 선배가 C++ 6.0 에서 짜놓은 코드로 visual studio 2015에서 해석을 돌리려 하는데 <div><br></div> <div>자꾸만 오류가 발생해서 일주일 째 골머리르 썩고 있습니다. 부디 지식 나눔을 해주시어 답을 내려주세요 ㅠㅠ</div> <div><br></div> <div>디버깅 에러문구 '<span style="color:#404040;font-family:gulim, sans-serif;font-size:13px;line-height:22.1px;">.exe : 0xC0000005: 0x00000000 위치를 읽는 동안 액세스 위반이 발생했습니다.'</span></div> <div><span style="color:#404040;font-family:gulim, sans-serif;font-size:13px;line-height:22.1px;">릴리즈 에러문구 '</span><span style="color:#404040;font-family:gulim, sans-serif;font-size:13px;line-height:22.1px;">0xC0000005: 0xFFFFFFF8 위치를 읽는 동안 액세스 위반이 발생했습니다.'</span></div> <div><span style="color:#404040;font-family:gulim, sans-serif;font-size:13px;line-height:22.1px;"><br></span></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">코드는 조금 길 수 있지만 고수분들의 실력이라면 가능할 거라 믿겠습니다...</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">F10으로 디버깅 실시했을떄 빨간색 부분부터 에러가 발생하였습니다.</span></font></div> <div><span style="color:#404040;font-family:gulim, sans-serif;font-size:13px;line-height:22.1px;"><br></span></div> <div><span style="color:#404040;font-family:gulim, sans-serif;font-size:13px;line-height:22.1px;"><br></span></div> <div><div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">// AB_MAP.cpp : Defines the entry point for the console application.</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">//</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">#include "stdafx.h"</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">#include <stdio.h></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">#include <stdlib.h></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">#include <iostream></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">#include <string.h></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">#include <fstream></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">#include <math.h></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">#include <time.h></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">int num_lon, num_lat, num_M;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">int control_b = 10;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">double W,E,S,N,dD;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">double Ml,Mu,dM;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">double aa[91][101];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">double bb[91][101];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">double probability(double[5], double, double, double, double, double*);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">int return_p[4];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">int iteration(int);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">void main()</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>int i, j, m, ii;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>int iter;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double del_b;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double lat, lon, b, RMSE;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double a10, v_old, v_new;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double RMSE_a, RMSE_b;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double* original1;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double* original2;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double* original3;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double* result1;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double* result2;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double* result3;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double* sm_b;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>char temp[80];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>return_p[1] = 500; return_p[2] =1000; return_p[3] = 2400;</span></font></div> <div><span class="Apple-tab-span" style="white-space:pre;font-size:13px;line-height:22.1px;"><font color="#404040" face="gulim, sans-serif"> </font></span></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>W = 123; E = 131; S = 33; N=43;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>dD = 0.1;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>num_lon = int((E-W)/dD); num_lat = int((N-S)/dD);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>Ml=4.5; Mu=7; dM=0.5;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>num_M = int((Mu-Ml)/dM);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>FILE* fin;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>FILE* fin_o1;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>FILE* fin_o2;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>FILE* fin_o3;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>FILE* fin_r1;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>FILE* fin_r2;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>FILE* fin_r3;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>FILE* finput;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>FILE* frmse;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>FILE* fsmoothb;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>///////////////////////////////////////////////////////////////////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>///////////////////////////////////////// Read Input data /////////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fin = fopen("input.txt","r");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fscanf(fin,"%24c\n",temp);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fscanf(fin,"%d<span class="Apple-tab-span" style="white-space:pre;"> </span>%lf<span class="Apple-tab-span" style="white-space:pre;"> </span>%lf<span class="Apple-tab-span" style="white-space:pre;"> </span>%lf\n",&iter,&del_b,&RMSE_a,&RMSE_b);</span></font></div> <div><span class="Apple-tab-span" style="white-space:pre;font-size:13px;line-height:22.1px;"><font color="#404040" face="gulim, sans-serif"> </font></span></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>printf("Iteration Number = %d\n",iter);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>///////////////////////////////////////////////////////////////////////</span></font></div> <div><span class="Apple-tab-span" style="white-space:pre;font-size:13px;line-height:22.1px;"><font color="#404040" face="gulim, sans-serif"> </font></span></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>///////////////////////////////////////////////////////////////////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>///////////////////////////////////////// Obtain Original Data ////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>sprintf(temp,"seismic_map_%d_original.dat",return_p[1]);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fin_o1 = fopen(temp,"r");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>sprintf(temp,"seismic_map_%d_original.dat",return_p[2]);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fin_o2 = fopen(temp,"r");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>sprintf(temp,"seismic_map_%d_original.dat",return_p[3]);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fin_o3 = fopen(temp,"r");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>b = double(control_b);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if(control_b == 10)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fsmoothb = fopen("smooth_b.txt","r");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>sm_b = new double[9001];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>original1 = new double[9001];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>original2 = new double[9001];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>original3 = new double[9001];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>for(j=1;j<=num_lat;j++)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>for(i=1;i<=num_lon;i++)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>m = (j-1)*num_lon + i;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fscanf(fin_o1,"%lf",&original1[m]);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fscanf(fin_o2,"%lf",&original2[m]);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fscanf(fin_o3,"%lf",&original3[m]);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if(control_b == 10)<span class="Apple-tab-span" style="white-space:pre;"> </span>fscanf(fsmoothb,"%lf\n",&sm_b[m]);<span class="Apple-tab-span" style="white-space:pre;"> </span></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fclose(fin_o1);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fclose(fin_o2);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fclose(fin_o3);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if(control_b == 10) fclose(fsmoothb);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>///////////////////////////////////////////////////////////////////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>///////////////////////////////////////////////////////////////////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>//////////////////////////////// Scale by Original / New data /////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>sprintf(temp,"RMSE_%.2f.txt",b);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>frmse = fopen(temp,"w");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>sprintf(temp,"copy ab_data_ceus_%.2f.dat ab_data_ceus.dat",b);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>system(temp);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if(iter == 0)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>sprintf(temp,"copy ab_data_ceus.dat ab_data_ceus_0.dat");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>system(temp);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>iteration(0);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>sprintf(temp,"copy seismic_map_500.dat seismic_map_500_0.dat");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>system(temp);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>sprintf(temp,"copy seismic_map_1000.dat seismic_map_1000_0.dat");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>system(temp);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>sprintf(temp,"copy seismic_map_2400.dat seismic_map_2400_0.dat");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>system(temp);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fin_r1 = fopen("seismic_map_500.dat","r");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fin_r2 = fopen("seismic_map_1000.dat","r");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fin_r3 = fopen("seismic_map_2400.dat","r");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>result1 = new double[9001];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>result2 = new double[9001];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>result3 = new double[9001];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>RMSE = 0;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>for(j=1;j<=num_lat;j++)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>for(i=1;i<=num_lon;i++)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>m = (j-1)*num_lon + i;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fscanf(fin_r1,"%lf",&result1[m]);<span class="Apple-tab-span" style="white-space:pre;"> </span></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fscanf(fin_r2,"%lf",&result2[m]);<span class="Apple-tab-span" style="white-space:pre;"> </span></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fscanf(fin_r3,"%lf",&result3[m]);<span class="Apple-tab-span" style="white-space:pre;"> </span></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>RMSE = RMSE + fabs(original1[m]-result1[m])+fabs(original2[m]-result2[m])+fabs(original3[m]-result3[m]);<span class="Apple-tab-span" style="white-space:pre;"> </span>// RMSE check</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>printf("RMSE_a = %lf\n",RMSE);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fprintf(frmse,"%.6f\n",RMSE);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fclose(fin_r1);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fclose(fin_r2);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fclose(fin_r3);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}<span class="Apple-tab-span" style="white-space:pre;"> </span></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>for(ii=1;ii<=iter;ii++)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if( ii != 1) system("copy ab_data_ceus_temp.dat ab_data_ceus.dat");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>sprintf(temp,"copy ab_data_ceus.dat ab_data_ceus_%d.dat",ii);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>system(temp);<span class="Apple-tab-span" style="white-space:pre;"> </span></span></font></div> <div><span class="Apple-tab-span" style="white-space:pre;font-size:13px;line-height:22.1px;"><font color="#404040" face="gulim, sans-serif"> </font></span></div> <div><font face="gulim, sans-serif" color="#ff0000"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>iteration(ii);<span class="Apple-tab-span" style="white-space:pre;"> </span>// Gain New Data</span></font></div> <div><font face="gulim, sans-serif" color="#ff0000"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font face="gulim, sans-serif" color="#ff0000"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>sprintf(temp,"copy seismic_map_500.dat seismic_map_500_%d.dat",ii);</span></font></div> <div><font face="gulim, sans-serif" color="#ff0000"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>system(temp);</span></font></div> <div><font face="gulim, sans-serif" color="#ff0000"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>sprintf(temp,"copy seismic_map_1000.dat seismic_map_1000_%d.dat",ii);</span></font></div> <div><font face="gulim, sans-serif" color="#ff0000"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>system(temp);</span></font></div> <div><font face="gulim, sans-serif" color="#ff0000"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>sprintf(temp,"copy seismic_map_2400.dat seismic_map_2400_%d.dat",ii);</span></font></div> <div><font face="gulim, sans-serif" color="#ff0000"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>system(temp);</span></font></div> <div><font face="gulim, sans-serif" color="#ff0000"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>RMSE = 0;</span></font></div> <div><font face="gulim, sans-serif" color="#ff0000"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fin_r1 = fopen("seismic_map_500.dat","r");</span></font></div> <div><font face="gulim, sans-serif" color="#ff0000"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fin_r2 = fopen("seismic_map_1000.dat","r");</span></font></div> <div><font face="gulim, sans-serif" color="#ff0000"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fin_r3 = fopen("seismic_map_2400.dat","r");</span></font></div> <div><font face="gulim, sans-serif" color="#ff0000"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>result1 = new double[9001];</span></font></div> <div><font face="gulim, sans-serif" color="#ff0000"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>result2 = new double[9001];</span></font></div> <div><font face="gulim, sans-serif" color="#ff0000"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>result3 = new double[9001];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>for(j=1;j<=num_lat;j++)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>for(i=1;i<=num_lon;i++)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>m = (j-1)*num_lon + i;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fscanf(fin_r1,"%lf",&result1[m]);<span class="Apple-tab-span" style="white-space:pre;"> </span></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fscanf(fin_r2,"%lf",&result2[m]);<span class="Apple-tab-span" style="white-space:pre;"> </span></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fscanf(fin_r3,"%lf",&result3[m]);<span class="Apple-tab-span" style="white-space:pre;"> </span></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>RMSE = RMSE + fabs(original1[m]-result1[m])+fabs(original2[m]-result2[m])+fabs(original3[m]-result3[m]);<span class="Apple-tab-span" style="white-space:pre;"> </span>// RMSE check</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>printf("RMSE_a = %lf\n",RMSE);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fprintf(frmse,"%.6f\n",RMSE);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fclose(fin_r1);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fclose(fin_r2);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fclose(fin_r3);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if(RMSE < RMSE_a || ii == iter) break;</span></font></div> <div><span class="Apple-tab-span" style="white-space:pre;font-size:13px;line-height:22.1px;"><font color="#404040" face="gulim, sans-serif"> </font></span></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>finput = fopen("ab_data_ceus_temp.dat","w");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>for(j=1;j<=num_lat;j++)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{<span class="Apple-tab-span" style="white-space:pre;"> </span></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>lat = 33.0+(j-1)*dD;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>for(i=1;i<=num_lon;i++)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>lon = 123.0+(i-1)*dD;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>m = (j-1)*num_lon + i;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if( aa[i][j] != 0 ) v_old = pow(10,log10(aa[i][j])-bb[i][j]*Ml);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else v_old = 0;<span class="Apple-tab-span" style="white-space:pre;"> </span></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if(original1[m] < 0.001) a10 = 0;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if( v_old < 1e-10 ) a10 = 0;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else if( result1[m] == 0) a10 = aa[i][j];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if( original1[m] >= result1[m] && original2[m] >= result2[m] && original3[m] >= result3[m] )</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>v_new = v_old * original2[m] / result2[m];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else if( original1[m] < result1[m] && original2[m] < result2[m] && original3[m] < result3[m] )</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>v_new = v_old * original2[m] / result2[m];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else v_new = v_old;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if( original1[m] >= result1[m] && original3[m] < result3[m] )</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>sm_b[m] = sm_b[m] + del_b;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else if( original1[m] < result1[m] && original3[m] >= result3[m] )</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>sm_b[m] = sm_b[m] - del_b;</span></font></div> <div><span class="Apple-tab-span" style="white-space:pre;font-size:13px;line-height:22.1px;"><font color="#404040" face="gulim, sans-serif"> </font></span></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if( sm_b[m] >= 2.0) sm_b[m] = 2.0;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else if( sm_b[m] <= 0.2 ) sm_b[m] = 0.2;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if(control_b == 10) a10 = pow(10,log10(v_new)+sm_b[m]*Ml);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else a10 = pow(10,log10(v_new)+b*Ml);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><span class="Apple-tab-span" style="white-space:pre;font-size:13px;line-height:22.1px;"><font color="#404040" face="gulim, sans-serif"> </font></span></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if(a10 == 0)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if(control_b == 10)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fprintf(finput,"%.2f<span class="Apple-tab-span" style="white-space:pre;"> </span>%.2f<span class="Apple-tab-span" style="white-space:pre;"> </span>%.8f<span class="Apple-tab-span" style="white-space:pre;"> </span>%.8f<span class="Apple-tab-span" style="white-space:pre;"> </span>0.00\n",lon,lat,a10,sm_b[m]);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fprintf(finput,"%.2f<span class="Apple-tab-span" style="white-space:pre;"> </span>%.2f<span class="Apple-tab-span" style="white-space:pre;"> </span>%.8f<span class="Apple-tab-span" style="white-space:pre;"> </span>%.8f<span class="Apple-tab-span" style="white-space:pre;"> </span>0.00\n",lon,lat,a10,b);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if(control_b == 10)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fprintf(finput,"%.2f<span class="Apple-tab-span" style="white-space:pre;"> </span>%.2f<span class="Apple-tab-span" style="white-space:pre;"> </span>%.8f<span class="Apple-tab-span" style="white-space:pre;"> </span>%.8f<span class="Apple-tab-span" style="white-space:pre;"> </span>%.2f\n",lon,lat,a10,sm_b[m],Mu);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fprintf(finput,"%.2f<span class="Apple-tab-span" style="white-space:pre;"> </span>%.2f<span class="Apple-tab-span" style="white-space:pre;"> </span>%.8f<span class="Apple-tab-span" style="white-space:pre;"> </span>%.8f<span class="Apple-tab-span" style="white-space:pre;"> </span>%.2f\n",lon,lat,a10,b,Mu);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fclose(finput);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}<span class="Apple-tab-span" style="white-space:pre;"> </span></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fclose(frmse);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>//////////////////////////////////////////////////////////////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">int iteration(int iter)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>int i,ii,j,k1,k2,k3,m,kk;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>int num_acc;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>int atten;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>int low1, mid1, high1, low2, mid2, high2, low3, mid3, high3;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double lon, lat, mw;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double rec_p[4];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double sigma;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double t_lon,t_lat,a10,b;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double a_est1, a_est2, a_est3;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double* a;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double c[5];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double* prob;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double* pro;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double result;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>clock_t start,end;</span></font></div> <div><span class="Apple-tab-span" style="white-space:pre;font-size:13px;line-height:22.1px;"><font color="#404040" face="gulim, sans-serif"> </font></span></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>rec_p[1] = 1./return_p[1];<span class="Apple-tab-span" style="white-space:pre;"> </span>// Return Period</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>rec_p[2] = 1./return_p[2];<span class="Apple-tab-span" style="white-space:pre;"> </span>// Return Period</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>rec_p[3] = 1./return_p[3];<span class="Apple-tab-span" style="white-space:pre;"> </span>// Return Period</span></font></div> <div><span class="Apple-tab-span" style="white-space:pre;font-size:13px;line-height:22.1px;"><font color="#404040" face="gulim, sans-serif"> </font></span></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>///////////////////////////////////////////////////////////////////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>//<span class="Apple-tab-span" style="white-space:pre;"> </span>Part 1 : Produce Seismic Hazard Map<span class="Apple-tab-span" style="white-space:pre;"> </span>///////////////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>///////////////////////////////////////////////////////////////////////</span></font></div> <div><span class="Apple-tab-span" style="white-space:pre;font-size:13px;line-height:22.1px;"><font color="#404040" face="gulim, sans-serif"> </font></span></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>atten = 1;<span class="Apple-tab-span" style="white-space:pre;"> </span>// Select Attenuation Relationship</span></font></div> <div><span class="Apple-tab-span" style="white-space:pre;font-size:13px;line-height:22.1px;"><font color="#404040" face="gulim, sans-serif"> </font></span></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>///////////////////////////////////////////////////////////////////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>////////////////////////////////// Determine Attenuation Relationship /</span></font></div> <div><span class="Apple-tab-span" style="white-space:pre;font-size:13px;line-height:22.1px;"><font color="#404040" face="gulim, sans-serif"> </font></span></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if(atten == 1)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>c[1] = 0.4854; c[2] = 1.2; c[3] = -0.8416; c[4] = -0.0061;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>sigma = 0.8036;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else if(atten == 2)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>c[1] = 0.5577; c[2] = 1.2; c[3] = -0.8587; c[4] = -0.0062;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>sigma = 0.7629;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else if(atten == 3)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>c[1] = 5.0244; c[2] = 0.5442; c[3] = -1.0020; c[4] = 0.0;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>sigma = 0.1;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>///////////////////////////////////////////////////////////////////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>///////////////////////////////////////////////////////////////////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>////////////////////////////////// Set Target Acceleration ////////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>a = new double[33];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>a[1]=0.001;a[2]=0.01;a[3]=0.02;a[4]=0.03;a[5]=0.04;a[6]=0.05;a[7]=0.06;a[8]=0.07;a[9]=0.08;a[10]=0.09;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>a[11]=0.10;a[12]=0.11;a[13]=0.12;a[14]=0.13;a[15]=0.14;a[16]=0.15;a[17]=0.16;a[18]=0.17;a[19]=0.18;a[20]=0.19;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>a[21]=0.20;a[22]=0.21;a[23]=0.22;a[24]=0.23;a[25]=0.24;a[26]=0.25;a[27]=0.26;a[28]=0.27;a[29]=0.28;a[30]=0.29;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>a[31]=0.30;a[32]=0.31;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>num_acc = 32;<span class="Apple-tab-span" style="white-space:pre;"> </span></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>///////////////////////////////////////////////////////////////////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>////////////////////////////////////////////////////////////////////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>/////////////////////////////////// Read Gutenberg-Richter Param. a, b /</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>FILE* fin;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>FILE* fout1;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>FILE* fout2;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>FILE* fout3;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>FILE* fpro;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fin = fopen("ab_data_ceus.dat","r");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fout1 = fopen("seismic_map_500.dat","w");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fout2 = fopen("seismic_map_1000.dat","w");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fout3 = fopen("seismic_map_2400.dat","w");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>for(j=1;j<=num_lat;j++)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>for(i=1;i<=num_lon;i++)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>m = (j-1)*num_lon + i;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fscanf(fin," %lf %lf %lf %lf %lf\n",&t_lon,&t_lat,&a10,&b,&mw);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>aa[i][j] = a10;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>bb[i][j] = b;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fclose(fin);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>////////////////////////////////////////////////////////////////////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>////////////////////////////////////////////////////////////////////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>/////////////////////////////////// Read Normaly Distributed Function //</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>pro = new double[400];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fpro = fopen("CDF.txt","r");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>for(i=0;i<=399;i++)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fscanf(fpro,"%lf",&pro[i]);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fclose(fpro);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>////////////////////////////////////////////////////////////////////////</span></font></div> <div><span class="Apple-tab-span" style="white-space:pre;font-size:13px;line-height:22.1px;"><font color="#404040" face="gulim, sans-serif"> </font></span></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>//////////////////////////////////////////////////////////////////////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>//////////////////// Determine Acceleration of Target Return Period //////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>prob = new double[17];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>lat = S;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>for(j=1;j<=num_lat;j++)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>start = clock();<span class="Apple-tab-span" style="white-space:pre;"> </span>// check time </span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>lon = W;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>for(i=1;i<=num_lon;i++)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>m = (j-1)*num_lon + i;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>a_est1 = 0.;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>a_est2 = 0.;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>a_est3 = 0.;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>kk = int(log(num_acc)/log(2));</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>low1 = low2 = low3 = 0;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>high1 = high2 = high3 = num_acc;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>mid1 = mid2 = mid3 = num_acc/2;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>for(ii=1;ii<=num_acc;ii++)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>prob[ii] = -1;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>for(ii=1;ii<=kk;ii++)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>k1 = mid1;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if( prob[k1] < 0 )</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>prob[k1] = probability(c,a[k1],sigma,lat,lon,pro);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if( prob[k1] < rec_p[1] ) {<span class="Apple-tab-span" style="white-space:pre;"> </span>high1 = mid1;<span class="Apple-tab-span" style="white-space:pre;"> </span>mid1 = (low1 + high1)/2;<span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else if( prob[k1] > rec_p[1] )<span class="Apple-tab-span" style="white-space:pre;"> </span>{<span class="Apple-tab-span" style="white-space:pre;"> </span>low1 = mid1;<span class="Apple-tab-span" style="white-space:pre;"> </span>mid1 = (low1 + high1)/2;<span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else<span class="Apple-tab-span" style="white-space:pre;"> </span>{<span class="Apple-tab-span" style="white-space:pre;"> </span>a_est1 = a[k1];<span class="Apple-tab-span" style="white-space:pre;"> </span>break;<span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>k2 = mid2;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if( prob[k2] < 0 )</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>prob[k2] = probability(c,a[k2],sigma,lat,lon,pro);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if( prob[k2] < rec_p[2] ) {<span class="Apple-tab-span" style="white-space:pre;"> </span>high2 = mid2;<span class="Apple-tab-span" style="white-space:pre;"> </span>mid2 = (low2 + high2)/2;<span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else if( prob[k2] > rec_p[2] )<span class="Apple-tab-span" style="white-space:pre;"> </span>{<span class="Apple-tab-span" style="white-space:pre;"> </span>low2 = mid2;<span class="Apple-tab-span" style="white-space:pre;"> </span>mid2 = (low2 + high2)/2;<span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else<span class="Apple-tab-span" style="white-space:pre;"> </span>{<span class="Apple-tab-span" style="white-space:pre;"> </span>a_est2 = a[k2];<span class="Apple-tab-span" style="white-space:pre;"> </span>break;<span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>k3 = mid3;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if( prob[k3] < 0 )</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>prob[k3] = probability(c,a[k3],sigma,lat,lon,pro);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if( prob[k3] < rec_p[3] ) {<span class="Apple-tab-span" style="white-space:pre;"> </span>high3 = mid3;<span class="Apple-tab-span" style="white-space:pre;"> </span>mid3 = (low3 + high3)/2;<span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else if( prob[k3] > rec_p[3] )<span class="Apple-tab-span" style="white-space:pre;"> </span>{<span class="Apple-tab-span" style="white-space:pre;"> </span>low3 = mid3;<span class="Apple-tab-span" style="white-space:pre;"> </span>mid3 = (low3 + high3)/2;<span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else<span class="Apple-tab-span" style="white-space:pre;"> </span>{<span class="Apple-tab-span" style="white-space:pre;"> </span>a_est3 = a[k3];<span class="Apple-tab-span" style="white-space:pre;"> </span>break;<span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if(rec_p[1] > prob[k1])</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{<span class="Apple-tab-span" style="white-space:pre;"> </span>if(k1 == 1) a_est1 = 0;<span class="Apple-tab-span" style="white-space:pre;"> </span>a_est1 = a[k1] - (a[k1]-a[k1-1])*(prob[k1]-rec_p[1])/(prob[k1]-prob[k1-1]);<span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else if(rec_p[1] < prob[k1])</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>a_est1 = a[k1] - (a[k1+1]-a[k1])*(prob[k1]-rec_p[1])/(prob[k1+1]-prob[k1]);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fprintf(fout1,"%.8f ",a_est1);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if(rec_p[2] > prob[k2])</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{<span class="Apple-tab-span" style="white-space:pre;"> </span>if(k2 == 1) a_est2 = 0;<span class="Apple-tab-span" style="white-space:pre;"> </span>a_est2 = a[k2] - (a[k2]-a[k2-1])*(prob[k2]-rec_p[2])/(prob[k2]-prob[k2-1]);<span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else if(rec_p[2] < prob[k2])</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>a_est2 = a[k2] - (a[k2+1]-a[k2])*(prob[k2]-rec_p[2])/(prob[k2+1]-prob[k2]);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fprintf(fout2,"%.8f ",a_est2);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if(rec_p[3] > prob[k3])</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{<span class="Apple-tab-span" style="white-space:pre;"> </span>if(k3 == 1) a_est3 = 0;<span class="Apple-tab-span" style="white-space:pre;"> </span>a_est3 = a[k3] - (a[k3]-a[k3-1])*(prob[k3]-rec_p[3])/(prob[k3]-prob[k3-1]);<span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else if(rec_p[3] < prob[k3])</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>a_est3 = a[k3] - (a[k3+1]-a[k3])*(prob[k3]-rec_p[3])/(prob[k3+1]-prob[k3]);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fprintf(fout3,"%.8f ",a_est3);</span></font></div> <div><span class="Apple-tab-span" style="white-space:pre;font-size:13px;line-height:22.1px;"><font color="#404040" face="gulim, sans-serif"> </font></span></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>lon = lon + dD;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>end = clock();<span class="Apple-tab-span" style="white-space:pre;"> </span>// check time</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>result=(double)(end-start) / CLOCKS_PER_SEC;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>printf("%3.2f<span class="Apple-tab-span" style="white-space:pre;"> </span>%2.1f sec<span class="Apple-tab-span" style="white-space:pre;"> </span>%d\n",lat,result,iter);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fprintf(fout1,"\n");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fprintf(fout2,"\n");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fprintf(fout3,"\n");</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>lat = lat + dD;<span class="Apple-tab-span" style="white-space:pre;"> </span></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fclose(fout1);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fclose(fout2);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>fclose(fout3);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>//////////////////////////////////////////////////////////////////////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>return 0;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">//////////////////////////////////////////////////////////////////////////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">//<span class="Apple-tab-span" style="white-space:pre;"> </span> Yearly generate probability above u0 acceleration<span class="Apple-tab-span" style="white-space:pre;"> </span>//</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">//////////////////////////////////////////////////////////////////////////////</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">double probability(double c[5],double u,double sigma,double lat2,double lon2,double* pro)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>int i, j, k, l;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>int f_M_check;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double lat, lon, M;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double alpha, beta;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double pi, dv, dh, z;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double CDF, p_M, f_MM;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double a, R;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double prob;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>double* f_M;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>f_M = new double[num_M+1];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><br></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>pi = 3.141592;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>dv = 2*pi*6400/360;<span class="Apple-tab-span" style="white-space:pre;"> </span>// unit distance of latitude</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>prob = 0;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>f_M_check = 1;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>lon = W;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>for(i=1;i<=num_lon;i++)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>lat = S;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>for(j=1;j<=num_lat;j++)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>beta = 2.303*bb[i][j];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if(aa[i][j] != 0)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>alpha = 2.303*log10(aa[i][j]);</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>p_M = exp(alpha-beta*Ml);<span class="Apple-tab-span" style="white-space:pre;"> </span>// V</span></font></div> <div><span class="Apple-tab-span" style="white-space:pre;font-size:13px;line-height:22.1px;"><font color="#404040" face="gulim, sans-serif"> </font></span></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>dh = dv*cos((lat+lat2+dD)/2*pi/180);<span class="Apple-tab-span" style="white-space:pre;"> </span>// unit distance of longitude</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>R = sqrt(pow(((lon2-lon)*dh),2)+pow(((lat2-lat)*dv),2));</span></font></div> <div><span class="Apple-tab-span" style="white-space:pre;font-size:13px;line-height:22.1px;"><font color="#404040" face="gulim, sans-serif"> </font></span></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>M = Ml;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>for(k=1;k<=num_M;k++)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if( control_b == 10 )</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>f_MM = beta*exp(-beta*(M+dM/2-Ml))/(1-exp(-beta*(Mu-Ml)))*dM;<span class="Apple-tab-span" style="white-space:pre;"> </span>// P[M]</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>f_M[k] = f_MM;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if( f_M_check != 0 )</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>f_MM = beta*exp(-beta*(M+dM/2-Ml))/(1-exp(-beta*(Mu-Ml)))*dM;<span class="Apple-tab-span" style="white-space:pre;"> </span>// P[M]</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>f_M[k] = f_MM;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if( k == num_M ) f_M_check = 0;<span class="Apple-tab-span" style="white-space:pre;"> </span></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>a = exp(c[1]+c[2]*(M+dM/2)+c[3]*log(sqrt(R*R+100))+c[4]*(sqrt(R*R+100)))/981;<span class="Apple-tab-span" style="white-space:pre;"> </span>// Acc. by ARship</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if( a > 0.001 )</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>{<span class="Apple-tab-span" style="white-space:pre;"> </span></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>z = (log(u*981)-log(a*981))/sigma;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>l = int(ceil(fabs(z)*100-0.5));</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if(l > 399) l=399;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>if(z < 0)</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>CDF = 1-pro[l];</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>else</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>CDF = pro[l];<span class="Apple-tab-span" style="white-space:pre;"> </span></span></font></div> <div><span class="Apple-tab-span" style="white-space:pre;font-size:13px;line-height:22.1px;"><font color="#404040" face="gulim, sans-serif"> </font></span></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>prob = prob + f_M[k]*CDF*p_M;<span class="Apple-tab-span" style="white-space:pre;"> </span>// Annual Occurence Rate</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>M=M+dM;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}<span class="Apple-tab-span" style="white-space:pre;"> </span></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}<span class="Apple-tab-span" style="white-space:pre;"> </span></span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>lat = lat + dD;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>lon = lon + dD;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>}</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;"><span class="Apple-tab-span" style="white-space:pre;"> </span>return prob;</span></font></div> <div><font color="#404040" face="gulim, sans-serif"><span style="font-size:13px;line-height:22.1px;">}</span></font></div> <div style="color:#404040;font-family:gulim, sans-serif;font-size:13px;line-height:22.1px;"><br></div></div>