プログラムのエラー（ファイルの出力がうまくいかない）

Question

###前提・実現したいこと
2次元拡散方程式のプログラムを作成しています。

計算結果をコンター図で表すため，それぞれの時刻と位置における値を出力したいです。

###発生している問題・エラーメッセージ
csvファイルは作成されているが、データが全く入っていない。

###該当のソースコード
Public Class Form1
Dim b, h, wcr, ic, bc, dx, dy, dt, t, tmax, rx, ry As Double
Dim i, j, l, m, n, p, q, r, s, x, y As Integer
Dim c0(b, tmax), c1(x, y), c2(x, y) As Double
Dim FileNum As Short
Dim FileOut As String

Private Sub Button1_Click(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles Button1.Click
    Call Cl_ion_Calculate()
End Sub

Sub Cl_ion_Calculate()

    Dim i, j, l, m, n, p, q, r, s, x, y As Integer

    b = Val(TextBox1.Text)
    h = Val(TextBox2.Text)
    wcr = Val(TextBox3.Text)
    dx = Val(TextBox4.Text)
    dy = Val(TextBox5.Text)
    dt = Val(TextBox6.Text)
    ic = Val(TextBox7.Text)
    bc = Val(TextBox8.Text)
    tmax = Val(TextBox9.Text)

    t = tmax / dt
    x = b / dx
    y = h / dy
    rx = dt / dx ^ 2
    ry = dt / dy ^ 2

    ReDim c0(b, tmax), c1(x, y), c2(x, y)


    For i = 0 To x
        For j = 0 To y
            c1(i, j) = ic
        Next j
    Next i

    For l = 1 To t
        For m = 1 To x
            c2(m, 0) = bc
        Next m
        For n = 1 To y
            c2(0, n) = bc
            c2(x, n) = bc
        Next n

        For p = 1 To x - 1
            For q = 1 To y - 1
                c2(p, q) = (1 - 2 * rx) * c2(p, q) + rx * c2(p - 1, q) + rx * c2(p + 1, q) + (1 - 2 * ry) * c2(p, q) + ry * c2(p, q - 1) + ry * c2(p, q + 1)
            Next
        Next

        For r = 0 To x
            For s = 0 To y
                c1(r, s) = c2(r, s)
            Next
        Next
    Next
    Call FOpen()
    Call FSave()

    MsgBox("Finish")

End Sub
Sub FOpen()

    FileOut = "D:result.csv"
    FileNum = FreeFile()
    FileOpen(FileNum, FileOut, OpenMode.Output)

End Sub
Sub FSave()

    Dim c2 As Object

    For r = 1 To x
        For s = 1 To y
            WriteLine(FileNum, l, r, s, c2(r, s))
        Next
    Next

End Sub

End Class

Answer 1

あまり内容を理解せずに回答しますので、間違ってたらすみません。
FSave関数の中でc2配列の中身を出力しようとしていますね。
ループの回数はxとyですが、これはCl_ion_Calculate関数の中で算出されています。
が、変数xとyはグローバル変数で宣言されているにも関わらず、Cl_ion_Calculate関数内でもローカル変数として宣言されてしまっています。
Cl_ion_Calculate関数内ではローカル変数の方が優先度は高いですから、グローバル変数のxとyに値が設定されていない可能性ああります。
Cl_ion_Calculate関数内の変数宣言を削除して実行してみてください。
この行。

Sub Cl_ion_Calculate()
'Dim i, j, l, m, n, p, q, r, s, x, y As Integer   コメントアウト

Answer 2

以下のようなイメージではないかと思うのですが、いかがでしょうか？
一旦ans(,)に入れているのは私の好みですので、直接l,r,s,c2(r,s)をWritelineで書き込んでも問題ありません。

Sub Cl_ion_Calculate()

        Dim i, j, l, m, n, p, q, r, s, x, y As Integer
        Dim ans(,) As Double
        FileOut = "D:result.csv"
        FileNum = FreeFile() 
        FileOpen(FileNum, FileOut, OpenMode.Output)


        b = Val(TextBox1.Text)
        h = Val(TextBox2.Text)
        wcr = Val(TextBox3.Text)
        dx = Val(TextBox4.Text)
        dy = Val(TextBox5.Text)
        dt = Val(TextBox6.Text)
        ic = Val(TextBox7.Text)
        bc = Val(TextBox8.Text)
        tmax = Val(TextBox9.Text)

        t = tmax / dt
        x = b / dx
        y = h / dy
        rx = dt / dx ^ 2
        ry = dt / dy ^ 2

        ReDim ans(3, t)

        ReDim c0(b, tmax), c1(x, y), c2(x, y)

        't=0
        For i = 0 To x
            For j = 0 To y
                c1(i, j) = ic '初期値
            Next j
        Next i

        For l = 1 To t
            ans(0, l) = l
            '境界条件
            For m = 1 To x
                c2(m, 0) = bc
            Next m
            For n = 1 To y
                c2(0, n) = bc
                c2(x, n) = bc
            Next n

            For p = 1 To x - 1
                For q = 1 To y - 1
                    c2(p, q) = (1 - 2 * rx) * c2(p, q) + rx * c2(p - 1, q) + rx * c2(p + 1, q) + (1 - 2 * ry) * c2(p, q) + ry * c2(p, q - 1) + ry * c2(p, q + 1)
                Next
            Next

            For r = 0 To x
                ans(1, l) = r
                For s = 0 To y
                    ans(2, l) = s
                    c1(r, s) = c2(r, s)
                    ans(3, l) = c2(r, s)
                    WriteLine(FileNum, ans(0, l), ans(1, l), ans(2, l), ans(3, l))
                Next s
            Next r
        Next l
        'Call FOpen()
        'Call FSave()

        MsgBox("Finish")

    End Sub

なお、初期値として恐らく適当ではないと思いますが、以下のような初期値、境界条件で実行すると、計算開始1time step、2time step後に、以下のようなプロットが得られます。