exm_5.bas

来自「VB版本_SIMPLE算法和11各个算例」· BAS 代码 · 共 83 行

Attribute VB_Name = "EXM_5"
'*****************************************************************'
'
'      DIMENSION U(22, 22), V(22, 22), PC(22, 22)
'      EQUIVALENCE (F(1,1,1),U(1,1)),(F(1,1,2),V(1,1)),(F(1,1,3),PC(1,1))
'      DIMENSION TH(22), THU(22), THDIF(22), THCV(22), THCVS(22)
 '     EQUIVALENCE(X,TH),(XU,THU),(XDIF,THDIF),(XCV,THCV)
'     &  ,(XCVS,THCVS),(XL,THL)
'*****************************************************************
'*
'*----------------- EXAMPLE FIVE -----------------------------
'*'
'*****************************************************************
'      DIMENSION T(22, 22)
'      EQUIVALENCE (F(1,1,4),T(1,1))

Option Explicit
Public AMU As Double
Public DEN As Double, WBAR As Double, COND As Double, CP As Double, RHOCP As Double
Public DPDZ As Double, DTBDZ As Double
Public ASUM#, TSUM#, WSUM#, TW#, AR#, WR#, TB#, DH#, RE#, FRE#, QW#, ANU#
Dim j As Long, i As Long
Dim THEAT#, DTDZ#, PI#, WP#, TBULK#, HTP#, HTC#
Sub GRID()
      LSolve(4) = True
      Title(1) = ".VEL  U."
      Title(2) = ".VEL  V."
      Title(3) = ".STR FN."
      Title(4) = ". TEMP ."
      LPrint(1) = True
      LPrint(2) = True
      LPrint(3) = True
      LPrint(4) = True
      XL = 1#
      YL = 1#
      L1 = 7
      M1 = 7
      Call UGRID
 End Sub

      Sub START()
      For j = 1 To M1
          For i = 1 To L1
               F(i, j, 1) = 10# * CDbl(i - 2) / CDbl(L1 - 2)
               F(i, j, 2) = -10# * CDbl(j - 2) / CDbl(M1 - 2)
               F(i, j, 4) = 50#
               F(1, j, 4) = 100#
          Next
      Next
      End Sub

      Sub DENSE()
      End Sub

      Sub BOUND()
      For j = 2 To M2
           F(L1, j, 4) = F(L2, j, 4)
      Next j
      For i = 2 To L2
           F(i, 1, 4) = F(i, 2, 4)
      Next i
      End Sub
      Sub OUTPUT()
      If (ITER = 0) Then
        Print #1,
        Print #1, "   ITER"; Space(8); "T(4,5)"; Space(8); "T(5,3)"
      End If
       Print #1,
       Print #1, Format(ITER, "0#     "); Format(F(4, 5, 4), "##.00000"), Format(F(5, 3, 4), "##.00000")
      If (ITER = LAST) Then
          Call ToPrint
      End If
      End Sub
      Sub GAMSOR()
        For j = 1 To M1
              For i = 1 To L1
                   GAM(i, j) = 1#
                   GAM(i, 1) = 0#
                   GAM(L1, j) = 0#
               Next i
      Next j
      End Sub