get_ngfld.F File Reference

#include "cppdefs.h"

Include dependency graph for get_ngfld.F:

Go to the source code of this file.

Functions/Subroutines
subroutine	get_ngfld (ng, model, ifield, ncid, piofile, nfiles, s, recordless, update, lbi, ubi, ubj, ubk, istr, iend, jrec, fout)
subroutine	get_ngfld_nf90 (ng, model, ifield, ncid, nfiles, s, recordless, update, lbi, ubi, ubj, ubk, istr, iend, jrec, fout)
subroutine	get_ngfld_pio (ng, model, ifield, piofile, nfiles, s, recordless, update, lbi, ubi, ubj, ubk, istr, iend, jrec, fout)

Function/Subroutine Documentation

get_ngfld()

subroutine get_ngfld	(	integer, intent(in)	ng,
		integer, intent(in)	model,
		integer, intent(in)	ifield,
		integer	ncid,
		type (file_desc_t), intent(inout)	piofile,
		integer, intent(in)	nfiles,
		type(t_io), dimension(nfiles), intent(inout)	s,
		logical, intent(in)	recordless,
		logical, intent(out)	update,
		integer, intent(in)	lbi,
		integer, intent(in)	ubi,
		integer, intent(in)	ubj,
		integer, intent(in)	ubk,
		integer, intent(in)	istr,
		integer, intent(in)	iend,
		integer, intent(in)	jrec,
		real(r8), dimension(lbi:ubi,ubj,ubk), intent(inout)	fout )

Definition at line 2 of file get_ngfld.F.

!
!git $Id$
!================================================== Hernan G. Arango ===
!  Copyright (c) 2002-2025 The ROMS Group                              !
!    Licensed under a MIT/X style license                              !
!    See License_ROMS.md                                               !
!=======================================================================
!                                                                      !
!  This routine reads in requested non-grided field from specified     !
!  NetCDF file.  A non-grided field has different dimensions  than     !
!  model spatial dimensions. Forward time processing.                  !
!                                                                      !
!  On Input:                                                           !
!                                                                      !
!     ng         Nested grid number.                                   !
!     model      Calling model identifier.                             !
!     ifield     Field ID.                                             !
!     ncid       NetCDF file ID.                                       !
# if defined PIO_LIB && defined DISTRIBUTE
!     pioFile    PIO file descriptor structure, TYPE(file_desc_t)      !
!                  pioFile%fh         file handler                     !
!                  pioFile%iosystem   IO system descriptor (struct)    !
# endif
!     nfiles     Number of input NetCDF files.                         !
!     S          I/O derived type structure, TYPE(T_IO).               !
!     recordless Switch for time invariant field (logical).            !
!     LBi        "Fout" 1st dimension lower-bound value.               !
!     UBi        "Fout" 1st dimension upper-bound value.               !
!     UBj        "Fout" 2nd dimension upper-bound value, if any.       !
!                  Otherwise, a value of one is expected.              !
!     UBk        "Fout" time dimension upper-bound value, if any.      !
!                  Otherwise, a value of one is expected.              !
!     Istr       Starting location of read data in the 1st dimension.  !
!     Iend       Ending location of read data in the 1st dimension;    !
!                  Number of records read is: Iend-Istr+1.             !
!     Jrec       Number of records read in the 2st dimension.          !
!                                                                      !
!  On Output:                                                          !
!                                                                      !
!     Fout       Read field.                                           !
!     update     Switch indicating reading of the requested field      !
!                  the current time step.                              !
!                                                                      !
!=======================================================================
!
      USE mod_param
      USE mod_parallel
      USE mod_iounits
      USE mod_ncparam
      USE mod_scalars
!
      USE strings_mod, ONLY : founderror
!
      implicit none
!
!  Imported variable declarations.
!
      logical, intent(in)  :: recordless
      logical, intent(out) :: update
!
      integer, intent(in) :: ng, model, ifield, nfiles
      integer, intent(in) :: LBi, UBi, UBj, UBk, Istr, Iend, Jrec
      integer :: ncid
