Run on geomag in /home/purucker/MF3 directory

./crustal_thickness_iterate_tf.x

 MF3 RMS values are:   2.33946157  2.47065043 nT

Initial update from 3SMAC to MF3
 Rem reduced RMS values are:   2.63242435  2.89068604 nT
 Calculating z and total field from crustal thickness
 3000 out of 11562 complete.
 6000 out of 11562 complete.
 9000 out of 11562 complete.
 Converting field to spherical harmonics
 Converting field to spherical harmonics
 Removing degrees 1-15 of field
 Removing degrees 1-15 of field
 Spherical Harmonic Filtering complete
 Obs-model RMS values are:   2.72527385  3.14401102 nT
 Observations read in
 Information matrix accumulation begins
 3000 of 11562 A and b
 6000 of 11562 A and b
 9000 of 11562 A and b

   1637825 Length of ijc and sc
 Conjugate gradient inversion beginning
 Observations read in
Initial SD on fit =    0.273E+01
 Conjugate gradient complete after  26 iterations
  with an error of   0.0293997973

 LATITUDE   LONGITUDE     Del T (km)
      40.59    280.29         6.1216
      17.65     15.88        -8.3050
      -1.76    302.82        -9.6696
     -22.94    230.82         1.9410
     -45.88    220.32        -7.0357

Iteration 2
 Rem reduced RMS values are:   2.63242435  2.89068604 nT
 Calculating z and total field from crustal thickness
 3000 out of 11562 complete.
 6000 out of 11562 complete.
 9000 out of 11562 complete.
 Converting field to spherical harmonics
 Converting field to spherical harmonics
 Removing degrees 1-15 of field
 Removing degrees 1-15 of field
 Spherical Harmonic Filtering complete
 Obs-model RMS values are:   0.63499701  1.49830317 nT
 Observations read in
 Information matrix accumulation begins
 3000 of 11562 A and b
 6000 of 11562 A and b
 9000 of 11562 A and b

   1637825 Length of ijc and sc
 Conjugate gradient inversion beginning
 Observations read in
Initial SD on fit =    0.635E+00
 Conjugate gradient complete after  26 iterations
  with an error of   0.0234158039

 LATITUDE   LONGITUDE     Del T (km)
      40.59    280.29        -1.1197
      17.65     15.88        -3.7027
      -1.76    302.82         9.2609
     -22.94    230.82        -0.3795
     -45.88    220.32        -1.4868


Iteration 3
 Rem reduced RMS values are:   2.63242435  2.89068604 nT
 Calculating z and total field from crustal thickness
 3000 out of 11562 complete.
 6000 out of 11562 complete.
 9000 out of 11562 complete.
 Converting field to spherical harmonics
 Converting field to spherical harmonics
 Removing degrees 1-15 of field
 Removing degrees 1-15 of field
 Spherical Harmonic Filtering complete
 Obs-model RMS values are:   0.36540395  1.39289045 nT
 Observations read in
 Information matrix accumulation begins
 3000 of 11562 A and b
 6000 of 11562 A and b
 9000 of 11562 A and b

 A description of the production and contents of 
    working_model_MF3.xyz
----
   1637825 Length of ijc and sc
 Conjugate gradient inversion beginning
 Observations read in
Initial SD on fit =    0.365E+00
 Conjugate gradient complete after  26 iterations
  with an error of   0.0206605848


 LATITUDE   LONGITUDE     Del T (km)
      40.59    280.29        -0.9298
      17.65     15.88        -1.7781
      -1.76    302.82         0.0135
     -22.94    230.82        -0.5086
     -45.88    220.32        -0.2244

---------------
The individual elements of working_model_MF3.xyz are

Column 1: Geocentric latitude
Column 2: Geocentric longitude
Column 3: Direction cosine 1 (ss)
Column 4: Direction cosine 2 (cc)
Column 5: Direction cosine 3 (cs)
Column 6: Main field strength (nT)
Column 7: Magnetic crustal thickness (km)
Column 8: MF-3 total field at 400 km (nT)
Column 9: MF-3 Z field at 400 km (nT)

This file is working_model_MF3.readme in /home/purucker/MF3 on geomag.gsfc.nasa.gov
13 February 2004