#!/bin/sh
#
# 表題: 時間平均 flux 図作成
# 履歴: 2003.05.04 yukiko@ep.sci.hokudai.ac.jp
# 
# tako:/work/aqua3d/yukiko/sh/mkfig-flux.sh
#

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

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

# -----------------------------------------------------
# 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 rain evap olr slr sens
  do
    gtsel ${data}con/adj700con1000/${list}.sum x=0 && gtext out:${data}/con/adj700con1000/${list}-zonal.sum x:GLON128
    gtsub ${dir}${list}.sum ${data}con/adj700con1000/${list}-zonal.sum out:${dir}${list}-zonal.anm
  done

}

# -----------------------------------------------------
# データ item の指定
# -----------------------------------------------------

sub_item(){

  for list in rain evap olr slr sens
  do
    gtset ${dir}${list}.sum out:${dir}${list}.sum item=${list}
#    gtset ${dir}${list}-zonal.anm out:${dir}${list}-zonal.anm item=${list}
  done

}

# -----------------------------------------------------
# main: A0[フラックス] の図示
# -----------------------------------------------------

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

con pa

# for name in pa

for exp in paN paE paC paCC
do 
  for alph in a b c d
  do

    name=${exp}${alph}
    dir=${data}${exp}/tmean/adj700${name}1000/
#    sub_gtsub
    sub_item

    # 子午面 各フラックス図
    gtcurv-y ${dir}rain.sum ${dir}evap.sum ${dir}olr.sum ${dir}slr.sum ${dir}sens.sum  lidx=23,33,43,53,13 ltype=1,1,1,1,1 x=0 range=0,300  div=20,100  DSET="${name}" title:"heat flux" tick=F lay=2 -sg:lcorner=.false. -print
    file=${figs}flux-zonal-${name}
    sub_figs


    # 子午面 各フラックス図 - アノマリ
#    gtcurv-y ${dir}rain.anm ${dir}evap.anm ${dir}olr.anm ${dir}slr.anm ${dir}sens.anm  lidx=23,33,43,53,13 ltype=1,1,1,1,1 x=0 range=-10,20  div=1,10  DSET="${name}-A0" title:"heat flux" -print
#    file=${figs}flux-zonal-${name}-anm
#    sub_figs

  done
done

# 子午面 各フラックス図 con, pa
for name in con pa 
do

  dir=${data}${name}/adj700${name}1000/
#    sub_gtsub
  sub_item

  gtcurv-y ${dir}rain.sum ${dir}evap.sum ${dir}olr.sum ${dir}slr.sum ${dir}sens.sum  lidx=23,33,43,53,13 ltype=1,1,1,1,1 x=0 range=0,300  div=20,100  DSET="${name}" title:"heat flux" tick=F lay=2 -sg:lcorner=.false. -print
  file=${figs}flux-zonal-${name}
  sub_figs

    # 子午面 各フラックス図 - アノマリ
#    gtcurv-y ${dir}rain.anm ${dir}evap.anm ${dir}olr.anm ${dir}slr.anm ${dir}sens.anm  lidx=23,33,43,53,13 ltype=1,1,1,1,1 x=0 range=-10,20  div=1,10  DSET="${name}-con" title:"heat flux" -print
#    file=${figs}flux-zonal-${name}-anm
#    sub_figs

done