!
#if defined PIO_LIB && defined DISTRIBUTE
      TYPE (File_desc_t), intent(inout) :: pioFile
#endif
      TYPE(T_IO), intent(inout) :: S(nfiles)
!
      real(r8), intent(inout) :: Fout(LBi:UBi,UBj,UBk)
!
!  Local variable declarations.
!
      character (len=*), parameter :: MyFile =                          &
     &  __FILE__
!
!-----------------------------------------------------------------------
!  Read in requested 2D field according to IO type.
!-----------------------------------------------------------------------
!
      SELECT CASE (s(1)%IOtype)
        CASE (io_nf90)
          CALL get_ngfld_nf90 (ng, model, ifield, ncid,                 &
     &                         nfiles, s, recordless, update,           &
     &                         lbi, ubi, ubj, ubk, istr, iend, jrec,    &
     &                         fout)
 
#if defined PIO_LIB && defined DISTRIBUTE
        CASE (io_pio)
          CALL get_ngfld_pio (ng, model, ifield, piofile,               &
     &                        nfiles, s, recordless, update,            &
     &                        lbi, ubi, ubj, ubk, istr, iend, jrec,     &
     &                        fout)
#endif
        CASE DEFAULT
          IF (master) WRITE (stdout,10) s(1)%IOtype
          exit_flag=2
      END SELECT
      IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
!
  10  FORMAT (' GET_NGFLD - Illegal input file type, io_type = ',i0,    &
     &        /,13x,'Check KeyWord ''INP_LIB'' in ''roms.in''.')
!
      RETURN

References mod_scalars::exit_flag, strings_mod::founderror(), get_ngfld_nf90(), get_ngfld_pio(), mod_ncparam::io_nf90, mod_ncparam::io_pio, mod_parallel::master, mod_scalars::noerror, and mod_iounits::stdout.

Referenced by get_data(), get_idata(), rp_get_data(), rp_get_idata(), tl_get_data(), and tl_get_idata().

Here is the call graph for this function:

Here is the caller graph for this function:

get_ngfld_nf90()

subroutine get_ngfld_nf90	(	integer, intent(in)	ng,
		integer, intent(in)	model,
		integer, intent(in)	ifield,
		integer	ncid,
		integer, intent(in)	nfiles,
		type(t_io), dimension(nfiles), intent(inout)	s,
		logical, intent(in)	recordless,
		logical, intent(out)	update,
		integer, intent(in)	lbi,
		integer, intent(in)	ubi,
		integer, intent(in)	ubj,
		integer, intent(in)	ubk,
		integer, intent(in)	istr,
		integer, intent(in)	iend,
		integer, intent(in)	jrec,
		real(r8), dimension(lbi:ubi,ubj,ubk), intent(inout)	fout )

Definition at line 116 of file get_ngfld.F.

!***********************************************************************
!
      USE mod_param
      USE mod_parallel
      USE mod_iounits
      USE mod_ncparam
      USE mod_netcdf
      USE mod_scalars
!
      USE dateclock_mod, ONLY : time_string
#ifdef CHECKSUM
      USE get_hash_mod,  ONLY : get_hash
#endif
      USE inquiry_mod,   ONLY : inquiry
      USE strings_mod,   ONLY : founderror
!
      implicit none
!
!  Imported variable declarations.
!
      logical, intent(in)  :: recordless
      logical, intent(out) :: update
!
      integer, intent(in) :: ng, model, ifield, nfiles
      integer, intent(in) :: LBi, UBi, UBj, UBk, Istr, Iend, Jrec
      integer :: ncid
!
      TYPE(T_IO), intent(inout) :: S(nfiles)
!
      real(r8), intent(inout) :: Fout(LBi:UBi,UBj,UBk)
!
!  Local variable declarations.
!
      logical :: Linquire, Liocycle, Lmulti, Lonerec
!
      integer :: Nrec, Tid, Tindex, Trec, Vid, Vtype
      integer :: i, ic, j, job, lend, lstr, npts, nvdim, status
#ifdef CHECKSUM
      integer(i8b) :: Fhash
