Add Pack.py.test.in

configure.ac

@@ -247,6 +247,7 @@ AC_CONFIG_FILES([test/pytest/cdoTest.py
                 test/pytest/Monarith.py.test
                 test/pytest/Multiyearstat.py.test
                 test/pytest/Ninfo.py.test
+                test/pytest/Pack.py.test
                 test/pytest/Percentile.py.test
                 test/pytest/Read_grib.py.test
                 test/pytest/Read_netcdf.py.test

test/data/Makefile.am

@@ -4,7 +4,7 @@ INPUTDATA = ts_1d_5years ts_ym_5years ts_mm_5years ts_mm_1year ts_mm_1991 ts_6h_
             grib_testfile01.grb grib_testfile02.grb grib_testfile03.grb netcdf_testfile01.nc netcdf_testfile02.nc netcdf_testfile03.nc testfile01c.nc \
             datar.nc datac.nc datau.nc datag.nc arith1.srv expr1.srv arithmask.srv psl_DJF_anom.grb tsurf_spain.grb spain.grid \
             topo_eu5.grb vars_data.grb math_data tsurf_1d_1year tsurf_runpctl_1d_1year mpiom_tho_sao.srv topo5.srv \
-            temp_hp32_nest.nc temp_hp16_nest.nc temp_global_2.nc
+            temp_hp32_nest.nc temp_hp16_nest.nc temp_global_2.nc tsurf_5steps.nc tsurf_5steps_land.nc
 
 FILE         = file_F32_srv_ref
 GRIB         = grib_testfile01_sinfo_ref grib_testfile01_info_ref grib_testfile02_sinfo_ref grib_testfile02_info_ref grib_testfile03_sinfo_ref grib_testfile03_info_ref
@@ -180,6 +180,7 @@ VERTFILLMISS = vertfillmiss_1.srv vertfillmiss_2.srv vertfillmiss_3.srv vertfill
                vertfillmiss_forward_1_ref vertfillmiss_forward_2_ref vertfillmiss_forward_3_ref vertfillmiss_forward_4_ref \
                vertfillmiss_linear_1_ref vertfillmiss_linear_2_ref vertfillmiss_linear_3_ref vertfillmiss_linear_4_ref \
                vertfillmiss_nearest_1_ref vertfillmiss_nearest_2_ref vertfillmiss_nearest_3_ref vertfillmiss_nearest_4_ref
+PACK         = pack_i16_ref pack_i8_ref pack_u16_ref pack_u8_ref
 
 EXTRA_DIST = $(INPUTDATA) $(FILE) $(GRIB) $(NETCDF) $(ADISIT) $(EOF) $(COMP) $(YDRUNSTAT) $(YDAYSTAT) $(YDAYSTATM) \
              $(YMONSTAT) $(YMONSTATM) $(YSEASSTAT) $(YSEASSTATM) $(SEASSTAT) $(SEASSTATM) $(TIMPCTL) $(TIMSELPCTL) \
@@ -190,4 +191,5 @@ EXTRA_DIST = $(INPUTDATA) $(FILE) $(GRIB) $(NETCDF) $(ADISIT) $(EOF) $(COMP) $(Y
              $(MERSTAT) $(ZONSTAT) $(ENSSTAT) $(ENSSTATM) $(ENSPCTL) $(SPECTRAL) $(WIND) $(INTTIME) $(VERTINT) \
              $(REMAPHEALPIX) $(REMAPGRID) $(REMAPSTAT) $(REMAP) $(REMAP2) $(REMAP3) $(REMAP4) $(SELECT) $(DETREND) \
              $(THREAD) $(EXPR) $(GRADSDES) $(ARITH) $(DAYARITH) $(MONARITH) $(YEARARITH) $(YDAYARITH) $(YMONARITH) \
-             $(MAPREDUCE) $(MERGETIME) $(REMAPETA) $(SMOOTH) $(SETMISS) $(FILTER) $(PERCENTILE) $(TIMFILLMISS) $(VERTFILLMISS)
+             $(MAPREDUCE) $(MERGETIME) $(REMAPETA) $(SMOOTH) $(SETMISS) $(FILTER) $(PERCENTILE) $(TIMFILLMISS) $(VERTFILLMISS) \
+             $(PACK)

test/data/gen_refdata.sh

@@ -6,6 +6,19 @@ FORMAT="-f srv -b F32"
 #
 ########################################################################
 #
+# PACK
+#
+set -x
+IFILE=tsurf_5steps_land.nc
+OPERATOR=pack
+for NBITS in i16 i8 u16 u8; do
+    $CDO -f nc4 -b $NBITS $OPERATOR $IFILE ${OPERATOR}_${NBITS}_ref
+    $CDO -f nc unpack ${OPERATOR}_${NBITS}_ref unpack_${NBITS}_ref
+done
+exit
+#
+########################################################################
+#
 # Timfill
 #
 IFILE=timfilldata

test/pytest/Pack.py.test.in

+#! @PYTHON@
+
+from cdoTest import *
+
+HAS_NETCDF=cdo_check_req("has-nc")
+
+FORMAT="-f nc4"
+
+ABSLIMMAX=0.001
+
+IFILE=f'{DATAPATH}/tsurf_5steps_land.nc'
+
+OPERATOR="pack"
+
+test_module = TestModule()
+
+for NBITS in ["i16", "i8", "u16", "u8"]:
+    RFILE=f'{DATAPATH}{OPERATOR}_{NBITS}_ref'
+    OFILE=f'{OPERATOR}_{NBITS}_res'
+    if (HAS_NETCDF):
+        t=TAPTest(f'{OPERATOR} {NBITS}')
+        t.add(f'{CDO}  {FORMAT} -b {NBITS} {OPERATOR} {IFILE} {OFILE}')
+        t.add(f'{CDO}  diff,abslim={ABSLIMMAX} {RFILE} {OFILE}')
+        t.clean(OFILE)     
+        test_module.add(t)
+    else:
+        test_module.add_skip("NetCDF not enabled")
+
+OPERATOR="unpack"
+
+for NBITS in ["i16", "i8", "u16", "u8"]:
+    RFILE=f'{DATAPATH}pack_{NBITS}_ref'
+    OFILE=f'{OPERATOR}_{NBITS}_res'
+    if (HAS_NETCDF):
+        t=TAPTest(f'{OPERATOR} {NBITS}')
+        t.add(f'{CDO}  -f nc {OPERATOR} {RFILE} {OFILE}')
+        t.add(f'{CDO}  diff,abslim={ABSLIMMAX} {RFILE} {OFILE}')
+        t.clean(OFILE)     
+        test_module.add(t)
+    else:
+        test_module.add_skip("NetCDF not enabled")
+
+test_module.run()