mod_netcdf::netcdf_get_fvar Interface Reference

Public Member Functions
subroutine	netcdf_get_fvar_0dp (ng, model, ncname, myvarname, a, ncid, start, total, broadcast, min_val, max_val)
subroutine	netcdf_get_fvar_1dp (ng, model, ncname, myvarname, a, ncid, start, total, broadcast, min_val, max_val)
subroutine	netcdf_get_fvar_2dp (ng, model, ncname, myvarname, a, ncid, start, total, broadcast, min_val, max_val)
subroutine	netcdf_get_fvar_3dp (ng, model, ncname, myvarname, a, ncid, start, total, broadcast, min_val, max_val)
subroutine	netcdf_get_fvar_0d (ng, model, ncname, myvarname, a, ncid, start, total, broadcast, min_val, max_val)
subroutine	netcdf_get_fvar_1d (ng, model, ncname, myvarname, a, ncid, start, total, broadcast, min_val, max_val)
subroutine	netcdf_get_fvar_2d (ng, model, ncname, myvarname, a, ncid, start, total, broadcast, min_val, max_val)
subroutine	netcdf_get_fvar_3d (ng, model, ncname, myvarname, a, ncid, start, total, broadcast, min_val, max_val)
subroutine	netcdf_get_fvar_4d (ng, model, ncname, myvarname, a, ncid, start, total, broadcast, min_val, max_val)

Detailed Description

Definition at line 45 of file mod_netcdf.F.

Member Function/Subroutine Documentation

netcdf_get_fvar_0d()

subroutine mod_netcdf::netcdf_get_fvar::netcdf_get_fvar_0d	(	integer, intent(in)	ng,
		integer, intent(in)	model,
		character (len=*), intent(in)	ncname,
		character (len=*), intent(in)	myvarname,
		real(r8), intent(out)	a,
		integer, intent(in), optional	ncid,
		integer, dimension(:), intent(in), optional	start,
		integer, dimension(:), intent(in), optional	total,
		logical, intent(in), optional	broadcast,
		real(r8), intent(out), optional	min_val,
		real(r8), intent(out), optional	max_val )

Definition at line 3430 of file mod_netcdf.F.

!
!=======================================================================
!                                                                      !
!  This routine reads requested floating-point scalar variable from    !
!  specified NetCDF file.                                              !
!                                                                      !
!  On Input:                                                           !
!                                                                      !
!     ng           Nested grid number (integer)                        !
!     model        Calling model identifier (integer)                  !
!     ncname       NetCDF file name (string)                           !
!     myVarName    Variable name (string)                              !
!     ncid         NetCDF file ID (integer, OPTIONAL)                  !
!     start        Starting index where the first of the data values   !
!                    will be read along each dimension (integer,       !
!                    OPTIONAL)                                         !
!     total        Number of data values to be read along each         !
!                    dimension (integer, OPTIONAL)                     !
!     broadcast    Switch to broadcast read values from root to all    !
!                    members of the communicator in distributed-       !
!                    memory applications (logical, OPTIONAL,           !
!                    default=TRUE)                                     !
!                                                                      !
!  On Ouput:                                                           !
!                                                                      !
!     A            Read scalar variable (real)                         !
!     min_val      Read data minimum value (real, OPTIONAL)            !
!     max_val      Read data maximum value (real, OPTIONAL)            !
!                                                                      !
!  Examples:                                                           !
!                                                                      !
!    CALL netcdf_get_fvar (ng, iNLM, 'file.nc', 'VarName', fvar)       !
!    CALL netcdf_get_fvar (ng, iNLM, 'file.nc', 'VarName', fvar(1))    !
!                                                                      !
!=======================================================================
!
!  Imported variable declarations.
!
      logical, intent(in), optional :: broadcast
!
      integer, intent(in) :: ng, model
 
      integer, intent(in), optional :: ncid
      integer, intent(in), optional :: start(:)
      integer, intent(in), optional :: total(:)
!
      character (len=*), intent(in) :: ncname
      character (len=*), intent(in) :: myVarName
!
      real(r8), intent(out), optional :: min_val
      real(r8), intent(out), optional :: max_val
 
      real(r8), intent(out) :: A
!
!  Local variable declarations.
!
#if !defined PARALLEL_IO && defined DISTRIBUTE
      logical :: DoBroadcast
!
#endif
 
      integer :: my_ncid, status, varid
 
#if !defined PARALLEL_IO && defined DISTRIBUTE
      integer, dimension(2) :: ibuffer
#endif
!
      real(r8), dimension(1) :: my_A
!
      character (len=*), parameter :: MyFile =                          &
     &  __FILE__//", netcdf_get_fvar_0d"
!
!-----------------------------------------------------------------------
!  Read in a floating-point scalar variable.
!-----------------------------------------------------------------------
!
!  If NetCDF file ID is not provided, open NetCDF for reading.
!
      IF (.not.PRESENT(ncid)) THEN
        CALL netcdf_open (ng, model, trim(ncname), 0, my_ncid)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
      ELSE
        my_ncid=ncid
      END IF
 
#if !defined PARALLEL_IO && defined DISTRIBUTE
!
!  Determine if the read data is only needed and allocated by the master
!  node ins distribute-memory applications.
!
      IF (.not.PRESENT(broadcast)) THEN
        dobroadcast=.true.
      ELSE
        dobroadcast=broadcast
      END IF
#endif
!
!  Read in variable.
!
      IF (inpthread) THEN
        status=nf90_inq_varid(my_ncid, trim(myvarname), varid)
        IF (status.eq.nf90_noerr) THEN
          IF (PRESENT(start).and.PRESENT(total)) THEN
            status=nf90_get_var(my_ncid, varid, my_a, start, total)
            a=my_a(1)
          ELSE
            status=nf90_get_var(my_ncid, varid, a)
          END IF
          IF (founderror(status, nf90_noerr, __line__, myfile)) THEN
            WRITE (stdout,10) trim(myvarname), trim(ncname),            &
     &                        trim(sourcefile), nf90_strerror(status)
            exit_flag=2
            ioerror=status
          END IF
        ELSE
          WRITE (stdout,20) trim(myvarname), trim(ncname),              &
     &                      trim(sourcefile), nf90_strerror(status)
          exit_flag=2
          ioerror=status
        END IF
      END IF
 
#if !defined PARALLEL_IO && defined DISTRIBUTE
!
!  Broadcast read variable to all processors in the group.
!
      ibuffer(1)=exit_flag
      ibuffer(2)=ioerror
      CALL mp_bcasti (ng, model, ibuffer)
      exit_flag=ibuffer(1)
      ioerror=ibuffer(2)
      IF (dobroadcast.and.(exit_flag.eq.noerror)) THEN
        CALL mp_bcastf (ng, model, a)
      END IF
#endif
!
!  Compute minimum and maximum values of read variable. Notice that
!  the same read value is assigned since a scalar variable was
!  processed.
!
      IF (PRESENT(min_val)) THEN
        min_val=a
      END IF
      IF (PRESENT(max_val)) THEN
        max_val=a
      END IF
!
!  If NetCDF file ID is not provided, close input NetCDF file.
!
      IF (.not.PRESENT(ncid)) THEN
        CALL netcdf_close (ng, model, my_ncid, ncname, .false.)
      END IF
!
  10  FORMAT (/,' NETCDF_GET_FVAR_0D - error while reading variable:',  &
     &        2x,a,/,22x,'in input file:',2x,a,/,22x,'call from:',2x,a, &
     &        /,22x,a)
  20  FORMAT (/,' NETCDF_GET_FVAR_0D - error while inquiring ID for ',  &
     &        'variable:',2x,a,/,22x,'in input file:',2x,a,/,22x,       &
     &        'call from:',2x,a,/,22x,a)
!
      RETURN

References mod_scalars::exit_flag, mod_parallel::inpthread, mod_iounits::ioerror, mod_netcdf::netcdf_close(), mod_netcdf::netcdf_open(), mod_scalars::noerror, mod_iounits::sourcefile, and mod_iounits::stdout.

Here is the call graph for this function:

netcdf_get_fvar_0dp()

subroutine mod_netcdf::netcdf_get_fvar::netcdf_get_fvar_0dp	(	integer, intent(in)	ng,
		integer, intent(in)	model,
		character (len=*), intent(in)	ncname,
		character (len=*), intent(in)	myvarname,
		real(dp), intent(out)	a,
		integer, intent(in), optional	ncid,
		integer, dimension(:), intent(in), optional	start,
		integer, dimension(:), intent(in), optional	total,
		logical, intent(in), optional	broadcast,
		real(dp), intent(out), optional	min_val,
		real(dp), intent(out), optional	max_val )

Definition at line 2521 of file mod_netcdf.F.

!
!=======================================================================
!                                                                      !
!  This routine reads requested double-precision scalar variable from  !
!  specified NetCDF file.                                              !
!                                                                      !
!  On Input:                                                           !
!                                                                      !
!     ng           Nested grid number (integer)                        !
!     model        Calling model identifier (integer)                  !
!     ncname       NetCDF file name (string)                           !
!     myVarName    Variable name (string)                              !
!     ncid         NetCDF file ID (integer, OPTIONAL)                  !
!     start        Starting index where the first of the data values   !
!                    will be read along each dimension (integer,       !
!                    OPTIONAL)                                         !
!     total        Number of data values to be read along each         !
!                    dimension (integer, OPTIONAL)                     !
!     broadcast    Switch to broadcast read values from root to all    !
!                    members of the communicator in distributed-       !
!                    memory applications (logical, OPTIONAL,           !
!                    default=TRUE)                                     !
!                                                                      !
!  On Ouput:                                                           !
!                                                                      !
!     A            Read scalar variable (double precision)             !
!     min_val      Read data minimum value (double precision, OPTIONAL)!
!     max_val      Read data maximum value (double precision, OPTIONAL)!
!                                                                      !
!  Examples:                                                           !
!                                                                      !
!    CALL netcdf_get_fvar (ng, iNLM, 'file.nc', 'VarName', fvar)       !
!    CALL netcdf_get_fvar (ng, iNLM, 'file.nc', 'VarName', fvar(1))    !
!                                                                      !
!=======================================================================
!
!  Imported variable declarations.
!
      logical, intent(in), optional :: broadcast
