Added test case for perturbation data.

.gitignore

@@ -12,3 +12,6 @@ src/sampling
 src/spec2tab
 src/variance
 src/volcano
+tests/data/SNDR*nc
+tests/data/cris_l1b_full.zip
+tests/pert_test/data/

src/Makefile

@@ -5,6 +5,9 @@
 # Executables...
 EXC = day2doy doy2day map_pert map_rad noise_pert perturbation sampling spec2tab variance volcano
 
+# Tests...
+TESTS = pert_test
+
 # Installation directory...
 DESTDIR ?= ../bin
 
@@ -26,6 +29,9 @@ INFO ?= 0
 # Compile for profiling...
 PROF ?= 0
 
+# Compile for coverage testing...
+COV ?= 0
+
 # -----------------------------------------------------------------------------
 # Set flags for GNU compiler...
 # -----------------------------------------------------------------------------
@@ -53,6 +59,11 @@ ifeq ($(PROF),1)
   CFLAGS += -pg
 endif
 
+# Compile for coverage testing...
+ifeq ($(COV),1)
+  CFLAGS += --coverage
+endif
+
 # Static compilation...
 ifeq ($(STATIC),1)
   CFLAGS += -static
@@ -76,6 +87,14 @@ libcris.o: libcris.c libcris.h Makefile
 jurassic.o: jurassic.c jurassic.h Makefile
 	$(CC) $(CFLAGS) -c -o jurassic.o jurassic.c
 
+check: $(TESTS)
+
+$(TESTS): all
+	@(echo "\n===== Running \"$@\" ... =====") ; \
+	  cd ../tests/$@ ; ./run.sh \
+	  && (echo "\n===== Test \"$@\" passed! =====") \
+	  || (echo "\n===== Test \"$@\" failed! =====" ; exit 1)
+
 clean:
 	rm -rf $(EXC) *.o *.gcda *.gcno *.gcov coverage* *~
 
@@ -89,7 +108,7 @@ cppcheck:
 
 dist:
 	cd ..  && tar cvfz cris_$(VERSION).tgz bin/* \
-          docs/Doxyfile example/* libs/build.sh libs/*gz src/*
+          docs/Doxyfile example/* libs/build.sh libs/*gz src/* tests/*
 
 doxygen:
 	cd ../docs && doxygen

tests/pert_test/data.ref/noise_15mu_high.tab

@@ -0,0 +1,6 @@
+# $1 = longitude [deg]
+# $2 = latitude [deg]
+# $3 = mean brightness temperature [K]
+# $4 = noise estimate [K]
+
+6.83922 -56.997 262.576 0.294744

tests/pert_test/data.ref/noise_15mu_low.tab

@@ -0,0 +1,6 @@
+# $1 = longitude [deg]
+# $2 = latitude [deg]
+# $3 = mean brightness temperature [K]
+# $4 = noise estimate [K]
+
+6.83922 -56.997 239.119 0.0764632

tests/pert_test/data.ref/noise_4mu.tab

@@ -0,0 +1,6 @@
+# $1 = longitude [deg]
+# $2 = latitude [deg]
+# $3 = mean brightness temperature [K]
+# $4 = noise estimate [K]
+
+6.83922 -56.997 252.844 0.0665062

tests/pert_test/run.sh

@@ -0,0 +1,43 @@
+#! /bin/bash
+
+# Set environment...
+export LD_LIBRARY_PATH=../../libs/build/lib:$LD_LIBRARY_PATH
+export OMP_NUM_THREADS=4
+
+# Setup...
+cris=../../src
+
+# Uncompress test data...
+cd ../data \
+    && zip -s 0 cris_l1b.zip --out cris_l1b_full.zip \
+    && unzip -o cris_l1b_full.zip \
+    && cd - || exit
+
+# Create directory...
+rm -rf data && mkdir -p data || exit
+
+# Create perturbation file...
+$cris/perturbation - data/pert.nc \
+		   ../data/SNDR.SNPP.CRIS.20220115T0000.m06.g001.L1B.std.v03_08.G.220115064421.nc
+
+# Loop over channel sets...
+for pert in 4mu 15mu_low 15mu_high ; do
+
+    # Get map data...
+    $cris/map_pert - data/pert.nc data/map_$pert.tab PERTNAME $pert
+
+    # Estimate noise...
+    $cris/noise_pert - data/pert.nc data/noise_$pert.tab PERTNAME $pert
+
+    # Get variance...
+    $cris/variance - data/var_$pert.tab data/pert.nc PERTNAME $pert NX 60 NY 30
+
+done
+
+# Compare files...
+echo -e "\nCompare results..."
+error=0
+for f in $(ls data.ref/*.nc data.ref/*.tab) ; do
+    diff -q -s data/"$(basename "$f")" "$f" || error=1
+done
+exit $error