#!/bin/csh

# c-shell script for MJO life cycle composite

# !!!!!! What should be modified !!!!!!
# HHH  : home directory
# var  : name of data
# TTT  : total numer of time (daily)
# LNX  : whether the machine is linux or not
#  This is for the record length problem.
#  (e.g. 4 : linux machine, 1 : other machine)
# MMM  : missing value
# PPP  : period of data

# !!!!!! What should be modified !!!!!!

# setenv HHH /jdata4/cdl/kim/msm_release/msd
#enver setenv HHH /jdata4/cdl/kim/msd
 setenv HHH ~/work/programs/MJOWG/msd


 foreach var ( ceof )

 setenv TTT 9855
 setenv LNX 1
 setenv MMM -999.

 setenv SN1 1
 setenv SN2 -1

 mkdir -p $HHH/level_2/comp/data

 mkdir -p $HHH/level_2/comp/$var/src
 cd $HHH/level_2/comp/$var/src

 cp -f $HHH/level_2/sample/comp/norm_pc_devide.f90.sample .
 cp -f $HHH/level_2/sample/comp/amp_pha.ctl.sample .

 sed "s/num_t/$TTT/g"   amp_pha.ctl.sample > tmp1
 sed "s/missing/$MMM/g"               tmp1 > amp_pha.ctl

 sed "s#homedir#$HHH#g" norm_pc_devide.f90.sample > tmp1
 sed "s/num_t/$TTT/g"                        tmp1 > tmp2
 sed "s/sign1/$SN1/g"                        tmp2 > tmp1
 sed "s/sign2/$SN2/g"                        tmp1 > tmp2
 sed "s/linux_recl/$LNX/g"                   tmp2 > norm_pc_devide.f90

f95 -ffree-line-length-none norm_pc_devide.f90 
a.out
rm -f a.out

 end