#endif
!
      real(r8) :: Aval, Fmax, Fmin
 
      real(dp) :: Clength, Tdelta, Tend
      real(dp) :: Tmax, Tmin, Tmono, Tscale, Tstr
      real(dp) :: Tsec, Tval
 
      real(r8), dimension((UBi-LBi+1)*UBj) :: Awrk
!
      character (len=22) :: t_code
 
      character (len=*), parameter :: MyFile =                          &
     &  __FILE__//", get_ngfld_nf90"
!
      sourcefile=myfile
!
!-----------------------------------------------------------------------
!  Initialize.
!-----------------------------------------------------------------------
!
      IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
!
!  Determine if inquiring about the field to process in input NetCDF
!  file(s).  This usually happens on first call or when the field
!  time records are split (saved) in several multi-files.
!
      linquire=.false.
      lmulti=.false.
      IF (iic(ng).eq.0) linquire=.true.
      IF (.not.linquire.and.                                            &
     &    ((iinfo(10,ifield,ng).gt.1).and.                              &
     &     (linfo( 6,ifield,ng).or.                                     &
     &     (finfo( 2,ifield,ng)*day2sec.lt.time(ng))))) THEN
        linquire=.true.
        lmulti=.true.
      END IF
!
!  If appropriate, inquire about the contents of input NetCDF file and
!  fill information arrays.
!
!  Also, if appropriate, deactivate the Linfo(6,ifield,ng) switch after
!  the query for the UPPER snapshot interpolant from previous multifile
!  in the list. The switch was activated previously to indicate the
!  processing of the FIRST record of the file for the LOWER snapshot
!  interpolant.
!
      IF (linquire) THEN
        job=1
        CALL inquiry (ng, model, job, ubk, iend, ubi, ifield, ncid,     &
     &                lmulti, nfiles, s)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        IF (linfo(6,ifield,ng)) THEN
          linfo(6,ifield,ng)=.false.
        END IF
      END IF
!
!-----------------------------------------------------------------------
!  If appropriate, read in new data.
!-----------------------------------------------------------------------
!
      update=.false.
      tmono=finfo(7,ifield,ng)
!
      IF ((tmono.lt.time(ng)).or.(iic(ng).eq.0).or.                     &
     &    (iic(ng).eq.ntstart(ng))) THEN
!
!  Load properties for requested field from information arrays.
!
        liocycle=linfo( 2,ifield,ng)
        lonerec =linfo( 3,ifield,ng)
        vtype   =iinfo( 1,ifield,ng)
        vid     =iinfo( 2,ifield,ng)
        tid     =iinfo( 3,ifield,ng)
        nrec    =iinfo( 4,ifield,ng)
        tindex  =iinfo( 8,ifield,ng)
        trec    =iinfo( 9,ifield,ng)
        nvdim   =iinfo(11,ifield,ng)
        tmin    =finfo( 1,ifield,ng)
        tmax    =finfo( 2,ifield,ng)
        clength =finfo( 5,ifield,ng)
        tscale  =finfo( 6,ifield,ng)
        ncfile  =cinfo(ifield,ng)
!
        IF (liocycle) THEN
          trec=mod(trec,nrec)+1
        ELSE
          trec=trec+1
        END IF
        iinfo(9,ifield,ng)=trec
!
        IF (trec.le.nrec) THEN
!
!  Set rolling index for two-time record storage of input data.  If
!  "UBk" is unity, input data is stored in recordless array by the
!  calling program.
!
          IF (recordless) THEN
            tindex=1
          ELSE
            tindex=3-tindex
          END IF
          iinfo(8,ifield,ng)=tindex
!
!  Read in time coordinate.
!
          IF (.not.recordless.and.(tid.ge.0).and.(tid.ne.vid)) THEN
            CALL netcdf_get_time (ng, model, ncfile, tname(ifield),     &
     &                            rclock%DateNumber, tval,              &
     &                            ncid = ncid,                          &
     &                            start = (/trec/),                     &
     &                            total = (/1/))
            IF (founderror(exit_flag, noerror, __line__, myfile)) THEN
              IF (master) WRITE (stdout,50) trim(tname(ifield)), trec
              RETURN
            END IF
            tval=tval*tscale
            vtime(tindex,ifield,ng)=tval
