MSMデータのバーチャル・サウンディング (現地観測と同じ地点・時刻のデータを抽出する)

注意

入力データがnctCDF形式なので、netCDFライブラリをインストールしておく必要があります。

プログラムの所在

プログラム一式の圧縮ファイルは、このページの下にあります（MSM.profiles.at.Stations.zip）。解凍して利用してください。

使用例

[2014年 11月 6日 木曜日 16:12:35 JST]

[~/to_manda_sensei/Tools/MSM.profiles.at.Stations]

[kunoki@aofd165]

$ cd src

$ ls

bilinear.f90 print_header.f90

cal_vpt.f90 print_ioinfo.f90

interpolate_space.f90 read_msmp.f90

interpolate_time.f90 read_msms.f90

makefile set_out_range.f90

module_inout.f90 srcbackup/

module_main.f90 virtual_potential_temperature.f90

module_msmp_var.f90 write_msmp.f90

module_msms_var.f90 write_msms.f90

msmstn.f90 write_stn.f90

$ make clean; make

if [ ! -d ../obj ]; then \

mkdir -p ../obj ; \

fi

ifort -c -CB -traceback -fpe0 -module ../obj -c -o ../obj/module_main.o module_main.f90

ifort -c -CB -traceback -fpe0 -module ../obj -c -o ../obj/module_msms_var.o module_msms_var.f90

ifort -c -CB -traceback -fpe0 -module ../obj -c -o ../obj/module_msmp_var.o module_msmp_var.f90

ifort -c -CB -traceback -fpe0 -module ../obj -c -o ../obj/module_inout.o module_inout.f90

ifort -c -CB -traceback -fpe0 -module ../obj -c -o ../obj/msmstn.o msmstn.f90

ifort -c -CB -traceback -fpe0 -module ../obj -c -o ../obj/read_msms.o read_msms.f90

ifort -c -CB -traceback -fpe0 -module ../obj -c -o ../obj/read_msmp.o read_msmp.f90

ifort -c -CB -traceback -fpe0 -module ../obj -c -o ../obj/interpolate_space.o interpolate_space.f90

ifort -c -CB -traceback -fpe0 -module ../obj -c -o ../obj/bilinear.o bilinear.f90

ifort -c -CB -traceback -fpe0 -module ../obj -c -o ../obj/interpolate_time.o interpolate_time.f90

ifort -c -CB -traceback -fpe0 -module ../obj -c -o ../obj/write_stn.o write_stn.f90

ifort -c -CB -traceback -fpe0 -module ../obj -c -o ../obj/cal_vpt.o cal_vpt.f90

ifort -c -CB -traceback -fpe0 -module ../obj -c -o ../obj/virtual_potential_temperature.o virtual_potential_temperature.f90

ifort -o ../msmstn ../obj/module_main.o ../obj/module_msms_var.o ../obj/module_msmp_var.o ../obj/module_inout.o ../obj/msmstn.o ../obj/read_msms.o ../obj/read_msmp.o ../obj/interpolate_space.o ../obj/bilinear.o ../obj/interpolate_time.o ../obj/write_stn.o ../obj/cal_vpt.o ../obj/virtual_potential_temperature.o -module ../obj -L/usr/local/lib -lnetcdf

$ cd ..

$ msmstn.all.run.sh

..................

..................

write_stn:

Output: output/msmstn.MandA2012Leg3.19.txt

Done program msmstn.

付録

鉛直プロファイルを高度方向に内挿するスクリプト

$ interpolate.in.z.run.sh

....................

....................

output/msmstn.MandA2012Leg3.19.txt

output/msmstn.MandA2012Leg3.19.intp.txt

風速の鉛直プロファイルの作図

$ pl.uv.run.sh

....................

....................

Input : output/msmstn.MandA2012Leg3.19.txt

Output : ./fig/msmstn.MandA2012Leg3.19.ps