# makefile for dcpam SyncRot
#
# 2009-10-01 S. Noda
# 2009-08-21 S. Noda
# 2009-08-18 S. Noda
# 2009-08-17 S. Noda
# 2009-07-09 S. Noda
# 2009-06-08 S. Noda
# 2009-05-30 S. Noda
# 2008-11-24 Y. Morikawa
# 2007-06-14 M. Ishiwatari
# 2007-06-13 M. Ishiwatari
# 2007-05-30 M. Ishiwatari
#

#GPVIEW=gpview
GPVIEW=./gpview_custom-ym06

.SUFFIXES : .ps .png .obj .eps

.ps.eps :
#	dclpsrmcm $< | dclpsrot | sed 's/2\.0 W/13.0 W/g' | sed 's/6\.0 W/20.0 W/g' > tmp.ps
	dclpsrot $< | sed 's/2\.0 W/13.0 W/g' | sed 's/6\.0 W/20.0 W/g' > tmp.ps
	eps2eps tmp.ps $@
	rm tmp.ps

.ps.png :
	convert -rotate 90 $< $@

STRTIME = 0
ENDTIME = 2000
AVRSTRTIME = 1000
AVRENDTIME = 2000
AVRDELTATIME = 10   # 動画に使う画像の時間平均幅
ANIMOPTDELAY = 50   # 動画のコマ送り間隔. 1/100 秒単位.
ANIMSTRTIME = $(STRTIME)
ANIMENDTIME = $(ENDTIME)

ANIMS = \
Tg_anim-xy.gif \
Rain_anim-xy.gif \
Temp_anim-xy-sigma0.5.gif \
U_anim-xy-sigma0.2.gif \
Ps_anim-xy.gif


IMGS = \
Tg_timeseries-xymean \
ps_timeseries-xymean \
OLR_timeseries-xymean \
SLR_timeseries-xymean \
Rain_timeseries-xymean \
EVap_timeseries-xymean \
Sens_timeseries-xymean \
Tg-x0-180_timeseries-xymean \
ps-x0-180_timeseries-xymean \
OLR-x0-180_timeseries-xymean \
SLR-x0-180_timeseries-xymean \
Rain-x0-180_timeseries-xymean \
EVap-x0-180_timeseries-xymean \
Sens-x0-180_timeseries-xymean \
Tg-x180-360_timeseries-xymean \
ps-x180-360_timeseries-xymean \
OLR-x180-360_timeseries-xymean \
SLR-x180-360_timeseries-xymean \
Rain-x180-360_timeseries-xymean \
EVap-x180-360_timeseries-xymean \
Sens-x180-360_timeseries-xymean \
flux_timeseries-xymean \
flux-x0-180_timeseries-xymean \
flux-x180-360_timeseries-xymean \
Tg_xy \
ps_xy \
SSR_xy \
SLR_xy \
Rain_xy \
OLR_xy \
Sens_xy \
EVap_xy \
u_xz-y0 \
u_yz-xmean \
u_yz-x90 \
u_yz-x270 \
v_yz-xmean \
v_yz-x270 \
v_yz-x90 \
T_xy-z0 \
T_xy-sigma0.7 \
T_yz-xmean \
T_yz-x270 \
T_yz-x90 \
T_xz-y0 \
SigDot_xz-y0 \
q_xy-z0 \
q_yz-xmean \
q_yz-x270 \
q_yz-x90 \
SigDot_xz-y0 \
u_xy-sigma0.5 \
v_xy-sigma0.5 \
v_xy-sigma0.2 \
surface-flux_y0 \
surface-flux_y-xmean \
surface-flux_y-x0-180 \
surface-flux_y-x180-360

PSFIGS = $(IMGS:%=%.ps)
EPSFIGS = $(IMGS:%=%.eps)
PNGFIGS = $(IMGS:%=%.png)

all:ps eps png anim

ps: $(PSFIGS)
eps: $(EPSFIGS)
png: $(PNGFIGS)
anim: $(ANIMS)

#ps_timeseries-xymean.ps \

# SSR_timeseries-xymean.ps \
# SSR-x0-180_timeseries-xymean.ps \
# SSR-x180-360_timeseries-xymean.ps \



# 熱収支, 水収支の表を別 html に出力 (出力先のパスに注意)
budget.txt: diagnosis_budget.rb
	ruby diagnosis_budget.rb > ../budget.htm
	echo '* ((<"熱収支と水収支"|URL:./budget.htm>))' > $@




# 動画
Tg_anim-xy.gif: SurfTemp.nc
	anim_str=$(ANIMSTRTIME);                                     \
	anim_end=`expr $(ANIMENDTIME) - $(AVRDELTATIME)`;            \
	timelist=`seq -w $${anim_str} $(AVRDELTATIME) $${anim_end}`; \
	for t in $${timelist}; do                                    \
          t_str=`expr $${t} + 0`;                                    \
	  t_end=`expr $${t} + $(AVRDELTATIME)`;                      \
	  echo time: $${t_str} - $${t_end};                          \
	  $(GPVIEW) SurfTemp.nc@SurfTemp,time=$${t_str}:$${t_end}    \
	    --mean time --wsn 2 --range 210:330                      \
	  && convert -rotate 90 dcl.ps Tg$${t}.gif;                  \
	done;                                                        \
	convert -delay $(ANIMOPTDELAY) Tg[0-9]*.gif $@               \
	&& rm Tg[0-9]*.gif

