diff --git a/model/atmosphere/subgrid_scale_physics/muphys/tests/muphys/integration_tests/run_graupel_only.py b/model/atmosphere/subgrid_scale_physics/muphys/tests/muphys/integration_tests/run_graupel_only.py
index 1bda067112..be6d97a9dd 100755
--- a/model/atmosphere/subgrid_scale_physics/muphys/tests/muphys/integration_tests/run_graupel_only.py
+++ b/model/atmosphere/subgrid_scale_physics/muphys/tests/muphys/integration_tests/run_graupel_only.py
@@ -477,3 +477,12 @@ def write_fields(
     pflx=np.transpose(pflx_out.asnumpy()),
     pre_gsp=data.pre_gsp,
 )
+
+from cdo import Cdo
+
+cdo = Cdo()
+
+for varname in ['ta','hus','clw','cli','qr','qs','qg']:
+    cdo.diffn('abslim=3e-13',
+              input=f' -selname,{varname} junk.nc -selname,{varname}',
+              output='../../MUPHYS/muphys/build_apple/new_fields.nc')