!
      integer, intent(in) :: ng, model
 
      integer, intent(in), optional :: ncid
      integer, intent(in), optional :: start(:)
      integer, intent(in), optional :: total(:)
!
      character (len=*), intent(in) :: ncname
      character (len=*), intent(in) :: myVarName
!
      real(dp), intent(out), optional :: min_val
      real(dp), intent(out), optional :: max_val
 
      real(dp), intent(out) :: A
!
!  Local variable declarations.
!
# if !defined PARALLEL_IO && defined DISTRIBUTE
      logical :: DoBroadcast
!
# endif
 
      integer :: my_ncid, status, varid
 
# if !defined PARALLEL_IO && defined DISTRIBUTE
      integer, dimension(2) :: ibuffer
# endif
!
      real(dp), dimension(1) :: my_A
!
      character (len=*), parameter :: MyFile =                          &
     &  __FILE__//", netcdf_get_fvar_0dp"
!
!-----------------------------------------------------------------------
!  Read in a double-precision scalar variable.
!-----------------------------------------------------------------------
!
!  If NetCDF file ID is not provided, open NetCDF for reading.
!
      IF (.not.PRESENT(ncid)) THEN
        CALL netcdf_open (ng, model, trim(ncname), 0, my_ncid)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
      ELSE
        my_ncid=ncid
      END IF
 
# if !defined PARALLEL_IO && defined DISTRIBUTE
!
!  Determine if the read data is only needed and allocated by the master
!  node ins distribute-memory applications.
!
      IF (.not.PRESENT(broadcast)) THEN
        dobroadcast=.true.
      ELSE
        dobroadcast=broadcast
      END IF
# endif
!
!  Read in variable.
!
      IF (inpthread) THEN
        status=nf90_inq_varid(my_ncid, trim(myvarname), varid)
        IF (status.eq.nf90_noerr) THEN
          IF (PRESENT(start).and.PRESENT(total)) THEN
            status=nf90_get_var(my_ncid, varid, my_a, start, total)
            a=my_a(1)
          ELSE
            status=nf90_get_var(my_ncid, varid, a)
          END IF
          IF (founderror(status, nf90_noerr, __line__, myfile)) THEN
            WRITE (stdout,10) trim(myvarname), trim(ncname),            &
     &                        trim(sourcefile), nf90_strerror(status)
            exit_flag=2
            ioerror=status
          END IF
        ELSE
          WRITE (stdout,20) trim(myvarname), trim(ncname),              &
     &                      trim(sourcefile), nf90_strerror(status)
          exit_flag=2
          ioerror=status
        END IF
      END IF
 
# if !defined PARALLEL_IO && defined DISTRIBUTE
!
!  Broadcast read variable to all processors in the group.
!
      ibuffer(1)=exit_flag
      ibuffer(2)=ioerror
      CALL mp_bcasti (ng, model, ibuffer)
      exit_flag=ibuffer(1)
      ioerror=ibuffer(2)
      IF (dobroadcast.and.(exit_flag.eq.noerror)) THEN
        CALL mp_bcastf (ng, model, a)
      END IF
# endif
!
!  Compute minimum and maximum values of read variable. Notice that
!  the same read value is assigned since a scalar variable was
!  processed.
!
      IF (PRESENT(min_val)) THEN
        min_val=a
      END IF
      IF (PRESENT(max_val)) THEN
        max_val=a
      END IF
!
!  If NetCDF file ID is not provided, close input NetCDF file.
!
      IF (.not.PRESENT(ncid)) THEN
        CALL netcdf_close (ng, model, my_ncid, ncname, .false.)
      END IF
!
  10  FORMAT (/,' NETCDF_GET_FVAR_0DP - error while reading variable:', &
     &        2x,a,/,23x,'in input file:',2x,a,/,23x,'call from:',2x,a, &
     &        /,23x,a)
  20  FORMAT (/,' NETCDF_GET_FVAR_0DP - error while inquiring ID for ', &
     &        'variable:',2x,a,/,23x,'in input file:',2x,a,/,23x,       &
     &        'call from:',2x,a,/,23x,a)
!
      RETURN

References mod_scalars::exit_flag, mod_parallel::inpthread, mod_iounits::ioerror, mod_netcdf::netcdf_close(), mod_netcdf::netcdf_open(), mod_scalars::noerror, mod_iounits::sourcefile, and mod_iounits::stdout.

Here is the call graph for this function:

netcdf_get_fvar_1d()

subroutine mod_netcdf::netcdf_get_fvar::netcdf_get_fvar_1d	(	integer, intent(in)	ng,
		integer, intent(in)	model,
		character (len=*), intent(in)	ncname,
		character (len=*), intent(in)	myvarname,
		real(r8), dimension(:), intent(out)	a,
		integer, intent(in), optional	ncid,
		integer, dimension(:), intent(in), optional	start,
		integer, dimension(:), intent(in), optional	total,
		logical, intent(in), optional	broadcast,
		real(r8), intent(out), optional	min_val,
		real(r8), intent(out), optional	max_val )

Definition at line 3596 of file mod_netcdf.F.

!
!=======================================================================
!                                                                      !
!  This routine reads requested floating-point 1D-array variable from  !
!  specified NetCDF file.                                              !
!                                                                      !
!  On Input:                                                           !
!                                                                      !
!     ng           Nested grid number (integer)                        !
!     model        Calling model identifier (integer)                  !
!     ncname       NetCDF file name (string)                           !
!     myVarName    Variable name (string)                              !
!     ncid         NetCDF file ID (integer, OPTIONAL)                  !
!     start        Starting index where the first of the data values   !
!                    will be read along each dimension (integer,       !
!                    OPTIONAL)                                         !
!     total        Number of data values to be read along each         !
!                    dimension (integer, OPTIONAL)                     !
!     broadcast    Switch to broadcast read values from root to all    !
!                    members of the communicator in distributed-       !
!                    memory applications (logical, OPTIONAL,           !
!                    default=TRUE)                                     !
!                                                                      !
!  On Ouput:                                                           !
!                                                                      !
!     A            Read 1D-array variable (real)                       !
!     min_val      Read data minimum value (real, OPTIONAL)            !
!     max_val      Read data maximum value (real, OPTIONAL)            !
!                                                                      !
!  Examples:                                                           !
!                                                                      !
!    CALL netcdf_get_fvar (ng, iNLM, 'file.nc', 'VarName', fvar)       !
!    CALL netcdf_get_fvar (ng, iNLM, 'file.nc', 'VarName', fvar(0:))   !
!    CALL netcdf_get_fvar (ng, iNLM, 'file.nc', 'VarName', fvar(:,1))  !
!                                                                      !
!=======================================================================
!
!  Imported variable declarations.
!
      logical, intent(in), optional :: broadcast
!
      integer, intent(in) :: ng, model
 
      integer, intent(in), optional :: ncid
      integer, intent(in), optional :: start(:)
      integer, intent(in), optional :: total(:)
!
      character (len=*), intent(in) :: ncname
      character (len=*), intent(in) :: myVarName
!
      real(r8), intent(out), optional :: min_val
      real(r8), intent(out), optional :: max_val
 
      real(r8), intent(out) :: A(:)
!
!  Local variable declarations.
!
#if !defined PARALLEL_IO && defined DISTRIBUTE
      logical :: DoBroadcast
#endif
      logical, dimension(3) :: foundit
!
      integer :: i, my_ncid, status, varid
 
      integer, dimension(1) :: Asize
#if !defined PARALLEL_IO && defined DISTRIBUTE
      integer, dimension(2) :: ibuffer
#endif
!
      real(r8) :: Afactor, Aoffset, Aspval
 
      real(r8), parameter :: Aepsilon = 1.0e-8_r8
 
      real(r8), dimension(3) :: AttValue
!
      character (len=12), dimension(3) :: AttName
 
      character (len=*), parameter :: MyFile =                          &
     &  __FILE__//", netcdf_get_fvar_1d"
!
!-----------------------------------------------------------------------
!  Read in a floating-point 1D-array variable.
!-----------------------------------------------------------------------
!
      IF (PRESENT(start).and.PRESENT(total)) THEN
        asize(1)=1
        DO i=1,SIZE(total)              ! this logic is for the case
          asize(1)=asize(1)*total(i)    ! of reading multidimensional
        END DO                          ! data into a compact 1D array
      ELSE
        asize(1)=ubound(a, dim=1)
      END IF
!
!  If NetCDF file ID is not provided, open NetCDF for reading.
!
      IF (.not.PRESENT(ncid)) THEN
        CALL netcdf_open (ng, model, trim(ncname), 0, my_ncid)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
      ELSE
        my_ncid=ncid
      END IF
 
#if !defined PARALLEL_IO && defined DISTRIBUTE
!
!  Determine if the read data is only needed and allocated by the master
!  node ins distribute-memory applications.
!
      IF (.not.PRESENT(broadcast)) THEN
        dobroadcast=.true.
      ELSE
        dobroadcast=broadcast
      END IF
#endif
!
!  Read in variable.
!
      IF (inpthread) THEN
        status=nf90_inq_varid(my_ncid, trim(myvarname), varid)
        IF (status.eq.nf90_noerr) THEN
          IF (PRESENT(start).and.PRESENT(total)) THEN
            status=nf90_get_var(my_ncid, varid, a, start, total)
          ELSE
            status=nf90_get_var(my_ncid, varid, a)
          END IF
          IF (founderror(status, nf90_noerr, __line__, myfile)) THEN
            WRITE (stdout,10) trim(myvarname), trim(ncname),            &
     &                        trim(sourcefile), nf90_strerror(status)
            exit_flag=2
            ioerror=status
          END IF
        ELSE
          WRITE (stdout,20) trim(myvarname), trim(ncname),              &
     &                      trim(sourcefile), nf90_strerror(status)
          exit_flag=2
          ioerror=status
        END IF
      END IF
 
