From b7bec590d6f67c37ac1afa7ca1e8ec428c08d0ba Mon Sep 17 00:00:00 2001
From: Lorenzo Figini <lorenzo.figini@istp.cnr.it>
Date: Tue, 1 Feb 2022 16:06:32 +0100
Subject: [PATCH] Close open units explicitly

---
 src/main.f90 | 16 ++++++++++++++++
 1 file changed, 16 insertions(+)

diff --git a/src/main.f90 b/src/main.f90
index 8dc5c7e..85ce893 100644
--- a/src/main.f90
+++ b/src/main.f90
@@ -127,6 +127,22 @@ program main_std
   print*,'Pabs (MW), Icd (kA) = ', pec,icd*1.0e3_wp_
 ! ======= control prints  END  =======
 
+! ======= flush open files BEGIN =======
+  close(ucenr)
+  close(usumm)
+  close(uprj0)
+  close(uprj0+1)
+  close(uwbm)
+  close(udisp)
+  close(ubres)
+  close(ucnt)
+  close(uoutr)
+  close(ueq)
+  close(uprfin)
+  close(uflx)
+  close(upec)
+! ======= flush open files  END  =======
+
 ! ======= free memory BEGIN =======
   if(allocated(psrad)) deallocate(psrad)
   if(allocated(terad)) deallocate(terad, derad, zfc)