module1.bas

VB版本_SIMPLE算法和11各个算例

Attribute VB_Name = "Exm_2"

'*****************************************************************
'*
'*-------------- EXAMPLE TWO --------------------------------
'*
'*****************************************************************
Option Explicit
Public Q As Double, GAM1 As Double, GY As Double, HTC As Double, TF As Double
Dim METHOD%, i%, j%, RES#, ARES#
 Sub GRID()
      LAST = 100
      Title(4) = "TEMP"
      Title(13) = "COND"
      LPrint(4) = True
      LSolve(4) = True
      LPrint(13) = True
      MODE = 2
      XL = 2#
      YL = 1#
      R(1) = 1#
      L1 = 7
      M1 = 7
      Call UGRID
End Sub
Sub START()
      METHOD = 1
      For j = 1 To M1
           For i = 1 To L1
                f(i, j, 4) = 200#
                If (i = 1) Then
                     f(i, j, 4) = 100# * (1# + Y(j))
                End If
           Next i
      Next j
      HTC = 5#
      Q = 50#
      TF = 100#
      GAM1 = 1#
End Sub

Sub DENSE()

End Sub
Sub BOUND()
      For i = 2 To L2
            f(i, M1, 4) = f(i, M2, 4)
            f(i, 1, 4) = f(i, 2, 4) + Q * YDIF(2) / GAM1
      Next i
      GY = GAM1 / XDIF(L1)
      For j = 2 To M2
            f(L1, j, 4) = (HTC * TF + GY * f(L2, j, 4)) / (HTC + GY)
      Next j
End Sub
Sub OUTPUT()
         If (ITER = 0) Then
             Print #1, "   Method="; METHOD
             Print #1,
             Print #1, "     ITER"; Space(11); "T(4,5)"; Space(14); "T(5,3)"
          End If
          Print #1,
          Print #1, Format(ITER, "######"), f(4, 5, 4), f(5, 3, 4)
          If (ITER = LAST) Then Call ToPrint
End Sub
Sub GAMSOR()
      For j = 1 To M1
          For i = 1 To L1
              GAM(i, j) = GAM1
          Next i
      Next j
      For j = 4 To 7
          For i = 4 To 5
               GAM(i, j) = 0.2 * (1# + f(i, j, 4) / 100#)
          Next i
      Next j
      For j = 2 To M2
           For i = 2 To L2
              CON(i, j) = 100#: AP(i, j) = -0.5
           Next i
      Next j
      If (METHOD = 1) Then Exit Sub
      For i = 2 To L2
          GAM(i, M1) = 0#
          GAM(i, 1) = 0#
          CON(i, 2) = CON(i, 2) + Q * R(1) / ARX(2)
      Next i
      RES = 1# / HTC + 1# / GY
      ARES = 1# / (RES * XCV(L2))
      For j = 2 To M2
          GAM(L1, j) = 0#
          CON(L2, j) = CON(L2, j) + ARES * TF
          AP(L2, j) = AP(L2, j) - ARES
      Next j
End Sub