workaround for gfortran 5.3 bug

_pygotm.pyx

@@ -10,7 +10,7 @@ cimport numpy
 import numpy
 
 cdef extern void initialize(int nlev, const char* nml_file, const char* yaml_file, double** ptke, double** peps, double** pL, double** pnum, double** pnuh) nogil
-cdef extern void finalize() nogil
+cdef extern void finalize2() nogil
 cdef extern void calculate(int nlev, double dt, const double* h, double D, double taus, double taub, double z0s, double z0b, const double* SS, const double* NN) nogil
 cdef extern void calculate_3d(int nx, int nz, int nz, int istart, int istop, int jstart, int jstop, double dt, const int* mask, const double* h3d, const double* D, const double* u_taus, const double* u_taub, const double* z0s, const double* z0b, const double* NN, const double* SS, double* tke, double* eps, double* L, double* num, double* nuh)
 cdef extern void diff(int nlev, double dt, double cnpar, int posconc, const double* h, int Bcup, int Bcdw, double Yup, double Ydw, const double* nuY, const double* Lsour, const double* Qsour, const double* Taur, const double* Yobs, double* Y)
@@ -25,7 +25,7 @@ stdout = open(stdout_path)
 stderr = open(stderr_path)
 
 def cleanup():
-    finalize()
+    finalize2()
     close_redirected_output()
     stdout.close()
     stderr.close()

pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "pygotm"
-version = "0.1.1"
+version = "0.1.2"
 authors = [
   {name = "Jorn Bruggeman", email = "jorn@bolding-bruggeman.com"},
   {name = "Karsten Bolding", email = "karsten@bolding-bruggeman.com"}

src/pygotm_wrapper.F90

@@ -38,7 +38,7 @@ subroutine initialize(nlev, pnml_file, pyaml_file, ptke, peps, pL, pnum, pnuh) b
       call c_f_pointer(c_loc(pnml_file), nml_file)
       call c_f_pointer(c_loc(pyaml_file), yaml_file)
 
-      call finalize()
+      call finalize2()
       allocate(store)
       pathlen = index(yaml_file, C_NULL_CHAR) - 1
       if (pathlen > 0) then
@@ -63,7 +63,9 @@ subroutine initialize(nlev, pnml_file, pyaml_file, ptke, peps, pL, pnum, pnuh) b
       flush(unit=6)
    end subroutine
 
-   subroutine finalize() bind(c)
+   ! Note: finalize was renamend to finalize2 to prevent gfortran confusing
+   ! it with store%finalize
+   subroutine finalize2() bind(c)
       if (associated(store)) then
          call clean_turbulence()
          call clean_tridiagonal()