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

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

}


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

sub_gtsub(){

  for list in qcnd qradl qvdf
  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[scnd, qradl, qvdf] の図示
# -----------------------------------------------------

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

for name in pa
do

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


  # 赤道上 経度-高さ図

  # 凝結加熱 qcnd
  gtcont-y ${dir}qcnd.anm y=32 DSET="${name}-A0" tone=,-25E-6,-15E-6,-5E-6,5E-6,15E-6,25E-6  pat=25999,35999,45999,55999,70999,75999,85999 cont=5E-6 -print
  file=${figs}qcnd-xzeq-${name}-anm
  sub_figs

  # 長波加熱 qradl
  gtcont-y ${dir}qradl.anm y=32 DSET="${name}-A0" tone=,-5E-6,-3E-6,-1E-6,1E-6,3E-6,5E-6  pat=25999,35999,45999,55999,70999,75999,85999 cont=1E-6 -print
  file=${figs}qradl-xzeq-${name}-anm
  sub_figs

  # 拡散による加熱 qvdf 
  gtcont-y ${dir}qvdf.anm y=32 DSET="${name}-A0" tone=,-25E-6,-15E-6,-5E-6,5E-6,15E-6,25E-6  pat=25999,35999,45999,55999,70999,75999,85999 cont=5E-6 -print
  file=${figs}qvdf-xzeq-${name}-anm
  sub_figs


  # 帯状平均図

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

  # 長波加熱 qradl
  gtcont-y ${dir}qradl.anm x=0 DSET="${name}-A0" tone=,-5E-7,-3E-7,-1E-7,1E-7,3E-7,5E-7  pat=25999,35999,45999,55999,70999,75999,85999 cont=1E-7 -print

###  gtcont-y ${dir}qradl.anm x=0 DSET="${name}-A0" tone=,-10E-7,-6E-7,-2E-7,2E-7,6E-7,10E-7  pat=25999,35999,45999,55999,70999,75999,85999 cont=2E-7 -print

  file=${figs}qradl-zonal-${name}-anm
  sub_figs

  # 拡散による加熱 qvdf 
  gtcont-y ${dir}qvdf.anm x=0 DSET="${name}-A0" tone=,-5E-7,-3E-7,-1E-7,1E-7,3E-7,5E-7  pat=25999,35999,45999,55999,70999,75999,85999 cont=1E-7 -print
  file=${figs}qvdf-zonal-${name}-anm
  sub_figs

done