リストを与えて、鉛直分布を得るためにgrads

・CCYYMMDDHHのリストを与える。
・時刻を調べる。ファイルは複数開く⇒日付ごとにOPENする
・set t ##
・PLOT
・xps

参考はjobs/2005/mrr/mk_grads

grads.scriptの中では引数をとってスクリプトに渡しても、結果を受け取ることはできない(たぶん)。外側のスクリプトで、結果をもらって'set t INUM'に入るべき数を引数で渡す。

2009/syaka/kdp/bb/vertical(sn3)
#!/usr/bin/perl
# script for make vertical profile 6/17'10
require "date2inum.pl";
if( $ARGV[0] eq "" ){die "Melting Layer list is necessary!\n";}
$mllst = $ARGV[0];#読み込みのファイルリストを与える(CCYYMMDDHHFFの並び)
$GSF = "verp2.gs"; # plot vertical profile
$PRG = "grads";
@ccyymmddhhff = `cat $mllst`;
chomp @ccyymmddhhff;
$n=0;
foreach $ccyymmddhhff(@ccyymmddhhff){
$n++;
$inum = date2inum($ccyymmddhhff);#サブルーチンで何時刻目かを調べる。
$ccyymmdd = substr($ccyymmddhhff,0,8);#ファイル名用の日付
$grads_argv = "\"run $GSF $ccyymmdd $inum\"";
print $grads_argv,"\n";
system("$PRG -lbc $grads_argv > /dev/null");
}