#!/bin/sh
#
# 表題: 時間平均実験 Fig6 図作成
# 履歴: 2003.05.05 yukiko@ep.sci.hokudai.ac.jp
# 
# tako:/work/aqua3d/yukiko/nies2003/sh/mkfig-fig6.sh
#

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

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

# -----------------------------------------------------
# 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 *

}

# -----------------------------------------------------
# 1000 日目のみのデータの切り出し, sigdot にファクターをかける
# -----------------------------------------------------

sub_gt(){

  cd ${data}
  mkdir snp1000
  for list in *.out *.sum
  do
    # データの切り出し
    gtsel ${list} str=400 end=400 out:snp1000/${list}
  done

  # sigdot にファクターをかける
  gtset ${dir}sigdot.out out:${dir}fact-sigdot.out fact=-1000000
  cd /work/aqua3d/yukiko/nies2003/tmp/

}

# -----------------------------------------------------
# main: A0 snapshot[ps, rain, T and [u,-sigdot] ] の図示
# -----------------------------------------------------

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

dir=${data}snp1000/
sub_gt

# ps
gtcont ${dir}ps.out tone=,960,970,980,990,1000,1010 pat=25999,35999,45999,55999,70999,75999,85999 DSET:"A0 snapshot" lay=2 -sg:lcorner=.false. -print 
file=${figs}ps-xy-con1000
sub_figs

# rain
gtcont ${dir}rain.out tone=50,250,450,650,850,1050,1250 pat=25999,35999,45999,55999,70999,75999,85999 -nocont DSET:"A0 snapshot" lay=2 -sg:lcorner=.false. -print 
file=${figs}rain-xy-con1000
sub_figs

# t, [u, -sigdot]
gtvecon2 ${dir}t.out ${dir}u.out ${dir}fact-sigdot.out  title:"t,(u,-sigdot*10**-6)" intv=3,1 -nocont  y=32 fact=3 tone=,230,240,250,260,270,280 pat=25999,35999,45999,55999,70999,75999,85999 DSET:"A0 snapshot" lay=2 -sg:lcorner=.false. -print
file=${figs}t-wind-xzeq-con1000
sub_figs