From 5a447e7e2767473ed4f90f1982bfacba375d266c Mon Sep 17 00:00:00 2001
From: "Igor S. Gerasimov" <i.s.ger>
Date: Mon, 7 Feb 2022 07:40:11 +0900
Subject: [PATCH 1/2] Fix compilation errors with gfortan

---
 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 13c5c73..455e30d 100644
--- a/fileIO.f90
+++ b/fileIO.f90
@@ -3455,7 +3455,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 9a8b0d1..066cc59 100644
--- a/grid.f90
+++ b/grid.f90
@@ -818,7 +818,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 4c1ba08..1e46863 100644
--- a/util.f90
+++ b/util.f90
@@ -2461,7 +2461,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.33.1