/[escript]/trunk/doc/examples/inversion/mag_netcdf.py
ViewVC logotype

Diff of /trunk/doc/examples/inversion/mag_netcdf.py

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 4286 by caltinay, Fri Mar 1 04:50:54 2013 UTC revision 4287 by caltinay, Thu Mar 7 05:26:44 2013 UTC
# Line 31  from esys.escript import saveDataCSV Line 31  from esys.escript import saveDataCSV
31  # Set parameters  # Set parameters
32  DATASET = 'data/MagneticSmall.nc'  DATASET = 'data/MagneticSmall.nc'
33  DATA_UNITS = U.Nano * U.V * U.sec / (U.m**2)  DATA_UNITS = U.Nano * U.V * U.sec / (U.m**2)
 latitude = -20.5  
34  PAD_X = 0.2  PAD_X = 0.2
35  PAD_Y = 0.2  PAD_Y = 0.2
36  thickness = 40. * U.km  thickness = 40. * U.km
37  l_air = 6. * U.km  l_air = 6. * U.km
38  n_cells_v = 25  n_cells_v = 25
39  MU = 0.1  MU = 0.1
40    # background magnetic field components (B_North, B_East, B_Vertical)
41    B_b = [31232.*U.Nano*U.Tesla, 2201.*U.Nano*U.Tesla, -41405.*U.Nano*U.Tesla]
42    
43  # Setup and run the inversion  # Setup and run the inversion
 B_b=simpleGeoMagneticFluxDensity(latitude=latitude)  
44  source=NetCdfData(NetCdfData.MAGNETIC, DATASET, scale_factor=DATA_UNITS)  source=NetCdfData(NetCdfData.MAGNETIC, DATASET, scale_factor=DATA_UNITS)
45  db=DomainBuilder(dim=3)  db=DomainBuilder(dim=3)
46  db.addSource(source)  db.addSource(source)

Legend:
Removed from v.4286  
changed lines
  Added in v.4287

  ViewVC Help
Powered by ViewVC 1.1.26