*"表題   BARO 共通 定数(座標値) GCM5(ACRGET)
*
*"履歴   90/10/15 沼口  敦
*"       92/04/02 竹広真一
*"       96/09/14 竹広真一
*
***********************************************************************
*"        << モデル座標定数の計算 >>
***********************************************************************
      SUBROUTINE ACRGET
     O         ( ALON  , DLON  ,
     O           ALAT  , DLAT   )
*
#if   SYS_IBMS
      INCLUDE   (ZCDIM)                  !" 格子点数, 波数
#else
#include        "zcdim.F"                    !" 格子点数, 波数
#endif
*
      REAL       ALON  ( IDIM )              !" 経度
      REAL       DLON  ( IDIM )              !" 経度荷重
      REAL       ALAT  ( JDIM )              !" 緯度
      REAL       DLAT  ( JDIM )              !" 緯度荷重
*
*"          < 1. 座標値のセット >
*
      CALL ACLON
     O         ( ALON , DLON )
*
      CALL ACLAT
     O         ( ALAT , DLAT )
*
      RETURN
      END
**********************************************************************
      SUBROUTINE ACLON
     O            ( ALON , DLON )
*
#if   SYS_IBMS
      INCLUDE   (ZCDIM)                  !" 格子点数, 波数
#else
#include        "zcdim.F"                    !" 格子点数, 波数
#endif
*
*  [OUTPUT]
      REAL       ALON  ( IDIM )              !" 経度
      REAL       DLON  ( IDIM )              !" 経度荷重
*
*  [INTERNAL SAVE]
      REAL       ALONZ ( IDIM )              !" 経度
      REAL       DLONZ ( IDIM )              !" 経度荷重
      LOGICAL    OFIRST
      DATA       OFIRST /.TRUE./
      SAVE
*
      IF ( OFIRST )THEN
         OFIRST = .FALSE.
*
         CALL SETLON
     O         ( ALONZ  , DLONZ  , 
     D           IMAX   , IDIM      )
      ENDIF
*
      CALL COPY ( ALON , ALONZ  , IDIM   )
      CALL COPY ( DLON , DLONZ  , IDIM   )
*
      RETURN
      END
**********************************************************************
      SUBROUTINE ACLAT
     O           ( ALAT , DLAT )
*
#if   SYS_IBMS
      INCLUDE   (ZCDIM)                  !" 格子点数, 波数
#else
#include        "zcdim.F"                    !" 格子点数, 波数
#endif
*
*  [OUTPUT]
      REAL       ALAT  ( JDIM )              !" 緯度
      REAL       DLAT  ( JDIM )              !" 緯度荷重
*
      REAL * 8   QSINLA( JMAX )              !" sin（緯度）：倍精度
      REAL * 8   QGW   ( JMAX )              !" ガウス荷重 ：倍精度
      REAL * 8   QPN( 0:JMAX+1 )            !" Ｐn (倍精度)
*
*  [INTERNAL SAVE]
      REAL       ALATZ ( JDIM )              !" 緯度
      REAL       DLATZ ( JDIM )              !" 緯度荷重
*
      LOGICAL    OFIRST
      DATA       OFIRST /.TRUE./
      SAVE
*
      IF ( OFIRST )THEN
         OFIRST = .FALSE.
*
         CALL SETLAT
     O         ( ALATZ  , DLATZ  , 
     D           JMAX   , JDIM   ,  
     W           QSINLA , QGW   , QPN   )
      ENDIF
*
      CALL COPY ( ALAT , ALATZ  , JDIM   )
      CALL COPY ( DLAT , DLATZ  , JDIM   )
*
      RETURN
      END