Rain_anim-xy.gif: Rain.nc
	anim_str=$(ANIMSTRTIME);                                     \
	anim_end=`expr $(ANIMENDTIME) - $(AVRDELTATIME)`;            \
	timelist=`seq -w $${anim_str} $(AVRDELTATIME) $${anim_end}`; \
	for t in $${timelist}; do                                    \
          t_str=`expr $${t} + 0`;                                    \
	  t_end=`expr $${t} + $(AVRDELTATIME)`;                      \
	  echo time: $${t_str} - $${t_end};                          \
	  $(GPVIEW) Rain.nc@Rain,time=$${t_str}:$${t_end}            \
	    --mean time --wsn 2 --range 0:2000                       \
	  && convert -rotate 90 dcl.ps Rain$${t}.gif;                \
	done;                                                        \
	convert -delay $(ANIMOPTDELAY) Rain[0-9]*.gif $@             \
	&& rm Rain[0-9]*.gif

Temp_anim-xy-sigma0.5.gif: Temp.nc
	anim_str=$(ANIMSTRTIME);                                     \
	anim_end=`expr $(ANIMENDTIME) - $(AVRDELTATIME)`;            \
	timelist=`seq -w $${anim_str} $(AVRDELTATIME) $${anim_end}`; \
	for t in $${timelist}; do                                    \
          t_str=`expr $${t} + 0`;                                    \
	  t_end=`expr $${t} + $(AVRDELTATIME)`;                      \
	  echo time: $${t_str} - $${t_end};                          \
	  $(GPVIEW) Temp.nc@Temp,sig=0.5,time=$${t_str}:$${t_end}    \
	    --mean time --wsn 2 --range 220:285                      \
	  && convert -rotate 90 dcl.ps Temp$${t}.gif;                \
	done;                                                        \
	convert -delay $(ANIMOPTDELAY) Temp[0-9]*.gif $@             \
	&& rm Temp[0-9]*.gif

U_anim-xy-sigma0.2.gif: U.nc
	anim_str=$(ANIMSTRTIME);                                     \
	anim_end=`expr $(ANIMENDTIME) - $(AVRDELTATIME)`;            \
	timelist=`seq -w $${anim_str} $(AVRDELTATIME) $${anim_end}`; \
	for t in $${timelist}; do                                    \
          t_str=`expr $${t} + 0`;                                    \
	  t_end=`expr $${t} + $(AVRDELTATIME)`;                      \
	  echo time: $${t_str} - $${t_end};                          \
	  $(GPVIEW) U.nc@U,sig=0.2,time=$${t_str}:$${t_end}          \
	    --mean time --wsn 2 --range -40:110                      \
	  && convert -rotate 90 dcl.ps U$${t}.gif;                   \
	done;                                                        \
	convert -delay $(ANIMOPTDELAY) U[0-9]*.gif $@                \
	&& rm U[0-9]*.gif

Ps_anim-xy.gif: Ps.nc
	anim_str=$(ANIMSTRTIME);                                     \
	anim_end=`expr $(ANIMENDTIME) - $(AVRDELTATIME)`;            \
	timelist=`seq -w $${anim_str} $(AVRDELTATIME) $${anim_end}`; \
	for t in $${timelist}; do                                    \
          t_str=`expr $${t} + 0`;                                    \
	  t_end=`expr $${t} + $(AVRDELTATIME)`;                      \
	  echo time: $${t_str} - $${t_end};                          \
	  $(GPVIEW) Ps.nc@Ps,time=$${t_str}:$${t_end}                \
	    --mean time --wsn 2 --range 88000:112000                 \
	  && convert -rotate 90 dcl.ps Ps$${t}.gif;                  \
	done;                                                        \
	convert -delay $(ANIMOPTDELAY) Ps[0-9]*.gif $@               \
	&& rm Ps[0-9]*.gif


# 時系列 (全球平均値)
Tg_timeseries-xymean.ps: SurfTemp.nc
	$(GPVIEW) SurfTemp.nc@SurfTemp --mean lon,lat --wsn=2 && mv dcl.ps $@
ps_timeseries-xymean.ps: Ps.nc
	$(GPVIEW) Ps.nc@Ps --mean lon,lat --wsn=2 && mv dcl.ps $@
OLR_timeseries-xymean.ps: OLR.nc
	$(GPVIEW) OLR.nc@OLR --mean lon,lat --wsn=2 && mv dcl.ps $@
SLR_timeseries-xymean.ps: SLR.nc
	$(GPVIEW) SLR.nc@SLR --mean lon,lat --wsn=2 && mv dcl.ps $@
Rain_timeseries-xymean.ps: Rain.nc
	$(GPVIEW) Rain.nc@Rain --mean lon,lat --wsn=2 && mv dcl.ps $@
EVap_timeseries-xymean.ps: EVap.nc
	$(GPVIEW) EVap.nc@EVap --mean lon,lat --wsn=2 && mv dcl.ps $@
Sens_timeseries-xymean.ps: Sens.nc
	$(GPVIEW) Sens.nc@Sens --mean lon,lat --wsn=2 && mv dcl.ps $@
SSR_timeseries-xymean.ps: SSR.nc
	$(GPVIEW) SSR.nc@SSR --mean lon,lat --wsn=2 && mv dcl.ps $@

# 時系列 (経度0-180 度)
Tg-x0-180_timeseries-xymean.ps: SurfTemp.nc
	$(GPVIEW) SurfTemp.nc@SurfTemp,lon=0:180 --mean=lon,lat --wsn=2 && mv dcl.ps $@
ps-x0-180_timeseries-xymean.ps: Ps.nc
	$(GPVIEW) Ps.nc@Ps,lon=0:180 --mean=lon,lat --wsn=2 && mv dcl.ps $@
OLR-x0-180_timeseries-xymean.ps: OLR.nc
	$(GPVIEW) OLR.nc@OLR,lon=0:180 --mean=lon,lat --wsn=2 && mv dcl.ps $@
SLR-x0-180_timeseries-xymean.ps: SLR.nc
	$(GPVIEW) SLR.nc@SLR,lon=0:180 --mean=lon,lat --wsn=2 && mv dcl.ps $@
