diff --git a/pyutils/src/icon4py/f2ser/parse.py b/pyutils/src/icon4py/f2ser/parse.py
index 609859485e..09755f92a4 100644
--- a/pyutils/src/icon4py/f2ser/parse.py
+++ b/pyutils/src/icon4py/f2ser/parse.py
@@ -159,11 +159,12 @@ def _parse_derived_types(self, derived_types: dict) -> dict:
             MissingDerivedTypeError: If the type definition for a derived type could not be found in any of the dependency files.
         """
         derived_type_defs = {}
-        for dep in self.dependencies:
-            parsed = crack(dep)
-            for block in parsed["body"]:
-                if block["block"] == "type":
-                    derived_type_defs[block["name"]] = block["vars"]
+        if self.dependencies:
+            for dep in self.dependencies:
+                parsed = crack(dep)
+                for block in parsed["body"]:
+                    if block["block"] == "type":
+                        derived_type_defs[block["name"]] = block["vars"]
 
         for _, subroutine_vars in derived_types.items():
             for _, intent_vars in subroutine_vars.items():
diff --git a/pyutils/tests/f2ser/test_f2ser_cli.py b/pyutils/tests/f2ser/test_f2ser_cli.py
index ac544211ea..b32d87e1d9 100644
--- a/pyutils/tests/f2ser/test_f2ser_cli.py
+++ b/pyutils/tests/f2ser/test_f2ser_cli.py
@@ -43,6 +43,15 @@ def test_cli_no_deps(no_deps_source_file, outfile, cli):
     assert result.exit_code == 0
 
 
+def test_cli_wrong_deps(diffusion_granule, samples_path, outfile, cli):
+    inp = str(diffusion_granule)
+    deps = [str(samples_path / "wrong_derived_types_example.f90")]
+    args = [inp, outfile, "-d", *deps]
+    result = cli.invoke(main, args)
+    assert result.exit_code == 2
+    assert "Invalid value for '--dependencies' / '-d'" in result.output
+
+
 def test_cli_missing_deps(diffusion_granule, outfile, cli):
     inp = str(diffusion_granule)
     args = [inp, outfile]
@@ -50,10 +59,16 @@ def test_cli_missing_deps(diffusion_granule, outfile, cli):
     assert isinstance(result.exception, MissingDerivedTypeError)
 
 
+def test_cli_wrong_source(outfile, cli):
+    inp = str("foo.90")
+    args = [inp, outfile]
+    result = cli.invoke(main, args)
+    assert "Invalid value for 'GRANULE_PATH'" in result.output
+
+
 def test_cli_missing_source(not_existing_diffusion_granule, outfile, cli):
     inp = str(not_existing_diffusion_granule)
     args = [inp, outfile]
     result = cli.invoke(main, args)
-    error_search = result.stdout.find("Invalid value for 'GRANULE_PATH'")
-    assert error_search != -1
     assert isinstance(result.exception, SystemExit)
+    assert "Invalid value for 'GRANULE_PATH'" in result.output
diff --git a/pyutils/tests/f2ser/test_f2ser_codegen.py b/pyutils/tests/f2ser/test_f2ser_codegen.py
index 9f5368373b..b6d0c38ed3 100644
--- a/pyutils/tests/f2ser/test_f2ser_codegen.py
+++ b/pyutils/tests/f2ser/test_f2ser_codegen.py
@@ -11,18 +11,100 @@
 #
 # SPDX-License-Identifier: GPL-3.0-or-later
 
+import pytest
+
 from icon4py.f2ser.deserialise import ParsedGranuleDeserialiser
 from icon4py.f2ser.parse import GranuleParser
 from icon4py.liskov.codegen.serialisation.generate import (
     SerialisationCodeGenerator,
 )
+from icon4py.liskov.codegen.shared.types import GeneratedCode
 
 
 def test_deserialiser_diffusion_codegen(diffusion_granule, diffusion_granule_deps):
-    parser = GranuleParser(diffusion_granule, diffusion_granule_deps)
-    parsed = parser()
-    deserialiser = ParsedGranuleDeserialiser(parsed, directory=".", prefix="test")
-    interface = deserialiser()
-    generator = SerialisationCodeGenerator(interface)
-    generated = generator()
+    parsed = GranuleParser(diffusion_granule, diffusion_granule_deps)()
+    interface = ParsedGranuleDeserialiser(parsed, directory=".", prefix="test")()
+    generated = SerialisationCodeGenerator(interface)()
     assert len(generated) == 3
+
+
+@pytest.fixture
+def expected_no_deps_serialization_directives():
+    serialization_directives = [
+        GeneratedCode(
+            startln=12,
+            source="\n"
+            '    !$ser init directory="." prefix="test"\n'
+            "\n"
+            "    !$ser savepoint no_deps_init_in\n"
+            "\n"
+            "    PRINT *, 'Serializing a=a'\n"
+            "\n"
+            "    !$ser data a=a\n"
+            "\n"
+            "    PRINT *, 'Serializing b=b'\n"
+            "\n"
+            "    !$ser data b=b",
+        ),
+        GeneratedCode(
+            startln=14,
+            source="\n"
+            "    !$ser savepoint no_deps_init_out\n"
+            "\n"
+            "    PRINT *, 'Serializing c=c'\n"
+            "\n"
+            "    !$ser data c=c\n"
+            "\n"
+            "    PRINT *, 'Serializing b=b'\n"
+            "\n"
+            "    !$ser data b=b",
+        ),
+        GeneratedCode(
+            startln=20,
+            source="\n"
+            "    !$ser savepoint no_deps_run_in\n"
+            "\n"
+            "    PRINT *, 'Serializing a=a'\n"
+            "\n"
+            "    !$ser data a=a\n"
+            "\n"
+            "    PRINT *, 'Serializing b=b'\n"
+            "\n"
+            "    !$ser data b=b",
+        ),
+        GeneratedCode(
+            startln=22,
+            source="\n"
+            "    !$ser savepoint no_deps_run_out\n"
+            "\n"
+            "    PRINT *, 'Serializing c=c'\n"
+            "\n"
+            "    !$ser data c=c\n"
+            "\n"
+            "    PRINT *, 'Serializing b=b'\n"
+            "\n"
+            "    !$ser data b=b",
+        ),
+    ]
+    return serialization_directives
+
+
+def test_deserialiser_directives_no_deps_codegen(
+    no_deps_source_file, expected_no_deps_serialization_directives
+):
+    parsed = GranuleParser(no_deps_source_file)()
+    interface = ParsedGranuleDeserialiser(parsed, directory=".", prefix="test")()
+    generated = SerialisationCodeGenerator(interface)()
+    assert generated == expected_no_deps_serialization_directives
+
+
+def test_deserialiser_directives_diffusion_codegen(
+    diffusion_granule, diffusion_granule_deps, samples_path
+):
+    parsed = GranuleParser(diffusion_granule, diffusion_granule_deps)()
+    interface = ParsedGranuleDeserialiser(parsed, directory=".", prefix="test")()
+    generated = SerialisationCodeGenerator(interface)()
+    reference_savepoint = (
+        samples_path / "expected_diffusion_granule_savepoint.f90"
+    ).read_text()
+    assert generated[0].source == reference_savepoint.rstrip()
diff --git a/testutils/src/icon4py/testutils/fortran/expected_diffusion_granule_savepoint.f90 b/testutils/src/icon4py/testutils/fortran/expected_diffusion_granule_savepoint.f90
new file mode 100644
index 0000000000..19a3653d97
--- /dev/null
+++ b/testutils/src/icon4py/testutils/fortran/expected_diffusion_granule_savepoint.f90
@@ -0,0 +1,656 @@
+
+    !$ser verbatim real, dimension(:,:,:), allocatable :: edges_primal_normal_vert_v1
+
+    !$ser verbatim real, dimension(:,:,:), allocatable :: edges_primal_normal_vert_v2
+
+    !$ser verbatim real, dimension(:,:,:), allocatable :: edges_dual_normal_vert_v1
+
+    !$ser verbatim real, dimension(:,:,:), allocatable :: edges_dual_normal_vert_v2
+
+    !$ser verbatim real, dimension(:,:,:), allocatable :: edges_primal_normal_cell_v1
+
+    !$ser verbatim real, dimension(:,:,:), allocatable :: edges_primal_normal_cell_v2
+
+    !$ser verbatim real, dimension(:,:,:), allocatable :: edges_dual_normal_cell_v1
+
+    !$ser verbatim real, dimension(:,:,:), allocatable :: edges_dual_normal_cell_v2
+
+    !$ser init directory="." prefix="test"
+
+    !$ser savepoint diffusion_init_in
+
+    !$ser verbatim allocate(edges_primal_normal_vert_v1(size(edges_primal_normal_vert, 1),size(edges_primal_normal_vert, 2),size(edges_primal_normal_vert, 3)))
+    !$ser data edges_primal_normal_vert_v1=edges_primal_normal_vert_v1(:,:,:)
+    !$ser verbatim deallocate(edges_primal_normal_vert_v1)
+
+    !$ser verbatim allocate(edges_primal_normal_vert_v2(size(edges_primal_normal_vert, 1),size(edges_primal_normal_vert, 2),size(edges_primal_normal_vert, 3)))
+    !$ser data edges_primal_normal_vert_v2=edges_primal_normal_vert_v2(:,:,:)
+    !$ser verbatim deallocate(edges_primal_normal_vert_v2)
+
+    !$ser verbatim allocate(edges_dual_normal_vert_v1(size(edges_dual_normal_vert, 1),size(edges_dual_normal_vert, 2),size(edges_dual_normal_vert, 3)))
+    !$ser data edges_dual_normal_vert_v1=edges_dual_normal_vert_v1(:,:,:)
+    !$ser verbatim deallocate(edges_dual_normal_vert_v1)
+
+    !$ser verbatim allocate(edges_dual_normal_vert_v2(size(edges_dual_normal_vert, 1),size(edges_dual_normal_vert, 2),size(edges_dual_normal_vert, 3)))
+    !$ser data edges_dual_normal_vert_v2=edges_dual_normal_vert_v2(:,:,:)
+    !$ser verbatim deallocate(edges_dual_normal_vert_v2)
+
+    !$ser verbatim allocate(edges_primal_normal_cell_v1(size(edges_primal_normal_cell, 1),size(edges_primal_normal_cell, 2),size(edges_primal_normal_cell, 3)))
+    !$ser data edges_primal_normal_cell_v1=edges_primal_normal_cell_v1(:,:,:)
+    !$ser verbatim deallocate(edges_primal_normal_cell_v1)
+
+    !$ser verbatim allocate(edges_primal_normal_cell_v2(size(edges_primal_normal_cell, 1),size(edges_primal_normal_cell, 2),size(edges_primal_normal_cell, 3)))
+    !$ser data edges_primal_normal_cell_v2=edges_primal_normal_cell_v2(:,:,:)
+    !$ser verbatim deallocate(edges_primal_normal_cell_v2)
+
+    !$ser verbatim allocate(edges_dual_normal_cell_v1(size(edges_dual_normal_cell, 1),size(edges_dual_normal_cell, 2),size(edges_dual_normal_cell, 3)))
+    !$ser data edges_dual_normal_cell_v1=edges_dual_normal_cell_v1(:,:,:)
+    !$ser verbatim deallocate(edges_dual_normal_cell_v1)
+
+    !$ser verbatim allocate(edges_dual_normal_cell_v2(size(edges_dual_normal_cell, 1),size(edges_dual_normal_cell, 2),size(edges_dual_normal_cell, 3)))
+    !$ser data edges_dual_normal_cell_v2=edges_dual_normal_cell_v2(:,:,:)
+    !$ser verbatim deallocate(edges_dual_normal_cell_v2)
+
+    PRINT *, 'Serializing cvd_o_rd=cvd_o_rd'
+
+    !$ser data cvd_o_rd=cvd_o_rd
+
+    PRINT *, 'Serializing grav=grav'
+
+    !$ser data grav=grav
+
+    PRINT *, 'Serializing jg=jg'
+
+    !$ser data jg=jg
+
+    PRINT *, 'Serializing nproma=nproma'
+
+    !$ser data nproma=nproma
+
+    PRINT *, 'Serializing nlev=nlev'
+
+    !$ser data nlev=nlev
+
+    PRINT *, 'Serializing nblks_e=nblks_e'
+
+    !$ser data nblks_e=nblks_e
+
+    PRINT *, 'Serializing nblks_v=nblks_v'
+
+    !$ser data nblks_v=nblks_v
+
+    PRINT *, 'Serializing nblks_c=nblks_c'
+
+    !$ser data nblks_c=nblks_c
+
+    PRINT *, 'Serializing nshift=nshift'
+
+    !$ser data nshift=nshift
+
+    PRINT *, 'Serializing nshift_total=nshift_total'
+
+    !$ser data nshift_total=nshift_total
+
+    PRINT *, 'Serializing nrdmax=nrdmax'
+
+    !$ser data nrdmax=nrdmax
+
+    PRINT *, 'Serializing ndyn_substeps=ndyn_substeps'
+
+    !$ser data ndyn_substeps=ndyn_substeps
+
+    PRINT *, 'Serializing hdiff_order=hdiff_order'
+
+    !$ser data hdiff_order=hdiff_order
+
+    PRINT *, 'Serializing itype_comm=itype_comm'
+
+    !$ser data itype_comm=itype_comm
+
+    PRINT *, 'Serializing itype_sher=itype_sher'
+
+    !$ser data itype_sher=itype_sher
+
+    PRINT *, 'Serializing itype_vn_diffu=itype_vn_diffu'
+
+    !$ser data itype_vn_diffu=itype_vn_diffu
+
+    PRINT *, 'Serializing itype_t_diffu=itype_t_diffu'
+
+    !$ser data itype_t_diffu=itype_t_diffu
+
+    PRINT *, 'Serializing hdiff_smag_z=hdiff_smag_z'
+
+    !$ser data hdiff_smag_z=hdiff_smag_z
+
+    PRINT *, 'Serializing hdiff_smag_z2=hdiff_smag_z2'
+
+    !$ser data hdiff_smag_z2=hdiff_smag_z2
+
+    PRINT *, 'Serializing hdiff_smag_z3=hdiff_smag_z3'
+
+    !$ser data hdiff_smag_z3=hdiff_smag_z3
+
+    PRINT *, 'Serializing hdiff_smag_z4=hdiff_smag_z4'
+
+    !$ser data hdiff_smag_z4=hdiff_smag_z4
+
+    PRINT *, 'Serializing hdiff_smag_fac=hdiff_smag_fac'
+
+    !$ser data hdiff_smag_fac=hdiff_smag_fac
+
+    PRINT *, 'Serializing hdiff_smag_fac2=hdiff_smag_fac2'
+
+    !$ser data hdiff_smag_fac2=hdiff_smag_fac2
+
+    PRINT *, 'Serializing hdiff_smag_fac3=hdiff_smag_fac3'
+
+    !$ser data hdiff_smag_fac3=hdiff_smag_fac3
+
+    PRINT *, 'Serializing hdiff_smag_fac4=hdiff_smag_fac4'
+
+    !$ser data hdiff_smag_fac4=hdiff_smag_fac4
+
+    PRINT *, 'Serializing hdiff_efdt_ratio=hdiff_efdt_ratio'
+
+    !$ser data hdiff_efdt_ratio=hdiff_efdt_ratio
+
+    PRINT *, 'Serializing k4=k4'
+
+    !$ser data k4=k4
+
+    PRINT *, 'Serializing k4w=k4w'
+
+    !$ser data k4w=k4w
+
+    PRINT *, 'Serializing nudge_max_coeff=nudge_max_coeff'
+
+    !$ser data nudge_max_coeff=nudge_max_coeff
+
+    PRINT *, 'Serializing denom_diffu_v=denom_diffu_v'
+
+    !$ser data denom_diffu_v=denom_diffu_v
+
+    PRINT *, 'Serializing p_test_run=p_test_run'
+
+    !$ser data p_test_run=p_test_run
+
+    PRINT *, 'Serializing lphys=lphys'
+
+    !$ser data lphys=lphys
+
+    PRINT *, 'Serializing lhdiff_rcf=lhdiff_rcf'
+
+    !$ser data lhdiff_rcf=lhdiff_rcf
+
+    PRINT *, 'Serializing lhdiff_w=lhdiff_w'
+
+    !$ser data lhdiff_w=lhdiff_w
+
+    PRINT *, 'Serializing lhdiff_temp=lhdiff_temp'
+
+    !$ser data lhdiff_temp=lhdiff_temp
+
+    PRINT *, 'Serializing l_zdiffu_t=l_zdiffu_t'
+
+    !$ser data l_zdiffu_t=l_zdiffu_t
+
+    PRINT *, 'Serializing l_limited_area=l_limited_area'
+
+    !$ser data l_limited_area=l_limited_area
+
+    PRINT *, 'Serializing lfeedback=lfeedback'
+
+    !$ser data lfeedback=lfeedback
+
+    PRINT *, 'Serializing ltkeshs=ltkeshs'
+
+    !$ser data ltkeshs=ltkeshs
+
+    PRINT *, 'Serializing lsmag_3d=lsmag_3d'
+
+    !$ser data lsmag_3d=lsmag_3d
+
+    PRINT *, 'Serializing lvert_nest=lvert_nest'
+
+    !$ser data lvert_nest=lvert_nest
+
+    PRINT *, 'Serializing ltimer=ltimer'
+
+    !$ser data ltimer=ltimer
+
+    PRINT *, 'Serializing ddt_vn_hdf_is_associated=ddt_vn_hdf_is_associated'
+
+    !$ser data ddt_vn_hdf_is_associated=ddt_vn_hdf_is_associated
+
+    PRINT *, 'Serializing ddt_vn_dyn_is_associated=ddt_vn_dyn_is_associated'
+
+    !$ser data ddt_vn_dyn_is_associated=ddt_vn_dyn_is_associated
+
+    PRINT *, 'Serializing vct_a=vct_a(:)'
+
+    IF (SIZE(vct_a) > 0) THEN
+       !$ser data vct_a=vct_a(:)
+    ELSE
+       PRINT *, 'Warning: Array vct_a has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing c_lin_e=c_lin_e(:,:,:)'
+
+    IF (SIZE(c_lin_e) > 0) THEN
+       !$ser data c_lin_e=c_lin_e(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array c_lin_e has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing e_bln_c_s=e_bln_c_s(:,:,:)'
+
+    IF (SIZE(e_bln_c_s) > 0) THEN
+       !$ser data e_bln_c_s=e_bln_c_s(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array e_bln_c_s has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing e_bln_c_u=e_bln_c_u(:,:,:)'
+
+    IF (SIZE(e_bln_c_u) > 0) THEN
+       !$ser data e_bln_c_u=e_bln_c_u(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array e_bln_c_u has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing e_bln_c_v=e_bln_c_v(:,:,:)'
+
+    IF (SIZE(e_bln_c_v) > 0) THEN
+       !$ser data e_bln_c_v=e_bln_c_v(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array e_bln_c_v has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing cells_aw_verts=cells_aw_verts(:,:,:)'
+
+    IF (SIZE(cells_aw_verts) > 0) THEN
+       !$ser data cells_aw_verts=cells_aw_verts(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array cells_aw_verts has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing geofac_div=geofac_div(:,:,:)'
+
+    IF (SIZE(geofac_div) > 0) THEN
+       !$ser data geofac_div=geofac_div(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array geofac_div has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing geofac_rot=geofac_rot(:,:,:)'
+
+    IF (SIZE(geofac_rot) > 0) THEN
+       !$ser data geofac_rot=geofac_rot(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array geofac_rot has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing geofac_n2s=geofac_n2s(:,:,:)'
+
+    IF (SIZE(geofac_n2s) > 0) THEN
+       !$ser data geofac_n2s=geofac_n2s(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array geofac_n2s has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing geofac_grg=geofac_grg(:,:,:,:)'
+
+    IF (SIZE(geofac_grg) > 0) THEN
+       !$ser data geofac_grg=geofac_grg(:,:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array geofac_grg has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing nudgecoeff_e=nudgecoeff_e(:,:)'
+
+    IF (SIZE(nudgecoeff_e) > 0) THEN
+       !$ser data nudgecoeff_e=nudgecoeff_e(:,:)
+    ELSE
+       PRINT *, 'Warning: Array nudgecoeff_e has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing rbf_vec_idx_v=rbf_vec_idx_v(:,:,:)'
+
+    IF (SIZE(rbf_vec_idx_v) > 0) THEN
+       !$ser data rbf_vec_idx_v=rbf_vec_idx_v(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array rbf_vec_idx_v has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing rbf_vec_blk_v=rbf_vec_blk_v(:,:,:)'
+
+    IF (SIZE(rbf_vec_blk_v) > 0) THEN
+       !$ser data rbf_vec_blk_v=rbf_vec_blk_v(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array rbf_vec_blk_v has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing rbf_vec_coeff_v=rbf_vec_coeff_v(:,:,:,:)'
+
+    IF (SIZE(rbf_vec_coeff_v) > 0) THEN
+       !$ser data rbf_vec_coeff_v=rbf_vec_coeff_v(:,:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array rbf_vec_coeff_v has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing enhfac_diffu=enhfac_diffu(:)'
+
+    IF (SIZE(enhfac_diffu) > 0) THEN
+       !$ser data enhfac_diffu=enhfac_diffu(:)
+    ELSE
+       PRINT *, 'Warning: Array enhfac_diffu has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing zd_intcoef=zd_intcoef(:,:)'
+
+    IF (SIZE(zd_intcoef) > 0) THEN
+       !$ser data zd_intcoef=zd_intcoef(:,:)
+    ELSE
+       PRINT *, 'Warning: Array zd_intcoef has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing zd_geofac=zd_geofac(:,:)'
+
+    IF (SIZE(zd_geofac) > 0) THEN
+       !$ser data zd_geofac=zd_geofac(:,:)
+    ELSE
+       PRINT *, 'Warning: Array zd_geofac has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing zd_diffcoef=zd_diffcoef(:)'
+
+    IF (SIZE(zd_diffcoef) > 0) THEN
+       !$ser data zd_diffcoef=zd_diffcoef(:)
+    ELSE
+       PRINT *, 'Warning: Array zd_diffcoef has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing wgtfac_c=wgtfac_c(:,:,:)'
+
+    IF (SIZE(wgtfac_c) > 0) THEN
+       !$ser data wgtfac_c=wgtfac_c(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array wgtfac_c has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing wgtfac_e=wgtfac_e(:,:,:)'
+
+    IF (SIZE(wgtfac_e) > 0) THEN
+       !$ser data wgtfac_e=wgtfac_e(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array wgtfac_e has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing wgtfacq_e=wgtfacq_e(:,:,:)'
+
+    IF (SIZE(wgtfacq_e) > 0) THEN
+       !$ser data wgtfacq_e=wgtfacq_e(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array wgtfacq_e has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing wgtfacq1_e=wgtfacq1_e(:,:,:)'
+
+    IF (SIZE(wgtfacq1_e) > 0) THEN
+       !$ser data wgtfacq1_e=wgtfacq1_e(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array wgtfacq1_e has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing ddqz_z_full_e=ddqz_z_full_e(:,:,:)'
+
+    IF (SIZE(ddqz_z_full_e) > 0) THEN
+       !$ser data ddqz_z_full_e=ddqz_z_full_e(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array ddqz_z_full_e has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing theta_ref_mc=theta_ref_mc(:,:,:)'
+
+    IF (SIZE(theta_ref_mc) > 0) THEN
+       !$ser data theta_ref_mc=theta_ref_mc(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array theta_ref_mc has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing zd_indlist=zd_indlist(:,:)'
+
+    IF (SIZE(zd_indlist) > 0) THEN
+       !$ser data zd_indlist=zd_indlist(:,:)
+    ELSE
+       PRINT *, 'Warning: Array zd_indlist has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing zd_blklist=zd_blklist(:,:)'
+
+    IF (SIZE(zd_blklist) > 0) THEN
+       !$ser data zd_blklist=zd_blklist(:,:)
+    ELSE
+       PRINT *, 'Warning: Array zd_blklist has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing zd_vertidx=zd_vertidx(:,:)'
+
+    IF (SIZE(zd_vertidx) > 0) THEN
+       !$ser data zd_vertidx=zd_vertidx(:,:)
+    ELSE
+       PRINT *, 'Warning: Array zd_vertidx has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing zd_listdim=zd_listdim'
+
+    !$ser data zd_listdim=zd_listdim
+
+    PRINT *, 'Serializing edges_start_block=edges_start_block(min_rledge:)'
+
+    IF (SIZE(edges_start_block) > 0) THEN
+       !$ser data edges_start_block=edges_start_block(min_rledge:)
+    ELSE
+       PRINT *, 'Warning: Array edges_start_block has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing edges_end_block=edges_end_block(min_rledge:)'
+
+    IF (SIZE(edges_end_block) > 0) THEN
+       !$ser data edges_end_block=edges_end_block(min_rledge:)
+    ELSE
+       PRINT *, 'Warning: Array edges_end_block has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing edges_start_index=edges_start_index(min_rledge:)'
+
+    IF (SIZE(edges_start_index) > 0) THEN
+       !$ser data edges_start_index=edges_start_index(min_rledge:)
+    ELSE
+       PRINT *, 'Warning: Array edges_start_index has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing edges_end_index=edges_end_index(min_rledge:)'
+
+    IF (SIZE(edges_end_index) > 0) THEN
+       !$ser data edges_end_index=edges_end_index(min_rledge:)
+    ELSE
+       PRINT *, 'Warning: Array edges_end_index has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing edges_vertex_idx=edges_vertex_idx(:,:,:)'
+
+    IF (SIZE(edges_vertex_idx) > 0) THEN
+       !$ser data edges_vertex_idx=edges_vertex_idx(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array edges_vertex_idx has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing edges_vertex_blk=edges_vertex_blk(:,:,:)'
+
+    IF (SIZE(edges_vertex_blk) > 0) THEN
+       !$ser data edges_vertex_blk=edges_vertex_blk(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array edges_vertex_blk has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing edges_cell_idx=edges_cell_idx(:,:,:)'
+
+    IF (SIZE(edges_cell_idx) > 0) THEN
+       !$ser data edges_cell_idx=edges_cell_idx(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array edges_cell_idx has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing edges_cell_blk=edges_cell_blk(:,:,:)'
+
+    IF (SIZE(edges_cell_blk) > 0) THEN
+       !$ser data edges_cell_blk=edges_cell_blk(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array edges_cell_blk has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing edges_tangent_orientation=edges_tangent_orientation(:,:)'
+
+    IF (SIZE(edges_tangent_orientation) > 0) THEN
+       !$ser data edges_tangent_orientation=edges_tangent_orientation(:,:)
+    ELSE
+       PRINT *, 'Warning: Array edges_tangent_orientation has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing edges_inv_vert_vert_length=edges_inv_vert_vert_length(:,:)'
+
+    IF (SIZE(edges_inv_vert_vert_length) > 0) THEN
+       !$ser data edges_inv_vert_vert_length=edges_inv_vert_vert_length(:,:)
+    ELSE
+       PRINT *, 'Warning: Array edges_inv_vert_vert_length has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing edges_inv_primal_edge_length=edges_inv_primal_edge_length(:,:)'
+
+    IF (SIZE(edges_inv_primal_edge_length) > 0) THEN
+       !$ser data edges_inv_primal_edge_length=edges_inv_primal_edge_length(:,:)
+    ELSE
+       PRINT *, 'Warning: Array edges_inv_primal_edge_length has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing edges_inv_dual_edge_length=edges_inv_dual_edge_length(:,:)'
+
+    IF (SIZE(edges_inv_dual_edge_length) > 0) THEN
+       !$ser data edges_inv_dual_edge_length=edges_inv_dual_edge_length(:,:)
+    ELSE
+       PRINT *, 'Warning: Array edges_inv_dual_edge_length has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing edges_area_edge=edges_area_edge(:,:)'
+
+    IF (SIZE(edges_area_edge) > 0) THEN
+       !$ser data edges_area_edge=edges_area_edge(:,:)
+    ELSE
+       PRINT *, 'Warning: Array edges_area_edge has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing cells_start_block=cells_start_block(min_rlcell:)'
+
+    IF (SIZE(cells_start_block) > 0) THEN
+       !$ser data cells_start_block=cells_start_block(min_rlcell:)
+    ELSE
+       PRINT *, 'Warning: Array cells_start_block has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing cells_end_block=cells_end_block(min_rlcell:)'
+
+    IF (SIZE(cells_end_block) > 0) THEN
+       !$ser data cells_end_block=cells_end_block(min_rlcell:)
+    ELSE
+       PRINT *, 'Warning: Array cells_end_block has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing cells_start_index=cells_start_index(min_rlcell:)'
+
+    IF (SIZE(cells_start_index) > 0) THEN
+       !$ser data cells_start_index=cells_start_index(min_rlcell:)
+    ELSE
+       PRINT *, 'Warning: Array cells_start_index has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing cells_end_index=cells_end_index(min_rlcell:)'
+
+    IF (SIZE(cells_end_index) > 0) THEN
+       !$ser data cells_end_index=cells_end_index(min_rlcell:)
+    ELSE
+       PRINT *, 'Warning: Array cells_end_index has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing cells_neighbor_idx=cells_neighbor_idx(:,:,:)'
+
+    IF (SIZE(cells_neighbor_idx) > 0) THEN
+       !$ser data cells_neighbor_idx=cells_neighbor_idx(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array cells_neighbor_idx has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing cells_neighbor_blk=cells_neighbor_blk(:,:,:)'
+
+    IF (SIZE(cells_neighbor_blk) > 0) THEN
+       !$ser data cells_neighbor_blk=cells_neighbor_blk(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array cells_neighbor_blk has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing cells_edge_idx=cells_edge_idx(:,:,:)'
+
+    IF (SIZE(cells_edge_idx) > 0) THEN
+       !$ser data cells_edge_idx=cells_edge_idx(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array cells_edge_idx has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing cells_edge_blk=cells_edge_blk(:,:,:)'
+
+    IF (SIZE(cells_edge_blk) > 0) THEN
+       !$ser data cells_edge_blk=cells_edge_blk(:,:,:)
+    ELSE
+       PRINT *, 'Warning: Array cells_edge_blk has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing cells_area=cells_area(:,:)'
+
+    IF (SIZE(cells_area) > 0) THEN
+       !$ser data cells_area=cells_area(:,:)
+    ELSE
+       PRINT *, 'Warning: Array cells_area has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing verts_start_block=verts_start_block(min_rlvert:)'
+
+    IF (SIZE(verts_start_block) > 0) THEN
+       !$ser data verts_start_block=verts_start_block(min_rlvert:)
+    ELSE
+       PRINT *, 'Warning: Array verts_start_block has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing verts_end_block=verts_end_block(min_rlvert:)'
+
+    IF (SIZE(verts_end_block) > 0) THEN
+       !$ser data verts_end_block=verts_end_block(min_rlvert:)
+    ELSE
+       PRINT *, 'Warning: Array verts_end_block has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing verts_start_index=verts_start_index(min_rlvert:)'
+
+    IF (SIZE(verts_start_index) > 0) THEN
+       !$ser data verts_start_index=verts_start_index(min_rlvert:)
+    ELSE
+       PRINT *, 'Warning: Array verts_start_index has size 0. Not serializing array.'
+    END IF
+
+    PRINT *, 'Serializing verts_end_index=verts_end_index(min_rlvert:)'
+
+    IF (SIZE(verts_end_index) > 0) THEN
+       !$ser data verts_end_index=verts_end_index(min_rlvert:)
+    ELSE
+       PRINT *, 'Warning: Array verts_end_index has size 0. Not serializing array.'
+    END IF