#!/bin/sh
#
# bsnsq_2d モデル実行サンプルスクリプト   96/03/01 S.Takehiro
#
#    注意 : make lib, make  を実行し, 実行ファイルを作っておく
#	    初期値ファイルを作成しておく(go-init)
#	    解析・表示のための軸ファイルを作っておく(go-axis)
#
##################### 以下環境に合わせて変更する #######################
#
#	bsnsq_2d インストールディレクトリ名
BSNDIR=$HOME/model/bsnsq_2d
#
#	モデル解像度
####RESOLUTION=nx32nz16
RESOLUTION=nx64nz8
#
# 	モデル実行用 NAMELIST 入出力ファイル名
CNTLDIR=$BSNDIR/cntl/sample
INPUT=$CNTLDIR/sample.in
OUTPUT=$CNTLDIR/sample.out
#
# 	出力データファイルを置くディレクトリ
#	(プログラム実行時のカレントパス)
DATADIR=$BSNDIR/data/test
#
# 	初期値データファイル名
INITDAT=sample.init
#
# 	出力データファイル名(リスタートファイル名)
RSTDAT=sample.rst
DFORMAT=UR4            # 単精度
### DFORMAT=UR8            # 倍精度
#
#	実行プログラム名
RUN=$BSNDIR/bin/bsn.$RESOLUTION
#
########################################################################
cat << END_OF_DATA > $INPUT
#
# Control File for bsnsq_2d
#
 &nmrun  run='sample', sign='$USER'              &end
#
# << 時間積分設定 >>
 &nmtime  start=  0, end=  5000, intrst=  500 &end
 &nmintg  delt=0.0001, scheme='RUNGE' &end
#
# << 計算領域 >>
 &nmcord  xlngth=  10.0, zlngth=    1.0  &end
# << 境界条件 >>
 &nmbndr  rigidb=  F, rigidt=  F, tfixb=  T, tfixt=  T &end
#
# << 初期条件 >>
 &nminit file='$INITDAT'    &end
#
# << 無次元数 >>
 &nmdyn  rayli=  10000., prndl=    1.00000, qintr=  0. &end
#
# << 出力 >>
#
 &nmrstr file='$RSTDAT' DFMT='$DFORMAT' &end
END_OF_DATA
##################################################################
#
        cd $DATADIR
	echo 'Program started at' `date` 'on' `hostname` 1> $OUTPUT
        time $RUN < $INPUT 1>> $OUTPUT 2>&1
	echo 'Program ended at' `date` 1>> $OUTPUT