ncid => file_id

Author: scivision

src/attributes.f90

@@ -11,18 +11,18 @@
 
 if(.not.self%is_open) error stop 'nc4fortran:write_attribute: file handle not open'
 
-ier = nf90_inq_varid(self%ncid, dname, varid)
+ier = nf90_inq_varid(self%file_id, dname, varid)
 
 if(ier == nf90_noerr) then
 select type(value)
 type is (character(*))
-   ier = nf90_put_att(self%ncid, varid, attrname, value)
+   ier = nf90_put_att(self%file_id, varid, attrname, value)
 type is (real(real32))
-  ier = nf90_put_att(self%ncid, varid, attrname, value)
+  ier = nf90_put_att(self%file_id, varid, attrname, value)
 type is (real(real64))
-  ier = nf90_put_att(self%ncid, varid, attrname, value)
+  ier = nf90_put_att(self%file_id, varid, attrname, value)
 type is (integer(int32))
-  ier = nf90_put_att(self%ncid, varid, attrname, value)
+  ier = nf90_put_att(self%file_id, varid, attrname, value)
 class default
   ier = NF90_EBADTYPE
 end select
@@ -38,18 +38,18 @@
 
 if(.not.self%is_open) error stop 'nc4fortran:read_attribute: file handle not open'
 
-ier = nf90_inq_varid(self%ncid, dname, varid)
+ier = nf90_inq_varid(self%file_id, dname, varid)
 
 if(ier == nf90_noerr) then
 select type (value)
 type is (character(*))
-  ier = nf90_get_att(self%ncid, varid, attrname, value)
+  ier = nf90_get_att(self%file_id, varid, attrname, value)
 type is (real(real32))
-  ier = nf90_get_att(self%ncid, varid, attrname, value)
+  ier = nf90_get_att(self%file_id, varid, attrname, value)
 type is (real(real64))
-  ier = nf90_get_att(self%ncid, varid, attrname, value)
+  ier = nf90_get_att(self%file_id, varid, attrname, value)
 type is (integer(int32))
-  ier = nf90_get_att(self%ncid, varid, attrname, value)
+  ier = nf90_get_att(self%file_id, varid, attrname, value)
 class default
   ier = NF90_EBADTYPE
 end select

src/interface.f90

@@ -21,7 +21,7 @@ module nc4fortran
 type :: netcdf_file
 
 character(:), allocatable  :: filename
-integer :: ncid   !< location identifier
+integer :: file_id   !< location identifier
 
 integer :: comp_lvl = 0 !< compression level (1-9)  0: disable compression
 logical :: verbose = .false.
@@ -201,46 +201,53 @@ module subroutine nc_read_scalar(self, dname, value)
 !! inout for character
 end subroutine
 
-module subroutine nc_read_1d(self, dname, value)
+module subroutine nc_read_1d(self, dname, value, istart, iend)
 class(netcdf_file), intent(in)     :: self
 character(*), intent(in)         :: dname
 class(*), intent(inout)      :: value(:)
+integer, intent(in), dimension(1), optional :: istart, iend
 end subroutine
 
-module subroutine nc_read_2d(self, dname, value)
+module subroutine nc_read_2d(self, dname, value, istart, iend)
 class(netcdf_file), intent(in)     :: self
 character(*), intent(in)         :: dname
 class(*), intent(inout)      :: value(:,:)
+integer, intent(in), dimension(2), optional :: istart, iend
 end subroutine
 
-module subroutine nc_read_3d(self, dname, value)
+module subroutine nc_read_3d(self, dname, value, istart, iend)
 class(netcdf_file), intent(in)     :: self
 character(*), intent(in)         :: dname
 class(*), intent(inout)      :: value(:,:,:)
+integer, intent(in), dimension(3), optional :: istart, iend
 end subroutine
 
-module subroutine nc_read_4d(self, dname, value)
+module subroutine nc_read_4d(self, dname, value, istart, iend)
 class(netcdf_file), intent(in)     :: self
 character(*), intent(in)         :: dname
 class(*), intent(inout)      :: value(:,:,:,:)
+integer, intent(in), dimension(4), optional :: istart, iend
 end subroutine
 
-module subroutine nc_read_5d(self, dname, value)
+module subroutine nc_read_5d(self, dname, value, istart, iend)
 class(netcdf_file), intent(in)     :: self
 character(*), intent(in)         :: dname
 class(*), intent(inout)      :: value(:,:,:,:,:)
+integer, intent(in), dimension(5), optional :: istart, iend
 end subroutine
 
