gdncfileopen.f90

Go to the documentation of this file.

 
 
subroutine gdncfileopen(fileid, filename, writable, overwrite, stat, err)
  use gtdata_netcdf_file_types, only: gd_nc_file_id_entry
  use gtdata_netcdf_file_internal, only: id_head, id_used
  use netcdf, only: &
    & nf90_write, &
    & nf90_nowrite, &
    & nf90_noerr, &
    & nf90_noclobber, &
    & nf90_clobber, &
    & nf90_64bit_offset, &
    & nf90_open, &
    & nf90_create
  use dc_message, only: messagenotify
  use dc_error, only: storeerror
  use dc_types, only: string
  use dc_trace, only: beginsub, endsub
  implicit none
  integer, intent(out):: fileid
  character(len = *), intent(in):: filename
  logical, intent(in), optional:: writable
  logical, intent(in), optional:: overwrite
  logical, intent(out), optional:: err
  integer, intent(out), optional:: stat
  logical:: writable_required
  logical:: overwrite_required
  type(gd_nc_file_id_entry), pointer:: identptr, prev
  integer:: mystat, mode
  character(len = 256):: real_filename
  character(len = STRING):: cause_c
  character(*), parameter:: subname = "GDNcFileOpen"
continue
  fileid = -1
  !
  ! オプションの解釈
  !
  writable_required = .false.
  overwrite_required = .false.
  if (present(writable)) writable_required = writable
  if (present(overwrite)) overwrite_required = overwrite
  call beginsub(subname, 'writable=%y overwrite=%y file=%c', &
    & l=(/writable_required, overwrite_required/), c1=trim(filename))
  !
  ! 同じ名前で書込み可能性も適合していれば nf90_open しないで済ませる
  !
  if (id_used) then
    identptr => id_head
    nullify(prev)
    do
      if ((identptr % filename == filename) &
        & .and. (identptr % writable .or. .not. writable_required)) then
        fileid = identptr % id
        identptr % count = identptr % count + 1
        if (present(err)) err = .false.
        if (present(stat)) stat = nf90_noerr
        mystat = nf90_noerr
        goto 999
      endif
      prev => identptr
      identptr => identptr % next
      if (.not. associated(identptr)) exit
    enddo
    allocate(identptr)
    prev%next => identptr
  else
    nullify(prev)
    allocate(id_head)
    identptr => id_head
    id_used = .true.
  endif
  nullify(identptr % next)
  identptr % filename = filename
  identptr % writable = writable_required
  identptr % count = 1
  !
  ! URL の部分的サポート
  !
  real_filename = filename
  if (real_filename(1:8) == 'file:///') then
    real_filename = real_filename(8: )
  else if (real_filename(1:5) == 'file:' .AND. real_filename(6:6) /= '/') then
    real_filename = real_filename(6: )
  endif
  !
  ! いざ nf90_open
  !
  mode = nf90_nowrite
  if (writable_required) mode = ior(mode, nf90_write)
  ! 既に nc ファイルがあると思って開けてみる
  mystat = nf90_open(real_filename, mode, identptr % id)
  !
  ! ファイルが既に存在する場合
  !
  if (mystat == nf90_noerr) then
    ! 書き込みモードの場合
    if (writable_required) then
      if (overwrite_required) then
        ! 上書きモードの場合
        mode = nf90_clobber
        call messagenotify('M', subname, &
          & '"%c" is overwritten.', c1=trim(filename), rank_mpi = -1)
      else
        ! 上書き禁止モードの場合
        mode = nf90_noclobber
        call messagenotify('W', subname, &
          & '"%c" is opened in write-protect mode.', c1=trim(filename), rank_mpi = -1)
      end if
      mode = ior(mode,nf90_64bit_offset)
      mystat = nf90_create(real_filename, mode, identptr % id)
      if (mystat /= nf90_noerr) then
        cause_c=filename
        if (present(stat)) stat = mystat
        goto 999
      end if
    endif
    ! 読み込みモードの場合は何もしない
  else
    !
    ! ファイルが無かった場合
    !
    if (.not. writable_required) then
      ! 読み込みモードの場合
      !
      ! 「無いよ」とエラーを吐いて終了
      if (mystat /= nf90_noerr) then
        cause_c=filename
        if (present(stat)) stat = mystat
        goto 999
      end if
    else
      ! 書き込みモードの場合
      mode = nf90_clobber
      ! ファイルを作成する
      mode = ior(mode,nf90_64bit_offset)
      mystat = nf90_create(real_filename, mode, identptr % id)
      if (mystat /= nf90_noerr) then
        cause_c=filename
        if (present(stat)) stat = mystat
        goto 999
      end if
    endif
  endif
 
  fileid = identptr % id
 
  ! 失敗したら消しておく
  if (mystat /= nf90_noerr) then
    if (associated(prev)) then
      prev%next => identptr % next
    else
      id_head => identptr % next
      if (.not. associated(id_head)) id_used = .false.
    endif
    deallocate(identptr)
    fileid = -1
  endif
 
  if (present(stat)) then
    stat = mystat
    if (present(err)) err = (stat /= nf90_noerr)
  else
    cause_c=filename
    goto 999
  endif
999 continue
  call storeerror(mystat, subname, err, cause_c)
  call endsub(subname, 'id=%d stat=%d', i=(/fileid, mystat/))
end subroutine gdncfileopen