WaveWatch3学习笔记——模拟实例（2）

2、预处理

mkdir tmp_case_Pacific work_case_Pacific
cd tmp_case_Pacific

# 如果环境变量发生变化，请重新编译
../../bin/w3_clean all
../../bin/w3_new

# 写switch文件
# 参考switch_ser: 
# echo "NC4 F90 NOGRB NOPA LRB4 SHRD PR3 UQ FLX2 LN1 ST2 STAB2 NL1 BT1 DB1 MLIM TR0 BS0 IC0 IS0 REF0 XX0 WNT1 WNX1 CRT1 CRX1 O0 O1 O2 O3 O4 O5 O6 O7 O11 O14" > ../../bin/switch
# 参考switch_mp: 
# echo "NC4 F90 NOGRB NOPA LRB4 MPI DIST PR3 UQ FLX2 LN1 ST2 STAB2 NL1 BT1 DB1 MLIM TR0 BS0 IC0 IS0 REF0 XX0 WNT1 WNX1 CRT1 CRX1 O0 O1 O2 O3 O4 O5 O6 O7 O11 O14" > ../../bin/switch

../../bin/w3_make ww3_grid
../../bin/w3_make ww3_strt
../../bin/w3_make ww3_prep
../../bin/w3_make ww3_multi     # 注意ww3_multi的编译为需要mpi，因此switch文件与其他不同
../../bin/w3_make gx_outf

2.1 网格处理

# 地形网格
cp ../mww3_data_Pacific/ww3_grid.inp.Pacific .
cp ww3_grid.inp.Pacific ww3_grid.inp
cp ../mww3_data_Pacific/Pacific.bot ../mww3_data_Pacific/Pacific.mask ../mww3_data_Pacific/Pacific.obs .

# 处理地形文件
../../exe/ww3_grid > ../work_case_Pacific/ww3_grid.out    # 得到文件mod_def.ww3
mv mod_def.ww3 mod_def.Pacific

# 风场网格
cp ../mww3_data_Pacific/ww3_grid.inp.gwind .
cp ww3_grid.inp.gwind ww3_grid.inp

# 处理地形文件
../../exe/ww3_grid > ../work_case_Pacific/ww3_grid.out    # 得到文件mod_def.ww3
mv mod_def.ww3 mod_def.gwind

2.2 初始条件

cp ../mww3_data_Pacific/ww3_strt.inp .
cp mod_def.Pacific mod_def.ww3 

# 处理初始条件
../../exe/ww3_strt > ../work_case_Pacific/ww3_strt.out    # 得到文件restart.ww3
mv restart.ww3 restart.Pacific

2.3 风场处理

cp ../mww3_data_Pacific/ww3_prep.inp .
cp ../mww3_data_Pacific/wind.raw .
cp mod_def.gwind mod_def.ww3 

# 处理风场文件
../../exe/ww3_prep > ../work_case_Pacific/ww3_prep.out    # 得到文件wind.ww3
mv wind.ww3 wind.gwind

3 运行计算

cp ../mww3_data_Pacific/ww3_multi.inp .
mpirun -np 2 ../../exe/ww3_multi > ../work_case_Pacific/ww3_multi.out

4.1 场输出的后处理（使用gx_outf）

cp ../mww3_data_Pacific/gx_outf.inp .
../../exe/gx_outf > ../work_case_Pacific/gx_outf.out

4.1 场输出的后处理（使用ww3_ounf）

cp ../mww3_data_Pacific/ww3_ounf.inp .
cp mod_def.Pacific mod_def.ww3
../../exe/ww3_ounf > ../work_case_Pacific/ww3_ounf.out