#!/bin/sh
#
# 表題: x-t ダイヤグラム 
# 履歴: 2003.06.07 yukiko@ep.sci.hokudai.ac.jp
# 
# tako:/work/aqua3d/yukiko/sh/mkfig-fig3.sh
#

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

# データの置場
exp=pa100
data=/work/aqua3d/yukiko/${exp}/

# -----------------------------------------------------
# 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 dcl.gif dcl001.ppm 

}


# -----------------------------------------------------
# ps, gph, u の東西平均からの偏差をとる...
# -----------------------------------------------------

sub_gt(){

  for list in gph u ps 
  do

    gtsel ${list}.out.seq x=0 && gtext out:${list}-zonal.out.seq x:GLON128
    gtsub ${list}.out.seq ${list}-zonal.out.seq out:${list}.out.seq.anm 

  done

}



# -----------------------------------------------------
# main: A0[ps, gph, u, v, rain] の図示
# -----------------------------------------------------

for num in 1 2 3 4 5 6 7 8 9 10

do 

for name in pa${num}00
#${exp}a${num}00 ${exp}b${num}00 ${exp}c${num}00 ${exp}d${num}00
do

  dir=${data}adj700${name}/
  cd ${dir}
  sub_gt

  # ps 
  gtcont ps.out.seq.anm  -nocont  tone=,-2.5,-1.5,-0.5,0.5,1.5,2.5 pat=25999,35999,45999,55999,70999,75999,85999 lay=2 -sg:lcorner=.false. -print
  file=${figs}ps-xteq-${name}-anm
  sub_figs

  # u y=1 
  gtcont u.out.seq.anm y=1 -nocont tone=,-5,-3,-1,1,3,5 pat=25999,35999,45999,55999,70999,75999,85999 lay=2 -sg:lcorner=.false. -print
  file=${figs}u-xteq-z1-${name}-anm
  sub_figs

  # u y=5 
  gtcont u.out.seq.anm y=5 -nocont tone=,-10,-6,-2,2,6,10 pat=25999,35999,45999,55999,70999,75999,85999 -sg:lcorner=.false. lay=2 -print
  file=${figs}u-xteq-z5-${name}-anm
  sub_figs

  # u y=13
  gtcont u.out.seq.anm y=13 -nocont  tone=,-10,-6,-2,2,6,10 pat=25999,35999,45999,55999,70999,75999,85999 -sg:lcorner=.false. lay=2 -print
  file=${figs}u-xteq-z13-${name}-anm
  sub_figs

  # gph z=5
  gtcont gph.out.seq.anm y=5 -nocont  tone=,-10,-6,-2,2,6,10 pat=25999,35999,45999,55999,70999,75999,85999 lay=2 -sg:lcorner=.false. -print
  file=${figs}gph-xteq-z5-${name}-anm
  sub_figs

  # gph z=13
  gtcont gph.out.seq.anm y=13 -nocont tone=,-50,-30,-10,10,30,50 pat=25999,35999,45999,55999,70999,75999,85999 lay=2 -sg:lcorner=.false. -print
  file=${figs}gph-xteq-z13-${name}-anm
  sub_figs

  # rain 
  gtcont rain.sum.seq  tone=50,250,450,650,850,1050,1250 pat=25999,35999,45999,55999,70999,75999,85999 -nocont lay=2 -sg:lcorner=.false. -print
  file=${figs}rain-xteq-${name}
  sub_figs

done

done