#if !defined PARALLEL_IO && defined DISTRIBUTE
!
!  Broadcast read variable to all processors in the group.
!
      ibuffer(1)=exit_flag
      ibuffer(2)=ioerror
      CALL mp_bcasti (ng, model, ibuffer)
      exit_flag=ibuffer(1)
      ioerror=ibuffer(2)
      IF (dobroadcast.and.(exit_flag.eq.noerror)) THEN
        CALL mp_bcastf (ng, model, a)
      END IF
#endif
!
!  Check if the following attributes: "scale_factor", "add_offset", and
!  "_FillValue" are present in the input NetCDF variable:
!
!  If the "scale_value" attribute is present, the data is multiplied by
!  this factor after reading.
!  If the "add_offset" attribute is present, this value is added to the
!  data after reading.
!  If both "scale_factor" and "add_offset" attributes are present, the
!  data are first scaled before the offset is added.
!  If the "_FillValue" attribute is present, the data having this value
!  is treated as missing and it is replaced with zero. This feature it
!  is usually related with the land/sea masking.
!
      attname(1)='scale_factor'
      attname(2)='add_offset  '
      attname(3)='_FillValue  '
 
      CALL netcdf_get_fatt (ng, model, ncname, varid, attname,          &
     &                      attvalue, foundit,                          &
     &                      ncid = my_ncid)
 
      IF (exit_flag.eq.noerror) THEN
        IF (.not.foundit(1)) THEN
          afactor=1.0_r8
        ELSE
          afactor=attvalue(1)
        END IF
 
        IF (.not.foundit(2)) THEN
          aoffset=0.0_r8
        ELSE
          aoffset=attvalue(2)
        END IF
 
        IF (.not.foundit(3)) THEN
          aspval=spval_check
        ELSE
          aspval=attvalue(3)
        END IF
 
        DO i=1,asize(1)                       ! zero out missing values
          IF ((foundit(3).and.(abs(a(i)-aspval).lt.aepsilon)).or.       &
     &        (.not.foundit(3).and.(abs(a(i)).ge.abs(aspval)))) THEN
            a(i)=0.0_r8
          END IF
        END DO
 
        IF (foundit(1)) THEN                  ! scale data
          DO i=1,asize(1)
            a(i)=afactor*a(i)
          END DO
        END IF
 
        IF (foundit(2)) THEN                  ! add data offset
          DO i=1,asize(1)
            a(i)=a(i)+aoffset
          END DO
        END IF
      END IF
!
!  Compute minimum and maximum values of read variable.
!
      IF (PRESENT(min_val)) THEN
        min_val=minval(a)
      END IF
      IF (PRESENT(max_val)) THEN
        max_val=maxval(a)
      END IF
!
!  If NetCDF file ID is not provided, close input NetCDF file.
!
      IF (.not.PRESENT(ncid)) THEN
        CALL netcdf_close (ng, model, my_ncid, ncname, .false.)
      END IF
!
  10  FORMAT (/,' NETCDF_GET_FVAR_1D - error while reading variable:',  &
     &        2x,a,/,22x,'in input file:',2x,a,/,22x,'call from:',2x,a, &
     &        /,22x,a)
  20  FORMAT (/,' NETCDF_GET_FVAR_1D - error while inquiring ID for ',  &
     &        'variable:',2x,a,/,22x,'in input file:',2x,a,/,22x,       &
     &        'call from:',2x,a,/,22x,a)
!
      RETURN

References mod_scalars::exit_flag, mod_parallel::inpthread, mod_iounits::ioerror, mod_netcdf::netcdf_close(), mod_netcdf::netcdf_open(), mod_scalars::noerror, mod_iounits::sourcefile, mod_scalars::spval_check, and mod_iounits::stdout.

Here is the call graph for this function:

netcdf_get_fvar_1dp()

subroutine mod_netcdf::netcdf_get_fvar::netcdf_get_fvar_1dp	(	integer, intent(in)	ng,
		integer, intent(in)	model,
		character (len=*), intent(in)	ncname,
		character (len=*), intent(in)	myvarname,
		real(dp), dimension(:), intent(out)	a,
		integer, intent(in), optional	ncid,
		integer, dimension(:), intent(in), optional	start,
		integer, dimension(:), intent(in), optional	total,
		logical, intent(in), optional	broadcast,
		real(dp), intent(out), optional	min_val,
		real(dp), intent(out), optional	max_val )

Definition at line 2687 of file mod_netcdf.F.

!
!=======================================================================
!                                                                      !
!  This routine reads requested double-precision 1D-array variable     !
!  froms specified NetCDF file.                                        !
!                                                                      !
!  On Input:                                                           !
!                                                                      !
!     ng           Nested grid number (integer)                        !
!     model        Calling model identifier (integer)                  !
!     ncname       NetCDF file name (string)                           !
!     myVarName    Variable name (string)                              !
!     ncid         NetCDF file ID (integer, OPTIONAL)                  !
!     start        Starting index where the first of the data values   !
!                    will be read along each dimension (integer,       !
!                    OPTIONAL)                                         !
!     total        Number of data values to be read along each         !
!                    dimension (integer, OPTIONAL)                     !
!     broadcast    Switch to broadcast read values from root to all    !
!                    members of the communicator in distributed-       !
!                    memory applications (logical, OPTIONAL,           !
!                    default=TRUE)                                     !
!                                                                      !
!  On Ouput:                                                           !
!                                                                      !
!     A            Read 1D-array variable (double precision)           !
!     min_val      Read data minimum value (double precision, OPTIONAL)!
!     max_val      Read data maximum value (double precision, OPTIONAL)!
!                                                                      !
!  Examples:                                                           !
!                                                                      !
!    CALL netcdf_get_fvar (ng, iNLM, 'file.nc', 'VarName', fvar)       !
!    CALL netcdf_get_fvar (ng, iNLM, 'file.nc', 'VarName', fvar(0:))   !
!    CALL netcdf_get_fvar (ng, iNLM, 'file.nc', 'VarName', fvar(:,1))  !
!                                                                      !
!=======================================================================
!
!  Imported variable declarations.
!
      logical, intent(in), optional :: broadcast
!
      integer, intent(in) :: ng, model
 
      integer, intent(in), optional :: ncid
      integer, intent(in), optional :: start(:)
      integer, intent(in), optional :: total(:)
!
      character (len=*), intent(in) :: ncname
      character (len=*), intent(in) :: myVarName
!
      real(dp), intent(out), optional :: min_val
      real(dp), intent(out), optional :: max_val
 
      real(dp), intent(out) :: A(:)
!
!  Local variable declarations.
!
# if !defined PARALLEL_IO && defined DISTRIBUTE
      logical :: DoBroadcast
# endif
      logical, dimension(3) :: foundit
!
      integer :: i, my_ncid, status, varid
 
      integer, dimension(1) :: Asize
# if !defined PARALLEL_IO && defined DISTRIBUTE
      integer, dimension(2) :: ibuffer
# endif
!
      real(dp) :: Afactor, Aoffset, Aspval
 
      real(dp), parameter :: Aepsilon = 1.0e-8_dp
 
      real(dp), dimension(3) :: AttValue
!
      character (len=12), dimension(3) :: AttName
 
      character (len=*), parameter :: MyFile =                          &
     &  __FILE__//", netcdf_get_fvar_1dp"
!
!-----------------------------------------------------------------------
!  Read in a double-precision 1D-array variable.
!-----------------------------------------------------------------------
!
      IF (PRESENT(start).and.PRESENT(total)) THEN
        asize(1)=1
        DO i=1,SIZE(total)              ! this logic is for the case
          asize(1)=asize(1)*total(i)    ! of reading multidimensional
        END DO                          ! data into a compact 1D array
      ELSE
        asize(1)=ubound(a, dim=1)
      END IF
!
!  If NetCDF file ID is not provided, open NetCDF for reading.
!
      IF (.not.PRESENT(ncid)) THEN
        CALL netcdf_open (ng, model, trim(ncname), 0, my_ncid)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
      ELSE
        my_ncid=ncid
      END IF
 
# if !defined PARALLEL_IO && defined DISTRIBUTE
!
!  Determine if the read data is only needed and allocated by the master
!  node ins distribute-memory applications.
!
      IF (.not.PRESENT(broadcast)) THEN
        dobroadcast=.true.
      ELSE
        dobroadcast=broadcast
      END IF
# endif
!
!  Read in variable.
!
      IF (inpthread) THEN
        status=nf90_inq_varid(my_ncid, trim(myvarname), varid)
        IF (status.eq.nf90_noerr) THEN
          IF (PRESENT(start).and.PRESENT(total)) THEN
            status=nf90_get_var(my_ncid, varid, a, start, total)
          ELSE
            status=nf90_get_var(my_ncid, varid, a)
          END IF
          IF (founderror(status, nf90_noerr, __line__, myfile)) THEN
            WRITE (stdout,10) trim(myvarname), trim(ncname),            &
     &                        trim(sourcefile), nf90_strerror(status)
            exit_flag=2
            ioerror=status
          END IF
        ELSE
          WRITE (stdout,20) trim(myvarname), trim(ncname),              &
     &                      trim(sourcefile), nf90_strerror(status)
          exit_flag=2
          ioerror=status
        END IF
      END IF
 
# if !defined PARALLEL_IO && defined DISTRIBUTE
!
!  Broadcast read variable to all processors in the group.
!
      ibuffer(1)=exit_flag
      ibuffer(2)=ioerror
      CALL mp_bcasti (ng, model, ibuffer)
      exit_flag=ibuffer(1)
      ioerror=ibuffer(2)
      IF (dobroadcast.and.(exit_flag.eq.noerror)) THEN
        CALL mp_bcastf (ng, model, a)
      END IF
