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

}


# -----------------------------------------------------
# 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=${exp}-${name}-${list}
#    gtset ${dir}${list}-zonal.anm out:${dir}${list}-zonal.anm item=${list}
  done

}

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

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

    dir=${data}${exp}700${name}1000/
    cd $dir
#    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=1 -sg:lcorner=.false. -print
    file=${figs}flux-zonal-${name}-${exp}
    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

# 各フラックス比較図

dir1=${data}adj700con1000/
dir2=${data}adj700nrads1000/

# 子午面 slr
gtcurv-y ${dir1}slr.sum ${dir2}slr.sum lidx=23,33 ltype=1,1 x=0 range=100,200  div=10,50  DSET="slr" tick=F lay=1 -sg:lcorner=.false. -print
    file=${figs}slr-zonal-nrads
    sub_figs

# 子午面 olr
gtcurv-y ${dir1}olr.sum ${dir2}olr.sum lidx=23,33 ltype=1,1 x=0 range=150,320  div=10,50  DSET="olr" tick=F lay=1 -sg:lcorner=.false. -print
    file=${figs}olr-zonal-nrads
    sub_figs

# 子午面 rain
gtcurv-y ${dir1}rain.sum ${dir2}rain.sum lidx=23,33 ltype=1,1 x=0 range=0,160  div=10,50  DSET="rain" tick=F lay=1 -sg:lcorner=.false. -print
    file=${figs}rain-zonal-nrads
    sub_figs

# 子午面 evap
gtcurv-y ${dir1}evap.sum ${dir2}evap.sum lidx=23,33 ltype=1,1 x=0 range=0,120  div=10,50  DSET="evap" tick=F lay=1 -sg:lcorner=.false. -print
    file=${figs}evap-zonal-nrads
    sub_figs

# 子午面 sens
gtcurv-y ${dir1}sens.sum ${dir2}sens.sum lidx=23,33 ltype=1,1 x=0 range=-10,30  div=2,10  DSET="sens" tick=F lay=1 -sg:lcorner=.false. -print
    file=${figs}sens-zonal-nrads
    sub_figs