!
!  Activate switch Linfo(6,ifield,ng) if processing the LAST record of
!  the file for the LOWER time snapshot. We need to get the UPPER time
!  snapshot from NEXT multifile.
!
            IF ((trec.eq.nrec).and.(tval*day2sec.le.time(ng))) THEN
              linfo(6,ifield,ng)=.true.
            END IF
          END IF
!
!  Read in non-grided data. The conditional statement on Jrec is to
!  differentiate between reading a 3D and 2D array.
!
          IF (vid.ge.0) THEN
            fmin=0.0_r8
            fmax=0.0_r8
            IF (nvdim.eq.1) THEN
              npts=iend-istr+1
              CALL netcdf_get_fvar (ng, model, ncfile,                  &
     &                              vname(1,ifield), awrk,              &
     &                              ncid = ncid,                        &
     &                              start = (/1/),                      &
     &                              total = (/iend-istr+1/))
            ELSE IF (nvdim.eq.2) THEN
              IF (recordless) THEN
                npts=(iend-istr+1)*jrec
                CALL netcdf_get_fvar (ng, model, ncfile,                &
     &                                vname(1,ifield), awrk,            &
     &                                ncid = ncid,                      &
     &                                start = (/1,1/),                  &
     &                                total = (/iend-istr+1,jrec/))
              ELSE
                npts=iend-istr+1
                CALL netcdf_get_fvar (ng, model, ncfile,                &
     &                                vname(1,ifield), awrk,            &
     &                                ncid = ncid,                      &
     &                                start = (/1,trec/),               &
     &                                total = (/iend-istr+1,1/))
              END IF
            ELSE IF (nvdim.eq.3) THEN
              npts=(iend-istr+1)*jrec
              CALL netcdf_get_fvar (ng, model, ncfile,                  &
     &                              vname(1,ifield), awrk,              &
     &                              ncid = ncid,                        &
     &                              start = (/1,1,trec/),               &
     &                              total = (/iend-istr+1,jrec,1/))
            END IF
            IF (founderror(exit_flag, noerror, __line__, myfile)) THEN
              IF (master) WRITE (stdout,50) trim(vname(1,ifield)), trec
              RETURN
            END IF
            fmin=awrk(1)*fscale(ifield,ng)
            fmax=awrk(1)*fscale(ifield,ng)
            ic=0
            DO j=1,jrec
              DO i=istr,iend
                ic=ic+1
                aval=awrk(ic)*fscale(ifield,ng)
                fmin=min(fmin,aval)
                fmax=max(fmax,aval)
                fout(i,j,tindex)=aval
              END DO
            END DO
            finfo(8,ifield,ng)=fmin
            finfo(9,ifield,ng)=fmax
#ifdef CHECKSUM
            CALL get_hash (awrk, npts, fhash)
#endif
            IF (master) THEN
              IF (recordless) THEN
                WRITE (stdout,60) trim(vname(2,ifield)), ng, fmin, fmax
              ELSE
                lstr=scan(ncfile,'/',back=.true.)+1
                lend=len_trim(ncfile)
                tsec=tval*day2sec
                CALL time_string (tsec, t_code)
                WRITE (stdout,70) trim(vname(2,ifield)), t_code,        &
     &                            ng, trec, tindex, ncfile(lstr:lend),  &
     &                            tmin, tmax, tval, fmin, fmax
              END IF
#ifdef CHECKSUM
              WRITE (stdout,80) fhash
#endif
            END IF
            update=.true.
          END IF
        END IF
