gdncvarattrtrue.f90

Go to the documentation of this file.

 
 
subroutine gdncvargetattrlogical(var, name, value, default)
 
  use gtdata_netcdf_types, only: gd_nc_variable, gd_nc_variable_entry
  use gtdata_netcdf_internal, only: vtable_lookup
  use gtdata_netcdf_generic, only: get_attr, inquireplus
  use dc_types, only: string
  use netcdf, only: nf90_max_name, nf90_noerr, nf90_char, nf90_inquire_attribute, nf90_get_att
  use dc_error
  use dc_string
  implicit none
  type(gd_nc_variable), intent(in):: var
  character(len = *), intent(in):: name
  logical, intent(out):: value
  logical, intent(in), optional:: default
  type(gd_nc_variable_entry):: ent
  character(len = STRING):: cbuffer
  character(len = 7):: c_default
  character(len = NF90_MAX_NAME):: aname
  real, allocatable:: rbuf(:)
  integer:: stat, xtype, attrlen
  integer:: varid
  stat = vtable_lookup(var, ent)
  if (stat /= nf90_noerr) goto 999
  ! 大域属性サポート
  call inquireplus(var, name, &
    & varid=varid, nf_attrname=aname)
  stat = nf90_inquire_attribute(ent%fileid, varid, aname, xtype=xtype, len=attrlen)
  if (stat /= nf90_noerr) goto 999
  if (xtype == nf90_char) then
    c_default = "0"
    if (present(default)) then
      if (default) c_default = "1"
    endif
    call get_attr(var, name, cbuffer, c_default, stat)
    ! もうちょっとましな方法があるべきだが。
    select case(cbuffer)
    case("", "0", "0.0", "0.", ".0", "FALSE", "false", ".FALSE.", &
      & ".false.", "F", "f", "0.0D0", "0.0d0")
      value = .false.
    case default
      value = .true.
    end select
  else
    allocate(rbuf(attrlen))
    stat = nf90_get_att(ent%fileid, varid, aname, rbuf)
    if (stat /= nf90_noerr) goto 999
    value = (abs(rbuf(1)) > tiny(0.0))
  endif
  return
 
999 continue
  value = .false.
  if (present(default)) value = default
end subroutine gdncvargetattrlogical