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

# ps, gif の置場
figs=/home/yukiko/tmp/agcm-data/figs/tmp/
mkdir ${figs}

# データの置場
data=/home/yukiko/tmp/agcm-data/nrads/


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

}


# -----------------------------------------------------
# main: A0[qcnd, qradl, qvdf] の図示
# -----------------------------------------------------


#for exp in kuo adj
for exp in adj
do 
  for name in nrads
  do

  dir=${data}${exp}700${name}1000/
  cd ${dir}
  sub_gt

  # 凝結加熱 qcnd
  gtcont ${dir}qcnd.sum x=0 DSET="${exp}-${name}" cont=2E-6  tone=,-2E-6,2E-6,6E-6,10E-6,14E-6,18E-6 pat=25999,35999,45999,55999,70999,75999,85999 lay=1 -sg:lcorner=.false. -print 
  file=${figs}qcnd-zonal-${name}-${exp}
  sub_figs

  # 長波加熱 qradl
  gtcont ${dir}qradl.sum x=0 DSET="${exp}-${name}" cont=2E-6 tone=,-20E-6,-16E-6,-12E-6,-8E-6,-4E-6,0E-6 pat=25999,35999,45999,55999,70999,75999,85999 lay=1 -sg:lcorner=.false. -print 
  file=${figs}qradl-zonal-${name}-${exp}
  sub_figs

  # 拡散による加熱 qvdf 
  gtcont ${dir}qvdf.sum x=0 cont=2.5E-6 tone=,-5E-6,0,5E-6,10E-6,15E-6,20E-6 pat=25999,35999,45999,55999,70999,75999,85999 DSET="${exp}-${name}"  lay=1 -sg:lcorner=.false. -print 
  file=${figs}qvdf-zonal-${name}-${exp}
  sub_figs

done
done