# endif
!
!  Check if the following attributes: "scale_factor", "add_offset", and
!  "_FillValue" are present in the input NetCDF variable:
!
!  If the "scale_value" attribute is present, the data is multiplied by
!  this factor after reading.
!  If the "add_offset" attribute is present, this value is added to the
!  data after reading.
!  If both "scale_factor" and "add_offset" attributes are present, the
!  data are first scaled before the offset is added.
!  If the "_FillValue" attribute is present, the data having this value
!  is treated as missing and it is replaced with zero. This feature it
!  is usually related with the land/sea masking.
!
      attname(1)='scale_factor'
      attname(2)='add_offset  '
      attname(3)='_FillValue  '
 
      CALL netcdf_get_fatt (ng, model, ncname, varid, attname,          &
     &                      attvalue, foundit,                          &
     &                      ncid = my_ncid)
 
      IF (exit_flag.eq.noerror) THEN
        IF (.not.foundit(1)) THEN
          afactor=1.0_dp
        ELSE
          afactor=attvalue(1)
        END IF
 
        IF (.not.foundit(2)) THEN
          aoffset=0.0_dp
        ELSE
          aoffset=attvalue(2)
        END IF
 
        IF (.not.foundit(3)) THEN
          aspval=spval_check
        ELSE
          aspval=attvalue(3)
        END IF
 
        DO i=1,asize(1)                       ! zero out missing values
          IF ((foundit(3).and.(abs(a(i)-aspval).lt.aepsilon)).or.       &
     &        (.not.foundit(3).and.(abs(a(i)).ge.abs(aspval)))) THEN
            a(i)=0.0_dp
          END IF
        END DO
 
        IF (foundit(1)) THEN                  ! scale data
          DO i=1,asize(1)
            a(i)=afactor*a(i)
          END DO
        END IF
 
        IF (foundit(2)) THEN                  ! add data offset
          DO i=1,asize(1)
            a(i)=a(i)+aoffset
          END DO
        END IF
      END IF
!
!  Compute minimum and maximum values of read variable.
!
      IF (PRESENT(min_val)) THEN
        min_val=minval(a)
      END IF
      IF (PRESENT(max_val)) THEN
        max_val=maxval(a)
      END IF
!
!  If NetCDF file ID is not provided, close input NetCDF file.
!
      IF (.not.PRESENT(ncid)) THEN
        CALL netcdf_close (ng, model, my_ncid, ncname, .false.)
      END IF
!
  10  FORMAT (/,' NETCDF_GET_FVAR_1DP - error while reading variable:', &
     &        2x,a,/,23x,'in input file:',2x,a,/,23x,'call from:',2x,a, &
     &        /,23x,a)
  20  FORMAT (/,' NETCDF_GET_FVAR_1DP - error while inquiring ID for ', &
     &        'variable:',2x,a,/,23x,'in input file:',2x,a,/,23x,       &
     &        'call from:',2x,a,/,23x,a)
!
      RETURN

References mod_scalars::exit_flag, mod_parallel::inpthread, mod_iounits::ioerror, mod_netcdf::netcdf_close(), mod_netcdf::netcdf_open(), mod_scalars::noerror, mod_iounits::sourcefile, mod_scalars::spval_check, and mod_iounits::stdout.

Here is the call graph for this function:

netcdf_get_fvar_2d()

subroutine mod_netcdf::netcdf_get_fvar::netcdf_get_fvar_2d	(	integer, intent(in)	ng,
		integer, intent(in)	model,
		character (len=*), intent(in)	ncname,
		character (len=*), intent(in)	myvarname,
		real(r8), dimension(:,:), intent(out)	a,
		integer, intent(in), optional	ncid,
		integer, dimension(:), intent(in), optional	start,
		integer, dimension(:), intent(in), optional	total,
		logical, intent(in), optional	broadcast,
		real(r8), intent(out), optional	min_val,
		real(r8), intent(out), optional	max_val )

Definition at line 3836 of file mod_netcdf.F.

!
!=======================================================================
!                                                                      !
!  This routine reads requested floating-point 2D-array variable from  !
!  specified NetCDF file.                                              !
!                                                                      !
!  On Input:                                                           !
!                                                                      !
!     ng           Nested grid number (integer)                        !
!     model        Calling model identifier (integer)                  !
!     ncname       NetCDF file name (string)                           !
!     myVarName    Variable name (string)                              !
!     ncid         NetCDF file ID (integer, OPTIONAL)                  !
!     start        Starting index where the first of the data values   !
!                    will be read along each dimension (integer,       !
!                    OPTIONAL)                                         !
!     total        Number of data values to be read along each         !
!                    dimension (integer, OPTIONAL)                     !
!     broadcast    Switch to broadcast read values from root to all    !
!                    members of the communicator in distributed-       !
!                    memory applications (logical, OPTIONAL,           !
!                    default=TRUE)                                     !
!                                                                      !
!  On Ouput:                                                           !
!                                                                      !
!     A            Read 2D-array variable (real)                       !
!     min_val      Read data minimum value (real, OPTIONAL)            !
!     max_val      Read data maximum value (real, OPTIONAL)            !
!                                                                      !
!  Examples:                                                           !
!                                                                      !
!    CALL netcdf_get_fvar (ng, iNLM, 'file.nc', 'VarName', fvar)       !
!    CALL netcdf_get_fvar (ng, iNLM, 'file.nc', 'VarName', fvar(0:,:)) !
!    CALL netcdf_get_fvar (ng, iNLM, 'file.nc', 'VarName', fvar(0:,0:))!
!    CALL netcdf_get_fvar (ng, iNLM, 'file.nc', 'VarName', fvar(:,:,1))!
!                                                                      !
!=======================================================================
!
!  Imported variable declarations.
!
      logical, intent(in), optional :: broadcast
!
      integer, intent(in) :: ng, model
 
      integer, intent(in), optional :: ncid
      integer, intent(in), optional :: start(:)
      integer, intent(in), optional :: total(:)
!
      character (len=*), intent(in) :: ncname
      character (len=*), intent(in) :: myVarName
!
      real(r8), intent(out), optional :: min_val
      real(r8), intent(out), optional :: max_val
 
      real(r8), intent(out) :: A(:,:)
!
!  Local variable declarations.
!
#if !defined PARALLEL_IO && defined DISTRIBUTE
      logical :: Dobroadcast
#endif
      logical, dimension(3) :: foundit
!
      integer :: i, j, my_ncid, status, varid
 
      integer, dimension(2) :: Asize
 
#if !defined PARALLEL_IO && defined DISTRIBUTE
      integer, dimension(2) :: ibuffer
#endif
!
      real(r8) :: Afactor, Aoffset, Aspval
 
      real(r8), parameter :: Aepsilon = 1.0e-8_r8
 
      real(r8), dimension(3) :: AttValue
!
      character (len=12), dimension(3) :: AttName
 
      character (len=*), parameter :: MyFile =                          &
     &  __FILE__//", netcdf_get_fvar_2d"
!
!-----------------------------------------------------------------------
!  Read in a floating-point 2D-array variable.
!-----------------------------------------------------------------------
!
      IF (PRESENT(start).and.PRESENT(total)) THEN
        asize(1)=total(1)
        asize(2)=total(2)
      ELSE
        asize(1)=ubound(a, dim=1)
        asize(2)=ubound(a, dim=2)
      END IF
!
!  If NetCDF file ID is not provided, open NetCDF for reading.
!
      IF (.not.PRESENT(ncid)) THEN
        CALL netcdf_open (ng, model, trim(ncname), 0, my_ncid)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
      ELSE
        my_ncid=ncid
      END IF
 
#if !defined PARALLEL_IO && defined DISTRIBUTE
!
!  Determine if the read data is only needed and allocated by the master
!  node ins distribute-memory applications.
!
      IF (.not.PRESENT(broadcast)) THEN
        dobroadcast=.true.
      ELSE
        dobroadcast=broadcast
      END IF
#endif
!
!  Read in variable.
!
      IF (inpthread) THEN
        status=nf90_inq_varid(my_ncid, trim(myvarname), varid)
        IF (status.eq.nf90_noerr) THEN
          IF (PRESENT(start).and.PRESENT(total)) THEN
            status=nf90_get_var(my_ncid, varid, a, start, total)
          ELSE
            status=nf90_get_var(my_ncid, varid, a)
          END IF
          IF (founderror(status, nf90_noerr, __line__, myfile)) THEN
            WRITE (stdout,10) trim(myvarname), trim(ncname),            &
     &                        trim(sourcefile), nf90_strerror(status)
            exit_flag=2
            ioerror=status
          END IF
        ELSE
          WRITE (stdout,20) trim(myvarname), trim(ncname),              &
     &                      trim(sourcefile), nf90_strerror(status),    &
     &                      nf90_strerror(status)
          exit_flag=2
          ioerror=status
        END IF
      END IF
 
#if !defined PARALLEL_IO && defined DISTRIBUTE
!
!  Broadcast read variable to all processors in the group.
!
      ibuffer(1)=exit_flag
      ibuffer(2)=ioerror
      CALL mp_bcasti (ng, model, ibuffer)
      exit_flag=ibuffer(1)
      ioerror=ibuffer(2)
      IF (dobroadcast.and.(exit_flag.eq.noerror)) THEN
        CALL mp_bcastf (ng, model, a)
      END IF
#endif
!
!  Check if the following attributes: "scale_factor", "add_offset", and
!  "_FillValue" are present in the input NetCDF variable:
!
!  If the "scale_value" attribute is present, the data is multiplied by
!  this factor after reading.
!  If the "add_offset" attribute is present, this value is added to the
!  data after reading.
!  If both "scale_factor" and "add_offset" attributes are present, the
!  data are first scaled before the offset is added.
!  If the "_FillValue" attribute is present, the data having this value
!  is treated as missing and it is replaced with zero. This feature it
!  is usually related with the land/sea masking.
!
      attname(1)='scale_factor'
      attname(2)='add_offset  '
      attname(3)='_FillValue  '
 
      CALL netcdf_get_fatt (ng, model, ncname, varid, attname,          &
     &                      attvalue, foundit,                          &
     &                      ncid = my_ncid)
 
      IF (exit_flag.eq.noerror) THEN
        IF (.not.foundit(1)) THEN
          afactor=1.0_r8
        ELSE
          afactor=attvalue(1)
        END IF
 
        IF (.not.foundit(2)) THEN
          aoffset=0.0_r8
        ELSE
          aoffset=attvalue(2)
        END IF
 
        IF (.not.foundit(3)) THEN
          aspval=spval_check
        ELSE
          aspval=attvalue(3)
        END IF
 
        DO j=1,asize(2)                       ! zero out missing values
          DO i=1,asize(1)
            IF ((foundit(3).and.(abs(a(i,j)-aspval).lt.aepsilon)).or.   &
     &          (.not.foundit(3).and.(abs(a(i,j)).ge.abs(aspval)))) THEN
              a(i,j)=0.0_r8
            END IF
          END DO
        END DO
 
        IF (foundit(1)) THEN                  ! scale data
          DO j=1,asize(2)
            DO i=1,asize(1)
              a(i,j)=afactor*a(i,j)
            END DO
          END DO
        END IF
 
        IF (foundit(2)) THEN                  ! add data offset
          DO j=1,asize(2)
            DO i=1,asize(1)
              a(i,j)=a(i,j)+aoffset
            END DO
          END DO
        END IF
      END IF
