calculate cutoff rigidity from Virtual Dipole Moment time series

calculate cutoff rigidity from Virtual Dipole Moment time series

vdm2rc(vdm, lat, model = "elsasser54")

Arguments

vdm

Virtual Dipole Moment (A.m^2)

lat

Latitude (deg)

model

Model used to compute Rc from the virtual dipole moment value (one of "elsasser54" or "lifton14")

Examples

vdm2rc
#> function (vdm, lat, model = "elsasser54") 
#> {
#>     M0 = 7.746e+22
#>     if (model == "elsasser54") {
#>         Rc = 14.31187 * vdm/M0 * (cos(d2r(lat)))^4
#>     }
#>     else if (model == "lifton14") {
#>         dd = c(6.89901, -103.241, 522.061, -1152.15, 1189.18, 
#>             -448.004)
#>         Rc = vdm/M0 * (dd[1] * cos(d2r(lat)) + dd[2] * (cos(d2r(lat)))^2 + 
#>             dd[3] * (cos(d2r(lat)))^3 + dd[4] * (cos(d2r(lat)))^4 + 
#>             dd[5] * (cos(d2r(lat)))^5 + dd[6] * (cos(d2r(lat)))^6)
#>     }
#>     else {
#>         stop("model must be one of elsasser54 or lifton14")
#>     }
#>     return(Rc)
#> }
#> <bytecode: 0x56376e9f3c10>
#> <environment: namespace:TCNtools>