Rain-x0-180_timeseries-xymean.ps: Rain.nc
	$(GPVIEW) Rain.nc@Rain,lon=0:180 --mean=lon,lat --wsn=2 && mv dcl.ps $@
EVap-x0-180_timeseries-xymean.ps: EVap.nc
	$(GPVIEW) EVap.nc@EVap,lon=0:180 --mean=lon,lat --wsn=2 && mv dcl.ps $@
Sens-x0-180_timeseries-xymean.ps: Sens.nc
	$(GPVIEW) Sens.nc@Sens,lon=0:180 --mean=lon,lat --wsn=2 && mv dcl.ps $@
SSR-x0-180_timeseries-xymean.ps: SSR.nc
	$(GPVIEW) SSR.nc@SSR,lon=0:180 --mean=lon,lat --wsn=2 && mv dcl.ps $@

# 時系列 (経度180-360 度)
Tg-x180-360_timeseries-xymean.ps: SurfTemp.nc
	$(GPVIEW) SurfTemp.nc@SurfTemp,lon=180:360 --mean=lon,lat --wsn=2 && mv dcl.ps $@
ps-x180-360_timeseries-xymean.ps: Ps.nc
	$(GPVIEW) Ps.nc@Ps,lon=180:360 --mean=lon,lat --wsn=2 && mv dcl.ps $@
OLR-x180-360_timeseries-xymean.ps: OLR.nc
	$(GPVIEW) OLR.nc@OLR,lon=180:360 --mean=lon,lat --wsn=2 && mv dcl.ps $@
SLR-x180-360_timeseries-xymean.ps: SLR.nc
	$(GPVIEW) SLR.nc@SLR,lon=180:360 --mean=lon,lat --wsn=2 && mv dcl.ps $@
Rain-x180-360_timeseries-xymean.ps: Rain.nc
	$(GPVIEW) Rain.nc@Rain,lon=180:360 --mean=lon,lat --wsn=2 && mv dcl.ps $@
EVap-x180-360_timeseries-xymean.ps: EVap.nc
	$(GPVIEW) EVap.nc@EVap,lon=180:360 --mean=lon,lat --wsn=2 && mv dcl.ps $@
Sens-x180-360_timeseries-xymean.ps: Sens.nc
	$(GPVIEW) Sens.nc@Sens,lon=180:360 --mean=lon,lat --wsn=2 && mv dcl.ps $@
SSR-x180-360_timeseries-xymean.ps: SSR.nc
	$(GPVIEW) SSR.nc@SSR,lon=180:360 --mean=lon,lat --wsn=2 && mv dcl.ps $@

# 時系列 (格子点)
Tg_timeseries_x0y0.ps: SurfTemp.nc
	$(GPVIEW) SurfTemp.nc@SurfTemp,lon=0,lat=0 --wsn=2 && mv dcl.ps $@
ps_timeseries_x0y0.ps: Ps.nc
	$(GPVIEW) Ps.nc@Ps,lon=0,lat=0     --wsn=2 && mv dcl.ps $@
OLR_timeseries_x0y0.ps: OLR.nc
	$(GPVIEW) OLR.nc@OLR,lon=0,lat=0   --wsn=2 && mv dcl.ps $@
SLR_timeseries_x0y0.ps: SLR.nc
	$(GPVIEW) SLR.nc@SLR,lon=0,lat=0   --wsn=2 && mv dcl.ps $@
Rain_timeseries_x0y0.ps: Rain.nc
	$(GPVIEW) Rain.nc@Rain,lon=0,lat=0 --wsn=2 && mv dcl.ps $@
EVap_timeseries_x0y0.ps: EVap.nc
	$(GPVIEW) EVap.nc@EVap,lon=0,lat=0 --wsn=2 && mv dcl.ps $@
Sens_timeseries_x0y0.ps: Sens.nc
	$(GPVIEW) Sens.nc@Sens,lon=0,lat=0 --wsn=2 && mv dcl.ps $@
SSR_timeseries_x0y0.ps: SSR.nc
	$(GPVIEW) SSR.nc@SSR,lon=0,lat=0   --wsn=2 && mv dcl.ps $@

Tg_timeseries_x90y0.ps: SurfTemp.nc
	$(GPVIEW) SurfTemp.nc@SurfTemp,lon=90,lat=0 --wsn=2 && mv dcl.ps $@
ps_timeseries_x90y0.ps: Ps.nc
	$(GPVIEW) Ps.nc@Ps,lon=90,lat=0     --wsn=2 && mv dcl.ps $@
OLR_timeseries_x90y0.ps: OLR.nc
	$(GPVIEW) OLR.nc@OLR,lon=90,lat=0   --wsn=2 && mv dcl.ps $@
SLR_timeseries_x90y0.ps: SLR.nc
	$(GPVIEW) SLR.nc@SLR,lon=90,lat=0   --wsn=2 && mv dcl.ps $@
Rain_timeseries_x90y0.ps: Rain.nc
	$(GPVIEW) Rain.nc@Rain,lon=90,lat=0 --wsn=2 && mv dcl.ps $@
EVap_timeseries_x90y0.ps: EVap.nc
	$(GPVIEW) EVap.nc@EVap,lon=90,lat=0 --wsn=2 && mv dcl.ps $@
Sens_timeseries_x90y0.ps: Sens.nc
	$(GPVIEW) Sens.nc@Sens,lon=90,lat=0 --wsn=2 && mv dcl.ps $@
SSR_timeseries_x90y0.ps: SSR.nc
	$(GPVIEW) SSR.nc@SSR,lon=90,lat=0   --wsn=2 && mv dcl.ps $@

Tg_timeseries_x180y0.ps: SurfTemp.nc
	$(GPVIEW) SurfTemp.nc@SurfTemp,lon=180,lat=0 --wsn=2 && mv dcl.ps $@