!
!  Compute minimum and maximum values of read variable.
!
      IF (PRESENT(min_val)) THEN
        min_val=minval(a)
      END IF
      IF (PRESENT(max_val)) THEN
        max_val=maxval(a)
      END IF
!
!  If NetCDF file ID is not provided, close input NetCDF file.
!
      IF (.not.PRESENT(ncid)) THEN
        CALL netcdf_close (ng, model, my_ncid, ncname, .false.)
      END IF
!
  10  FORMAT (/,' NETCDF_GET_FVAR_2D - error while reading variable:',  &
     &        2x,a,/,22x,'in input file:',2x,a,/,22x,'call from:',2x,a, &
     &        /,22x,a)
  20  FORMAT (/,' NETCDF_GET_FVAR_2D - error while inquiring ID for ',  &
     &        'variable:',2x,a,/,22x,'in input file:',2x,a,/,22x,       &
     &        'call from:',2x,a,/,22x,a)
!
      RETURN

References mod_scalars::exit_flag, mod_parallel::inpthread, mod_iounits::ioerror, mod_netcdf::netcdf_close(), mod_netcdf::netcdf_open(), mod_scalars::noerror, mod_iounits::sourcefile, mod_scalars::spval_check, and mod_iounits::stdout.

Here is the call graph for this function:

netcdf_get_fvar_2dp()

subroutine mod_netcdf::netcdf_get_fvar::netcdf_get_fvar_2dp	(	integer, intent(in)	ng,
		integer, intent(in)	model,
		character (len=*), intent(in)	ncname,
		character (len=*), intent(in)	myvarname,
		real(dp), dimension(:,:), intent(out)	a,
		integer, intent(in), optional	ncid,
		integer, dimension(:), intent(in), optional	start,
		integer, dimension(:), intent(in), optional	total,
		logical, intent(in), optional	broadcast,
		real(dp), intent(out), optional	min_val,
		real(dp), intent(out), optional	max_val )

Definition at line 2927 of file mod_netcdf.F.

!
!=======================================================================
!                                                                      !
!  This routine reads requested double-precision 2D-array variable     !
!  from specified NetCDF file.                                         !
!                                                                      !
!  On Input:                                                           !
!                                                                      !
!     ng           Nested grid number (integer)                        !
!     model        Calling model identifier (integer)                  !
!     ncname       NetCDF file name (string)                           !
!     myVarName    Variable name (string)                              !
!     ncid         NetCDF file ID (integer, OPTIONAL)                  !
!     start        Starting index where the first of the data values   !
!                    will be read along each dimension (integer,       !
!                    OPTIONAL)                                         !
!     total        Number of data values to be read along each         !
!                    dimension (integer, OPTIONAL)                     !
!     broadcast    Switch to broadcast read values from root to all    !
!                    members of the communicator in distributed-       !
!                    memory applications (logical, OPTIONAL,           !
!                    default=TRUE)                                     !
!                                                                      !
!  On Ouput:                                                           !
!                                                                      !
!     A            Read 2D-array variable (real)                       !
!     min_val      Read data minimum value (real, OPTIONAL)            !
!     max_val      Read data maximum value (real, OPTIONAL)            !
!                                                                      !
!  Examples:                                                           !
!                                                                      !
!    CALL netcdf_get_fvar (ng, iNLM, 'file.nc', 'VarName', fvar)       !
!    CALL netcdf_get_fvar (ng, iNLM, 'file.nc', 'VarName', fvar(0:,:)) !
!    CALL netcdf_get_fvar (ng, iNLM, 'file.nc', 'VarName', fvar(0:,0:))!
!    CALL netcdf_get_fvar (ng, iNLM, 'file.nc', 'VarName', fvar(:,:,1))!
!                                                                      !
!=======================================================================
!
!  Imported variable declarations.
!
      logical, intent(in), optional :: broadcast
!
      integer, intent(in) :: ng, model
 
      integer, intent(in), optional :: ncid
      integer, intent(in), optional :: start(:)
      integer, intent(in), optional :: total(:)
!
      character (len=*), intent(in) :: ncname
      character (len=*), intent(in) :: myVarName
!
      real(dp), intent(out), optional :: min_val
      real(dp), intent(out), optional :: max_val
 
      real(dp), intent(out) :: A(:,:)
!
!  Local variable declarations.
!
# if !defined PARALLEL_IO && defined DISTRIBUTE
      logical :: Dobroadcast
# endif
      logical, dimension(3) :: foundit
!
      integer :: i, j, my_ncid, status, varid
 
      integer, dimension(2) :: Asize
 
# if !defined PARALLEL_IO && defined DISTRIBUTE
      integer, dimension(2) :: ibuffer
# endif
!
      real(dp) :: Afactor, Aoffset, Aspval
 
      real(dp), parameter :: Aepsilon = 1.0e-8_r8
 
      real(dp), dimension(3) :: AttValue
!
      character (len=12), dimension(3) :: AttName
 
      character (len=*), parameter :: MyFile =                          &
     &  __FILE__//", netcdf_get_fvar_2dp"
!
!-----------------------------------------------------------------------
!  Read in a double-precision 2D-array variable.
!-----------------------------------------------------------------------
!
      IF (PRESENT(start).and.PRESENT(total)) THEN
        asize(1)=total(1)
        asize(2)=total(2)
      ELSE
        asize(1)=ubound(a, dim=1)
        asize(2)=ubound(a, dim=2)
      END IF
!
!  If NetCDF file ID is not provided, open NetCDF for reading.
!
      IF (.not.PRESENT(ncid)) THEN
        CALL netcdf_open (ng, model, trim(ncname), 0, my_ncid)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
      ELSE
        my_ncid=ncid
      END IF
 
# if !defined PARALLEL_IO && defined DISTRIBUTE
!
!  Determine if the read data is only needed and allocated by the master
!  node ins distribute-memory applications.
!
      IF (.not.PRESENT(broadcast)) THEN
        dobroadcast=.true.
      ELSE
        dobroadcast=broadcast
      END IF
# endif
!
!  Read in variable.
!
      IF (inpthread) THEN
        status=nf90_inq_varid(my_ncid, trim(myvarname), varid)
        IF (status.eq.nf90_noerr) THEN
          IF (PRESENT(start).and.PRESENT(total)) THEN
            status=nf90_get_var(my_ncid, varid, a, start, total)
          ELSE
            status=nf90_get_var(my_ncid, varid, a)
          END IF
          IF (founderror(status, nf90_noerr, __line__, myfile)) THEN
            WRITE (stdout,10) trim(myvarname), trim(ncname),            &
     &                        trim(sourcefile), nf90_strerror(status)
            exit_flag=2
            ioerror=status
          END IF
        ELSE
          WRITE (stdout,20) trim(myvarname), trim(ncname),              &
     &                      trim(sourcefile), nf90_strerror(status),    &
     &                      nf90_strerror(status)
          exit_flag=2
          ioerror=status
        END IF
      END IF
 
# if !defined PARALLEL_IO && defined DISTRIBUTE
!
!  Broadcast read variable to all processors in the group.
!
      ibuffer(1)=exit_flag
      ibuffer(2)=ioerror
      CALL mp_bcasti (ng, model, ibuffer)
      exit_flag=ibuffer(1)
      ioerror=ibuffer(2)
      IF (dobroadcast.and.(exit_flag.eq.noerror)) THEN
        CALL mp_bcastf (ng, model, a)
      END IF
# endif
!
!  Check if the following attributes: "scale_factor", "add_offset", and
!  "_FillValue" are present in the input NetCDF variable:
!
!  If the "scale_value" attribute is present, the data is multiplied by
!  this factor after reading.
!  If the "add_offset" attribute is present, this value is added to the
!  data after reading.
!  If both "scale_factor" and "add_offset" attributes are present, the
!  data are first scaled before the offset is added.
!  If the "_FillValue" attribute is present, the data having this value
!  is treated as missing and it is replaced with zero. This feature it
!  is usually related with the land/sea masking.
!
      attname(1)='scale_factor'
      attname(2)='add_offset  '
      attname(3)='_FillValue  '
 
      CALL netcdf_get_fatt (ng, model, ncname, varid, attname,          &
     &                      attvalue, foundit,                          &
     &                      ncid = my_ncid)
 
      IF (exit_flag.eq.noerror) THEN
        IF (.not.foundit(1)) THEN
          afactor=1.0_r8
        ELSE
          afactor=attvalue(1)
        END IF
 
        IF (.not.foundit(2)) THEN
          aoffset=0.0_r8
        ELSE
          aoffset=attvalue(2)
        END IF
 
        IF (.not.foundit(3)) THEN
          aspval=spval_check
        ELSE
          aspval=attvalue(3)
        END IF
 
        DO j=1,asize(2)                       ! zero out missing values
          DO i=1,asize(1)
            IF ((foundit(3).and.(abs(a(i,j)-aspval).lt.aepsilon)).or.   &
     &          (.not.foundit(3).and.(abs(a(i,j)).ge.abs(aspval)))) THEN
              a(i,j)=0.0_r8
            END IF
          END DO
        END DO
 
        IF (foundit(1)) THEN                  ! scale data
          DO j=1,asize(2)
            DO i=1,asize(1)
              a(i,j)=afactor*a(i,j)
            END DO
          END DO
        END IF
 
        IF (foundit(2)) THEN                  ! add data offset
          DO j=1,asize(2)
            DO i=1,asize(1)
              a(i,j)=a(i,j)+aoffset
            END DO
          END DO
        END IF
      END IF
