netcdf TGCM.data.magdat {
dimensions:
	zjmx = 36 ;
	zimxp1 = 73 ;
	zjmxp2 = 38 ;
	imaxmp = 81 ;
	jmaxm = 97 ;
	dim2 = 2 ;
	dim3 = 3 ;
	dim4 = 4 ;
variables:
	double ALATM(zjmxp2, zimxp1) ;
		ALATM:long_name = "geomagnetic latitude at each geographic grid point" ;
		ALATM:units = "radians" ;
	double ALONM(zjmxp2, zimxp1) ;
		ALONM:long_name = "geomagnetic longitude at each geographic grid point" ;
		ALONM:units = "radians" ;
	double XB(zjmxp2, zimxp1) ;
		XB:long_name = "northward component of magnetic field" ;
	double YB(zjmxp2, zimxp1) ;
		YB:long_name = "eastward component of magnetic field" ;
	double ZB(zjmxp2, zimxp1) ;
		ZB:long_name = "downward component of magnetic field" ;
	double BMOD(zjmxp2, zimxp1) ;
		BMOD:long_name = "magnitude of magnetic field" ;
		BMOD:units = "gauss" ;
	double DMLAT(zjmxp2, zimxp1) ;
		DMLAT:long_name = "dipole latitude corresponding to apex of field line" ;
	double RJAC(dim2, dim2, zjmxp2, zimxp1) ;
		RJAC:long_name = "scaled derivatives of geomagnetic coordinates w.r.t. geographic coordinates" ;
	double AV(dim2, dim3, zjmxp2, zimxp1) ;
		AV:long_name = "magnetic vectors A1 and A2" ;
	double P(zjmxp2, zimxp1) ;
		P:long_name = "scale factor between geographic and geomagnetic systems divided by sin(dip)*cos(thetas)" ;
	double RMAG11(zjmxp2, zimxp1) ;
		RMAG11:long_name = "(A1.A1)/P*SIN(I)*COS(THETAS)" ;
	double RMAGC(zjmxp2, zimxp1) ;
		RMAGC:long_name = "(A1.A2)/P*SIN(I)" ;
	double RMAG2(zjmxp2, zimxp1) ;
		RMAG2:long_name = "1./BMOD" ;
	double RMAG22(zjmxp2, zimxp1) ;
		RMAG22:long_name = "(A2.A2)/P*SIN(I)/COS(THETAS)" ;
	double RJACD(zjmxp2, zimxp1) ;
		RJACD:long_name = "determinant of RJAC" ;
	int IG(jmaxm, imaxmp) ;
		IG:long_name = "geographic latitude indices containing each geomagnetic latitude grid point" ;
	int JG(jmaxm, imaxmp) ;
		JG:long_name = "geographic longitude indices containing each geomagnetic longitude grid point" ;
	int IM(zjmxp2, zimxp1) ;
		IM:long_name = "geomagnetic latitude indices containing each geographic latitude grid point" ;
	int JM(zjmxp2, zimxp1) ;
		JM:long_name = "geomagnetic longitude indices containing each geographic longitude grid point" ;
	double WT(jmaxm, imaxmp, dim4) ;
		WT:long_name = "interpolation weights" ;
	double DIM(zjmxp2, zimxp1) ;
		DIM:long_name = "fractions in I direction for interpolation" ;
	double DJM(zjmxp2, zimxp1) ;
		DJM:long_name = "fractions in J direction for interpolation" ;
	double CSLATM(zjmx, zimxp1) ;
		CSLATM:long_name = "COS(THETAS)" ;
	double SNLATM(zjmx, zimxp1) ;
		SNLATM:long_name = "SIN(THETAS)" ;
	double CSLONM(zjmx, zimxp1) ;
		CSLONM:long_name = "COS(LAMDAS)" ;
	double SNLONM(zjmx, zimxp1) ;
		SNLONM:long_name = "SIN(LAMDAS)" ;
	double CSLATG(zjmx) ;
		CSLATG:long_name = "COS(THETA)" ;
	double SNLATG(zjmx) ;
		SNLATG:long_name = "SIN(THETA)" ;
	double CSLONG(zimxp1) ;
		CSLONG:long_name = "COS(LAMDA)" ;
	double SNLONG(zimxp1) ;
		SNLONG:long_name = "SIN(LAMDA)" ;

// global attributes:
		:title = "Magnetic field data file for use by tgcm" ;
		:msspath = "/TGCM/data/magdat.nc" ;
		:zjmx = "dimension: number of geographic latitudes" ;
		:zimxp = "dimension: number of geographic longitudes plus 1 periodic point" ;
		:zjmxp2 = "dimension: number of geographic latitudes plus 2 (model vars are dimensioned 0:zjmxp1)" ;
		:imaxmp = "dimension: number of apex longitudes plus 1 periodic point" ;
		:jmaxm = "dimension: number of apex latitudes" ;
data:
}