-module subroutine nc_read_6d(self, dname, value)
+module subroutine nc_read_6d(self, dname, value, istart, iend)
 class(netcdf_file), intent(in)     :: self
 character(*), intent(in)         :: dname
 class(*), intent(inout)      :: value(:,:,:,:,:,:)
+integer, intent(in), dimension(6), optional :: istart, iend
 end subroutine
 
-module subroutine nc_read_7d(self, dname, value)
+module subroutine nc_read_7d(self, dname, value, istart, iend)
 class(netcdf_file), intent(in)     :: self
 character(*), intent(in)         :: dname
 class(*), intent(inout)      :: value(:,:,:,:,:,:,:)
+integer, intent(in), dimension(7), optional :: istart, iend
 end subroutine
 
 end interface

src/read.f90

@@ -14,14 +14,14 @@
 
 chunk_size = -1
 
-i = nf90_inq_varid(self%ncid, dname, varid)
+i = nf90_inq_varid(self%file_id, dname, varid)
 if (i/=NF90_NOERR) error stop 'ERROR:nc4fortran:chunk: cannot find variable: ' // dname
 
-i = nf90_inquire_variable(self%ncid, varid, contiguous=contig)
+i = nf90_inquire_variable(self%file_id, varid, contiguous=contig)
 if (i/=NF90_NOERR) error stop 'ERROR:nc4fortran:chunk: cannot get variable properties' // dname
 
 if(contig) return
-i = nf90_inquire_variable(self%ncid, varid, chunksizes=chunk_size)
+i = nf90_inquire_variable(self%file_id, varid, chunksizes=chunk_size)
 if (i/=NF90_NOERR) error stop 'ERROR:nc4fortran:chunk: cannot get variable properties' // dname
 
 end procedure get_chunk
@@ -34,14 +34,14 @@
 
 get_deflate = .false.
 
-i = nf90_inq_varid(self%ncid, dname, varid)
+i = nf90_inq_varid(self%file_id, dname, varid)
 if (i/=NF90_NOERR) error stop 'ERROR:nc4fortran:get_deflate: cannot find variable: ' // dname
 
-i = nf90_inquire_variable(self%ncid, varid, contiguous=contig)
+i = nf90_inquire_variable(self%file_id, varid, contiguous=contig)
 if (i/=NF90_NOERR) error stop 'ERROR:nc4fortran:get_deflate: cannot get variable properties' // dname
 
 if(contig) return
-i = nf90_inquire_variable(self%ncid, varid, deflate_level=deflate_level)
+i = nf90_inquire_variable(self%file_id, varid, deflate_level=deflate_level)
 if (i/=NF90_NOERR) error stop 'ERROR:nc4fortran:get_deflate: cannot get variable properties' // dname
 
 get_deflate = deflate_level /= 0
@@ -56,10 +56,10 @@
 
 drank = -1
 
-ierr = nf90_inq_varid(self%ncid, dname, varid)
+ierr = nf90_inq_varid(self%file_id, dname, varid)
 if(ierr/=NF90_NOERR) error stop 'ERROR:nc4fortran:get_ndim: could not get variable ID for ' // dname
 
-ierr = nf90_inquire_variable(self%ncid, varid, ndims=drank)
+ierr = nf90_inquire_variable(self%file_id, varid, ndims=drank)
 if(ierr/=NF90_NOERR) error stop 'ERROR:nc4fortran:get_ndim: could not get rank for ' // dname
 
 end procedure get_ndim
@@ -75,19 +75,19 @@
 
 allocate(dimids(N), dims(N))
 
-ier = nf90_inq_varid(self%ncid, dname, varid)
+ier = nf90_inq_varid(self%file_id, dname, varid)
 if(check_error(ier, dname)) error stop 'ERROR:nc4fortran:get_shape: could not get variable ID for: ' // dname
 
-ier = nf90_inquire_variable(self%ncid, varid, dimids = dimids)
+ier = nf90_inquire_variable(self%file_id, varid, dimids = dimids)
 if(check_error(ier, dname)) error stop 'ERROR:nc4fortran:get_shape: could not get dimension IDs for: ' // dname
 
 if (present(dimnames)) allocate(tempnames(N))
 
 do i = 1,N
   if(present(dimnames)) then
-    ier = nf90_inquire_dimension(self%ncid, dimid=dimids(i), name=tempnames(i), len=dims(i))
+    ier = nf90_inquire_dimension(self%file_id, dimid=dimids(i), name=tempnames(i), len=dims(i))
   else
