#!/bin/sh
#
# 表題: SST データ作成 
# 履歴: 2003.05.04 yukiko@ep.sci.hokudai.ac.jp
# 
# tako:/work/aqua3d/yukiko/sh/mkfig-sst.sh
#

# ps, gif の置場
figs=/work/aqua3d/yukiko/figs
mkdir ${figs}

# データの置場
data=/work/aqua3d/yukiko/sst/

# -----------------------------------------------------
# ps -> gif ファイル, ps ファイル生成
# -----------------------------------------------------

sub_figs(){

  pstopnm dcl.ps 
  ppmtogif dcl001.ppm >  dcl.gif 
  gifsicle --rotate-90 dcl.gif > ${file}.gif

  mv dcl.ps ${file}.ps
  rm *

}


# -----------------------------------------------------
# データ item の指定
# -----------------------------------------------------

sub_item(){

   gtset ${data}sst${list}.init out:${data}sst${list}.init item=${list}

}


# -----------------------------------------------------
# main: A4[sstpa.init] の図示
# -----------------------------------------------------

cd /work/aqua3d/yukiko/tmp/
rm *

for name in con pa paNa paNb paNc paNd paEa paEb paEc paEd paCa paCb paCc paCd paCCa paCCb paCCc paCCd 
do

  sub_item

  # コンター図
  gtcont ${data}sst${name}.init cont=2 color=26 lay=2 -sg:lcorner=.false. -print 
  file=${figs}sst-xy-${name}
  sub_figs

  # 子午面図
  gtcurv-y ${data}sst${name}.init x=64 tick=F lay=2 -sg:lcorner=.false. -print ps:dcl.ps
  file=${figs}sst-y-${name}
  sub_figs

done


# 北ずらし子午面図
gtcurv-y ${data}sstpaNa.init ${data}sstpaNb.init ${data}sstpaNc.init ${data}sstpaNd.init ${data}sstpa.init lidx=21,31,41,61,11 ltype=1,1,1,1 x=64 DSET="paN" tick=F lay=2 -sg:lcorner=.false. -print 
file=${figs}sst-y-paN
sub_figs


# 赤道上経度面図
gtcurv-y ${data}sstpaEa.init ${data}sstpaEb.init ${data}sstpaEc.init ${data}sstpaEd.init lidx=11,21,31,41 ltype=1,1,1,1 y=32 DSET="paE" tick=F lay=2 -sg:lcorner=.false. -print 
file=${figs}sst-xeq-paE
sub_figs

# Cos(x) 波長1 赤道上経度面図
gtcurv-y ${data}sstpaCa.init ${data}sstpaCb.init ${data}sstpaCc.init ${data}sstpaCd.init ${data}sstpa.init lidx=21,31,41,61,11 ltype=1,1,1,1 y=32 DSET="paC" tick=F lay=2 -sg:lcorner=.false. -print 
file=${figs}sst-xeq-paC
sub_figs

# Cos(x) 波長2 赤道上経度面図
gtcurv-y ${data}sstpaCCa.init ${data}sstpaCCb.init ${data}sstpaCCc.init ${data}sstpaCCd.init ${data}sstpaCa.init lidx=21,31,41,61,11 ltype=1,1,1,1 y=32 DSET="paCC" tick=F lay=2 -sg:lcorner=.false. -print 
file=${figs}sst-xeq-paCC
sub_figs