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:

OUTPUT:

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.

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.