Close open units explicitly

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)