* PACKAGE XMINIT    !" $B=i4|CM=PNO(B
*
*"  [HIS] 96/03/13 (takepiro) 
***********************************************************************
      PROGRAM XMKINT
*
* [PARAM]
#include        "zcdim.F"            !" $B3J;RE@?t(B, $BGH?t(B
#include        "zhdim.F"            !" $BJ8;zNsJ8;z?t(B
*
* [VAR]
*" $B3J;RE@%G!<%?(B
      REAL PSI  ( 0:NX+1 , 0:NZ+1 ) !" $BN.@~4X?t(B
      REAL ZETA ( 0:NX+1 , 0:NZ+1 ) !" $B12EY(B
      REAL T    ( 0:NX+1 , 0:NZ+1 ) !" $B29EY(B
*
      REAL DX                   !" X $B<4>e$N3J;RE@4V3V(B
      REAL DZ                   !" Z $B<4>e$N3J;RE@4V3V(B
*
      CHARACTER  RUN  *(NCC)                 !" $B<B83L>(B
      CHARACTER  SIGN *(NCC)                 !" $B<B83<TL>(B
*
      INTEGER IX, IZ
      REAL XX, ZZ, PI
      INTEGER IT                !" $B8=:_%9%F%C%W?t(B
      REAL    DT                !" $B@QJ,;~4V4V3V(B($B%@%_!<(B)
*
      DATA   RUN, SIGN / 'TEST', 'momoko' /
      DATA IT, DT / 0, 1.0 /
*
*"    < 1. $B=i4|CM7W;;(B >
*
      PI = ATAN2( 0.0, -1.0 )
*
      CALL ATSTRT               !" $B<B83=i4|@_Dj(B($B;~4V(B, $B:BI8(B)
     O     ( IT , DT )
*
      CALL ACRGET               !" $B:BI8CM(B
     O     (  DX , DZ  )
*
      DO 1000 IX = 0, NX+1
         DO 1100 IZ = 0, NZ+1 
            XX = DX*(IX-1)
            ZZ = DZ*IZ
*
c$$$            PSI  ( IX , IZ ) = COS( XX ) * SIN( PI*ZZ )
c$$$     &           + ZZ
c$$$            PSI  ( IX , IZ ) = 0.0
c$$$            PSI  ( IX , IZ ) = 10*(1.0 - ZZ)
c$$$            ZETA  ( IX , IZ ) = 0.0
c$$$            T    ( IX , IZ ) = 1.0
            T    ( IX , IZ ) = 1.0 - ZZ
c$$$            T    ( IX , IZ ) = COS( XX ) * SIN( PI*ZZ )
c$$$     &           + ZZ
 1100    CONTINUE 
 1000 CONTINUE 
      T( NX/2+1, 1 ) = T( NX/2+1, 1 ) + 0.1 
*
      CALL DBNSET               !" $B6-3&>r7o$N@_Dj(B
     I     ( PSI ,              !" $BN.@~4X?t$KBP$7$FMQ$$$k(B
     I       ZETA, T  ,
     I       DX  , DZ   )
*
      CALL DLAPLA               !" PSI $B$rM?$($?;~(B -> ZETA $B$r7W;;(B
     O     ( ZETA ,
     I       PSI  ,
     C       DX  , DZ  )
*
      CALL DBNDRZ               !" $B1?F03X(B/$BNO3X6-3&>r7o$NE,MQ$HN.@~4X?t$N7W;;(B
     O       ( PSI  , 
     M         ZETA , 
     I         DX   , DZ  )
*
*"    < 2. $B%j%9%?!<%H%U%!%$%k=PNO(B >
*
      CALL WRRSTR               !" $B%j%9%?!<%H%U%!%$%k=q$-9~$_(B
     I     ( PSI  , ZETA , T   , 
     I        IT                    )
*
      STOP
      END