!
!  Compute minimum and maximum values of read variable.
!
      IF (PRESENT(min_val)) THEN
        min_val=minval(a)
      END IF
      IF (PRESENT(max_val)) THEN
        max_val=maxval(a)
      END IF
!
!  If NetCDF file ID is not provided, close input NetCDF file.
!
      IF (.not.PRESENT(ncid)) THEN
        CALL netcdf_close (ng, model, my_ncid, ncname, .false.)
      END IF
!
  10  FORMAT (/,' NETCDF_GET_FVAR_2DP - error while reading variable:', &
     &        2x,a,/,23x,'in input file:',2x,a,/,23x,'call from:',2x,a, &
     &        /,23x,a)
  20  FORMAT (/,' NETCDF_GET_FVAR_2DP - error while inquiring ID for ', &
     &        'variable:',2x,a,/,23x,'in input file:',2x,a,/,23x,       &
     &        'call from:',2x,a,/,23x,a)
!
      RETURN

References mod_scalars::exit_flag, mod_parallel::inpthread, mod_iounits::ioerror, mod_netcdf::netcdf_close(), mod_netcdf::netcdf_open(), mod_scalars::noerror, mod_iounits::sourcefile, mod_scalars::spval_check, and mod_iounits::stdout.

Here is the call graph for this function:

netcdf_get_fvar_3d()

subroutine mod_netcdf::netcdf_get_fvar::netcdf_get_fvar_3d	(	integer, intent(in)	ng,
		integer, intent(in)	model,
		character (len=*), intent(in)	ncname,
		character (len=*), intent(in)	myvarname,
		real(r8), dimension(:,:,:), intent(out)	a,
		integer, intent(in), optional	ncid,
		integer, dimension(:), intent(in), optional	start,
		integer, dimension(:), intent(in), optional	total,
		logical, intent(in), optional	broadcast,
		real(r8), intent(out), optional	min_val,
		real(r8), intent(out), optional	max_val )

Definition at line 4084 of file mod_netcdf.F.

!
!=======================================================================
!                                                                      !
!  This routine reads requested floating-point 3D-array variable from  !
!  specified NetCDF file.                                              !
!                                                                      !
!  On Input:                                                           !
!                                                                      !
!     ng           Nested grid number (integer)                        !
!     model        Calling model identifier (integer)                  !
!     ncname       NetCDF file name (string)                           !
!     myVarName    Variable name (string)                              !
!     ncid         NetCDF file ID (integer, OPTIONAL)                  !
!     start        Starting index where the first of the data values   !
!                    will be read along each dimension (integer,       !
!                    OPTIONAL)                                         !
!     total        Number of data values to be read along each         !
!                    dimension (integer, OPTIONAL)                     !
!     broadcast    Switch to broadcast read values from root to all    !
!                    members of the communicator in distributed-       !
!                    memory applications (logical, OPTIONAL,           !
!                    default=TRUE)                                     !
!                                                                      !
!  On Ouput:                                                           !
!                                                                      !
!     A            Read 3D-array variable (real)                       !
!     min_val      Read data minimum value (real, OPTIONAL)            !
!     max_val      Read data maximum value (real, OPTIONAL)            !
!                                                                      !
!  Examples:                                                           !
!                                                                      !
!    CALL netcdf_get_fvar (ng, iNLM, 'file.nc', 'VarName', fvar)       !
!                                                                      !
!=======================================================================
!
!  Imported variable declarations.
!
      logical, intent(in), optional :: broadcast
!
      integer, intent(in) :: ng, model
 
      integer, intent(in), optional :: ncid
      integer, intent(in), optional :: start(:)
      integer, intent(in), optional :: total(:)
!
      character (len=*), intent(in) :: ncname
      character (len=*), intent(in) :: myVarName
!
      real(r8), intent(out), optional :: min_val
      real(r8), intent(out), optional :: max_val
 
      real(r8), intent(out) :: A(:,:,:)
!
!  Local variable declarations.
!
#if !defined PARALLEL_IO && defined DISTRIBUTE
      logical :: DoBroadcast
#endif
      logical, dimension(3) :: foundit
!
      integer :: i, j, k, my_ncid, status, varid
 
      integer, dimension(3) :: Asize
 
#if !defined PARALLEL_IO && defined DISTRIBUTE
      integer, dimension(2) :: ibuffer
#endif
!
      real(r8) :: Afactor, Aoffset, Aspval
 
      real(r8), parameter :: Aepsilon = 1.0e-8_r8
 
      real(r8), dimension(3) :: AttValue
!
      character (len=12), dimension(3) :: AttName
 
      character (len=*), parameter :: MyFile =                          &
     &  __FILE__//", netcdf_get_fvar_3d"
!
!-----------------------------------------------------------------------
!  Read in a floating-point 2D-array variable.
!-----------------------------------------------------------------------
!
      IF (PRESENT(start).and.PRESENT(total)) THEN
        asize(1)=total(1)
        asize(2)=total(2)
        asize(3)=total(3)
      ELSE
        asize(1)=ubound(a, dim=1)
        asize(2)=ubound(a, dim=2)
        asize(3)=ubound(a, dim=3)
      END IF
!
!  If NetCDF file ID is not provided, open NetCDF for reading.
!
      IF (.not.PRESENT(ncid)) THEN
        CALL netcdf_open (ng, model, trim(ncname), 0, my_ncid)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
      ELSE
        my_ncid=ncid
      END IF
 
#if !defined PARALLEL_IO && defined DISTRIBUTE
!
!  Determine if the read data is only needed and allocated by the master
!  node ins distribute-memory applications.
!
      IF (.not.PRESENT(broadcast)) THEN
        dobroadcast=.true.
      ELSE
        dobroadcast=broadcast
      END IF
#endif
!
!  Read in variable.
!
      IF (inpthread) THEN
        status=nf90_inq_varid(my_ncid, trim(myvarname), varid)
        IF (status.eq.nf90_noerr) THEN
          IF (PRESENT(start).and.PRESENT(total)) THEN
            status=nf90_get_var(my_ncid, varid, a, start, total)
          ELSE
            status=nf90_get_var(my_ncid, varid, a)
          END IF
          IF (founderror(status, nf90_noerr, __line__, myfile)) THEN
            WRITE (stdout,10) trim(myvarname), trim(ncname),            &
     &                        trim(sourcefile), nf90_strerror(status)
            exit_flag=2
            ioerror=status
          END IF
        ELSE
          WRITE (stdout,20) trim(myvarname), trim(ncname),              &
     &                      trim(sourcefile), nf90_strerror(status)
          exit_flag=2
          ioerror=status
        END IF
      END IF
 
#if !defined PARALLEL_IO && defined DISTRIBUTE
!
!  Broadcast read variable to all processors in the group.
!
      ibuffer(1)=exit_flag
      ibuffer(2)=ioerror
      CALL mp_bcasti (ng, model, ibuffer)
      exit_flag=ibuffer(1)
      ioerror=ibuffer(2)
      IF (dobroadcast.and.(exit_flag.eq.noerror)) THEN
        CALL mp_bcastf (ng, model, a)
      END IF
#endif
!
!  Check if the following attributes: "scale_factor", "add_offset", and
!  "_FillValue" are present in the input NetCDF variable:
!
!  If the "scale_value" attribute is present, the data is multiplied by
!  this factor after reading.
!  If the "add_offset" attribute is present, this value is added to the
!  data after reading.
!  If both "scale_factor" and "add_offset" attributes are present, the
!  data are first scaled before the offset is added.
!  If the "_FillValue" attribute is present, the data having this value
!  is treated as missing and it is replaced with zero. This feature it
!  is usually related with the land/sea masking.
!
      attname(1)='scale_factor'
      attname(2)='add_offset  '
      attname(3)='_FillValue  '
 
      CALL netcdf_get_fatt (ng, model, ncname, varid, attname,          &
     &                      attvalue, foundit,                          &
     &                      ncid = my_ncid)
 
      IF (exit_flag.eq.noerror) THEN
        IF (.not.foundit(1)) THEN
          afactor=1.0_r8
        ELSE
          afactor=attvalue(1)
        END IF
 
        IF (.not.foundit(2)) THEN
          aoffset=0.0_r8
        ELSE
          aoffset=attvalue(2)
        END IF
 
        IF (.not.foundit(3)) THEN
          aspval=spval_check
        ELSE
          aspval=attvalue(3)
        END IF
 
        DO k=1,asize(3)                       ! zero out missing values
          DO j=1,asize(2)
            DO i=1,asize(1)
              IF ((foundit(3).and.                                      &
     &             (abs(a(i,j,k)-aspval).lt.aepsilon)).or.              &
     &            (.not.foundit(3).and.                                 &
     &             (abs(a(i,j,k)).ge.abs(aspval)))) THEN
                a(i,j,k)=0.0_r8
              END IF
            END DO
          END DO
        END DO
 
        IF (foundit(1)) THEN                  ! scale data
          DO k=1,asize(3)
            DO j=1,asize(2)
              DO i=1,asize(1)
                a(i,j,k)=afactor*a(i,j,k)
              END DO
            END DO
          END DO
        END IF
 
        IF (foundit(2)) THEN                  ! add data offset
          DO k=1,asize(3)
            DO j=1,asize(2)
              DO i=1,asize(1)
                a(i,j,k)=a(i,j,k)+aoffset
              END DO
            END DO
          END DO
        END IF
      END IF
!
!  Compute minimum and maximum values of read variable.
!
      IF (PRESENT(min_val)) THEN
        min_val=minval(a)
      END IF
      IF (PRESENT(max_val)) THEN
        max_val=maxval(a)
      END IF
