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

# ps, gif の置場
figs=/work/aqua3d/yukiko/nies2003/figs/tpa/
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東西平均データ
# -----------------------------------------------------

sub_gtzonal(){

  for list in t 
  do
    gtsel ${data}adj700con1000/${list}.sum x=0 && gtext out:${data}adj700con1000/${list}-zonal.sum x:GLON128

  done
}


# -----------------------------------------------------
# sigdot にファクターをかける
# -----------------------------------------------------

sub_gt(){

  # sigdot にファクターをかける
  gtset ${dir}sigdot.sum out:${dir}fact-sigdot.sum fact=-10000000

}


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

sub_gtsub(){

  for list in t u fact-sigdot
  do
    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: [中層 T] の図示
# -----------------------------------------------------

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

dir=${data}adj700con1000/
# sub_gtzonal
  sub_gt

for name in pa paNa paNb paNc paNd paEa paEb paEc paEd paCa paCb paCc paCd paCCa paCCb paCCc paCCd 
do

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

  # 中層 T sigma=0.549
  gtcont ${dir}t.anm z=8 DSET="${name}-A0" title:"t sigma=0.549" tone=,-1.5,-0.9,-0.3,0.3,0.9,1.5 pat=25999,35999,45999,55999,70999,75999,85999 cont=0.3 lay=2 -sg:lcorner=.false. -print  
  file=${figs}t-xymid-${name}-anm
  sub_figs

done


for name in pa paEa paEb paEc paEd paCa paCb paCc paCd paCCa paCCb paCCc paCCd 
do

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

  # T 赤道経度-高さ断面 
  gtvecon2 ${dir}t.anm ${dir}u.anm ${dir}fact-sigdot.anm DSET="${name}-A0" title:"t,(u,-sigdot*10**7)" intv=3,1 tone=,-1.0,-0.6,-0.2,0.2,0.6,1.0 pat=25999,35999,45999,55999,70999,75999,85999  cont=0.2 y=32 range=-5,5 lay=2 -sg:lcorner=.false. -print  
  file=${figs}t-wind-xzeq-${name}-anm
  sub_figs

done