Summary

This collection of files is a fortran wrapper of the magnetosphere field models of N. A. Tsyganenko. See his website for the original fortran codes and documentation. This project extends Tsyganenko's original codes with parallelized helper functions. These use the magnetosphere models to efficiently generate full-field data stored in standard netcdf format. Included with the original magnetosphere modelling functionallity, is a small program that can calculate the local reconnection rate as defined in this paper.

This directory contains:

• README.md
  • This file
• src/
  • main.f90
    • The main program that feeds the models data and manages in/out and program control.
  • inout.f90
    • A simple module that implements input/output through the net-CDF library.
  • TA16_RBF.f
    • A modification of the original TA16_RBF.f model code. The compute block is wrapped in a module for interfacing with the main program. Redundant code (generating the RBF centers grid and loading parameters) have been placed in their own subroutines.
  • geopack.f
    • As above, this is the original geopack.f file. The only modification was to wrap the included subroutines inside a Module interface.
  • recon.f90
    • A short routine that reads in the netcdf field files and calculates the local reconnection rate.
  • input_preparation/
    • Parameters.par
      • The parameter data for generating TS04c input data.
    • calculate_tilt.f
      • A helper script for calculating the geo-tilt from the geopack module.
    • format_data.py
      • Calculates the TA16 input parameters and outputs formatted input data for the TA16 model if given an OMNI dataset.
    • omni_data.fmt
      • The required OMNI data format for the format_data.py routine.

How to compile/install

This code requires the fortran NET-CDF library. Below are compilation commands for the Intel and GNU compilers.

Intel

ifort inout.f90 geopack.f TA16_RBF.f main.f90 -o main.out `nf-config --fflags --flibs`

Gnu Compilers

gfortran inout.f90 geopack.f TA16_RBF.f main.f90 -o main.out `nf-config --fflags --flibs`

Running the models require three files,

• input_parameters.lst -- The grid info for the models.
• input_data.lst -- The simulation input data prepared for the specific model being used (TS05 or TA16).
  • TA16_RBF.par -- The ta16 fitting parameters if the TA16 model is being used.

For an example of how to use this script, see the example directory, example/

To compile the reconnection code, the following work:

Intel

ifort inout.f90 recon.f90 -o recon.out `nf-config --fflags --flibs`

Gnu Compilers

gfortran inout.f90 recon.f90 -o recon.out `nf-config --fflags --flibs`