!
!  If NetCDF file ID is not provided, close input NetCDF file.
!
      IF (.not.PRESENT(ncid)) THEN
        CALL netcdf_close (ng, model, my_ncid, ncname, .false.)
      END IF
!
  10  FORMAT (/,' NETCDF_GET_FVAR_3D - error while reading variable:',  &
     &        2x,a,/,22x,'in input file:',2x,a,/,22x,'call from:',2x,a, &
     &        /,22x,a)
  20  FORMAT (/,' NETCDF_GET_FVAR_3D - error while inquiring ID for ',  &
     &        'variable:',2x,a,/,22x,'in input file:',2x,a,/,22x,       &
     &        'call from:',2x,a,/,22x,a)
!
      RETURN

References mod_scalars::exit_flag, mod_parallel::inpthread, mod_iounits::ioerror, mod_netcdf::netcdf_close(), mod_netcdf::netcdf_open(), mod_scalars::noerror, mod_iounits::sourcefile, mod_scalars::spval_check, and mod_iounits::stdout.

Here is the call graph for this function:

netcdf_get_fvar_3dp()

subroutine mod_netcdf::netcdf_get_fvar::netcdf_get_fvar_3dp	(	integer, intent(in)	ng,
		integer, intent(in)	model,
		character (len=*), intent(in)	ncname,
		character (len=*), intent(in)	myvarname,
		real(dp), dimension(:,:,:), intent(out)	a,
		integer, intent(in), optional	ncid,
		integer, dimension(:), intent(in), optional	start,
		integer, dimension(:), intent(in), optional	total,
		logical, intent(in), optional	broadcast,
		real(dp), intent(out), optional	min_val,
		real(dp), intent(out), optional	max_val )

Definition at line 3175 of file mod_netcdf.F.

!
!=======================================================================
!                                                                      !
!  This routine reads requested double-precision 3D-array variable     !
!  from specified NetCDF file.                                         !
!                                                                      !
!  On Input:                                                           !
!                                                                      !
!     ng           Nested grid number (integer)                        !
!     model        Calling model identifier (integer)                  !
!     ncname       NetCDF file name (string)                           !
!     myVarName    Variable name (string)                              !
!     ncid         NetCDF file ID (integer, OPTIONAL)                  !
!     start        Starting index where the first of the data values   !
!                    will be read along each dimension (integer,       !
!                    OPTIONAL)                                         !
!     total        Number of data values to be read along each         !
!                    dimension (integer, OPTIONAL)                     !
!     broadcast    Switch to broadcast read values from root to all    !
!                    members of the communicator in distributed-       !
!                    memory applications (logical, OPTIONAL,           !
!                    default=TRUE)                                     !
!                                                                      !
!  On Ouput:                                                           !
!                                                                      !
!     A            Read 3D-array variable (real)                       !
!     min_val      Read data minimum value (real, OPTIONAL)            !
!     max_val      Read data maximum value (real, OPTIONAL)            !
!                                                                      !
!  Examples:                                                           !
!                                                                      !
!    CALL netcdf_get_fvar (ng, iNLM, 'file.nc', 'VarName', fvar)       !
!                                                                      !
!=======================================================================
!
!  Imported variable declarations.
!
      logical, intent(in), optional :: broadcast
!
      integer, intent(in) :: ng, model
 
      integer, intent(in), optional :: ncid
      integer, intent(in), optional :: start(:)
      integer, intent(in), optional :: total(:)
!
      character (len=*), intent(in) :: ncname
      character (len=*), intent(in) :: myVarName
!
      real(dp), intent(out), optional :: min_val
      real(dp), intent(out), optional :: max_val
 
      real(dp), intent(out) :: A(:,:,:)
!
!  Local variable declarations.
!
# if !defined PARALLEL_IO && defined DISTRIBUTE
      logical :: DoBroadcast
# endif
      logical, dimension(3) :: foundit
!
      integer :: i, j, k, my_ncid, status, varid
 
      integer, dimension(3) :: Asize
 
# if !defined PARALLEL_IO && defined DISTRIBUTE
      integer, dimension(2) :: ibuffer
# endif
!
      real(dp) :: Afactor, Aoffset, Aspval
 
      real(dp), parameter :: Aepsilon = 1.0e-8_r8
 
      real(dp), dimension(3) :: AttValue
!
      character (len=12), dimension(3) :: AttName
 
      character (len=*), parameter :: MyFile =                          &
     &  __FILE__//", netcdf_get_fvar_3dp"
!
!-----------------------------------------------------------------------
!  Read in a double-precision 3D-array variable.
!-----------------------------------------------------------------------
!
      IF (PRESENT(start).and.PRESENT(total)) THEN
        asize(1)=total(1)
        asize(2)=total(2)
        asize(3)=total(3)
      ELSE
        asize(1)=ubound(a, dim=1)
        asize(2)=ubound(a, dim=2)
        asize(3)=ubound(a, dim=3)
      END IF
!
!  If NetCDF file ID is not provided, open NetCDF for reading.
!
      IF (.not.PRESENT(ncid)) THEN
        CALL netcdf_open (ng, model, trim(ncname), 0, my_ncid)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
      ELSE
        my_ncid=ncid
      END IF
 
# if !defined PARALLEL_IO && defined DISTRIBUTE
!
!  Determine if the read data is only needed and allocated by the master
!  node ins distribute-memory applications.
!
      IF (.not.PRESENT(broadcast)) THEN
        dobroadcast=.true.
      ELSE
        dobroadcast=broadcast
      END IF
# endif
!
!  Read in variable.
!
      IF (inpthread) THEN
        status=nf90_inq_varid(my_ncid, trim(myvarname), varid)
        IF (status.eq.nf90_noerr) THEN
          IF (PRESENT(start).and.PRESENT(total)) THEN
            status=nf90_get_var(my_ncid, varid, a, start, total)
          ELSE
            status=nf90_get_var(my_ncid, varid, a)
          END IF
          IF (founderror(status, nf90_noerr, __line__, myfile)) THEN
            WRITE (stdout,10) trim(myvarname), trim(ncname),            &
     &                        trim(sourcefile), nf90_strerror(status)
            exit_flag=2
            ioerror=status
          END IF
        ELSE
          WRITE (stdout,20) trim(myvarname), trim(ncname),              &
     &                      trim(sourcefile), nf90_strerror(status)
          exit_flag=2
          ioerror=status
        END IF
      END IF
 
# if !defined PARALLEL_IO && defined DISTRIBUTE
!
!  Broadcast read variable to all processors in the group.
!
      ibuffer(1)=exit_flag
      ibuffer(2)=ioerror
      CALL mp_bcasti (ng, model, ibuffer)
      exit_flag=ibuffer(1)
      ioerror=ibuffer(2)
      IF (dobroadcast.and.(exit_flag.eq.noerror)) THEN
        CALL mp_bcastf (ng, model, a)
      END IF
# endif
!
!  Check if the following attributes: "scale_factor", "add_offset", and
!  "_FillValue" are present in the input NetCDF variable:
!
!  If the "scale_value" attribute is present, the data is multiplied by
!  this factor after reading.
!  If the "add_offset" attribute is present, this value is added to the
!  data after reading.
!  If both "scale_factor" and "add_offset" attributes are present, the
!  data are first scaled before the offset is added.
!  If the "_FillValue" attribute is present, the data having this value
!  is treated as missing and it is replaced with zero. This feature it
!  is usually related with the land/sea masking.
!
      attname(1)='scale_factor'
      attname(2)='add_offset  '
      attname(3)='_FillValue  '
 
      CALL netcdf_get_fatt (ng, model, ncname, varid, attname,          &
     &                      attvalue, foundit,                          &
     &                      ncid = my_ncid)
 
      IF (exit_flag.eq.noerror) THEN
        IF (.not.foundit(1)) THEN
          afactor=1.0_r8
        ELSE
          afactor=attvalue(1)
        END IF
 
        IF (.not.foundit(2)) THEN
          aoffset=0.0_r8
        ELSE
          aoffset=attvalue(2)
        END IF
 
        IF (.not.foundit(3)) THEN
          aspval=spval_check
        ELSE
          aspval=attvalue(3)
        END IF
 
        DO k=1,asize(3)                       ! zero out missing values
          DO j=1,asize(2)
            DO i=1,asize(1)
              IF ((foundit(3).and.                                      &
     &             (abs(a(i,j,k)-aspval).lt.aepsilon)).or.              &
     &            (.not.foundit(3).and.                                 &
     &             (abs(a(i,j,k)).ge.abs(aspval)))) THEN
                a(i,j,k)=0.0_r8
              END IF
            END DO
          END DO
        END DO
 
        IF (foundit(1)) THEN                  ! scale data
          DO k=1,asize(3)
            DO j=1,asize(2)
              DO i=1,asize(1)
                a(i,j,k)=afactor*a(i,j,k)
              END DO
            END DO
          END DO
        END IF
 
        IF (foundit(2)) THEN                  ! add data offset
          DO k=1,asize(3)
            DO j=1,asize(2)
              DO i=1,asize(1)
                a(i,j,k)=a(i,j,k)+aoffset
              END DO
            END DO
          END DO
        END IF
      END IF
!
!  Compute minimum and maximum values of read variable.
!
      IF (PRESENT(min_val)) THEN
        min_val=minval(a)
      END IF
      IF (PRESENT(max_val)) THEN
        max_val=maxval(a)
      END IF
!
!  If NetCDF file ID is not provided, close input NetCDF file.
!
      IF (.not.PRESENT(ncid)) THEN
        CALL netcdf_close (ng, model, my_ncid, ncname, .false.)
      END IF
!
  10  FORMAT (/,' NETCDF_GET_FVAR_3DP - error while reading variable:', &
     &        2x,a,/,23x,'in input file:',2x,a,/,23x,'call from:',2x,a, &
     &        /,23x,a)
  20  FORMAT (/,' NETCDF_GET_FVAR_3DP - error while inquiring ID for ', &
     &        'variable:',2x,a,/,23x,'in input file:',2x,a,/,23x,       &
     &        'call from:',2x,a,/,23x,a)
!
      RETURN