!
!  Increment the local time variable "Tmono" by the interval between
!  snapshots. If the interval is negative, indicating cycling, add in
!  a cycle length.  Load time value (sec) into "Tintrp" which used
!  during interpolation between snapshots.
!
        IF (.not.lonerec.and.(.not.recordless)) THEN
          tdelta=vtime(tindex,ifield,ng)-vtime(3-tindex,ifield,ng)
          IF (liocycle.and.(tdelta.lt.0.0_r8)) THEN
            tdelta=tdelta+clength
          END IF
          tmono=tmono+tdelta*day2sec
          finfo(7,ifield,ng)=tmono
          tintrp(tindex,ifield,ng)=tmono
        END IF
      END IF
!
  10  FORMAT (/,' GET_NGFLD_NF90 - unable to find dimension ',a,        &
     &        /,18x,'for variable: ',a,/,18x,'in file: ',a,             &
     &        /,18x,'file is not CF compliant...')
  20  FORMAT (/,' GET_NGFLD_NF90 - too small dimension for variable ',  &
     &        a,': ',i0,2x,i0)
  30  FORMAT (/,' GET_NGFLD_NF90 - unable to find requested variable:', &
     &        1x,a,/,18x,'in file: ',a)
  40  FORMAT (/,' GET_NGFLD_NF90 - unable to open input NetCDF',        &
     &        ' file: ',a)
  50  FORMAT (/,' GET_NGFLD_NF90 - error while reading variable: ',a,   &
     &        2x,' at TIME index = ',i0)
  60  FORMAT (2x,'GET_NGFLD_NF90   - ',a,/,22x,'(Grid = ',i2.2,         &
     &        ', Min = ',1pe15.8,' Max = ', 1pe15.8,')')
  70  FORMAT (2x,'GET_NGFLD_NF90   - ',a,',',t75,a,/,22x,               &
     &        '(Grid= ',i2.2,', Rec=',i0,', Index=',i1,                 &
     &        ', File: ',a,')',/,22x,                                   &
     &        '(Tmin= ', f15.4, ' Tmax= ', f15.4,')',                   &
     &        t71, 't = ', f15.4 ,/,22x,                                &
     &        '(Min = ', 1p,e15.8,0p,' Max = ',1p,e15.8,0p,')')
#ifdef CHECKSUM
  80  FORMAT (22x,'(CheckSum = ',i0,')')
#endif
!
      RETURN

References mod_ncparam::cinfo, mod_scalars::day2sec, mod_scalars::exit_flag, mod_ncparam::finfo, strings_mod::founderror(), mod_ncparam::fscale, get_hash_mod::get_hash(), mod_scalars::iic, mod_ncparam::iinfo, mod_ncparam::linfo, mod_parallel::master, mod_iounits::ncfile, mod_scalars::noerror, mod_scalars::ntstart, mod_scalars::rclock, mod_iounits::sourcefile, mod_iounits::stdout, mod_scalars::time, dateclock_mod::time_string(), mod_ncparam::tintrp, mod_ncparam::tname, mod_ncparam::vname, and mod_ncparam::vtime.

Referenced by get_ngfld().

Here is the call graph for this function:

Here is the caller graph for this function:

get_ngfld_pio()

subroutine get_ngfld_pio	(	integer, intent(in)	ng,
		integer, intent(in)	model,
		integer, intent(in)	ifield,
		type (file_desc_t), intent(inout)	piofile,
		integer, intent(in)	nfiles,
		type(t_io), dimension(nfiles), intent(inout)	s,
		logical, intent(in)	recordless,
		logical, intent(out)	update,
		integer, intent(in)	lbi,
		integer, intent(in)	ubi,
		integer, intent(in)	ubj,
		integer, intent(in)	ubk,
		integer, intent(in)	istr,
		integer, intent(in)	iend,
		integer, intent(in)	jrec,
		real(r8), dimension(lbi:ubi,ubj,ubk), intent(inout)	fout )

Definition at line 410 of file get_ngfld.F.

!***********************************************************************
!
      USE mod_param
      USE mod_parallel
      USE mod_iounits
      USE mod_ncparam
      USE mod_pio_netcdf
      USE mod_scalars
!
      USE dateclock_mod, ONLY : time_string
# ifdef CHECKSUM
      USE get_hash_mod,  ONLY : get_hash
