Modeling planetary magnetospheres