* PACKAGE ASTRT    !" 初期値の設定(簡単版)
*
*"  [HIS] 95/08/14 (takepiro) 
***********************************************************************
      SUBROUTINE RDSTRT                      !" 初期値の設定
     O           ( PSI2  , ZETA2 , T2 ,
     O             PSI1  , ZETA1 , T1  )
*
* [PARAM]
#include        "zcdim.F"            !" 格子点数, 波数
*
* [OUTPUT]
      REAL PSI2  ( 0:NX+1 , 0:NZ+1 ) !" 流線関数(現在の値)
      REAL ZETA2 ( 0:NX+1 , 0:NZ+1 ) !" 渦度(現在の値)
      REAL T2    ( 0:NX+1 , 0:NZ+1 ) !" 温度(現在の値)
*
      REAL PSI1  ( 0:NX+1 , 0:NZ+1 ) !" 流線関数(1STEP前の値)
      REAL ZETA1 ( 0:NX+1 , 0:NZ+1 ) !" 渦度(1STEP前の値)
      REAL T1    ( 0:NX+1 , 0:NZ+1 ) !" 温度(1STEP前の値)
*
      INTEGER IX, IZ
*
      REAL    PI, XX, ZZ
*
      PI = ATAN2( 0.0, -1.0 )
*
      DO 1100 IZ = 0, NZ+1
         DO 1100 IX = 0, NX+1
           XX = 2*PI*(IX-1)/REAL(NX)
           ZZ = IZ / REAL(NZ+1)
           PSI2 ( IX,IZ ) = 0.0
           ZETA2( IX,IZ ) = 0.0
c$$$           T2   ( IX,IZ ) = COS( XX ) * SIN( PI*ZZ )
           T2   ( IX,IZ ) = - ZZ 
 1100 CONTINUE 
      T2( NX/2,1 ) =  T2( NX/2,1 ) + 0.001
*
      CALL COPY( PSI1  , PSI2  , (NX+2)*(NZ+2) )
      CALL COPY( ZETA1 , ZETA2 , (NX+2)*(NZ+2) )
      CALL COPY( T1    , T2    , (NX+2)*(NZ+2) )
*
      RETURN
      END