## data sets.
## Colatitude and longitude of geomagnetic observatories.
## Not run:
# coord=matrix(c(124.43, 19.23, 53.77,140.18,68.68,202.00,71.89,293.85),nrow=2,ncol=4)
#
#
# ## Estimation of SI index
# index.sample<- SAIndex(record, coord, wf="la8")
# si.v<-index.sample$SI
#
#
# ## generate datetime for one week period
# start.date="2001-3-1"
# end.date="2001-4-30"
#
# ## estimation of sq
# Sq<- SQ (datasq, si.v=si.v, wf = "la8")
#
# sqplot(Sq, Title="Sq variation", start = start.date, end=end.date, n.station=4,
# graphs.per.page=2, station.names=c("ABG","PHU","TUC","FRD"))
# ## End(Not run)
Run the code above in your browser using DataLab