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

Diff of /trunk/doc/examples/inversion/gravmag_nodriver.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 36  MAGNETIC_DATASET = 'data/MagneticSmall.n Line 36  MAGNETIC_DATASET = 'data/MagneticSmall.n
36  MAG_UNITS = U.Nano * U.V * U.sec / (U.m**2)  MAG_UNITS = U.Nano * U.V * U.sec / (U.m**2)
37  GRAVITY_DATASET = 'data/GravitySmall.nc'  GRAVITY_DATASET = 'data/GravitySmall.nc'
38  GRAV_UNITS = 1e-6 * U.m/(U.sec**2)  GRAV_UNITS = 1e-6 * U.m/(U.sec**2)
39  latitude=-20.5  # background magnetic field components (B_North, B_East, B_Vertical)
40    B_b = [31232.*U.Nano*U.Tesla, 2201.*U.Nano*U.Tesla, -41405.*U.Nano*U.Tesla]
41    
42  thickness = 40. * U.km # below surface  thickness = 40. * U.km # below surface
43  l_air = 6. * U.km      # above surface  l_air = 6. * U.km      # above surface
# Line 52  MU_MAGNETIC = 0.1 Line 53  MU_MAGNETIC = 0.1
53  # read data:  # read data:
54  source_g=NetCdfData(NetCdfData.GRAVITY, GRAVITY_DATASET, scale_factor=GRAV_UNITS)  source_g=NetCdfData(NetCdfData.GRAVITY, GRAVITY_DATASET, scale_factor=GRAV_UNITS)
55  source_m=NetCdfData(NetCdfData.MAGNETIC, MAGNETIC_DATASET, scale_factor=MAG_UNITS)  source_m=NetCdfData(NetCdfData.MAGNETIC, MAGNETIC_DATASET, scale_factor=MAG_UNITS)
 B_b=simpleGeoMagneticFluxDensity(latitude=latitude)  
56    
57  # create domain:  # create domain:
58  db=DomainBuilder(dim=3)  db=DomainBuilder(dim=3)

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

  ViewVC Help
Powered by ViewVC 1.1.26