#!/bin/sh
# SST パラメタ研究実験データ作成
# 2002.01.25 yukiko@ep.sci.hokudai.ac.jp
# tako 内で使う. con にかんしては真面目につくってない

cd /work/aqua3d/yukiko/timeavg/parameter/

for file in Aa Ab Ac Ad Ae Af Ag con
do
  # 100 日時系列データを用い,
  # u, t 東西平均データの作成
#  gtsel ../adj700${file}400/u.out x=0 >> err.log && gtext out:../adj700${file}400/u.zonal.out x:GLON128  >> err.log 
#  gtsel ../adj700${file}400/t.out x=0  >> err.log  && gtext out:../adj700${file}400/t.zonal.out x:GLON128  >> err.log
#  gtsel ../adj700${file}400/rain.sum x=0 >> err.log && gtext out:../adj700${file}400/rain.zonal.sum x:GLON128  >> err.log
#  echo u, t, rain zonal data criate 
#  date

  # u, t 東西平均からの差
#  gtsub ../adj700${file}400/u.out ../adj700${file}400/u.zonal.out out:$file/u.zonal.anm  >> err.log 
#  gtsub ../adj700${file}400/t.out ../adj700${file}400/t.zonal.out out:$file/t.zonal.anm  >> err.log 
#  gtsub ../adj700${file}400/rain.sum ../adj700${file}400/rain.zonal.sum out:$file/rain.zonal.sum.anm >> err.log
#  echo u, t, rain zonalsub data criate
#  date

  # rain, u, t 赤道上時系列データの生成
#  gtsel y=32 ../adj700${file}400/rain.out  >> err.log  && gtseq out:${file}/rain.out.seq  >> err.log 
#  gtsel y=32 ${file}/u.zonal.anm  >> err.log  && gtseq out:${file}/u.zonal.anm.seq >> err.log 
#  gtsel y=32 ${file}/t.zonal.anm  >> err.log && gtseq out:${file}/t.zonal.anm.seq >> err.log 
#  gtsel y=32 ${file}/rain.zonal.sum.anm  >> err.log && gtseq out:${file}/rain.zonal.sum.anm.seq >> err.log 
#  echo seq data end
#  date

  # 700 日時間平均を用い, 
  # 質量子午面流線関数生成
#  gtstrm ../adj700${file}1000/v.sum ps:../adj700${file}1000/ps.sum out:../adj700${file}1000/strm.sum >> err.log 
# gtset ../adj700${file}1000/strm.sum divs=1.5E+10 divl=7.5E+10 out:../adj700${file}1000/strm.sum >> err.log 

  # u, t, strm 東西平均データの作成
#  gtsel ../adj700${file}1000/u.sum x=0  >> err.log && gtext out:../adj700${file}1000/u.zonal.sum x:GLON128 >> err.log 
#  gtsel ../adj700${file}1000/t.sum x=0  >> err.log && gtext out:../adj700${file}1000/t.zonal.sum x:GLON128 >> err.log 
#  gtsel ../adj700${file}1000/strm.sum x=0 >> err.log &&   gtext out:../adj700${file}1000/strm.zonal.sum x:GLON128 >> err.log 

#  echo u, t, strm zonal data end
#  date
#  echo -----------

gtset ../adj700${file}1000/rain.sum out:../adj700${file}1000/rain.sum item=${file}rain
gtset ../adj700${file}1000/evap.sum out:../adj700${file}1000/evap.sum item=${file}evap
gtset ../adj700${file}1000/olr.sum out:../adj700${file}1000/olr.sum item=${file}olr
gtset ../adj700${file}1000/slr.sum out:../adj700${file}1000/slr.sum item=${file}slr
gtset ../adj700${file}1000/sens.sum out:../adj700${file}1000/sens.sum item=${file}sens

done

#gtsel ../adj1000con100/rain.sum x=0 && gtext out:../adj1000con100/rain.zonal.sum x:GLON128
#gtsub ../adj1000con100/rain.sum ../adj1000con100/rain.zonal.sum out:con/rain.zonal.sum.anm
#gtsel y=32 con/rain.zonal.sum.anm  >> err.log && gtseq out:con/rain.zonal.sum.anm.seq >> err.log