-    ier = nf90_inquire_dimension(self%ncid, dimid=dimids(i), len=dims(i))
+    ier = nf90_inquire_dimension(self%file_id, dimid=dimids(i), len=dims(i))
   endif
   if(ier/=NF90_NOERR) error stop 'ERROR:nc4fortran:get_shape: querying dimension size'
 enddo
@@ -107,7 +107,7 @@
 
 if(.not.self%is_open) error stop 'ERROR:nc4fortran:exist: file handle not open '
 
-ierr = nf90_inq_varid(self%ncid, dname, varid)
+ierr = nf90_inq_varid(self%file_id, dname, varid)
 
 select case (ierr)
 case (NF90_NOERR)

src/read_scalar.f90

@@ -12,32 +12,32 @@
 
 if(.not.self%is_open) error stop 'ERROR:nc4fortran:reader file handle not open'
 
-ier = nf90_inq_varid(self%ncid, dname, varid)
+ier = nf90_inq_varid(self%file_id, dname, varid)
 if (check_error(ier, dname)) error stop 'ERROR:nc4fortran:read inquire_id ' // dname // ' in ' // self%filename
 
 select type (value)
 type is (character(*))
   !! NetCDF4 requires knowing the exact character length as if it were an array without fill values
   !! HDF5 is not this strict; having a longer string variable than disk variable is OK in HDF5, but not NetCDF4
-  ier = nf90_inquire_variable(self%ncid, varid, dimids = dimids)
+  ier = nf90_inquire_variable(self%file_id, varid, dimids = dimids)
   if(check_error(ier, dname)) error stop 'ERROR:nc4fortran:read_scalar: could not get dimension IDs for: ' // dname
 
-  ier = nf90_inquire_dimension(self%ncid, dimid=dimids(1), len=dims)
+  ier = nf90_inquire_dimension(self%file_id, dimid=dimids(1), len=dims)
   if(ier/=NF90_NOERR) error stop 'ERROR:nc4fortran:read_scalar: querying dimension size'
 
   allocate(character(dims) :: buf)
-  ier = nf90_get_var(self%ncid, varid, buf)
+  ier = nf90_get_var(self%file_id, varid, buf)
   i = index(buf, c_null_char) - 1
   if(i == -1) i = len_trim(buf)
   value = buf(:i)
 type is (real(real64))
-  ier = nf90_get_var(self%ncid, varid, value)
+  ier = nf90_get_var(self%file_id, varid, value)
 type is (real(real32))
-  ier = nf90_get_var(self%ncid, varid, value)
+  ier = nf90_get_var(self%file_id, varid, value)
 type is (integer(int64))
-  ier = nf90_get_var(self%ncid, varid, value)
+  ier = nf90_get_var(self%file_id, varid, value)
 type is (integer(int32))
-  ier = nf90_get_var(self%ncid, varid, value)
+  ier = nf90_get_var(self%file_id, varid, value)
 class default
   ier = NF90_EBADTYPE
 end select

src/reader.inc

@@ -2,10 +2,10 @@ integer :: varid, ier, drank
 
 if(.not.self%is_open) error stop 'nc4fortran:reader file handle not open'
 
-ier = nf90_inq_varid(self%ncid, dname, varid)
+ier = nf90_inq_varid(self%file_id, dname, varid)
 if (check_error(ier, dname)) error stop 'nc4fortran:read inquire_id ' // dname // ' in ' // self%filename
 
-ier = nf90_inquire_variable(self%ncid, varid, ndims=drank)
+ier = nf90_inquire_variable(self%file_id, varid, ndims=drank)
 if (check_error(ier, dname)) error stop 'nc4fortran:read inquire ' // dname // ' in ' // self%filename
 if(drank /= rank(value)) then
   write(stderr,*) 'ERROR:nc4fortran:read ' // dname // ' rank ', drank, ' /= variable rank ',rank(value)
@@ -14,13 +14,13 @@ endif
 
 select type (value)
 type is (real(real64))
-  ier = nf90_get_var(self%ncid, varid, value)
+  ier = nf90_get_var(self%file_id, varid, value)
 type is (real(real32))
-  ier = nf90_get_var(self%ncid, varid, value)
+  ier = nf90_get_var(self%file_id, varid, value)
 type is (integer(int64))
-  ier = nf90_get_var(self%ncid, varid, value)
+  ier = nf90_get_var(self%file_id, varid, value)
 type is (integer(int32))