References mod_scalars::exit_flag, mod_parallel::inpthread, mod_iounits::ioerror, mod_netcdf::netcdf_close(), mod_netcdf::netcdf_open(), mod_scalars::noerror, mod_iounits::sourcefile, mod_scalars::spval_check, and mod_iounits::stdout.

Here is the call graph for this function:

netcdf_get_fvar_4d()

subroutine mod_netcdf::netcdf_get_fvar::netcdf_get_fvar_4d	(	integer, intent(in)	ng,
		integer, intent(in)	model,
		character (len=*), intent(in)	ncname,
		character (len=*), intent(in)	myvarname,
		real(r8), dimension(:,:,:,:), intent(out)	a,
		integer, intent(in), optional	ncid,
		integer, dimension(:), intent(in), optional	start,
		integer, dimension(:), intent(in), optional	total,
		logical, intent(in), optional	broadcast,
		real(r8), intent(out), optional	min_val,
		real(r8), intent(out), optional	max_val )

Definition at line 4338 of file mod_netcdf.F.

!
!=======================================================================
!                                                                      !
!  This routine reads requested floating-point 4D-array variable from  !
!  specified NetCDF file.                                              !
!                                                                      !
!  On Input:                                                           !
!                                                                      !
!     ng           Nested grid number (integer)                        !
!     model        Calling model identifier (integer)                  !
!     ncname       NetCDF file name (string)                           !
!     myVarName    Variable name (string)                              !
!     ncid         NetCDF file ID (integer, OPTIONAL)                  !
!     start        Starting index where the first of the data values   !
!                    will be read along each dimension (integer,       !
!                    OPTIONAL)                                         !
!     total        Number of data values to be read along each         !
!                    dimension (integer, OPTIONAL)                     !
!     broadcast    Switch to broadcast read values from root to all    !
!                    members of the communicator in distributed-       !
!                    memory applications (logical, OPTIONAL,           !
!                    default=TRUE)                                     !
!                                                                      !
!  On Ouput:                                                           !
!                                                                      !
!     A            Read 4D-array variable (real)                       !
!     min_val      Read data minimum value (real, OPTIONAL)            !
!     max_val      Read data maximum value (real, OPTIONAL)            !
!                                                                      !
!  Examples:                                                           !
!                                                                      !
!    CALL netcdf_get_fvar (ng, iNLM, 'file.nc', 'VarName', fvar)       !
!                                                                      !
!=======================================================================
!
!  Imported variable declarations.
!
      logical, intent(in), optional :: broadcast
!
      integer, intent(in) :: ng, model
 
      integer, intent(in), optional :: ncid
      integer, intent(in), optional :: start(:)
      integer, intent(in), optional :: total(:)
!
      character (len=*), intent(in) :: ncname
      character (len=*), intent(in) :: myVarName
!
      real(r8), intent(out), optional :: min_val
      real(r8), intent(out), optional :: max_val
 
      real(r8), intent(out) :: A(:,:,:,:)
!
!  Local variable declarations.
!
#if !defined PARALLEL_IO && defined DISTRIBUTE
      logical :: DoBroadcast
#endif
      logical, dimension(3) :: foundit
!
      integer :: i, j, k, l, my_ncid, status, varid
 
      integer, dimension(4) :: Asize
 
#if !defined PARALLEL_IO && defined DISTRIBUTE
      integer, dimension(2) :: ibuffer
#endif
!
      real(r8) :: Afactor, Aoffset, Aspval
 
      real(r8), parameter :: Aepsilon = 1.0e-8_r8
 
      real(r8), dimension(3) :: AttValue
!
      character (len=12), dimension(3) :: AttName
 
      character (len=*), parameter :: MyFile =                          &
     &  __FILE__//", netcdf_get_fvar_4d"
!
!-----------------------------------------------------------------------
!  Read in a floating-point 2D-array variable.
!-----------------------------------------------------------------------
!
      IF (PRESENT(start).and.PRESENT(total)) THEN
        asize(1)=total(1)
        asize(2)=total(2)
        asize(3)=total(3)
        asize(4)=total(4)
      ELSE
        asize(1)=ubound(a, dim=1)
        asize(2)=ubound(a, dim=2)
        asize(3)=ubound(a, dim=3)
        asize(4)=ubound(a, dim=4)
      END IF
!
!  If NetCDF file ID is not provided, open NetCDF for reading.
!
      IF (.not.PRESENT(ncid)) THEN
        CALL netcdf_open (ng, model, trim(ncname), 0, my_ncid)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
      ELSE
        my_ncid=ncid
      END IF
 
#if !defined PARALLEL_IO && defined DISTRIBUTE
!
!  Determine if the read data is only needed and allocated by the master
!  node ins distribute-memory applications.
!
      IF (.not.PRESENT(broadcast)) THEN
        dobroadcast=.true.
      ELSE
        dobroadcast=broadcast
      END IF
#endif
!
!  Read in variable.
!
      IF (inpthread) THEN
        status=nf90_inq_varid(my_ncid, trim(myvarname), varid)
        IF (status.eq.nf90_noerr) THEN
          IF (PRESENT(start).and.PRESENT(total)) THEN
            status=nf90_get_var(my_ncid, varid, a, start, total)
          ELSE
            status=nf90_get_var(my_ncid, varid, a)
          END IF
          IF (founderror(status, nf90_noerr, __line__, myfile)) THEN
            WRITE (stdout,10) trim(myvarname), trim(ncname),            &
     &                        trim(sourcefile), nf90_strerror(status)
            exit_flag=2
            ioerror=status
          END IF
        ELSE
          WRITE (stdout,20) trim(myvarname), trim(ncname),              &
     &                      trim(sourcefile), nf90_strerror(status)
          exit_flag=2
          ioerror=status
        END IF
      END IF
 
#if !defined PARALLEL_IO && defined DISTRIBUTE
!
!  Broadcast read variable to all processors in the group.
!
      ibuffer(1)=exit_flag
      ibuffer(2)=ioerror
      CALL mp_bcasti (ng, model, ibuffer)
      exit_flag=ibuffer(1)
      ioerror=ibuffer(2)
      IF (dobroadcast.and.(exit_flag.eq.noerror)) THEN
        CALL mp_bcastf (ng, model, a)
      END IF
#endif
!
!  Check if the following attributes: "scale_factor", "add_offset", and
!  "_FillValue" are present in the input NetCDF variable:
!
!  If the "scale_value" attribute is present, the data is multiplied by
!  this factor after reading.
!  If the "add_offset" attribute is present, this value is added to the
!  data after reading.
!  If both "scale_factor" and "add_offset" attributes are present, the
!  data are first scaled before the offset is added.
!  If the "_FillValue" attribute is present, the data having this value
!  is treated as missing and it is replaced with zero. This feature it
!  is usually related with the land/sea masking.
!
      attname(1)='scale_factor'
      attname(2)='add_offset  '
      attname(3)='_FillValue  '
 
      CALL netcdf_get_fatt (ng, model, ncname, varid, attname,          &
     &                      attvalue, foundit,                          &
     &                      ncid = my_ncid)
 
      IF (exit_flag.eq.noerror) THEN
        IF (.not.foundit(1)) THEN
          afactor=1.0_r8
        ELSE
          afactor=attvalue(1)
        END IF
 
        IF (.not.foundit(2)) THEN
          aoffset=0.0_r8
        ELSE
          aoffset=attvalue(2)
        END IF
 
        IF (.not.foundit(3)) THEN
          aspval=spval_check
        ELSE
          aspval=attvalue(3)
        END IF
 
        DO l=1,asize(4)                       ! zero out missing values
          DO k=1,asize(3)
            DO j=1,asize(2)
              DO i=1,asize(1)
                IF ((foundit(3).and.                                    &
     &               (abs(a(i,j,k,l)-aspval).lt.aepsilon)).or.          &
     &              (.not.foundit(3).and.                               &
     &               (abs(a(i,j,k,l)).ge.abs(aspval)))) THEN
                  a(i,j,k,l)=0.0_r8
                END IF
              END DO
            END DO
          END DO
        END DO
 
        IF (foundit(1)) THEN                  ! scale data
          DO l=1,asize(4)
            DO k=1,asize(3)
              DO j=1,asize(2)
                DO i=1,asize(1)
                  a(i,j,k,l)=afactor*a(i,j,k,l)
                END DO
              END DO
            END DO
          END DO
        END IF
 
        IF (foundit(2)) THEN                  ! add data offset
          DO l=1,asize(4)
            DO k=1,asize(3)
              DO j=1,asize(2)
                DO i=1,asize(1)
                  a(i,j,k,l)=a(i,j,k,l)+aoffset
                END DO
              END DO
            END DO
          END DO
        END IF
      END IF
!
!  Compute minimum and maximum values of read variable.
!
      IF (PRESENT(min_val)) THEN
        min_val=minval(a)
      END IF
      IF (PRESENT(max_val)) THEN
        max_val=maxval(a)
      END IF
!
!  If NetCDF file ID is not provided, close input NetCDF file.
!
      IF (.not.PRESENT(ncid)) THEN
        CALL netcdf_close (ng, model, my_ncid, ncname, .false.)
      END IF
!
  10  FORMAT (/,' NETCDF_GET_FVAR_4D - error while reading variable:',  &
     &        2x,a,/,22x,'in input file:',2x,a,/,22x,'call from:',2x,a, &
     &        /,22x,a)
  20  FORMAT (/,' NETCDF_GET_FVAR_4D - error while inquiring ID for ',  &
     &        'variable:',2x,a,/,22x,'in input file:',2x,a,/,22x,       &
     &        'call from:',2x,a,/,22x,a)
!
      RETURN

References mod_scalars::exit_flag, mod_parallel::inpthread, mod_iounits::ioerror, mod_netcdf::netcdf_close(), mod_netcdf::netcdf_open(), mod_scalars::noerror, mod_iounits::sourcefile, mod_scalars::spval_check, and mod_iounits::stdout.

Here is the call graph for this function:

---

The documentation for this interface was generated from the following file:

• ROMS/Modules/mod_netcdf.F