ps_timeseries_x180y0.ps: Ps.nc
	$(GPVIEW) Ps.nc@Ps,lon=180,lat=0     --wsn=2 && mv dcl.ps $@
OLR_timeseries_x180y0.ps: OLR.nc
	$(GPVIEW) OLR.nc@OLR,lon=180,lat=0   --wsn=2 && mv dcl.ps $@
SLR_timeseries_x180y0.ps: SLR.nc
	$(GPVIEW) SLR.nc@SLR,lon=180,lat=0   --wsn=2 && mv dcl.ps $@
Rain_timeseries_x180y0.ps: Rain.nc
	$(GPVIEW) Rain.nc@Rain,lon=180,lat=0 --wsn=2 && mv dcl.ps $@
EVap_timeseries_x180y0.ps: EVap.nc
	$(GPVIEW) EVap.nc@EVap,lon=180,lat=0 --wsn=2 && mv dcl.ps $@
Sens_timeseries_x180y0.ps: Sens.nc
	$(GPVIEW) Sens.nc@Sens,lon=180,lat=0 --wsn=2 && mv dcl.ps $@
SSR_timeseries_x180y0.ps: SSR.nc
	$(GPVIEW) SSR.nc@SSR,lon=180,lat=0   --wsn=2 && mv dcl.ps $@

Tg_timeseries_x270y0.ps: SurfTemp.nc
	$(GPVIEW) SurfTemp.nc@SurfTemp,lon=270,lat=0 --wsn=2 && mv dcl.ps $@
ps_timeseries_x270y0.ps: Ps.nc
	$(GPVIEW) Ps.nc@Ps,lon=270,lat=0     --wsn=2 && mv dcl.ps $@
OLR_timeseries_x270y0.ps: OLR.nc
	$(GPVIEW) OLR.nc@OLR,lon=270,lat=0   --wsn=2 && mv dcl.ps $@
SLR_timeseries_x270y0.ps: SLR.nc
	$(GPVIEW) SLR.nc@SLR,lon=270,lat=0   --wsn=2 && mv dcl.ps $@
Rain_timeseries_x270y0.ps: Rain.nc
	$(GPVIEW) Rain.nc@Rain,lon=270,lat=0 --wsn=2 && mv dcl.ps $@
EVap_timeseries_x270y0.ps: EVap.nc
	$(GPVIEW) EVap.nc@EVap,lon=270,lat=0 --wsn=2 && mv dcl.ps $@
Sens_timeseries_x270y0.ps: Sens.nc
	$(GPVIEW) Sens.nc@Sens,lon=270,lat=0 --wsn=2 && mv dcl.ps $@
SSR_timeseries_x270y0.ps: SSR.nc
	$(GPVIEW) SSR.nc@SSR,lon=270,lat=0   --wsn=2 && mv dcl.ps $@

Tg_timeseries_x90y30.ps: SurfTemp.nc
	$(GPVIEW) SurfTemp.nc@SurfTemp,lon=90,lat=30 --wsn=2 && mv dcl.ps $@
ps_timeseries_x90y30.ps: Ps.nc
	$(GPVIEW) Ps.nc@Ps,lon=90,lat=30     --wsn=2 && mv dcl.ps $@
OLR_timeseries_x90y30.ps: OLR.nc
	$(GPVIEW) OLR.nc@OLR,lon=90,lat=30   --wsn=2 && mv dcl.ps $@
SLR_timeseries_x90y30.ps: SLR.nc
	$(GPVIEW) SLR.nc@SLR,lon=90,lat=30   --wsn=2 && mv dcl.ps $@
Rain_timeseries_x90y30.ps: Rain.nc
	$(GPVIEW) Rain.nc@Rain,lon=90,lat=30 --wsn=2 && mv dcl.ps $@
EVap_timeseries_x90y30.ps: EVap.nc
	$(GPVIEW) EVap.nc@EVap,lon=90,lat=30 --wsn=2 && mv dcl.ps $@
Sens_timeseries_x90y30.ps: Sens.nc
	$(GPVIEW) Sens.nc@Sens,lon=90,lat=30 --wsn=2 && mv dcl.ps $@
SSR_timeseries_x90y30.ps: SSR.nc
	$(GPVIEW) SSR.nc@SSR,lon=90,lat=30   --wsn=2 && mv dcl.ps $@

Tg_timeseries_x270y30.ps: SurfTemp.nc
	$(GPVIEW) SurfTemp.nc@SurfTemp,lon=270,lat=30 --wsn=2 && mv dcl.ps $@
ps_timeseries_x270y30.ps: Ps.nc
	$(GPVIEW) Ps.nc@Ps,lon=270,lat=30     --wsn=2 && mv dcl.ps $@
OLR_timeseries_x270y30.ps: OLR.nc
	$(GPVIEW) OLR.nc@OLR,lon=270,lat=30   --wsn=2 && mv dcl.ps $@
SLR_timeseries_x270y30.ps: SLR.nc
	$(GPVIEW) SLR.nc@SLR,lon=270,lat=30   --wsn=2 && mv dcl.ps $@
Rain_timeseries_x270y30.ps: Rain.nc
	$(GPVIEW) Rain.nc@Rain,lon=270,lat=30 --wsn=2 && mv dcl.ps $@
EVap_timeseries_x270y30.ps: EVap.nc
	$(GPVIEW) EVap.nc@EVap,lon=270,lat=30 --wsn=2 && mv dcl.ps $@
Sens_timeseries_x270y30.ps: Sens.nc
	$(GPVIEW) Sens.nc@Sens,lon=270,lat=30 --wsn=2 && mv dcl.ps $@
SSR_timeseries_x270y30.ps: SSR.nc
	$(GPVIEW) SSR.nc@SSR,lon=270,lat=30   --wsn=2 && mv dcl.ps $@


