プログラムリスト
/* I--------------------------------------------I
I
I
I
海水溜への海水侵入シュミレーション I
I
第三回委員会資料用
I
I
I
I
1997- 1/27(月)
I
I
I
I--------------------------------------------I */
/**/
/**/
#include <stdio.h>
#include
/*---------- 定数設定 -------------*/
#define PI 3.14159265358
#define NAME "平潟港"
#define TO 15.0 /* 周 期 */
#define HH 7.6 /* 波 高 */
#define h 9.0 /* 水 深 */
#define HWL 1.5
/* */
#define h0 6.5-HWL /* 天端高 */
#define dc 1.6+HWL /* 没水部 */
/* */
#define ir_h 1.8 /* Hmax/H1/3 */
#define ir_w 2 /* 高調波/W */
/*------- 定数からの計算値 --------*/
#define w 2*PI/TO
/*------- 微分方程式係数 ----------*/
#define C1 12.2 /* 係数1 */
#define C2 -12.0 /* 係数2 */
#define NM 200 /* 分割数 */
/**/
/**/
/**/
/*-------- 加速度計算関数 --------*/
double alfa(doubule t, doube v, double C, double D) {
return (D*cos(w0*t)/v-1)
}
/**/
/**/
/*-------- 波数計算関数 ----------*/
double kn () {
double i,f,fl,x,w2 ;
/**/
w2 = w*w*h/9.086 ;
if (w2>0.672) x=sqrt(w2*w2+0.2) ;
else x=exp(0.532*log(w2/0.82))
for (i=0 ; i<=8 + i+="1)" {
f = x-w2/tanh(x) ;
fl= 1-w2+w2/(tanh(x)*tanh(x)) ;
x-= f/f1 ;
}
return x/h ; /* k を返却 */
}
/**/
/**/
前ページ 目次へ 次ページ