#!/bin/sh
############################################################
#to   prepare crop calendar calculation
#by   2010/03/31, hanasaki, NIES: H08 ver1.0
############################################################
# Edit here (basic)
############################################################
PRJ=L04_
RUN=____
YEAR=0000
MON=00
DAY=00
SUF=.one
############################################################
# Edit here (geography) 
############################################################
L=64800
XY="360 180"
L2X=../../map/dat/l2x_l2y_/l2x.one.txt
L2Y=../../map/dat/l2x_l2y_/l2y.one.txt
LONLAT="-180 180 -90 90"
############################################################
# Edit here (in)
############################################################
DIRORG=../../map/org/L04
DIRDAT=../../map/dat/hvs_ara_
HA01=L04_bar_.one
HA02=L04_cas_.one
HA03=L04_cot_.one
HA04=L04_grn_.one
HA05=L04_mai_.one
HA06=L04_mil_.one
HA07=L04_oil_.one
HA08=L04_oth_.one
HA09=L04_pot_.one
HA10=L04_pul_.one
HA11=L04_rap_.one
HA12=L04_ric_.one
HA13=L04_rye_.one
HA14=L04_sor_.one
HA15=L04_soy_.one
HA16=L04_sub_.one
HA17=L04_suc_.one
HA18=L04_sun_.one
HA19=L04_whe_.one
############################################################
# Edit here (out)
############################################################
DIRCRPTYP1ST=../../map/out/crp_typ1
DIRCRPTYP2ND=../../map/out/crp_typ2
DIRCRPTYP1STFIG=../../map/fig/crp_typ1
DIRCRPTYP2NDFIG=../../map/fig/crp_typ2
CRPTYP1ST=${DIRCRPTYP1ST}/${PRJ}${RUN}${YEAR}${MON}${DAY}${SUF}
CRPTYP2ND=${DIRCRPTYP2ND}/${PRJ}${RUN}${YEAR}${MON}${DAY}${SUF}
CRPTYP1STFIG=${DIRCRPTYP1STFIG}/${PRJ}${RUN}${YEAR}${MON}${DAY}.png
CRPTYP2NDFIG=${DIRCRPTYP2NDFIG}/${PRJ}${RUN}${YEAR}${MON}${DAY}.png
############################################################
# Make output directory
############################################################
if [ ! -d $DIRDAT ]; then
  mkdir -p $DIRDAT
fi
if [ ! -d $DIRCRPTYP1ST ]; then
  mkdir -p $DIRCRPTYP1ST
fi
if [ ! -d $DIRCRPTYP2ND ]; then
  mkdir -p $DIRCRPTYP2ND
fi
if [ ! -d $DIRCRPTYP1STFIG ]; then
  mkdir -p $DIRCRPTYP1STFIG
fi
if [ ! -d $DIRCRPTYP2NDFIG ]; then
  mkdir -p $DIRCRPTYP2NDFIG
fi
############################################################
# Binary
############################################################
htformat $L $XY $L2X $L2Y $LONLAT asciiu binary $DIRORG/$HA01.txt $DIRDAT/$HA01
htformat $L $XY $L2X $L2Y $LONLAT asciiu binary $DIRORG/$HA02.txt $DIRDAT/$HA02
htformat $L $XY $L2X $L2Y $LONLAT asciiu binary $DIRORG/$HA03.txt $DIRDAT/$HA03
htformat $L $XY $L2X $L2Y $LONLAT asciiu binary $DIRORG/$HA04.txt $DIRDAT/$HA04
htformat $L $XY $L2X $L2Y $LONLAT asciiu binary $DIRORG/$HA05.txt $DIRDAT/$HA05
htformat $L $XY $L2X $L2Y $LONLAT asciiu binary $DIRORG/$HA06.txt $DIRDAT/$HA06
htformat $L $XY $L2X $L2Y $LONLAT asciiu binary $DIRORG/$HA07.txt $DIRDAT/$HA07
htformat $L $XY $L2X $L2Y $LONLAT asciiu binary $DIRORG/$HA08.txt $DIRDAT/$HA08
htformat $L $XY $L2X $L2Y $LONLAT asciiu binary $DIRORG/$HA09.txt $DIRDAT/$HA09
htformat $L $XY $L2X $L2Y $LONLAT asciiu binary $DIRORG/$HA10.txt $DIRDAT/$HA10
htformat $L $XY $L2X $L2Y $LONLAT asciiu binary $DIRORG/$HA11.txt $DIRDAT/$HA11
htformat $L $XY $L2X $L2Y $LONLAT asciiu binary $DIRORG/$HA12.txt $DIRDAT/$HA12
htformat $L $XY $L2X $L2Y $LONLAT asciiu binary $DIRORG/$HA13.txt $DIRDAT/$HA13
htformat $L $XY $L2X $L2Y $LONLAT asciiu binary $DIRORG/$HA14.txt $DIRDAT/$HA14
htformat $L $XY $L2X $L2Y $LONLAT asciiu binary $DIRORG/$HA15.txt $DIRDAT/$HA15
htformat $L $XY $L2X $L2Y $LONLAT asciiu binary $DIRORG/$HA16.txt $DIRDAT/$HA16
htformat $L $XY $L2X $L2Y $LONLAT asciiu binary $DIRORG/$HA17.txt $DIRDAT/$HA17
htformat $L $XY $L2X $L2Y $LONLAT asciiu binary $DIRORG/$HA18.txt $DIRDAT/$HA18
htformat $L $XY $L2X $L2Y $LONLAT asciiu binary $DIRORG/$HA19.txt $DIRDAT/$HA19
############################################################
# Make set file
############################################################
SETFILE=temp.txt
if [ -f $SETFILE ]; then
  rm $SETFILE
fi
#
cat << EOF >> $SETFILE
&setcal
n0l=$L
c1hvsara(1)='$DIRDAT/$HA01'
c1hvsara(2)='$DIRDAT/$HA02'
c1hvsara(3)='$DIRDAT/$HA03'
c1hvsara(4)='$DIRDAT/$HA04'
c1hvsara(5)='$DIRDAT/$HA05'
c1hvsara(6)='$DIRDAT/$HA06'
c1hvsara(7)='$DIRDAT/$HA07'
c1hvsara(8)='$DIRDAT/$HA08'
c1hvsara(9)='$DIRDAT/$HA09'
c1hvsara(10)='$DIRDAT/$HA10'
c1hvsara(11)='$DIRDAT/$HA11'
c1hvsara(12)='$DIRDAT/$HA12'
c1hvsara(13)='$DIRDAT/$HA13'
c1hvsara(14)='$DIRDAT/$HA14'
c1hvsara(15)='$DIRDAT/$HA15'
c1hvsara(16)='$DIRDAT/$HA16'
c1hvsara(17)='$DIRDAT/$HA17'
c1hvsara(18)='$DIRDAT/$HA18'
c1hvsara(19)='$DIRDAT/$HA19'
c0crptyp1st='$CRPTYP1ST'
c0crptyp2nd='$CRPTYP2ND'
&end
EOF
############################################################
# Calculate
############################################################
prog_crptyp $SETFILE
############################################################
# Draw
############################################################
#CPT=../../cpt/crptyp.cpt
#htdraw $L $XY $L2X $L2Y $LONLAT $CRPTYP1ST  $CPT temp.eps
#htconv temp.eps $CRPTYP1STFIG rot
#htdraw $L $XY $L2X $L2Y $LONLAT $CRPTYP2ND  $CPT temp.eps
#htconv temp.eps $CRPTYP2NDFIG rot