#Dieses Script trennt Meeresregionen von Landregionen; Meeresregionen werden als "No Data"
#gekennzeichnet und fallen so bei der Landklassifikation nicht mehr ins Gewicht.
#Speziell für NOAA AVHRR DATEN 11-14, nicht NOAA 16, da dieser avhrr_ch3a hat!! 14.05.2002. 
#! /bin/csh -fx


cd /usr/data/maroc/masked_files/
foreach file (noaa-14_utm.wgs84_????.*)
cd /usr/data/maroc/masked_files/
emathp \
	expr_vars='avhrr_ch1 avhrr_ch2 avhrr_ch3 avhrr_ch4 avhrr_ch5 latitude longitude' \
	num_exprs='12' \
	y1_expr='x1' \
	y2_expr='x2' \
	y3_expr='x3' \
	y4_expr='x4' \
	y5_expr='x5' \
	y6_expr='x6' \
	y7_expr='x7' \
	y8_expr='land ? y1 : badval' \
	y9_expr='land ? y2 : badval' \
	y10_expr='land ? y3 : badval' \
	y11_expr='land ? y4 : badval' \
	y12_expr='land ? y5 : badval' \
	use_degrees='n' \
	poly_size='100' \
	save_exprs='8,9,10,11,12' \
	var_names='avhrr_ch1 avhrr_ch2 avhrr_ch3 avhrr_ch4 avhrr_ch5' \
	var_units='albedo*100% albedo*100% temp_deg_c temp_deg_c temp_deg_c' \
	var_types='sh sh sh sh sh' \
	scale_offsets='0 0 0 0 0' \
	scale_factors='1 1 1 1 1' \
$file $file.land

mv -f $file.land /usr/data/maroc/masked_files/ohne_meer/
cd /usr/data/maroc/masked_files/ohne_meer/
mv -f $file.land $file

end