output the difference of filtered to unfiltered

ecmwf · Nov 29, 2023 · 6240300 · 6240300
1 parent c5e0801
commit 6240300
Showing 1 changed file with 9 additions and 0 deletions.
diff --git a/src/sandbox/interpolation/atlas-filter.F90 b/src/sandbox/interpolation/atlas-filter.F90
@@ -249,6 +249,15 @@ program filtering
     print *, " filter.exe in seconds: ", end_time - start_time
 
     call sfield%rename("filtered")
+    call sfield%halo_exchange()
+    call gmsh%write(sfield)
+
+    ! output the difference of the unfiltered and the filtered field
+    do inode = 1, nb_nodes
+        sfield_v(inode) = sfield_v(inode) - MDPI_gulfstream(lonlat(1,inode), lonlat(2,inode))
+    end do
+    call sfield%rename("filt-unfilt")
+    call sfield%halo_exchange()
     call gmsh%write(sfield)
 
     call atlas_library%finalise()