NOAA/PMEL Ferret メモ　 2010年06月07日(更新）

NOAA/Pacific Marine Ecnvironmental Laboratory製の大気海洋系解析ソフト。

let a=$par　　　　　　　　　　　　　　　　　　　　　--->letで計算式の入力。この場合全く必要ない。list/x=208/y=71.5/format(f8.4)/file=$ofl2/clobber a　--->listコマンドを使う。Formatで出力形式を選択。
quit
EOF
ferret < fjnl
set num=`cat $ofl2 | wc -l`
awk 'FNR==10,FNR=='$num' {printf "%8.4f \n",$1}' $ofl2 >! tmp2

paste tmp1 tmp2 >! $ofl

@ yr++
end

#!/bin/csh
set jnl=godin.jnl
foreach ifile (dvel-M0104.dat)
set ofile=`echo $ifile|sed s/dat/godin.asc/g`
set fig=`echo $ifile|sed s/dat/godin.ps/g`
set fig2=`echo $ifile|sed s/dat/godin.jpg/g`
set dmax=`cat $ifile|wc -l`

cat <<EOF>! $jnl
set memory/size=80
set mode metafile
file/var="yr,dd,dd2,time,vvel1,vvel2,mvel" $ifile
sh d
let rm1=(vvel1[i=@shf:-12]+vvel1[i=@shf:-11]+vvel1[i=@shf:-10]+vvel1[i=@shf:-9]+vvel1[i=@shf:-8]+vvel1[i=@shf:-7]+vvel1[i=@shf:-6]+vvel1[i=@shf:-5]+vvel1[i=@shf:-4]+vvel1[i=@shf:-3]+vvel1[i=@shf:-2]+vvel1[i=@shf:-1]+vvel1[i=@shf:0]+vvel1[i=@shf:1]+vvel1[i=@shf:2]+vvel1[i=@shf:3]+vvel1[i=@shf:4]+vvel1[i=@shf:5]+vvel1[i=@shf:6]+vvel1[i=@shf:7]+vvel1[i=@shf:8]+vvel1[i=@shf:9]+vvel1[i=@shf:10]+vvel1[i=@shf:11])/24
let rm2=(rm1[i=@shf:-11]+rm1[i=@shf:-10]+rm1[i=@shf:-9]+rm1[i=@shf:-8]+rm1[i=@shf:-7]+rm1[i=@shf:-6]+rm1[i=@shf:-5]+rm1[i=@shf:-4]+rm1[i=@shf:-3]+rm1[i=@shf:-2]+rm1[i=@shf:-1]+rm1[i=@shf:0]+rm1[i=@shf:1]+rm1[i=@shf:2]+rm1[i=@shf:3]+rm1[i=@shf:4]+rm1[i=@shf:5]+rm1[i=@shf:6]+rm1[i=@shf:7]+rm1[i=@shf:8]+rm1[i=@shf:9]+rm1[i=@shf:10]+rm1[i=@shf:11]+rm1[i=@shf:12])/24
let/bad=-999 godin_vector1=rm2[l=@sbx:25]

let rm3=(vvel2[i=@shf:-12]+vvel2[i=@shf:-11]+vvel2[i=@shf:-10]+vvel2[i=@shf:-9]+vvel2[i=@shf:-8]+vvel2[i=@shf:-7]+vvel2[i=@shf:-6]+vvel2[i=@shf:-5]+vvel2[i=@shf:-4]+vvel2[i=@shf:-3]+vvel2[i=@shf:-2]+vvel2[i=@shf:-1]+vvel2[i=@shf:0]+vvel2[i=@shf:1]+vvel2[i=@shf:2]+vvel2[i=@shf:3]+vvel2[i=@shf:4]+vvel2[i=@shf:5]+vvel2[i=@shf:6]+vvel2[i=@shf:7]+vvel2[i=@shf:8]+vvel2[i=@shf:9]+vvel2[i=@shf:10]+vvel2[i=@shf:11])/24
let rm4=(rm3[i=@shf:-11]+rm3[i=@shf:-10]+rm3[i=@shf:-9]+rm3[i=@shf:-8]+rm3[i=@shf:-7]+rm3[i=@shf:-6]+rm3[i=@shf:-5]+rm3[i=@shf:-4]+rm3[i=@shf:-3]+rm3[i=@shf:-2]+rm3[i=@shf:-1]+rm3[i=@shf:0]+rm3[i=@shf:1]+rm3[i=@shf:2]+rm3[i=@shf:3]+rm3[i=@shf:4]+rm3[i=@shf:5]+rm3[i=@shf:6]+rm3[i=@shf:7]+rm3[i=@shf:8]+rm3[i=@shf:9]+rm3[i=@shf:10]+rm3[i=@shf:11]+rm3[i=@shf:12])/24
let/bad=-999 godin_vector2=rm4[l=@sbx:25]

let rm5=(mvel[i=@shf:-12]+mvel[i=@shf:-11]+mvel[i=@shf:-10]+mvel[i=@shf:-9]+mvel[i=@shf:-8]+mvel[i=@shf:-7]+mvel[i=@shf:-6]+mvel[i=@shf:-5]+mvel[i=@shf:-4]+mvel[i=@shf:-3]+mvel[i=@shf:-2]+mvel[i=@shf:-1]+mvel[i=@shf:0]+mvel[i=@shf:1]+mvel[i=@shf:2]+mvel[i=@shf:3]+mvel[i=@shf:4]+mvel[i=@shf:5]+mvel[i=@shf:6]+mvel[i=@shf:7]+mvel[i=@shf:8]+mvel[i=@shf:9]+mvel[i=@shf:10]+mvel[i=@shf:11])/24
let rm6=(rm5[i=@shf:-11]+rm5[i=@shf:-10]+rm5[i=@shf:-9]+rm5[i=@shf:-8]+rm5[i=@shf:-7]+rm5[i=@shf:-6]+rm5[i=@shf:-5]+rm5[i=@shf:-4]+rm5[i=@shf:-3]+rm5[i=@shf:-2]+rm5[i=@shf:-1]+rm5[i=@shf:0]+rm5[i=@shf:1]+rm5[i=@shf:2]+rm5[i=@shf:3]+rm5[i=@shf:4]+rm5[i=@shf:5]+rm5[i=@shf:6]+rm5[i=@shf:7]+rm5[i=@shf:8]+rm5[i=@shf:9]+rm5[i=@shf:10]+rm5[i=@shf:11]+rm5[i=@shf:12])/24
let/bad=-999 godin_vector3=rm6[l=@sbx:25]

list/x=1:$dmax/noheader/file=$ofile/clobber/format=(6f12.5) yr,dd,dd2,time,godin_vector1-godin_vector2,godin_vector3
plot/thickness=2/symbol=1/color=black/x=1:1000/title="Vector(Godin_filter applied)" vvel1
plot/thickness=1/line/symbol=12/over/color=black/x=1:1000 godin_vector1
plot/thickness=2/symbol=1/color=red/over/x=1:1000/title="Vector(Godin_filter applied)" vvel2
plot/thickness=1/line/symbol=12/over/color=red/x=1:1000 godin_vector2
plot/thickness=2/symbol=1/color=blue/over/x=1:1000/title="Vector(Godin_filter applied)" mvel
plot/thickness=1/line/symbol=12/over/color=blue/x=1:1000 godin_vector3
sh d
quit
EOF

cat <<EOF>! gofile
go godin.jnl
quit
EOF
ls $ifile

ferret < gofile
Fprint metafile*plt -l cps -p portrait -o $fig ---> PSファイルの保存、cpsオプションでカラーPSになる
convert -density 150 $fig $fig2　　　　　　--> convertコマンドでJPGへ変換
rm -fr metafile*plt $fig

echo $ifile, $ofile
rm -fr ferret.jnl* godin.jnl gofile
endif
end
exit