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

# ps, gif の置場
figs=/work/aqua3d/yukiko/nies2003/figs/fig4pa/
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 *

}

# -----------------------------------------------------
# 温位(gttheta), 相当温位(gtthetae)の生成
# -----------------------------------------------------

sub_gt(){

  # 温位生成
  gttheta ${dir}t.sum ps:${dir}ps.sum out:${dir}theta.sum
  gtset ${dir}theta.sum out:${dir}theta.sum  divs=10 divl=50

  # 相当温位生成
  gtthetae t:${dir}t.sum q:${dir}q.sum ps:${dir}ps.sum out:${dir}thetae.sum 

  # 飽和比湿生成
  gtqsat t:${dir}t.sum ps:${dir}ps.sum out:${dir}qsat.sum 

  # 飽和相当温位生成
  gtthetae t:${dir}t.sum q:${dir}qsat.sum ps:${dir}ps.sum out:${dir}thetaesat.sum 
  gtset ${dir}thetaesat.sum out:${dir}thetaesat.sum  title:"Sat. Eqv. Pot. Temp."

  # ジオポテンシャル生成
  gthydro t:${dir}t.sum q:${dir}q.sum zs:${data}zs.given out:${dir}/gph.sum 

  # 安定度生成
  gtbrunt ${dir}theta.sum GPH:${dir}gph.sum out:${dir}brunt.sum

}


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

sub_gtsub(){

  for list in theta thetae thetaesat brunt
  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

}


# -----------------------------------------------------
# main: A0[theta, thetae, theatesat, n2] の図示
# -----------------------------------------------------

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

for name in pa
do

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

  # 赤道上 経度-高さ図

  # theta
  gtcont-y ${dir}theta.anm y=32 DSET="${name}-A0" tone=,-10E-1,-6E-1,-2E-1,2E-1,6E-1,10E-1  pat=25999,35999,45999,55999,70999,75999,85999 cont=2E-1 -print
  file=${figs}theta-xzeq-${name}-anm
  sub_figs

  # thetae
  gtcont-y ${dir}thetae.anm y=32 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}thetae-xzeq-${name}-anm
  sub_figs

  # thetaesat
  gtcont-y ${dir}thetaesat.anm y=32 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}thetaesat-xzeq-${name}-anm
  sub_figs

  # stability N^2
  gtcont-y ${dir}brunt.anm y=32 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}brunt-xzeq-${name}-anm
  sub_figs


  # 帯状平均図 (いる!?)

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

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

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

  # stability N^2
  gtcont-y ${dir}brunt.anm x=0 DSET="${name}-A0" -print 
  file=${figs}brunt-zonal-${name}-anm
  sub_figs


done