Scribd
Upload
6 views1 page

Grid to Magnetic Field Adjustment

This document defines parameters for a rectangular grid of stations to compute gravitational and magnetic fields. It specifies variables like the number of corners, their coordinates, density, and includes code to add a water depth offset. It also defines the faces of blocks in the grid and starting values for the station grid like the south end coordinate, station spacing, north end, west end, profile spacing and east end.

Uploaded by

Phyphi Lovito
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
6 views1 page

Grid to Magnetic Field Adjustment

This document defines parameters for a rectangular grid of stations to compute gravitational and magnetic fields. It specifies variables like the number of corners, their coordinates, density, and includes code to add a water depth offset. It also defines the faces of blocks in the grid and starting values for the station grid like the south end coordinate, station spacing, north end, west end, profile spacing and east end.

Uploaded by

Phyphi Lovito
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd

PROGRAM

calgrv=1; % change to zero if gravity field is not required


calmag=1; %change to zero if magnetic field is not required
Ncor=26;
Hintn=50000;
Hincl=50;
Decl=0;
Susc=0.01
Mstrength=0;
Mincl=0;
Mdecl=0;
Corner=[100 0 0;75 -75 0;0 -100 0;-75 -75 0;100 0 0;...
-75 75 0;0 100 0;75 75 0;75 0 -75;60 -60 -60;0 -75 -75;...
-60 -60 -60;-75 0 -75;-60 60 -60;0 75 -75;60 60 -60;...
0 0 -100;75 0 75;60 -60 60;0 -75 75;-60 -60 60;-75 0 75;...
-60 60 60;0 75 75;60 60 60;0 0 100];
fht=200;dens=10;
%Add fht to depths of all corners
Corner(:,3)=Corner(:,3)+fht;
Face=zeros([50,9]); %initialize a sufficiently large array
Face(1,1:5)=[4 1 2 10 9];Face(2,1:5)=[4 2 3 11 10];
Face(3,1:5)=[4 2 4 12 11];Face(4,1:5)=[4 4 5 13 12];
Face(5,1:5)=[4 5 6 14 13];Face(6,1:5)=[4 6 7 15 14];
Face(7,1:5)=[4 7 8 16 15];Face(8,1:5)=[4 8 1 9 16];
Face(9,1:5)=[4 9 10 11 17];Face(10,1:5)=[4 11 12 13 17];
Face(11,1:5)=[4 13 14 15 17];Face(12,1:5)=[4 15 16 9 17];
Face(13,1:5)=[4 1 18 19 2];Face(14,1:5)=[4 2 19 20 3];
Face(15,1:5)=[4 3 20 21 4];Face(16,1:5)=[4 4 21 22 5];
Face(17,1:5)=[4 5 22 23 6];Face(18,1:5)=[4 6 23 24 7];
Face(19,1:5)=[4 7 24 25 8];Face(20,1:5)=[4 8 25 18 1];
Face(21,1:5)=[4 20 19 18 26];Face(22,1:5)=[4 22 21 20 26];
Face(23,1:5)=[4 24 23 22 26];Face(24,1:5)=[4 18 25 24 26];
%Comments:Rectangular grid of stations for computing fields. profiles are
%not long the x-axis (north-south direction). All values are in meters.
s_end=-320; %stating value of x;south end of profiles
stn_spcng=40; %Stepsize in north direction;stn interval
n_end=320; %Last x;maximum north coordinate
w_end=0; %y value for westernmost profile
prof_spcng=1; %profile spacing
e_end=0; %y value for easternmost profile

You might also like