Be exposure age calculator MATLAB function reference Version 2

More documents

Recommendations

Info

2.10 r2d.m Syntax: degrees = r2d(radians) This function converts angular measurements in radians to angular measurements in degrees. 40
2.11 skyline.m Syntax: [out,data,ver] = skyline(az,el,strike,dip) This function calculates the topographic shielding correction to the cosmic-ray flux at a site whose horizon is obstructed either by a dipping surface or by the surrounding topography. The input arguments strike and dip are the strike and dip of the sampled surface in degrees. The input arguments az and el are vectors containing the azimuths (degrees; 0 < az < 360; north = 0) and horizon angles (degrees; 0 < el < 360; zenith = 90) of points on the horizon. out is the shielding correction for the sample (nondimensional; the ratio of the nuclide production rate at the shielded site to the nuclide production rate at an unshielded site at the same location). data is a vector containing the interpolated horizon angle visible to the sample, in degrees, for 1-degree increments of azimuth. It can be used to plot the horizon, e.g., plot(0:360,data) ver is a string variable containing the version number of the function. The calculation goes as folllows: First, we divide the range of azimuths 0 < φ < 2π into 1-degree increments. Second, we calculate the horizon angle θ of the dipping surface as a function of azimuth: θ = arctan [tan θ d cos (φ − φ s )] (74) where θ d is the dip and φ d is the strike of the surface. Third, we interpolate the supplied horizon points to give the horizon angle of the topography as a function of azimuth. Fourth, at each increment of azimuth, we take the higher of the horizon angle due to the dipping surface or the horizon angle due to the surrounding topography to be the horizon angle visible to the sample. This results in a vector θ(φ) which defines the visible horizon angle θ at each increment of azimuth φ. Finally, the fraction of the cosmic-ray flux f(φ) that lies below the horizon angle in the increment of azimuth at φ is: f(φ) = δφ 2π (sin θ)3.3 (75) where δφ is the width of the azimuth increment in radians (π/180). The total fraction of the cosmic-ray flux that lies below the horizon f is the sum of the values of f(φ) for each increment of φ. The topographic shielding correction is then (1 − f). 41
Page 1 and 2: CRONUS-Earth 26 Al- 10 Be exposure
Page 3 and 4: 1 Wrapper scripts and control funct
Page 5 and 6: 1.2 al be erosion many v2.m Syntax:
Page 7 and 8: 1.3 get al be age.m results = get_a
Page 9 and 10: If sample.pressure is not set, calc
Page 11 and 12: the geographic latitude of the samp
Page 13 and 14: The external uncertainty in the exp
Page 15 and 16: Non-scaling-scheme-specific informa
Page 17 and 18: In like manner, we have already cal
Page 19 and 20: al be consts.Fsp10 fraction of refe
Page 21 and 22: 1.6 makeCplot.m Syntax: filename =
Page 23 and 24: 2 Subsidiary calculation functions
Page 25 and 26: 2.3 d2r.m Syntax: radians = d2r(deg
Page 27 and 28: Obtains the altitude scaling factor
Page 29 and 30: Finally, calculates the scaling fac
Page 31 and 32: where z is depth (g · cm −2 ),
Page 33 and 34: 2.8 NCEPatm 2.m Syntax: pressure =
Page 35 and 36: out.phi vert slhl Flux of verticall
Page 37 and 38: R v (Z) = R v,0 (Z)e δz Λµ(Z) (5
Page 39: f 10,C = 0.704 f 26,C = 0.296 (65)
Page 43 and 44: 2.13 stone2000.m Syntax: scalingfac
Page 45 and 46: 2.14 stone2000Rcsp.m Syntax: scalin
Page 47 and 48: 2.15 thickness.m Syntax: Sthick = t

Be exposure age calculator MATLAB function reference Version 2

You also want an ePaper? Increase the reach of your titles

Delete template?

Save as template?