# endif
      USE inquiry_mod,   ONLY : inquiry
      USE strings_mod,   ONLY : founderror
!
      implicit none
!
!  Imported variable declarations.
!
      logical, intent(in)  :: recordless
      logical, intent(out) :: update
!
      integer, intent(in) :: ng, model, ifield, nfiles
      integer, intent(in) :: LBi, UBi, UBj, UBk, Istr, Iend, Jrec
!
      TYPE (File_desc_t), intent(inout) :: pioFile
      TYPE(T_IO), intent(inout) :: S(nfiles)
!
      real(r8), intent(inout) :: Fout(LBi:UBi,UBj,UBk)
!
!  Local variable declarations.
!
      logical :: Linquire, Liocycle, Lmulti, Lonerec
!
      integer :: Nrec, Tindex, Trec, Vtype
      integer :: i, ic, j, job, lend, lstr, npts, nvdim, status
# ifdef CHECKSUM
      integer(i8b) :: Fhash
# endif
!
      real(r8) :: Aval, Fmax, Fmin
 
      real(dp) :: Clength, Tdelta, Tend
      real(dp) :: Tmax, Tmin, Tmono, Tscale, Tstr
      real(dp) :: Tsec, Tval
 
      real(r8), dimension((UBi-LBi+1)*UBj) :: Awrk
!
      character (len=22) :: t_code
 
      character (len=*), parameter :: MyFile =                          &
     &  __FILE__//", get_ngfld_pio"
!
      TYPE (My_VarDesc) :: TpioVar, VpioVar
!
      sourcefile=myfile
!
!-----------------------------------------------------------------------
!  Initialize.
!-----------------------------------------------------------------------
!
      IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
!
!  Determine if inquiring about the field to process in input NetCDF
!  file(s).  This usually happens on first call or when the field
!  time records are split (saved) in several multi-files.
!
      linquire=.false.
      lmulti=.false.
      IF (iic(ng).eq.0) linquire=.true.
      IF (.not.linquire.and.                                            &
     &    ((iinfo(10,ifield,ng).gt.1).and.                              &
     &     (linfo( 6,ifield,ng).or.                                     &
     &     (finfo( 2,ifield,ng)*day2sec.lt.time(ng))))) THEN
        linquire=.true.
        lmulti=.true.
      END IF
!
!  If appropriate, inquire about the contents of input NetCDF file and
!  fill information arrays.
!
!  Also, if appropriate, deactivate the Linfo(6,ifield,ng) switch after
!  the query for the UPPER snapshot interpolant from previous multifile
!  in the list. The switch was activated previously to indicate the
!  processing of the FIRST record of the file for the LOWER snapshot
!  interpolant.
!
      IF (linquire) THEN
        job=1
        CALL inquiry (ng, model, job, ubk, iend, ubi, ifield, piofile,  &
     &                lmulti, nfiles, s)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        IF (linfo(6,ifield,ng)) THEN
          linfo(6,ifield,ng)=.false.
        END IF
      END IF
!
!-----------------------------------------------------------------------
!  If appropriate, read in new data.
!-----------------------------------------------------------------------
!
      update=.false.
      tmono=finfo(7,ifield,ng)
!
      IF ((tmono.lt.time(ng)).or.(iic(ng).eq.0).or.                     &
     &    (iic(ng).eq.ntstart(ng))) THEN
!
!  Load properties for requested field from information arrays.
!
        liocycle=linfo( 2,ifield,ng)
        lonerec =linfo( 3,ifield,ng)
        vtype   =iinfo( 1,ifield,ng)
        vpiovar =dinfo( 1,ifield,ng)
        tpiovar =dinfo( 2,ifield,ng)
        nrec    =iinfo( 4,ifield,ng)
        tindex  =iinfo( 8,ifield,ng)
        trec    =iinfo( 9,ifield,ng)
        nvdim   =iinfo(11,ifield,ng)
        tmin    =finfo( 1,ifield,ng)
        tmax    =finfo( 2,ifield,ng)
        clength =finfo( 5,ifield,ng)
        tscale  =finfo( 6,ifield,ng)
        ncfile  =cinfo(ifield,ng)
