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

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

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

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

}

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

sub_gt(){

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

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

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

}


# -----------------------------------------------------
# con東西平均データ, con東西平均からのアノマリデータ作成
# -----------------------------------------------------

sub_gtsub(){

  for list in t u v ps q relh
  do
    gtsel ${data}adj700con1000/${list}.sum x=0 && gtext out:${data}adj700con1000/${list}-zonal.sum x:GLON128
    gtsub ${dir}${list}.sum ${data}adj700con1000/${list}.sum out:${dir}${list}.anm
    gtsub ${dir}${list}.sum ${data}adj700con1000/${list}-zonal.sum out:${dir}${list}-zonal.anm
  done
  # 質量流線関数のアノマリデータ
  gtsub ${dir}strm.sum ${data}adj700con1000/strm.sum out:${dir}strm.anm

}


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

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

for name in pa
do

  dir=${data}adj700${name}1000/
  sub_gt
  sub_gtsub

  # t
  gtcont-y ${dir}t.anm x=0 -nocont tone=,-1.0,-0.6,-0.2,0.2,0.6,1.0 pat=25999,35999,45999,55999,70999,75999,85999 cont=0.2 DSET="${name}-A0" unit:"K, kg/s" TITLE="T and PSI" -print
  file=${figs}t-zonal-${name}-anm
  sub_figs

  # strm
  gtcont-y ${dir}strm.anm DSET="${name}-A0" unit:"K, kg/s" TITLE="T and PSI" -print
  file=${figs}strm-zonal-${name}-anm
  sub_figs

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

  # u
  gtcont-y ${dir}u.anm x=0  DSET="${name}-A0" tone=,-2.5,-1.5,-0.5,0.5,1.5,2.5 pat=25999,35999,45999,55999,70999,75999,85999 cont=0.5 -print
  file=${figs}u-zonal-${name}-anm
  sub_figs

  # v
  gtcont-y ${dir}v.anm x=0  DSET="${name}-A0" -print
  file=${figs}v-zonal-${name}-anm
  sub_figs

  # ps
  gtcurv-y ${dir}ps.anm x=0  DSET="${name}-A0" -print
  file=${figs}ps-zonal-${name}-anm
  sub_figs

  # q
  gtcont-y ${dir}q.anm x=0  DSET="${name}-A0" tone=,-25E-5,-15E-5,-5E-5,5E-5,15E-5,25E-5  pat=25999,35999,45999,55999,70999,75999,85999 cont=5E-5 -print
  file=${figs}q-zonal-${name}-anm
  sub_figs

  # relh
  gtcont-y ${dir}relh.anm x=0  DSET="${name}-A0" tone=,-25E-5,-15E-5,-5E-5,5E-5,15E-5,25E-5  pat=25999,35999      -print
#gtcont-y ${dir}relh.anm x=0  DSET="${name}-A0"  tone=,-25E-3,-15E-3,-5E-3,5E-3,15E-3,25E-3  pat=25999,35999,45999,55999,70999,75999,85999 cont=5E-3 lay=2 -print
  file=${figs}relh-zonal-${name}-anm
  sub_figs

done