* PACKAGE GTINTFC  !" GTOOL3 実数精度変換インターフェース(ZFREAD)
*
*"  [HIS] 96/03/08(takepiro) 自動倍精度コンパイル禁止
*
**********************************************************************
*"      <<  GTOOL3標準フォーマットデータ読み込み >>
**********************************************************************
      SUBROUTINE ZFREAD
     O            ( HHEAD , GDATA , IEOD  ,
     I              IFILE , MODE               )
*
      CHARACTER  HHEAD ( * )*(*)         !" ヘッダー
#ifdef AUTO_DBL
      REAL*8     GDATA ( * )             !" データ
#else
      REAL*4     GDATA ( * )             !" データ
#endif
      INTEGER    IEOD                    !" データ存在フラグ
      INTEGER    IFILE                   !" ファイル番号
      INTEGER    MODE                    !" モード(1なら表示ON)
*
      CALL  ZFRDSL
     O         ( HHEAD , GDATA , IEOD  ,
     I           '###' , '###' , -1    ,
     I           IFILE , MODE           )
*
      RETURN
      END
**********************************************************************
*"      <<  GTOOL3標準フォーマットデータ読み込み: 限定つき
**********************************************************************
      SUBROUTINE ZFRDSL
     O         ( HHEAD , GDATA , IEOD  ,
     I           HPSEL , HSEL  , ISEL  ,
     I           IFILE , MODE           )
*
      CHARACTER  HHEAD ( * )*(*)         !" ヘッダー
#ifdef AUTO_DBL
      REAL*8     GDATA ( * )             !" データ
#else
      REAL*4     GDATA ( * )             !" データ
#endif
      INTEGER    IEOD                    !" データ存在フラグ
      CHARACTER  HPSEL *(*)              !" 選択する欄
      CHARACTER  HSEL  *(*)              !" 選択するデータ(文字)
*"   (INTEGER)   ISEL                    !  選択するデータ(数)
      INTEGER    ISEL                    !  選択するデータ(数)
      INTEGER    IFILE                   !" ファイル番号
      INTEGER    MODE                    !" モード(1なら表示ON)
*
* [PARAM]
#include        "zcdim.F"            !" 格子点数, 波数
#include        "zhdim.F"            !" 文字列文字数
*
      INTEGER    IJKDIM
        PARAMETER  ( IJKDIM = (NX+2)*(NZ+2) )
      INTEGER    IJK                         !" ループカウンター
      INTEGER    IJKSIZ                      !" 配列の大きさ
*
      CHARACTER  HDFMTI  *(NCC)          !" データフォーマット
*
      REAL*8     DDATA ( IJKDIM )
      REAL*4     SDATA ( IJKDIM )
      EQUIVALENCE (DDATA,SDATA)
*
      CALL GFRDSL
     O     ( HHEAD , DDATA , IEOD  ,
     I       HPSEL , HSEL  , ISEL  ,
     I       IFILE , MODE           )
*
      IF ( IEOD .EQ. 0 ) THEN
*
*" << CNVDBL : 精度変換 >>
*
         CALL GHCGET ( HHEAD, 'DFMT',  HDFMTI )
         CALL GUQTSZ ( HHEAD, IJKSIZ )
*
         IF ( HDFMTI .EQ. 'UR8' ) THEN
            DO 1110 IJK = 1, IJKSIZ
                  GDATA( IJK ) = DDATA( IJK )
 1110       CONTINUE
         ELSE
            DO 1120 IJK = 1, IJKSIZ
                  GDATA( IJK ) = SDATA( IJK )
*                            equiv DDATA
 1120       CONTINUE
         ENDIF
      ENDIF
*
      RETURN
      END