# 熱フラックス時間変化 全球平均値 重ね描き
flux_timeseries-xymean.ps: OLR.nc SLR.nc Rain.nc EVap.nc Sens.nc
	$(GPVIEW) OLR.nc@OLR --mean=lon,lat SLR.nc@SLR --mean=lon,lat Rain.nc@Rain --mean=lon,lat EVap.nc@EVap --mean=lon,lat Sens.nc@Sens --mean=lon,lat --overplot=5 --range=0:500 --wsn=2 && mv dcl.ps $@

# 熱フラックス時間変化 経度 0-180 度 重ね描き
flux-x0-180_timeseries-xymean.ps:
	$(GPVIEW) OLR.nc@OLR,lon=0:180 --mean=lon,lat SLR.nc@SLR,lon=0:180 --mean=lon,lat Rain.nc@Rain,lon=0:180 --mean=lon,lat EVap.nc@EVap,lon=0:180 --mean=lon,lat Sens.nc@Sens,lon=0:180 --mean=lon,lat --overplot=5 --range=0:500 --wsn=2 && mv dcl.ps $@

# 熱フラックス時間変化 経度 180-360 度 重ね描き
flux-x180-360_timeseries-xymean.ps:
	$(GPVIEW) OLR.nc@OLR,lon=180:360 --mean=lon,lat SLR.nc@SLR,lon=180:360 --mean=lon,lat Rain.nc@Rain,lon=180:360 --mean=lon,lat EVap.nc@EVap,lon=180:360 --mean=lon,lat Sens.nc@Sens,lon=180:360 --mean=lon,lat --overplot=5 --range=0:500 --wsn=2 && mv dcl.ps $@

# 熱フラックス時間変化 1 格子点 重ね描き
flux-x0y0_timeseries.ps:
	$(GPVIEW) OLR.nc@OLR,lon=0,lat=0 SLR.nc@SLR,lon=0,lat=0 Rain.nc@Rain,lon=0,lat=0 EVap.nc@EVap,lon=0,lat=0 Sens.nc@Sens,lon=0,lat=0 --overplot=5 --range=0:500 --wsn=2 && mv dcl.ps $@

flux-x90y0_timeseries.ps:
	$(GPVIEW) OLR.nc@OLR,lon=90,lat=0 SLR.nc@SLR,lon=90,lat=0 Rain.nc@Rain,lon=90,lat=0 EVap.nc@EVap,lon=90,lat=0 Sens.nc@Sens,lon=90,lat=0 --overplot=5 --range=0:500 --wsn=2 && mv dcl.ps $@

flux-x180y0_timeseries.ps:
	$(GPVIEW) OLR.nc@OLR,lon=180,lat=0 SLR.nc@SLR,lon=180,lat=0 Rain.nc@Rain,lon=180,lat=0 EVap.nc@EVap,lon=180,lat=0 Sens.nc@Sens,lon=180,lat=0 --overplot=5 --range=0:500 --wsn=2 && mv dcl.ps $@

flux-x270y0_timeseries.ps:
	$(GPVIEW) OLR.nc@OLR,lon=270,lat=0 SLR.nc@SLR,lon=270,lat=0 Rain.nc@Rain,lon=270,lat=0 EVap.nc@EVap,lon=270,lat=0 Sens.nc@Sens,lon=270,lat=0 --overplot=5 --range=0:500 --wsn=2 && mv dcl.ps $@

flux-x90y30_timeseries.ps:
	$(GPVIEW) OLR.nc@OLR,lon=90,lat=30 SLR.nc@SLR,lon=90,lat=30 Rain.nc@Rain,lon=90,lat=30 EVap.nc@EVap,lon=90,lat=30 Sens.nc@Sens,lon=90,lat=30 --overplot=5 --range=0:500 --wsn=2 && mv dcl.ps $@

flux-x270y30_timeseries.ps:
	$(GPVIEW) OLR.nc@OLR,lon=270,lat=30 SLR.nc@SLR,lon=270,lat=30 Rain.nc@Rain,lon=270,lat=30 EVap.nc@EVap,lon=270,lat=30 Sens.nc@Sens,lon=270,lat=30 --overplot=5 --range=0:500 --wsn=2 && mv dcl.ps $@


# 平面プロット
Tg_xy.ps: SurfTemp.nc
	$(GPVIEW) SurfTemp.nc@SurfTemp,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time --cint=5 --range=170:330 --wsn=2 && mv dcl.ps $@
ps_xy.ps: Ps.nc
	$(GPVIEW) Ps.nc@Ps,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time --cint=500 --range=92000:108000 --wsn=2 && mv dcl.ps $@
Tg_xy_mono.ps: SurfTemp.nc
	$(GPVIEW) SurfTemp.nc@SurfTemp,time=$(AVRSTRTIME):$(AVRENDTIME) --noshade --mean=time --cint=5 --wsn=2 && mv dcl.ps $@
SSR_xy.ps: SSR.nc
	$(GPVIEW) SSR.nc@SSR,time=$(ENDTIME) --wsn=2 && mv dcl.ps $@
SLR_xy.ps:
	$(GPVIEW) SLR.nc@SLR,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time --wsn=2 && mv dcl.ps $@
Rain_xy.ps:
	$(GPVIEW) Rain.nc@Rain,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time --range=0:2000 --cint=300 --wsn=2 && mv dcl.ps $@
OLR_xy.ps:
	$(GPVIEW) OLR.nc@OLR,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time --cint=10 --wsn=2 && mv dcl.ps $@
Sens_xy.ps:
	$(GPVIEW) Sens.nc@Sens,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time --cint=5 --wsn=2 && mv dcl.ps $@
EVap_xy.ps:
	$(GPVIEW) EVap.nc@EVap,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time --wsn=2 && mv dcl.ps $@
u_xy-sigma0.5.ps: U.nc
	$(GPVIEW) U.nc@U,sig=0.5,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time --range=-80:80 --cint=10  --wsn=2 && mv dcl.ps $@
u_xy-sigma0.2.ps: U.nc
	$(GPVIEW) U.nc@U,sig=0.2,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time --range=-80:100 --cint=10  --wsn=2 && mv dcl.ps $@
v_xy-sigma0.5.ps: V.nc
	$(GPVIEW) V.nc@V,sig=0.5,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time --range=-20:20 --cint=2  --wsn=2 && mv dcl.ps $@
v_xy-sigma0.2.ps: V.nc
	$(GPVIEW) V.nc@V,sig=0.2,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time --range=-50:50 --cint=5  --wsn=2 && mv dcl.ps $@
T_xy-z0.ps: Temp.nc
	$(GPVIEW) Temp.nc@Temp,sig=1,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time  --wsn=2 && mv dcl.ps $@
T_xy-sigma0.7.ps: Temp.nc
	$(GPVIEW) Temp.nc@Temp,sig=0.7,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time --range=250:320 --cint=5 --wsn=2 && mv dcl.ps $@

# 赤道断面
u_xz-y0.ps: U.nc
	$(GPVIEW) U.nc@U,lat=0,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time --range=-60:90 --cint=5 --wsn=2 && mv dcl.ps $@
SigDot_xz-y0.ps: SigDot.nc
#	$(GPVIEW) SigDot.nc@SigmaDot,lat=0,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time --range=-5e-6:5e-6 --cint=5e-7 --wsn=2 && mv dcl.ps $@
	$(GPVIEW) SigDot.nc@SigDot,lat=0,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time --range=-5e-6:5e-6 --cint=5e-7 --wsn=2 && mv dcl.ps $@

# 南北断面
u_yz-xmean.ps: U.nc
	$(GPVIEW) U.nc@U,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time,lon --wsn=2 && mv dcl.ps $@
u-x0-180_yz-xmean.ps: U.nc
	$(GPVIEW) U.nc@U,lon=0:180,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time,lon --wsn=2 && mv dcl.ps $@
u-x180-360_yz-xmean.ps: U.nc
	$(GPVIEW) U.nc@U,lon=180:360,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time,lon --wsn=2 && mv dcl.ps $@

u_yz-x270.ps: U.nc
	$(GPVIEW) U.nc@U,lon=270,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time  --wsn=2 && mv dcl.ps $@
u_yz-x90.ps: U.nc
	$(GPVIEW) U.nc@U,lon=90,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time  --wsn=2 && mv dcl.ps $@

v_yz-xmean.ps: V.nc
	$(GPVIEW) V.nc@V,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time,lon --wsn=2 && mv dcl.ps $@
v-x0-180_yz-xmean.ps: V.nc
	$(GPVIEW) V.nc@V,lon=0:180,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time,lon --wsn=2 && mv dcl.ps $@
v-x180-360_yz-xmean.ps: V.nc
	$(GPVIEW) V.nc@V,lon=180:360,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time,lon --wsn=2 && mv dcl.ps $@

v_yz-x270.ps: V.nc
	$(GPVIEW) V.nc@V,lon=270,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time  --wsn=2 && mv dcl.ps $@
v_yz-x90.ps: V.nc
	$(GPVIEW) V.nc@V,lon=90,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time  --wsn=2 && mv dcl.ps $@

T_yz-xmean.ps: Temp.nc
	$(GPVIEW) Temp.nc@Temp,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time,lon --wsn=2 && mv dcl.ps $@
T-x0-180_yz-xmean.ps: Temp.nc
	$(GPVIEW) Temp.nc@Temp,lon=0:180,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time,lon --wsn=2 && mv dcl.ps $@
T-x180-360_yz-xmean.ps: Temp.nc
	$(GPVIEW) Temp.nc@Temp,lon=180:360,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time,lon --wsn=2 && mv dcl.ps $@

T_yz-x270.ps: Temp.nc
	$(GPVIEW) Temp.nc@Temp,lon=270,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time  --wsn=2 && mv dcl.ps $@
T_yz-x90.ps: Temp.nc
	$(GPVIEW) Temp.nc@Temp,lon=90,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time --wsn=2 && mv dcl.ps $@

T_xz-y0.ps: Temp.nc
	$(GPVIEW) Temp.nc@Temp,lat=0,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time --cint=5 --range=220:320 --wsn=2 && mv dcl.ps $@

q_xy-z0.ps: QVap.nc
	$(GPVIEW) QVap.nc@QVap,sig=1,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time --wsn=2 && mv dcl.ps $@

q_yz-xmean.ps: QVap.nc
	$(GPVIEW) QVap.nc@QVap,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time,lon --wsn=2 && mv dcl.ps $@
q-x0-180_yz-xmean.ps: QVap.nc
	$(GPVIEW) QVap.nc@QVap,lon=0:180,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time,lon --wsn=2 && mv dcl.ps $@
q-x180-360_yz-xmean.ps: QVap.nc
	$(GPVIEW) QVap.nc@QVap,lon=180:360,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time,lon --wsn=2 && mv dcl.ps $@

q_yz-x270.ps: QVap.nc
	$(GPVIEW) QVap.nc@QVap,lon=270,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time --wsn=2 && mv dcl.ps $@
q_yz-x90.ps: QVap.nc
	$(GPVIEW) QVap.nc@QVap,lon=90,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time --wsn=2 && mv dcl.ps $@


surface-flux_y0.ps:
	$(GPVIEW) --overplot=5 OLR.nc@OLR,lat=0,time=$(AVRSTRTIME):$(AVRENDTIME) SLR.nc@SLR,lat=0,time=$(AVRSTRTIME):$(AVRENDTIME) Rain.nc@Rain,lat=0,time=$(AVRSTRTIME):$(AVRENDTIME) EVap.nc@EVap,lat=0,time=$(AVRSTRTIME):$(AVRENDTIME) Sens.nc@Sens,lat=0,time=$(AVRSTRTIME):$(AVRENDTIME) --mean=time --index=888 --title='Surface Heat Flux at the Equator' --range=0:500 --wsn=2 && mv dcl.ps $@

surface-flux_y-x0-180.ps:
	$(GPVIEW) --overplot=5 OLR.nc@OLR,lon=0:180,time=$(AVRSTRTIME):$(AVRENDTIME) SLR.nc@SLR,lon=0:180,time=$(AVRSTRTIME):$(AVRENDTIME) Rain.nc@Rain,lon=0:180,time=$(AVRSTRTIME):$(AVRENDTIME) EVap.nc@EVap,lon=0:180,time=$(AVRSTRTIME):$(AVRENDTIME) Sens.nc@Sens,lon=0:180,time=$(AVRSTRTIME):$(AVRENDTIME)  --mean=time,lon --index=888 --title='Surface Heat Flux lon0-180' --range=0:500 --wsn=2 && mv dcl.ps $@

surface-flux_y-x180-360.ps:
	$(GPVIEW) --overplot=5 OLR.nc@OLR,lon=180:360,time=$(AVRSTRTIME):$(AVRENDTIME) SLR.nc@SLR,lon=180:360,time=$(AVRSTRTIME):$(AVRENDTIME) Rain.nc@Rain,lon=180:360,time=$(AVRSTRTIME):$(AVRENDTIME)  EVap.nc@EVap,lon=180:360,time=$(AVRSTRTIME):$(AVRENDTIME) Sens.nc@Sens,lon=180:360,time=$(AVRSTRTIME):$(AVRENDTIME)  --mean=time,lon --index=888 --title='Surface Heat Flux lon180-360' --range=0:500  --wsn=2 && mv dcl.ps $@

surface-flux_y-xmean.ps:
	$(GPVIEW) --overplot=5 OLR.nc@OLR,time=$(AVRSTRTIME):$(AVRENDTIME) SLR.nc@SLR,time=$(AVRSTRTIME):$(AVRENDTIME)  Rain.nc@Rain,time=$(AVRSTRTIME):$(AVRENDTIME) EVap.nc@EVap,time=$(AVRSTRTIME):$(AVRENDTIME) Sens.nc@Sens,time=$(AVRSTRTIME):$(AVRENDTIME)  --mean=time,lon --index=888 --title='Surface Heat Flux zonal mean' --range=0:500 --wsn=2 && mv dcl.ps $@

# 帯状平均の時系列
ps_timeseries-xmean.ps: Ps.nc
	$(GPVIEW) Ps.nc@Ps,time=$(STRTIME):$(ENDTIME):5 --mean lon --wsn 2 --nocont --ex --range 90000:110000 && mv dcl.ps $@

Tg_timeseries-xmean.ps: SurfTemp.nc
	$(GPVIEW) SurfTemp.nc@SurfTemp,time=$(STRTIME):$(ENDTIME):5 --mean lon --wsn 2 --nocont --ex && mv dcl.ps $@

u_timeseries-xmean-sigma0.2.ps: U.nc
	$(GPVIEW) U.nc@U,sig=0.2,time=$(STRTIME):$(ENDTIME):5 --mean lon --wsn 2 --nocont --exch && mv dcl.ps $@

Rain_timeseries-xmean.ps: Rain.nc
	$(GPVIEW) Rain.nc@Rain,time=$(STRTIME):$(ENDTIME):5 --mean lon --wsn 2 --nocont --exch --range 0:1000 && mv dcl.ps $@


# その他
temp_eddy_variance.ps: temp_eddy_variance.rb Temp.nc
	ruby temp_eddy_variance.rb Temp.nc Temp && mv dcl.ps $@

u_eddy_var_wave-zmean.ps: u_eddy_var_wave-zmean.rb
	ruby u_eddy_var_wave-zmean.rb U.nc U && mv dcl.ps $@



# parameters
path_src_tree.txt: ../../arch_data/src_backup
	echo -n '  * ((<ソースツリー|URL:../arch_data/src_backup/>))' > tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
path_src_tgz.txt: ../../arch_data/src_backup.tar.gz
	echo -n '  * ((<[tar.gz]|URL:../arch_data/src_backup.tar.gz>))' > tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
path_namelist.txt: ../../arch_data
	echo -n '* ((<NAMELIST ファイル|URL:../arch_data/' > tmpfile \
	&& basename `ls -1 ../../arch_data/dcpam_*.nml` >> tmpfile \
	&& echo '>))' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
path_log.txt: ../../arch_data
	echo -n '* ((<標準出力|URL:../arch_data/' > tmpfile \
	&& basename `ls -1 ../../arch_data/dcpam_*.log` >> tmpfile \
	&& echo '>))' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile

FilterParameter.txt: ../../arch_data ez_param_extract_from_log.sh
	echo -n '      フィルター係数は ' > tmpfile \
	&& ./ez_param_extract_from_log.sh $@ ../../arch_data/*.log >> tmpfile \
	&& echo '.' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
imax.txt: ../../arch_data ez_param_extract_from_log.sh 
	echo -n '  * 経度格子点数: ' > tmpfile \
	&& ./ez_param_extract_from_log.sh $@ ../../arch_data/*.log >> tmpfile \
	&& imax=`./ez_param_extract_from_log.sh $@ ../../arch_data/*.log` \
	&& echo ' (' >> tmpfile \
	&& echo "scale=1; 360/`./ez_param_extract_from_log.sh $@ ../../arch_data/*.log`" | bc >> tmpfile \
	&& echo '度)' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
jmax.txt: ../../arch_data ez_param_extract_from_log.sh
	echo -n '  * 緯度格子点数: ' > tmpfile \
	&& ./ez_param_extract_from_log.sh $@ ../../arch_data/*.log >> tmpfile \
	&& jmax=`./ez_param_extract_from_log.sh $@ ../../arch_data/*.log` \
	&& echo ' (' >> tmpfile \
	&& echo "scale=1; 180/`./ez_param_extract_from_log.sh $@ ../../arch_data/*.log`" | bc >> tmpfile \
	&& echo '度)' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
kmax.txt: ../../arch_data ez_param_extract_from_log.sh
	echo -n '  * 鉛直層数: ' > tmpfile \
	&& ./ez_param_extract_from_log.sh $@ ../../arch_data/*.log >> tmpfile \
	&& echo '' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
nmax.txt: ../../arch_data ez_param_extract_from_log.sh
	echo -n '* 最大波数: ' > tmpfile \
	&& ./ez_param_extract_from_log.sh $@ ../../arch_data/*.log >> tmpfile \
	&& echo '' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
EndTime.txt: ../../arch_data ez_param_extract_from_log.sh
	echo -n '* 積分時間: ' > tmpfile \
	&& ./ez_param_extract_from_log.sh $@ ../../arch_data/*.log >> tmpfile \
	&& echo ' [day]' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
DelTime.txt: ../../arch_data ez_param_extract_from_log.sh
	echo -n '* 時間間隔: ' > tmpfile \
	&& ./ez_param_extract_from_log.sh $@ ../../arch_data/*.log >> tmpfile \
	&& echo ' 分' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
RPlanet.txt: ../../arch_data ez_param_extract_from_log.sh
	echo -n '* 惑星半径: ' > tmpfile \
	&& ./ez_param_extract_from_log.sh $@ ../../arch_data/*.log >> tmpfile \
	&& echo ' [m]' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
Omega.txt: ../../arch_data ez_param_extract_from_log.sh
	echo -n '* 自転角速度: ' > tmpfile \
	&& ./ez_param_extract_from_log.sh $@ ../../arch_data/*.log >> tmpfile \
	&& echo ' [s-1]' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
Grav.txt: ../../arch_data ez_param_extract_from_log.sh
	echo -n '* 重力加速度: ' > tmpfile \
	&& ./ez_param_extract_from_log.sh $@ ../../arch_data/*.log >> tmpfile \
	&& echo ' [m s-2]' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
CpDry.txt: ../../arch_data ez_param_extract_from_log.sh
	echo -n '* 乾燥大気定圧比熱: ' > tmpfile \
	&& ./ez_param_extract_from_log.sh $@ ../../arch_data/*.log >> tmpfile \
	&& echo ' [J K-1 kg-1]' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
MolWtDry.txt: ../../arch_data ez_param_extract_from_log.sh
	echo -n '* 乾燥大気平均分子量: ' > tmpfile \
	&& ./ez_param_extract_from_log.sh $@ ../../arch_data/*.log >> tmpfile \
	&& echo ' [kg mol-1]' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
GasRDry.txt: ../../arch_data ez_param_extract_from_log.sh
	echo -n '* 乾燥大気気体定数: ' > tmpfile \
	&& ./ez_param_extract_from_log.sh $@ ../../arch_data/*.log >> tmpfile \
	&& echo ' [J K-1 kg-1]' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
EFoldTime.txt: ../../arch_data ez_param_extract_from_log.sh
	echo -n '* 最大波数に対する e-folding time: ' > tmpfile \
	&& ./ez_param_extract_from_log.sh $@ ../../arch_data/*.log >> tmpfile \
	&& echo ' [hours]' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
VisOrder.txt: ../../arch_data ez_param_extract_from_log.sh
	echo -n '* 超粘性の次数: ' > tmpfile \
	&& ./ez_param_extract_from_log.sh $@ ../../arch_data/*.log >> tmpfile \
	&& echo '' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
SolarConst.txt: ../../arch_data ez_param_extract_from_log.sh
	echo -n '    * 太陽定数: ' > tmpfile \
	&& ./ez_param_extract_from_log.sh $@ ../../arch_data/*.log >> tmpfile \
	&& echo ' [W m-2]' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
AtmosAlbedo.txt: ../../arch_data ez_param_extract_from_log.sh
	echo -n '    * 大気アルベド: ' > tmpfile \
	&& ./ez_param_extract_from_log.sh $@ ../../arch_data/*.log >> tmpfile \
	&& echo '' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
ShortBandNum.txt: ../../arch_data ez_param_extract_from_log.sh
	echo -n '    * 短波放射: ' > tmpfile \
	&& ./ez_param_extract_from_log.sh $@ ../../arch_data/*.log >> tmpfile \
	&& echo ' 色' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
LongBandNum.txt: ../../arch_data ez_param_extract_from_log.sh
	echo -n '    * 長波放射: ' > tmpfile \
	&& ./ez_param_extract_from_log.sh $@ ../../arch_data/*.log >> tmpfile \
	&& echo ' 色' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
Albedo.txt: ../../arch_data ez_param_extract_from_log.sh
	echo -n '    * アルベド: ' > tmpfile \
	&& ./ez_param_extract_from_log.sh $@ ../../arch_data/*.log >> tmpfile \
	&& echo '' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile
HeatCapacity.txt: ../../arch_data ez_param_extract_from_log.sh
	echo -n '    * 熱容量: ' > tmpfile \
	&& ./ez_param_extract_from_log.sh $@ ../../arch_data/*.log >> tmpfile \
	&& echo '' >> tmpfile \
	&& cat tmpfile | tr -d '\012' > $@ \
	&& rm tmpfile







clean.all: clean clean.fig

clean.fig:
	-rm -f *.ps *.eps *.png

clean:
	-rm -f *~
