grids.inc

来自「shpf 1.9一个并行编译器」· INC 代码 · 共 45 行

!----------------------------------------------------------
! Used in 'mg00d ()', 'init ()', 'mgi ()', getdef ()' & 'getdif ()'.
!----------------------------------------------------------
  INTEGER, PARAMETER :: n1 = 2**1
  INTEGER, PARAMETER :: n2 = 2**2
  INTEGER, PARAMETER :: n3 = 2**3
  INTEGER, PARAMETER :: n4 = 2**4
  INTEGER, PARAMETER :: n5 = 2**5
  INTEGER, PARAMETER :: n6 = 2**6
  INTEGER, PARAMETER :: n7 = 2**7
  INTEGER, PARAMETER :: n8 = 2**8
  INTEGER, PARAMETER :: n9 = 2**9

  INCLUDE "tmplg.inc"    ! declares templates 'tmpl1/2/3...'
  
  DOUBLE PRECISION, DIMENSION (n1+1, n1+1) :: u1, f1
!HPF$ ALIGN WITH tmpl1 :: u1, f1

  DOUBLE PRECISION, DIMENSION (n2+1, n2+1) :: u2, f2
!HPF$ ALIGN WITH tmpl2 :: u2, f2

  DOUBLE PRECISION, DIMENSION (n3+1, n3+1) :: u3, f3
!HPF$ ALIGN WITH tmpl3 :: u3, f3

  DOUBLE PRECISION, DIMENSION (n4+1, n4+1) :: u4, f4
!HPF$ ALIGN WITH tmpl4 :: u4, f4

  DOUBLE PRECISION, DIMENSION (n5+1, n5+1) :: u5, f5
!HPF$ ALIGN WITH tmpl5 :: u5, f5

  DOUBLE PRECISION, DIMENSION (n6+1, n6+1) :: u6, f6
!HPF$ ALIGN WITH tmpl6 :: u6, f6

  DOUBLE PRECISION, DIMENSION (n7+1, n7+1) :: u7, f7
!HPF$ ALIGN WITH tmpl7 :: u7, f7

  DOUBLE PRECISION, DIMENSION (n8+1, n8+1) :: u8, f8
!HPF$ ALIGN WITH tmpl8 :: u8, f8

  DOUBLE PRECISION, DIMENSION (n9+1, n9+1) :: u9, f9
!HPF$ ALIGN WITH tmpl9 :: u9, f9

  COMMON /grids/ u1, u2, u3, u4, u5, u6, u7, u8, u9,        &
                 f1, f2, f3, f4, f5, f6, f7, f8, f9