!
        IF (liocycle) THEN
          trec=mod(trec,nrec)+1
        ELSE
          trec=trec+1
        END IF
        iinfo(9,ifield,ng)=trec
!
        IF (trec.le.nrec) THEN
!
!  Set rolling index for two-time record storage of input data.  If
!  "UBk" is unity, input data is stored in recordless array by the
!  calling program.
!
          IF (recordless) THEN
            tindex=1
          ELSE
            tindex=3-tindex
          END IF
          iinfo(8,ifield,ng)=tindex
!
!  Read in time coordinate.
!
          IF (.not.recordless.and.(tpiovar%vd%varID.ge.0).and.          &
     &        (tpiovar%vd%varID.ne.vpiovar%vd%varID)) THEN
            CALL pio_netcdf_get_time (ng, model, ncfile, tname(ifield), &
     &                                rclock%DateNumber, tval,          &
     &                                piofile = piofile,                &
     &                                start = (/trec/),                 &
     &                                total = (/1/))
            IF (founderror(exit_flag, noerror, __line__, myfile)) THEN
              IF (master) WRITE (stdout,50) trim(tname(ifield)), trec
              RETURN
            END IF
            tval=tval*tscale
            vtime(tindex,ifield,ng)=tval
!
!  Activate switch Linfo(6,ifield,ng) if processing the LAST record of
!  the file for the LOWER time snapshot. We need to get the UPPER time
!  snapshot from NEXT multifile.
!
            IF ((trec.eq.nrec).and.(tval*day2sec.le.time(ng))) THEN
              linfo(6,ifield,ng)=.true.
            END IF
          END IF
!
!  Read in non-grided data. The conditional statement on Jrec is to
!  differentiate between reading a 3D and 2D array.
!
          IF (vpiovar%vd%varID.ge.0) THEN
            fmin=0.0_r8
            fmax=0.0_r8
            IF (nvdim.eq.1) THEN
              npts=iend-istr+1
              CALL pio_netcdf_get_fvar (ng, model, ncfile,              &
     &                                  vname(1,ifield), awrk,          &
     &                                  piofile = piofile,              &
     &                                  start = (/1/),                  &
     &                                  total = (/iend-istr+1/))
            ELSE IF (nvdim.eq.2) THEN
              IF (recordless) THEN
                npts=(iend-istr+1)*jrec
                CALL pio_netcdf_get_fvar (ng, model, ncfile,            &
     &                                    vname(1,ifield), awrk,        &
     &                                    piofile = piofile,            &
     &                                    start = (/1,1/),              &
     &                                    total = (/iend-istr+1,jrec/))
 
              ELSE
                npts=iend-istr+1
                CALL pio_netcdf_get_fvar (ng, model, ncfile,            &
     &                                    vname(1,ifield), awrk,        &
     &                                    piofile = piofile,            &
     &                                    start = (/1,trec/),           &
     &                                    total = (/iend-istr+1,1/))
              END IF
            ELSE IF (nvdim.eq.3) THEN
              npts=(iend-istr+1)*jrec
              CALL pio_netcdf_get_fvar (ng, model, ncfile,              &
     &                                  vname(1,ifield), awrk,          &
     &                                  piofile = piofile,              &
     &                                  start = (/1,1,trec/),           &
     &                                  total = (/iend-istr+1,jrec,1/))
            END IF
            IF (founderror(exit_flag, noerror, __line__, myfile)) THEN
              IF (master) WRITE (stdout,50) trim(vname(1,ifield)), trec
              RETURN
            END IF
            fmin=awrk(1)*fscale(ifield,ng)
            fmax=awrk(1)*fscale(ifield,ng)
            ic=0
            DO j=1,jrec
              DO i=istr,iend
                ic=ic+1
                aval=awrk(ic)*fscale(ifield,ng)
                fmin=min(fmin,aval)
                fmax=max(fmax,aval)
                fout(i,j,tindex)=aval
              END DO
            END DO
            finfo(8,ifield,ng)=fmin
            finfo(9,ifield,ng)=fmax
