Please note: Save your scripts to a directory named Homework-5 in your home account on the server computer in SCA 219 (131.247.211.166); do not email your answers. A short report should be typed and submitted in PDF format. A flowchart is not necessary for this assignment. Multiple pages should be stapled together. Some additional explanation is helpful and often necessary for the person reading your assignment to understand exactly what you are doing and why. Assume that the person grading your assignment does not have a copy of the assignment. What is important here is that (1) you understand what you are doing, and (2) that you can communicate what you are doing to others.

Introduction

In this homework assignment we begin to explore and create maps using Generic Mapping Tools (GMT). Be sure to read the handout about maps (maps.pdf) and the introduction to GMT (gmt.pdf) before trying to start this assignment. Work through the examples provided in these handouts.

In assignment 4, you calculated the volume of water in Crater Lake using a script and a data file, filtered crater lake data.xyz. Now look more carefully at the contents of this file. Each line of the file contains two geographic coordinates and an elevation. For example, the first line of the file shows:

568000 4758500 2041

The first two columns contain UTM coordinates of the elevation estimate. The last column is the elevation. Crater Lake is in UTM zone 10. The USGS saved the data using datum NAD27.

Problem 1

Use the Proj.4 Cartographic Library (see the maps.pdf handout) to convert this file from UTM coordinates to:

longitude latitude elevation

using the map datum WGS84. If you specify the name of your data file (e.g., filtered crater lake data.xyz) on the command line, after the parameter list, you can then redirect the conversions to a file using > filename, the same syntax you used when executing your Perl or Python script. Note that when converting to longitudes and latitudes it is necessary to use at least 6 decimal places, for accuracy.

For this problem, be sure to include in your report 1) the command line you use to make this file conversion and 2) the first 2 lines of your converted data file.

Problem 2

Now that you know the approximate latitude and longitude of Crater Lake:

Make a location map of the volcano using the GMT program pscoast. Be sure to include sufficient area so it is easy to visualize the location of the volcano with respect to geographic features of the western US (e.g., the coastline) and neaby US state boundaries.

Plot the location of Crater Lake on your location map and label it. Use the coordinate from Problem 1 and the GMT program psxy to plot the location of Crater Lake as a large red triangle. Label your red triangle using the GMT command pstext. Also label 2 US states, and any large, nearby bodies of water. This means you will have at least 4 pstext commands.

1

Explain the meaning of each GMT command and parameters used with each of your commands. Use the GMT manual pages, technical reference and cookbook, or the tutorial to find this information.

Report your results. The report should include 1) your problem statement(s), 2) a labeled location map, 3) a figure caption for your map which explains what and where and any symbols you used, 4)the GMT command lines that you used to create the map, and a description of what each GMT command comtributes to your map and how each parameter contributes to the command. Commands used multiple times only need to be described once.

2