Running the code with different initial polar fields

The code was run with different initial polar fields.

1. Initial field of 10 Gauss

OUTPUT:


Peak initial field for the 3rd cycle is -7.55G in the northern hem which occurs at 57.1233 degree lat.@25.3yrs.
The lowest that occurs on the same latitude is -0.86G.@14.2yrs.
Difference=6.7G

So, the cycles go like:

    Year                         Polar field (Gauss)
    14.2                               -0.86
    25.3                               -7.55
    36.4                               -0.86
    47.5                               -7.55


2. Initial field of 9 Gauss

OUTPUT:


Peak initial field for 3rd cycle is -6.7827G in the northern hem which occurs at 57.1233 degree lat.@25.2133yrs.
The lowest that occurs on the same latitude is -0.10699G.@14.1321yrs. and -0.15G.@35.8679yrs.
Difference=6.68G

So, the cycle goes like:

  Year                         Polar field (Gauss)

    14.1                               -0.11

    25.2                               -6.78

    35.9                               -0.15

    47.5                               -6.82

3. Initial field of 8 Gauss

OUTPUT:

Peak initial field for 3rd cycle is -6.06 G in the northern hem which occurs at 57.1233 degree lat.@24.9yrs.

The lowest that occurs on the same latitude is +0.61 G.@14.1yrs. and -0.61G.@35.9yrs.

Difference=6.67G

So, the cycle goes like:

  Year                         Polar field (Gauss)

    14.1                               +0.61

    24.9                               -6.06

    35.9                               +0.61

    46.7                               -6.06

4. Initial field of 7 Gauss:

OUTPUT:

  Year                         Polar field (Gauss)

    14.1                               +1.41

    25.3                               -5.26

    35.9                               +1.38

    47.3                               -5.28

    58.0                              +1.36

    68.8                               -5.31

Difference=6.67G

5. Initial field of 6 Gauss:

OUTPUT:

Year                         Polar field (Gauss)

    14.1                               +2.16

    25.3                               -4.5

    35.9                               +2.16

    47.3                               -4.5

Difference=6.66G


For realistic result, the difference should be around 20 Gauss.

The input parameters that can be varied to achieve this are:
(1) No. of sunspots per cycle (These were modeled vaguely and hence remain doubtful).
(2) Separation between sunspots within a BMR(taken to be a constant=2R)
(3) The longitudes of eruption are totally random. The degree of randomness can be optimized.

Running the code with decreased randomization for calibration

To relatively stabilize the peak polar magnetic field, the input of one sunspot cycle was repeatedly fed into the surface flux transport code with changing polarity of magnetic field. The initial polar field is 4.5 Gauss above 60 degrees latitude.

The butterfly diagram of the simulation is :

The peak polar magnetic field values are:

Cycle          Peak polar magnetic field in Northern Hem        Peak polar magnetic field in Southern Hem
   1                                 +3.30                                                                         -3.15
   2                                  -3.35                                                                        +3.39
   3                                 +3.30                                                                         -3.15
   4                                  -3.35                                                                        +3.39
   5                                 +3.30                                                                         -3.15
   6                                  -3.35                                                                        +3.39
   7                                 +3.30                                                                         -3.15