日本財団 図書館


 

プログラムリスト

/* I--------------------------------------------I

     I                                                                                       I

     I           海水溜への海水侵入シュミレーション             I

     I                              第三回委員会資料用                    I

     I                                                                                       I

     I                                  1997- 1/27(月)                           I

     I                                                                                       I  

     I--------------------------------------------I */

/**/

/**/

#include <stdio.h>

#include <math .h>

/*---------- 定数設定 -------------*/

#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 を返却 */

}

/**/

/**/

 

 

 

前ページ   目次へ   次ページ

 






日本財団図書館は、日本財団が運営しています。

  • 日本財団 THE NIPPON FOUNDATION