-  ier = nf90_get_var(self%ncid, varid, value)
+  ier = nf90_get_var(self%file_id, varid, value)
 class default
   ier = NF90_EBADTYPE
 end select

src/utils.f90

@@ -26,18 +26,18 @@
 
 select case(laction)
 case('r')
-  ier = nf90_open(self%filename, NF90_NOWRITE, self%ncid)
+  ier = nf90_open(self%filename, NF90_NOWRITE, self%file_id)
 case('r+')
-  ier = nf90_open(self%filename, NF90_NETCDF4, self%ncid)
+  ier = nf90_open(self%filename, NF90_NETCDF4, self%file_id)
 case('rw', 'a')
   if(is_netcdf(filename)) then
     !! NF90_WRITE is necessary to be in true read/write mode
-    ier = nf90_open(self%filename, ior(NF90_WRITE, NF90_NETCDF4), self%ncid)
+    ier = nf90_open(self%filename, ior(NF90_WRITE, NF90_NETCDF4), self%file_id)
   else
-    ier = nf90_create(self%filename, ior(NF90_CLOBBER, NF90_NETCDF4), self%ncid)
+    ier = nf90_create(self%filename, ior(NF90_CLOBBER, NF90_NETCDF4), self%file_id)
   endif
 case('w')
-  ier = nf90_create(self%filename, ior(NF90_CLOBBER, NF90_NETCDF4), self%ncid)
+  ier = nf90_create(self%filename, ior(NF90_CLOBBER, NF90_NETCDF4), self%file_id)
 case default
   error stop 'nc4fortran: Unsupported action -> ' // laction
 end select
@@ -68,7 +68,7 @@
   return
 endif
 
-ier = nf90_close(self%ncid)
+ier = nf90_close(self%file_id)
 if (ier /= NF90_NOERR) error stop 'ERROR:close: ' // self%filename
 
 self%is_open = .false.
@@ -84,10 +84,10 @@
 module procedure is_contig
 integer :: ier, varid
 
-ier = nf90_inq_varid(self%ncid, dname, varid)
+ier = nf90_inq_varid(self%file_id, dname, varid)
 if (ier/=NF90_NOERR) error stop 'nc4fortran:is_contig: cannot find variable: ' // dname
 
-ier = nf90_inquire_variable(self%ncid, varid, contiguous=is_contig)
+ier = nf90_inquire_variable(self%file_id, varid, contiguous=is_contig)
 if (ier/=NF90_NOERR) error stop 'nc4fortran:is_contig: cannot get variable properties' // dname
 
 end procedure is_contig
@@ -96,27 +96,27 @@
 module procedure is_chunked
 integer :: ier, varid
 
-ier = nf90_inq_varid(self%ncid, dname, varid)
+ier = nf90_inq_varid(self%file_id, dname, varid)
 if (ier/=NF90_NOERR) error stop 'nc4fortran:is_chunked: cannot find variable: ' // dname
 
-ier = nf90_inquire_variable(self%ncid, varid, contiguous=is_chunked)
+ier = nf90_inquire_variable(self%file_id, varid, contiguous=is_chunked)
 if (ier/=NF90_NOERR) error stop 'nc4fortran:is_chunked: cannot get variable properties' // dname
 
 is_chunked = .not.is_chunked
 end procedure is_chunked
 
 
 module procedure is_netcdf
-integer :: ierr, ncid
+integer :: ierr, file_id
 
 inquire(file=filename, exist=is_netcdf)
 !! avoid warning/error messages
 if (.not. is_netcdf) return
 
-ierr = nf90_open(filename, NF90_NOWRITE, ncid)
+ierr = nf90_open(filename, NF90_NOWRITE, file_id)
 is_netcdf = ierr == 0
 
-ierr = nf90_close(ncid)
+ierr = nf90_close(file_id)
 
 end procedure is_netcdf
 
@@ -140,9 +140,9 @@
 case (NF90_EBADDIM)
   m = 'ERROR: ' // dname // ' invalid dimension ID or Name'
 case (NF90_EBADGRPID)
-  m = 'ERROR: ' // dname // ' bad group ID in ncid'
+  m = 'ERROR: ' // dname // ' bad group ID in file_id'
 case (NF90_EBADID)
-  m = 'ERROR: ' // dname // ' Bad group id or ncid invalid'
+  m = 'ERROR: ' // dname // ' Bad group id or file_id invalid'
 case (NF90_ENOTVAR)
   m = 'ERROR: ' // dname // ' variable not found'
 case (NF90_ENOTNC)