From 3b31b7f02e74647ca0e6f106163cdcba71fadfaf Mon Sep 17 00:00:00 2001
From: "Igor S. Gerasimov" <i.s.ger>
Date: Mon, 14 Feb 2022 14:13:57 +0900
Subject: [PATCH] Fix gfortran build

---
 fileIO.f90 | 4 +++-
 grid.f90   | 2 +-
 util.f90   | 2 +-
 3 files changed, 5 insertions(+), 3 deletions(-)

diff --git a/fileIO.f90 b/fileIO.f90
index dd3fd99..26b2ece 100644
--- a/fileIO.f90
+++ b/fileIO.f90
@@ -3456,7 +3456,9 @@ read(10,*)
 read(10,*) a%charge
 call loclabel(10,"<Nuclear Cartesian Coordinates>")
 read(10,*)
-read(10,*) ((a(i)%x,a(i)%y,a(i)%z),i=1,ncenter)
+do i = 1, ncenter
+  read(10,*) a(i)%x,a(i)%y,a(i)%z
+end do
 call loclabel(10,"<Number of Electrons>")
 read(10,*)
 read(10,*) nelec
diff --git a/grid.f90 b/grid.f90
index d449659..6f1841c 100644
--- a/grid.f90
+++ b/grid.f90
@@ -827,7 +827,7 @@ end subroutine
 
 
 !!----------- If three vectors of the grid data are parallel to X/Y/Z axis, respectively. Return 1 means yes, 0 means no
-integer function ifgridortho
+integer function ifgridortho()
 use defvar
 if (abs(gridv1(2))<1E-10.and.abs(gridv1(3))<1E-10.and.abs(gridv2(1))<1E-10.and.abs(gridv2(3))<1E-10.and.abs(gridv3(1))<1E-10.and.abs(gridv3(2))<1E-10) then
     ifgridortho=1
diff --git a/util.f90 b/util.f90
index 66e8c8b..8ab006a 100644
--- a/util.f90
+++ b/util.f90
@@ -2513,7 +2513,7 @@ if (itmp==4) then
 	c80tmp(79:79)='/'
 	itmp=0
 end if
-write(*,"(2a\)") trim(c80tmp),char(13)
+write(*,"(2a)", advance="no") trim(c80tmp),char(13)
 if (inow>=nall) write(*,*)
 end subroutine
 
-- 
2.25.1