# ifdef CHECKSUM
            CALL get_hash (awrk, npts, fhash)
# endif
            IF (master) THEN
              IF (recordless) THEN
                WRITE (stdout,60) trim(vname(2,ifield)), ng, fmin, fmax
              ELSE
                lstr=scan(ncfile,'/',back=.true.)+1
                lend=len_trim(ncfile)
                tsec=tval*day2sec
                CALL time_string (tsec, t_code)
                WRITE (stdout,70) trim(vname(2,ifield)), t_code,        &
     &                            ng, trec, tindex, ncfile(lstr:lend),  &
     &                            tmin, tmax, tval, fmin, fmax
              END IF
# ifdef CHECKSUM
              WRITE (stdout,80) fhash
# endif
            END IF
            update=.true.
          END IF
        END IF
!
!  Increment the local time variable "Tmono" by the interval between
!  snapshots. If the interval is negative, indicating cycling, add in
!  a cycle length.  Load time value (sec) into "Tintrp" which used
!  during interpolation between snapshots.
!
        IF (.not.lonerec.and.(.not.recordless)) THEN
          tdelta=vtime(tindex,ifield,ng)-vtime(3-tindex,ifield,ng)
          IF (liocycle.and.(tdelta.lt.0.0_r8)) THEN
            tdelta=tdelta+clength
          END IF
          tmono=tmono+tdelta*day2sec
          finfo(7,ifield,ng)=tmono
          tintrp(tindex,ifield,ng)=tmono
        END IF
      END IF
!
  10  FORMAT (/,' GET_NGFLD_PIO - unable to find dimension ',a,         &
     &        /,17x,'for variable: ',a,/,17x,'in file: ',a,             &
     &        /,17x,'file is not CF compliant...')
  20  FORMAT (/,' GET_NGFLD_PIO - too small dimension for variable ',   &
     &        a,': ',i0,2x,i0)
  30  FORMAT (/,' GET_NGFLD_PIO - unable to find requested variable:',  &
     &        1x,a,/,18x,'in file: ',a)
  40  FORMAT (/,' GET_NGFLD_PIO - unable to open input NetCDF',         &
     &        ' file: ',a)
  50  FORMAT (/,' GET_NGFLD_PIO - error while reading variable: ',a,    &
     &        2x,' at TIME index = ',i0)
  60  FORMAT (2x,'GET_NGFLD_PIO    - ',a,/,22x,'(Grid = ',i2.2,         &
     &        ', Min = ',1pe15.8,' Max = ', 1pe15.8,')')
  70  FORMAT (2x,'GET_NGFLD_PIO    - ',a,',',t75,a,/,22x,               &
     &        '(Grid= ',i2.2,', Rec=',i0,', Index=',i1,                 &
     &        ', File: ',a,')',/,22x,                                   &
     &        '(Tmin= ', f15.4, ' Tmax= ', f15.4,')',                   &
     &        t71, 't = ', f15.4 ,/,22x,                                &
     &        '(Min = ', 1p,e15.8,0p,' Max = ',1p,e15.8,0p,')')
# ifdef CHECKSUM
  80  FORMAT (22x,'(CheckSum = ',i0,')')
# endif
!
      RETURN

References mod_ncparam::cinfo, mod_scalars::day2sec, mod_ncparam::dinfo, mod_scalars::exit_flag, mod_ncparam::finfo, strings_mod::founderror(), mod_ncparam::fscale, get_hash_mod::get_hash(), mod_scalars::iic, mod_ncparam::iinfo, mod_ncparam::linfo, mod_parallel::master, mod_iounits::ncfile, mod_scalars::noerror, mod_scalars::ntstart, mod_scalars::rclock, mod_iounits::sourcefile, mod_iounits::stdout, mod_scalars::time, dateclock_mod::time_string(), mod_ncparam::tintrp, mod_ncparam::tname, mod_ncparam::vname, and mod_ncparam::vtime.

Referenced by get_ngfld().

Here is the call graph for this function:

Here is the caller graph for this function: