lengds.f

Go to the documentation of this file.

C> @file
C> @brief GIven a grid description section (in w3fi63 format),
C> return its size in terms of number of data points.
C> @author Mark Iredell @date 1996-07-19
 
C> Program history log:
C> - Mark Iredell 1996-07-19
C>
C> @param[in] KGDS Integer (200) gds parameters in w3fi63 format.
C> @return LENGDS Integer size of grid.
C>
C> @author Mark Iredell @date 1996-07-19
C-----------------------------------------------------------------------
      FUNCTION lengds(KGDS)
      INTEGER kgds(200)
C - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C  SPECIAL CASE OF STAGGERED ETA
      IF(kgds(1).EQ.201) THEN
        lengds=kgds(7)*kgds(8)-kgds(8)/2
C  SPECIAL CASE OF FILLED ETA
      ELSEIF(kgds(1).EQ.202) THEN
        lengds=kgds(7)*kgds(8)
C  SPECIAL CASE OF THINNED WAFS
      ELSEIF(kgds(19).EQ.0.AND.kgds(20).NE.255) THEN
        lengds=kgds(21)
C  GENERAL CASE
      ELSE
        lengds=kgds(2)*kgds(3)
      ENDIF
C - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
      END