mola2gmt - A program that runs various GMT programs to convert raw MOLA data into a format usable by GMT.


Copyright (C) 2005,2006 Ross A. Beyer,

This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version.

This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.

You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA or visit their website at


mola2gmt [--label_file <lbl file>] [--output_file <filename>] [--title <title string>] [--remark <remark string>] [--region <w/e/s/n> [--delete]] [--swap] [--no_execute] [--version] [--help] [--man] <MOLA gridded data file>


This program takes in a MOLA gridded data file and converts it into a netCDF file that GMT can work with. It is mostly a wrapper around the GMT program xyz2grd.


--label_file <lbl file>
This program needs to read the PDS label file for the gridded data file, to get information for the conversion. By default, the program will look for a file of the same name in the same directory as the input gridded data file but with an ``.lbl'' extension instead of an ``.img'' extension. If your data isn't set up like that, you can use this argument to give the location of the label file.

--output_file <filename>
You can specify the output netCDF filename with this argument. If not specified, the output netCDF file will be named similarly to the input file, but with the extension ``.grd''.

--title <title string>
You can specify the title that will be placed into the netCDF file. If not specified, it defaults to a sensible value read from the label file.

--remark <remark string>
You can specify the remark that will be placed into the netCDF file. If not specified, it defaults to ``n pixels per degree'' where the `n' is replaced with the appropriate resolution from the label file.

--region <w/e/s/n>
If this option is specified, in addition to converting the gridded data file into a netCDF file, GMT grdcut will also be run to make the size of the output gridfile restricted to the region of interest specified. The values in the region of interest will also be used to automatically create a filename for the results of GMT grdcut.

The underlying xyz2grd program needs to know how to read the binary data of the raw MOLA gridded data file. The MOLA label file indicates that the data (SAMPLE_TYPE) are MSB_INTEGER. If you are running this program on a big-endian computer (like a Sun system or a pre-Intel Mac for example) then you can safely ignore this option. However, if you are running this program on a small-endian computer (like most commodity PC hardware that runs Linux, FreeBSD, etc.), then you will need to use this option. If you're not sure, ask your system administrator or just grab a gridded data file of an area that you know, and try it one way. If the topography looks weird, try it the other way.

Print the commands that would be executed, but do not execute them.

If the --region argument is given, this will delete the output of the xyz2grd program after it has been cut down, leaving only the output of GMT grdcut.

Display the version number and information, and then exit.

Display a short usage statement and exit.

Display manual page and exit.


Ross A. Beyer,