Listing of a GIC run

 Calculate GICs in a given power grid.
 
 Note: if you have modified a data file "fyle" of ASCII format
 you must remove the corresponding "fyle.mat" file.
 
 
----- Choose the network data file -----
 
      1) Finland 1978-1979 (400 kV)       
      2) Finland 1979-1983 (400 kV)       
      3) Finland 1991 (400 kV)            
      4) Finland 1991a (220 kV and 400 kV)
      5) Finland 1991b (220 kV and 400 kV)
      6) Finland 1991c (220 kV and 400 kV)
      7) Finland 1992 (220 kV and 400 kV) 
      8) Finland 1993 (220 kV and 400 kV) 
      9) Fennoscandia 1984 (400 kV)       
      10) basic situation for Ph.D.        
      11) your own choice                  
 

Select a menu number: 1
 
 
----- Choose the model of the electric field -----
 
      1) constant field (1 V/km northwards or eastwards)
      2) other possibilities
 

Select a menu number: 2

 
----- Select the model -----
 
      1) constant field northwards (1 V/km)
      2) constant field eastwards (1 V/km)
      3) image line current
      4) image sheet current
      5) random field
      6) E=0 in Finland, constant elsewhere
      7) Laatokka_line
      8) fictitious E-W anomaly
      9) "realistic" model of Finland
      10) half-circle anomaly
      11) "potential" line current field
 

Select a menu number: 4
 
----- Give the size of one grid rectangular (kmxkm) -----
 
      1) DX=10, DY=10
      2) DX=20, DY=20
      3) DX=30, DY=30
      4) DX=40, DY=40
      5) DX=50, DY=50
      6) Your own choice
 

Select a menu number: 3
 

heigth of the sheet current (km) = 110

width of the sheet (km) = 400

depth of the perfect conductor (km) = 500

x-coordinate of the centre of the current (km) at y=0 = 700

direction of the current positive from E to S (degrees) = 0
 
 Calculating the electric field...
 
  = -0 V/km;  = 1 V/km
 
  = 1 V/km
 
 
----- Interpolation -----
 
      1) first order (4 points)
      2) third order (16 points)
 

Select a menu number: 1

 Calculating voltages...

 Calculating GICs...

 Elapsed time of the last step = 1.7 s.
 
 
Earthing currents (A); E given by yourself
 
   1 Inkoo                1.6
   2 Loviisa             32.5
   3 NurmijŠrvi          -0.7
   4 Lieto              -13.4
   5 HyvinkŠŠ            -0.3
   6 Koria               17.6
   7 Olkiluoto          -68.1
   8 Ulvila             -16.4
   9 Kangasala          -19.1
  10 Huutokoski          67.2
  11 AlajŠrvi           -18.7
  12 AlapitkŠ            30.1
  13 Pikkarala           69.6
  14 PetŠjŠskoski        52.0
  15 Pirttikoski         97.2
  16 Letsi             -151.4
  17 Messaure           -79.8
 
Mean (absolute) value      43
 
 
 
Transmission line currents (A); E as above
 
   1   4  INK-LIE    -46.4
   1   5  INK-HYV     44.8
   2   3  LOV-NUR    -55.8
   2   6  LOV-KOR     23.4
   3   5  NUR-HYV    -55.2
   4   7  LIE-OLK    -33.1
   5   9  HYV-KAN    -10.1
   6  10  KOR-HUU      5.8
   7   8  OLK-ULV     35.0
   8  11  ULV-ALJ     51.4
   9  11  KAN-ALJ      9.0
  10  11  HUU-ALJ    -53.7
  10  12  HUU-ALP     -7.7
  11  13  ALJ-PIK     25.5
  12  14  ALP-PET    -37.9
  13  15  PIK-PIR     35.7
  13  17  PIK-MES    -79.8
  14  15  PET-PIR     61.5
  14  16  PET-LET   -151.4
 
Mean (absolute) value   43
 
 

 |GIC| in ascending order at stations
 Letsi              151.4
 Pirttikoski         97.2
 Messaure            79.8
 Pikkarala           69.6
 Olkiluoto           68.1
 Huutokoski          67.2
 PetŠjŠskoski        52.0
 Loviisa             32.5
 AlapitkŠ            30.1
 Kangasala           19.1
 AlajŠrvi            18.7
 Koria               17.6
 Ulvila              16.4
 Lieto               13.4
 Inkoo                1.6
 NurmijŠrvi           0.7
 HyvinkŠŠ             0.3

 Results are in the file GIC_result.
Č