#! /bin/csh -fx

# Input consists of data in two directories: maskdir and datadir.
# Script cuts data out of files in datadir and places in outputdir
# using coresseponding avcloud (mask) file.

# Important: Files in maskdir must have same name as corresponding file in
# datadir

set maskdir=/usr/data/maroc/maskdir
set datadir=/usr/data/maroc/datadir
set outputdir=/usr/data/maroc/masked_files
set maskvar=cloud
if ( ! -d $outputdir ) mkdir $outputdir \
emaskp \
	file1_vars='cloud' \
	file2_vars='' \
	num_exprs='1' \
	y1_expr='x1==0 ? 1: badval' \
	mask_vars='avhrr_ch1 avhrr_ch2 avhrr_ch3a avhrr_ch4 avhrr_ch5' \
	constant_files='1' \
	maskfile \
	datafile \
	masked_files/
foreach datafile ( $datadir/* )
echo "Masking $datafile"
emaskp \
	file1_vars=$maskvar \
	file2_vars='' \
	num_exprs='1' \
	y1_expr='x1==0 ? 1: badval' \
	mask_vars='avhrr_ch1 avhrr_ch2 avhrr_ch3a avhrr_ch4 avhrr_ch5' \
	constant_files='1' \
	maskdir/$datafile:t \
	$datafile \
	$outputdir/
end
	
	

# Notes:	
# emaskp identifies variable to mask with
# x1==0 ? 1: badval: If cloud is zero then extract data variables to mask
# maskfile= File with mask variable
# datafile= File with variables to mask






	
# Insert code here for computing NDVI on files in $outputdir
# Insert code here for compositing NDVI files into single dataset