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

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

# データの置場
data=/work/aqua3d/yukiko/kuo/kuo-baseSST/

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

}

# -----------------------------------------------------
# 質量流線関数(strm), 相対比湿(relh)の生成, 比湿(q) の設定
# -----------------------------------------------------

sub_gt(){

  # 質量子午面流線関数生成
  gtstrm v.sum ps:ps.sum out:strm.sum 
  gtset strm.sum divs=1.5E+10 divl=7.5E+10 out:strm.sum

  # 相対比湿生成
  gtrelh q:q.sum t:t.sum ps:ps.sum out:relh.sum

  # 比湿のコンター設定
  gtset q.sum divs=0.002 divl=0.01 out:q.sum


}

# -----------------------------------------------------
# main: A0[t, strm, u, v, ps, q, relh] の図示
# -----------------------------------------------------

for name in Aa Ab Ac Ad con Ae Af Ag
do

  dir=${data}kuo700${name}1000/
  cd ${dir}
  sub_gt

  # t
  gtcont t.sum x=0 -nocont color=25 range=190,300 cont=5 DSET="kuo ${name}" unit:"K, kg/s" TITLE="T and PSI" lay=1 -sg:lcorner=.false. -print 
  file=${figs}t-zonal-${name}-kuo
  sub_figs

  # strm
  gtcont strm.sum DSET="kuo ${name}" unit:"K, kg/s" TITLE="T and PSI" lay=1 -sg:lcorner=.false. -print 
  file=${figs}strm-zonal-${name}-kuo
  sub_figs

  # T and strm
  comdclps.pl ${figs}t-zonal-${name}-kuo.ps ${figs}strm-zonal-${name}-kuo.ps > dcl.ps
  file=${figs}t-strm-zonal-${name}-kuo
  sub_figs

  # u
  gtcont u.sum x=0 tone=,-25,-15,-5,5,15,25 pat=25999,35999,45999,55999,70999,75999,85999 DSET="kuo ${name}" lay=1 -sg:lcorner=.false. -print 
  file=${figs}u-zonal-${name}-kuo
  sub_figs

  # v
  gtcont v.sum x=0 cont=0.5 DSET="kuo ${name}"  tone=,-2.5,-1.5,-0.5,0.5,1.5,2.5 pat=25999,35999,45999,55999,70999,75999,85999 lay=1 -sg:lcorner=.false. -print 
  file=${figs}v-zonal-${name}-kuo
  sub_figs

  # ps
  gtcurv-y ps.sum x=0  DSET="kuo ${name}" tick=F lay=1 -sg:lcorner=.false. -print 
  file=${figs}ps-zonal-${name}-kuo
  sub_figs

  # q
  gtcont q.sum x=0  DSET="kuo ${name}" cont=1E-3 tone=,2E-3,4E-3,6E-3,8E-3,10E-3,12E-3 pat=25999,35999,45999,55999,70999,75999,85999 lay=1 -sg:lcorner=.false. -print 
  file=${figs}q-zonal-${name}-kuo
  sub_figs

  # relh
  gtcont relh.sum x=0 cont=5E-02 DSET="kuo ${name}" tone=,10E-2,20E-2,30E-2,40E-2,50E-2,60E-2 pat=25999,35999,45999,55999,70999,75999,85999 lay=1 -sg:lcorner=.false. -print 
  file=${figs}relh-zonal-${name}-kuo
  sub_figs

done