eval '(exit $?0)' && eval 'exec perl -S $0 ${1+"$@"}' & eval 'exec perl -S $0 $argv:q' if 0; #-------------------------------------------------------------------- # The MB-system: mbm_xyplot.perl 8/6/95 # $Id: mbm_xyplot.pl 1867 2010-09-15 05:30:50Z caress $ # # Copyright (c) 1993-2009 by # D. W. Caress (caress@mbari.org) # Monterey Bay Aquarium Research Institute # Moss Landing, CA # and D. N. Chayes (dale@lamont.ldgo.columbia.edu) # Lamont-Doherty Earth Observatory # Palisades, NY 10964 # # See README file for copying and redistribution conditions. #-------------------------------------------------------------------- # # Command: # mbm_xyplot # # Purpose: # Macro to generate a shellscript of GMT commands which, # when executed, will generate a Postscript plot of xy data. # Axes may be linear, log, or any of several geographic # projections. Data may be plotted as symbols or lines. # The plot will be scaled to fit on the specified page size # or, if the scale is user defined, the page size will be # chosen in accordance with the plot size. The primary purpose # of this macro is to allow the simple, semi-automated # production of nice looking plots with a few command line # arguments. For users seeking more control over the plot # appearance, a number of additional optional arguments are # provided. Truly ambitious users may edit the plot shellscript # to take advantage of GMT capabilites not supported by this # macro. # # Basic Usage: # mbm_xyplot -I[filepars:]file [ -Gfill -H -Oroot # -Ppagesize -Ssymbol/size -Uorientation # -V -Wpen ] # # Additional Options: # [-Btickinfo -Jprojection[/scale | width] # -Ltitle[:scale_label] -Mmisc -Q -Rw/e/s/n -X] # # Miscellaneous Options: # [-MGDgmtdef/value -MGL[f][x]lon0/lat0/slat/length[m] # -MGTx/y/size/angle/font/just/text # -MGU[/dx/dy/][label] ] # # Author: # David W. Caress # Lamont-Doherty Earth Observatory # Palisades, NY 10964 # August 9, 1995 # # Version: # $Id: mbm_xyplot.pl 1867 2010-09-15 05:30:50Z caress $ # # Revisions: # $Log: mbm_xyplot.perl,v $ # Revision 5.15 2008/09/11 20:06:46 caress # Checking in updates made during cruise AT15-36. # # Revision 5.14 2007/10/08 04:29:06 caress # Added some large page definitions. # # Revision 5.13 2006/07/05 19:50:21 caress # Working towards 5.1.0beta # # Revision 5.12 2006/06/16 19:30:58 caress # Check in after the Santa Monica Basin Mapping AUV Expedition. # # Revision 5.11 2006/01/06 18:26:26 caress # Working towards 5.0.8 # # Revision 5.10 2005/11/05 01:34:20 caress # Much work over the past two months. # # Revision 5.9 2005/03/25 04:05:40 caress # Fixed handling of tickinfo string. # For mbm_plot only, added control on filename annotation direction. # # Revision 5.8 2004/09/16 19:12:55 caress # Supports postscript viewer ggv. # # Revision 5.7 2003/11/25 21:25:43 caress # Implemented Val Schmidt's changes. # # Revision 5.6 2003/04/17 20:42:48 caress # Release 5.0.beta30 # # Revision 5.5 2002/04/06 02:51:54 caress # Release 5.0.beta16 # # Revision 5.4 2001/12/18 04:26:12 caress # Version 5.0.beta11. # # Revision 5.3 2001/11/02 21:07:40 caress # Adjusted handling of segmented xy files. # # Revision 5.2 2001/10/10 23:56:01 dcaress # Regrettably, I don't remember what I changed... # # Revision 5.1 2001-06-30 10:36:53-07 caress # Release 5.0.beta02 # # Revision 5.0 2000/12/01 22:58:01 caress # First cut at Version 5.0. # # Revision 4.10 2000/10/03 21:42:17 caress # Snapshot for Dale. # # Revision 4.9 1999/06/25 17:55:47 caress # I must have changed something! # # Revision 4.8 1999/05/06 23:46:32 caress # Release 4.6a # # Revision 4.7 1999/04/16 01:25:51 caress # Version 4.6 final release? # # Revision 4.6 1999/04/15 19:28:52 caress # Fixed sprintf statements. # # Revision 4.5 1999/02/04 23:39:54 caress # MB-System version 4.6beta7 # # Revision 4.4 1998/10/05 17:00:15 caress # MB-System version 4.6beta # # Revision 4.3 1997/04/21 16:54:41 caress # MB-System 4.5 Beta Release. # # Revision 4.3 1997/04/17 15:06:49 caress # MB-System 4.5 Beta Release # # Revision 4.2 1996/03/12 17:28:19 caress # Check-in after flail with format 63. # # Revision 4.1 1995/09/28 18:05:43 caress # Various bug fixes working toward release 4.3. # # Revision 4.0 1995/08/17 14:53:25 caress # Revision for release 4.3. # # # # $program_name = "mbm_xyplot"; # set page size database @page_size_names = ( "a", "b", "c", "d", "e", "f", "e1", "a0", "a1", "a2", "a3", "a4", "a5", "a6", "a7", "a8", "a9", "a10", "b0", "b1", "b2", "b3", "b4", "b5", "b6", "b7", "b8", "b9", "b10", "c0", "c1", "c2", "c3", "c4", "c5", "c6", "c7", "m1", "m2", "m3", "m4", "m5", "m6"); %page_width_in = ( "a", 8.50, "b", 11.00, "c", 17.00, "d", 22.00, "e", 34.00, "f", 28.00, "e1", 44.00, "a0", 33.11, "a1", 23.39, "a2", 16.54, "a3", 11.69, "a4", 8.27, "a5", 5.83, "a6", 4.13, "a7", 2.91, "a8", 2.05, "a9", 1.46, "a10", 1.02, "b0", 39.37, "b1", 27.83, "b2", 19.68, "b3", 13.90, "b4", 9.84, "b5", 6.93, "b6", 4.92, "b7", 3.46, "b8", 2.44, "b9", 1.73, "b10", 1.22, "c0", 36.00, "c1", 25.60, "c2", 18.00, "c3", 12.80, "c4", 9.00, "c5", 6.40, "c6", 4.50, "c7", 3.20, "m1", 54.00, "m2", 54.00, "m3", 54.00, "m4", 60.00, "m5", 60.00, "m6", 60.00); %page_height_in = ( "a", 11.00, "b", 17.00, "c", 22.00, "d", 34.00, "e", 44.00, "f", 40.00, "e1", 68.00, "a0", 46.81, "a1", 33.11, "a2", 23.39, "a3", 16.54, "a4", 11.69, "a5", 8.27, "a6", 5.83, "a7", 4.13, "a8", 2.91, "a9", 2.05, "a10", 1.46, "b0", 56.67, "b1", 39.37, "b2", 27.83, "b3", 19.68, "b4", 13.90, "b5", 9.84, "b6", 6.93, "b7", 4.92, "b8", 3.46, "b9", 2.44, "b10", 1.73, "c0", 51.20, "c1", 36.00, "c2", 25.60, "c3", 18.00, "c4", 12.80, "c5", 9.00, "c6", 6.40, "c7", 4.50, "m1", 72.00, "m2", 84.00, "m3", 96.00, "m4", 72.00, "m5", 84.00, "m6", 96.00); %page_anot_font = ( "a", 8, "b", 12, "c", 16, "d", 24, "e", 24, "f", 24, "e1", 24, "a0", 24, "a1", 24, "a2", 16, "a3", 12, "a4", 8, "a5", 6, "a6", 6, "a7", 6, "a8", 4, "a9", 4, "a10", 4, "b0", 24, "b1", 24, "b2", 16, "b3", 16, "b4", 12, "b5", 8, "b6", 6, "b7", 4, "b8", 4, "b9", 4, "b10", 4, "c0", 24, "c1", 24, "c2", 16, "c3", 12, "c4", 8, "c5", 6, "c6", 6, "c7", 6, "m1", 24, "m2", 24, "m3", 24, "m4", 24, "m5", 24, "m6", 24); %page_header_font =( "a", 10, "b", 15, "c", 20, "d", 30, "e", 30, "f", 30, "e1", 30, "a0", 30, "a1", 30, "a2", 20, "a3", 15, "a4", 10, "a5", 8, "a6", 8, "a7", 8, "a8", 5, "a9", 5, "a10", 5, "b0", 30, "b1", 30, "b2", 20, "b3", 20, "b4", 15, "b5", 10, "b6", 8, "b7", 5, "b8", 5, "b9", 5, "b10", 5, "c0", 30, "c1", 30, "c2", 20, "c3", 15, "c4", 10, "c5", 8, "c6", 8, "c7", 8, "m1", 30, "m2", 30, "m3", 30, "m4", 30, "m5", 30, "m6", 30); %page_gmt_name = ( "a", "archA", "b", "archB", "c", "archC", "d", "archD", "e", "archE", "f", "B0", "e1", "B0", "a0", "A0", "a1", "A1", "a2", "A2", "a3", "A3", "a4", "A4", "a5", "A5", "a6", "A6", "a7", "A7", "a8", "A8", "a9", "A9", "a10", "A10", "b0", "B0", "b1", "B1", "b2", "B2", "b3", "B3", "b4", "B4", "b5", "B5", "b6", "A6", "b7", "A7", "b8", "A8", "b9", "A9", "b10", "A10", "c0", "B0", "c1", "B1", "c2", "B2", "c3", "B3", "c4", "B4", "c5", "B5", "c6", "B6", "c7", "B7", "m1", "Custom_4241x5655", "m2", "Custom_4241x6578", "m3", "Custom_4241x7540", "m4", "Custom_4712x5655", "m5", "Custom_4712x6578", "m6", "Custom_4712x7540"); %xpsview_mem = ( "a", "4m", "b", "6m", "c", "8m", "d", "12m", "e", "16m", "f", "16m", "e1", "16m", "a0", "16m", "a1", "12m", "a2", "8m", "a3", "6m", "a4", "4m", "a5", "4m", "a6", "4m", "a7", "4m", "a8", "4m", "a9", "4m", "a10", "4m", "b0", "16m", "b1", "12m", "b2", "8m", "b3", "8m", "b4", "6m", "b5", "4m", "b6", "4m", "b7", "4m", "b8", "4m", "b9", "4m", "b10", "4m", "c0", "16m", "c1", "12m", "c2", "8m", "c3", "6m", "c4", "4m", "c5", "4m", "c6", "4m", "c7", "4m", "m1", "16m", "m2", "16m", "m3", "16m", "m4", "16m", "m5", "16m", "m6", "16m"); # Determine the GMT version @grdinfo = `grdinfo 2>&1`; while (@grdinfo) { $line = shift @grdinfo; if ($line =~ /^grdinfo\s+(\S+)\s+\S+/) { ($gmt_version) = $line =~ /^grdinfo\s+(\S+)\s+\S+/; } } # Get the current working directory. $current_working_dir = `pwd`; # Deal with command line arguments $command_line = "@ARGV"; &MBGetopts('B:b:G:g:HhI+i+J:j:L:l:M+m+NnO:o:P:p:QqR:r:S:s:TtU:u:VvW:w:Xx:Zz'); $tick_info = ($opt_B || $opt_b); $xyfill = ($opt_G || $opt_g); $help = ($opt_H || $opt_h); $file_data = ($opt_I || $opt_i); $map_scale = ($opt_J || $opt_j); $labels = ($opt_L || $opt_l); $misc = ($opt_M || $opt_m); $xysegment = ($opt_N || $opt_n); $root = ($opt_O || $opt_o); $pagesize = ($opt_P || $opt_p); $no_view_ps = ($opt_Q || $opt_q); $bounds = ($opt_R || $opt_r); $xysymbol = ($opt_S || $opt_s); $coast_control = ($opt_T || $opt_t); $orientation = ($opt_U || $opt_u); $verbose = ($opt_V || $opt_v); $xypen = ($opt_W || $opt_w); $execute = ($opt_X || $opt_x); $delete_temp_files = ($opt_Z || $opt_z); # print out help message if required if ($help) { print "\n$program_name:\n"; print "\nVersion: $Id: mbm_xyplot.pl 1867 2010-09-15 05:30:50Z caress $\n"; print "\nMacro to generate a shellscript of GMT commands which, \n"; print "when executed, will generate a Postscript plot of xy data. \n"; print "Axes may be linear, log, or any of several geographic \n"; print "projections. Data may be plotted as symbols or lines. \n"; print "The plot will be scaled to fit on the specified page size \n"; print "or, if the scale is user defined, the page size will be \n"; print "chosen in accordance with the plot size. The primary purpose \n"; print "of this macro is to allow the simple, semi-automated \n"; print "production of nice looking plots with a few command line \n"; print "arguments. For users seeking more control over the plot \n"; print "appearance, a number of additional optional arguments are \n"; print "provided. Truly ambitious users may edit the plot shellscript \n"; print "to take advantage of GMT capabilites not supported by this \n"; print "macro. \n"; print "\nBasic Usage: \n"; print "\t$program_name -I[filepars:]file [ -Gfill -H -Oroot \n"; print "\t\t-Ppagesize -Ssymbol/size -Uorientation \n"; print "\t\t-V -Wpen ]\n"; print "Additional Options:\n"; print "\t\t[-Btickinfo -Jprojection[/scale | width]\n"; print "\t\t-Ltitle[:scale_label] -Mmisc -Q -Rw/e/s/n -X]\n"; print "Miscellaneous Options:\n"; print "\t\t[-MGDgmtdef/value -MGL[f][x]lon0/lat0/slat/length[m]\n"; print "\t\t-MGTx/y/size/angle/font/just/text\n"; print "\t\t-MGU[/dx/dy/][label] ]\n"; exit 0; } # check for input file if (!$file_data) { print "\a"; die "\nNo input file specified!\n$program_name aborted\n"; } # parse misc commands if ($misc) { @misc_cmd = split(/:::::::/, $misc); foreach $cmd (@misc_cmd) { # deal with general options ############################## # set GMT default values if ($cmd =~ /^[Gg][Dd]./) { ($gmt_def) = $cmd =~ /^[Gg][Dd](\S+)/; push(@gmt_defs, $gmt_def); } # set map scale if ($cmd =~ /^[Gg][Ll]./) { ($length_scale) = $cmd =~ /^[Gg][Ll](\S+)/; } # set text labels if ($cmd =~ /^[Gg][Tt]./) { ($tx, $ty, $tsize, $tangle, $font, $just, $txt) = $cmd =~ /^[Gg][Pp](\S+)\/(\S+)\/(\S+)\/(\S+)\/(\S+)\/(\S+)\/(.+)/; if ($txt) { ($text_info) = $cmd =~ /^[Gg][Tt](.*)/; push(@text, $text_info); } else { print "\nInvalid text label ignored: $cmd\n"; } } # set unix time stamp if ($cmd =~ /^[Gg][Uu]./) { ($unix_stamp) = $cmd =~ /^[Gg][Uu](\S+)/; $unix_stamp_on = 1; } elsif ($cmd =~ /^[Gg][Uu]/) { $unix_stamp_on = 1; } # deal with pscoast options ############################## # set pscoast lake fill if ($cmd =~ /^[Tt][Cc]./) { ($coast_lakefill) = $cmd =~ /^[Tt][Cc](.+)/; $coast_control = 1; } # set pscoast resolution if ($cmd =~ /^[Tt][Dd]./) { ($coast_resolution) = $cmd =~ /^[Tt][Dd](.+)/; } # set pscoast dry fill if ($cmd =~ /^[Tt][Gg]./) { ($coast_dryfill) = $cmd =~ /^[Tt][Gg](.+)/; $coast_control = 1; } # set pscoast rivers if ($cmd =~ /^[Tt][Ii]./) { ($coast_river) = $cmd =~ /^[Tt][Ii](.+)/; $coast_control = 1; } # set pscoast national boundaries if ($cmd =~ /^[Tt][Nn]./) { ($coast_boundary) = $cmd =~ /^[Tt][Nn](.+)/; push(@coast_boundaries, $coast_boundary); $coast_control = 1; } # set pscoast wet fill if ($cmd =~ /^[Tt][Ss]./) { ($coast_wetfill) = $cmd =~ /^[Tt][Ss](.+)/; $coast_control = 1; } # set pscoast coastline pen if ($cmd =~ /^[Tt][Ww]./) { ($coast_pen) = $cmd =~ /^[Tt][Ww](.+)/; $coast_control = 1; } } } # set page size if (!$pagesize) { $pagesize = "a"; } else { $pagesize =~ tr/A-Z/a-z/; if (!$page_width_in{$pagesize}) { $pagesize = "a"; } } # get postscript viewer # check environment variable if ($ENV{"MB_PS_VIEWER"}) { $ps_viewer = $ENV{"MB_PS_VIEWER"}; } # check for .mbio_defaults file if (!$ps_viewer) { $home = $ENV{"HOME"}; $mbdef = "$home/.mbio_defaults"; if (open(MBDEF,"<$mbdef")) { while () { if (/ps viewer:\s+(\S+)/) { ($ps_viewer) = /ps viewer:\s+(\S+)/; } } close MBDEF; } } # else just set it to ghostview if (!$ps_viewer) { $ps_viewer = "ghostview"; } # set default xy control values if (!$xysymbol) { $xysymbol = "N"; } if (!$xyfill) { $xyfill = "N"; } if (!$xypen) { $xypen = "N"; } if (!$xysegment) { $xysegment = "N"; } else { $xysegment = "Y"; } # Some notes about how data files are processed: # mbm_xyplot will plot multiple data sets on a single plot using # either of two mechanisms. The first method is to specify multiple # input file flags and associated arguments ( e.g. mbm_xyplot # -i[fileargs]:file -i[fileargs]:file ...). The second method is to # use the utilize the multi-segemnt feature psxy. The -M flag allows # multiple "segments" of x,y pairs to be listed in a single file each # separated by a single line having a special character. The default # character is ">", and indeed that is the ONLY character supported # here. # When multiple flags are specified as described above, MBGetopts will # return their arguments concatinated by ":::::::". The first line # below splits those arguments so they may be processed separately for # each file. # Each file may have any of several optional arguments. Most are quite # straightforward, however the column flag (C) warrents some # explaination. The column flag allows the user to specify the columns # to be plotted by utilizing a simple syntax in which a column is # specified by c[N] where N deones the column number. The full syntax # is then c#[xcol]_c#[ycol] where the underbar delimits the syntax for # the x column and y columns respectively. @data_file_list = split(/:::::::/, $file_data); foreach $xyfile_raw (@data_file_list) { # set default parameters $xysymbol_f = $xysymbol; $xyfill_f = $xyfill; $xypen_f = $xypen; $xysegment_f = $xysegment; $delimiter_f = '\s+'; $xmath_f='$line[0]'; $ymath_f='$line[1]'; # figure out how many parameters to set for each file @filearg = split(/:/, $xyfile_raw); for ($i = 0; $i < $#filearg; $i++) { if ($filearg[$i] =~ /S(\S+)/) { ($xysymbol_f) = $filearg[$i] =~ /S(\S+)/; } if ($filearg[$i] =~ /G(\S+)/) { ($xyfill_f) = $filearg[$i] =~ /G(\S+)/; } if ($filearg[$i] =~ /W(\S+)/) { ($xypen_f) = $filearg[$i] =~ /W(\S+)/; } if ($filearg[$i] =~ /M/) { $xysegment_f = "Y"; } if ($filearg[$i]=~ /D/) { ($delimiter_f) = $filearg[$i] =~ /D(\S+)/; } if ($filearg[$i]=~ /C/) { ($xmath_f,$ymath_f) = $filearg[$i] =~ /C(\S+)_(\S+)/; # Turn the c[] expresssions used in the command line # arguments to ones that perl can evaluate, by # substituing '$line' for c and subtracting one from # the column number. $xmath_f =~ s/c\[(\d+)\]/\$line[$1-1]/g; $ymath_f =~ s/c\[(\d+)\]/\$line[$1-1]/g; $xmath_f =~ s/\#/\$linecnt/g; $ymath_f =~ s/\#/\$linecnt/g; } } # Push each argument into parallel arrays. push(@xysymbols, $xysymbol_f); push(@xyfills, $xyfill_f); push(@xypens, $xypen_f); push(@xysegments, $xysegment_f); push(@xyfiles, $filearg[$#filearg]); push(@delimiters, $delimiter_f); push(@xmath, $xmath_f); push(@ymath, $ymath_f); } # check that files are ok if (scalar(@xyfiles) <= 0) { print "\a"; die "\nNo input file specified!\n$program_name aborted\n"; } foreach $xyfile (@xyfiles) { if (! -e $xyfile) { print "\a"; die "\nInput file $xyfile cannot be opened!\n$program_name aborted\n"; } } # set output root if needed if (!$root) { $root = $xyfiles[0]; } # Manually extract the data from a larger file. $i=0; foreach $xyfile(@xyfiles){ $INFILE="<$xyfile"; $OUTFILE=">$xyfile$i.$$"; # Strip off the leading path and substitute the current working directory. $OUTFILE =~ s/^\>\//.*\/\>$current_working_dir/g; open INFILE or die "Cannot open $xyfile: $'"; open OUTFILE or die "Cannot open temporary file: $'"; $linecnt=1; while () { chomp; # Here we check to see if a line has the standard psxy segment # separater. If it does, we print print it out directly to the # output file, and continue processing lines. if ( $_ =~ m/\>/ ) { print(OUTFILE "$_\n"); }else{ # Here we check to see if the line starts with a # delimiter. If it does, all the first column will be # empty. For files with white-space delimiters, this is # counter intuitive. So we remove the leading delimiter # when it occurs by substituting nothing for it. if ( $_ =~ /^$delimiters[$i]/ ) { $_ =~ s/^$delimiters[$i]//; } @line = split /$delimiters[$i]/, $_; # Evaluate the perl/math expressions. If there's an error, # don't quite but warn to STDIO. (should probably go to # STDERR). Divide-by-zero errors will be caught by this. my $xval = eval $xmath[$i]; print "WARNING!!! NON-NUMERIC RESULT DETECTED: $@" if $@; my $yval = eval $ymath[$i]; print "WARNING!!! NON-NUMERIC RESULT DETECTED: $@" if $@; # Verify that we got numbers... if (! ($xval =~ m/-?\d*\.?\d*|-?\.\d+/ && $yval =~ m/-?\d*\.?\d*|-?\.\d+/) ) { print "WARNING!!! NON-NUMERIC RESULT DETECTED! X: $xval Y: $yval Skipping...\n"; } push(@xvalues, $xval); push(@yvalues, $yval); print(OUTFILE "$xval\t$yval\n"); } $linecnt+=1; } $i+=1; close INFILE; close OUTFILE; # Reset the file names to the new temporarily files. $OUTFILE =~ s/^>//; $xyfile="$OUTFILE"; } # get limits of files by sorting all the data. if (!$bounds) { @sorted_xvalues = sort {$a <=> $b} @xvalues; @sorted_yvalues = sort {$a <=> $b} @yvalues; $xmin = $sorted_xvalues[0]; $xmax = $sorted_xvalues[$#sorted_xvalues]; $ymin = $sorted_yvalues[0]; $ymax = $sorted_yvalues[$#sorted_yvalues]; } # use user defined bounds if ($bounds) { if ($bounds =~ /^\S+\/\S+\/\S+\/\S+r$/) { ($xmin_raw,$ymin_raw,$xmax_raw,$ymax_raw) = $bounds =~ /^(\S+)\/(\S+)\/(\S+)\/(\S+)r$/; $xmin = &GetDecimalDegrees($xmin_raw); $xmax = &GetDecimalDegrees($xmax_raw); $ymin = &GetDecimalDegrees($ymin_raw); $ymax = &GetDecimalDegrees($ymax_raw); $use_corner_points = 1; $bounds_plot = $bounds; } elsif ($bounds =~ /^\S+\/\S+\/\S+\/\S+$/) { ($xmin_raw,$xmax_raw,$ymin_raw,$ymax_raw) = $bounds =~ /(\S+)\/(\S+)\/(\S+)\/(\S+)/; $xmin = &GetDecimalDegrees($xmin_raw); $xmax = &GetDecimalDegrees($xmax_raw); $ymin = &GetDecimalDegrees($ymin_raw); $ymax = &GetDecimalDegrees($ymax_raw); $bounds_plot = $bounds; } elsif ($bounds =~ /^r$/) { $use_corner_points = 1; $bounds_plot = sprintf ("%1.8g/%1.8g/%1.8g/%1.8gr", $xmin, $ymin, $xmax, $ymax); } } # set bounds string for plotting if not already set if (!$bounds_plot) { $bounds_plot = sprintf ("%1.8g/%1.8g/%1.8g/%1.8g", $xmin, $xmax, $ymin, $ymax); } ## This checking no longer makes sense when the data could be for an xy plot ## with one variable not changing. ie min=max. # check that there is data #if ($xmin >= $xmax || $ymin >= $ymax) # { # print "\a"; # die "Improper data limits: x: $xmin $xmax y: $ymin $ymax\n$program_name aborted.\n" # } # set the relevent page width and height &GetPageSize; # get user constraints on map scale if ($map_scale) { # sets $plot_scale or $plot_width if possible &GetProjection; } # set up for mapproject if (($use_scale && $plot_scale) || ($use_width && $plot_width)) { ($projection) = $map_scale =~ /^(\w)/; ($projection_pars) = $map_scale =~ /^$projection(\S+)/; } elsif ($use_scale || $use_width) { ($projection) = $map_scale =~ /^(\w)/; ($projection_pars) = $map_scale =~ /^$projection(\S+)/; $projection_pars = "$projection_pars" . "$separator" . "$trial_value"; } else { $projection = "x"; $projection_pars = "1/1"; $use_scale = 1; } # now find out the apparent size of the plot `echo $xmin $ymin > tmp$$.dat`; `echo $xmax $ymin >> tmp$$.dat`; `echo $xmax $ymax >> tmp$$.dat`; `echo $xmin $ymax >> tmp$$.dat`; @projected = `mapproject tmp$$.dat -J$projection$projection_pars -R$bounds_plot 2>&1 `; `/bin/rm -f tmp$$.dat`; print "mapproject tmp$$.dat -J$projection$projection_pars -R$bounds_plot\n"; while (@projected) { $line = shift @projected; print $line; if(!$xxmin){ ($xxmin,$yymin) = $line =~ /(\S+)\s+(\S+)/; $xxmax = $xxmin; $yymax = $yymin; } else { ($xx,$yy) = $line =~ /(\S+)\s+(\S+)/; $xxmin = ($xx < $xxmin ? $xx : $xxmin); $xxmax = ($xx > $xxmax ? $xx : $xxmax); $yymin = ($yy < $yymin ? $yy : $yymin); $yymax = ($yy > $yymax ? $yy : $yymax); } } $dxx = $xxmax - $xxmin; $dyy = $yymax - $yymin; # check for valid scaling if ($dxx == 0.0 && $dyy == 0.0) { print "\a"; die "Invalid projection specified - $program_name aborted\n"; } # figure out scaling issues if (($use_scale && $plot_scale) || ($use_width && $plot_width)) { $plot_width = $dxx; $plot_height = $dyy; # decide which plot orientation to use if ($orientation == 1) { $portrait = 1; $width = $page_width_in{$pagesize}; $height = $page_height_in{$pagesize}; $width_max = $width_max_portrait; $height_max = $height_max_portrait; } elsif ($orientation == 2) { $landscape = 1; $width = $page_height_in{$pagesize}; $height = $page_width_in{$pagesize}; $width_max = $width_max_landscape; $height_max = $height_max_landscape; } elsif ($dxx > $dyy) { $landscape = 1; $width = $page_height_in{$pagesize}; $height = $page_width_in{$pagesize}; $width_max = $width_max_landscape; $height_max = $height_max_landscape; } else { $portrait = 1; $width = $page_width_in{$pagesize}; $height = $page_height_in{$pagesize}; $width_max = $width_max_portrait; $height_max = $height_max_portrait; } # check if plot fits on page if ($plot_width > $width_max || $plot_height > $height_max) { # try to find a sufficiently large pagesize $pagesize_save = $pagesize; foreach $elem (@page_size_names) { $pagesize = "$elem"; &GetPageSize; if ($portrait) { $width_max = $width_max_portrait; $height_max = $height_max_portrait; } else { $width_max = $width_max_landscape; $height_max = $height_max_landscape; } if (!$good_page && $plot_width <= $width_max && $plot_height <= $height_max) { $good_page = 1; $pagesize_save = $pagesize; } } # print out warning if ($pagesize eq $pagesize_save) { print "\nWarning: Unable to fit plot on any available page size!\n"; print "\tThis plot will not be particularly useful!\n"; print "\tTry using a different scale or allow the program to set the scale!\n"; } # reset the page size $pagesize = $pagesize_save; &GetPageSize; if ($portrait) { $width = $page_width_in{$pagesize}; $height = $page_height_in{$pagesize}; } else { $width = $page_height_in{$pagesize}; $height = $page_width_in{$pagesize}; } } } elsif ($use_scale && $projection =~ /^x.*/ && !$geographic) { # set size of plot according to orientation if ($orientation == 1) { $portrait = 1; $plot_scale_x = $width_max_portrait/$dxx; $plot_scale_y = $height_max_portrait/$dyy; $plot_scale = "$plot_scale_x/$plot_scale_y"; $width = $page_width_in{$pagesize}; $height = $page_height_in{$pagesize}; } else { $landscape = 1; $plot_scale_x = $width_max_landscape/$dxx; $plot_scale_y = $height_max_landscape/$dyy; $plot_scale = "$plot_scale_x/$plot_scale_y"; $width = $page_height_in{$pagesize}; $height = $page_width_in{$pagesize}; } # set plot width $plot_width = $dxx * $plot_scale_x; $plot_height = $dyy * $plot_scale_y; # construct plot scale parameters ## These lines break things. $projection_pars was already properly ## defined and $map_scale may not yet be defined. Moreover, ## $projection_pars should now contain JUST the scale values ## with no leading projection type. The single line added ## is all that is required. ## Val Schmidt LDEO/Columbia OCTOBER 2003 #($projection_pars) = $map_scale =~ /^$projection(\S+)/; # $projection_pars = sprintf ("$projection_pars$separator%1.5g", # $plot_scale); $projection_pars="$plot_scale"; } elsif ($use_width && $projection =~ /^X.*/ && !$geographic) { # set size of plot according to orientation if ($orientation == 1) { $portrait = 1; $plot_width_x = $width_max_portrait; $plot_width_y = $height_max_portrait; $plot_width = "$plot_width_x/$plot_width_y"; $width = $page_width_in{$pagesize}; $height = $page_height_in{$pagesize}; } else { $landscape = 1; $plot_width_x = $width_max_landscape; $plot_width_y = $height_max_landscape; $plot_width = "$plot_width_x/$plot_width_y"; $width = $page_height_in{$pagesize}; $height = $page_width_in{$pagesize}; } # set plot width $plot_scale_x = $plot_width / $dxx; $plot_scale_y = $plot_height / $dyy; # construct plot scale parameters ## These lines break things. $projection_pars was already properly ## defined and $map_scale may not yet be defined. Moreover, ## $projection_pars should now contain JUST the scale values ## with no leading projection type. The single line added ## is all that is required. ## Val Schmidt LDEO/Columbia OCTOBER 2003 #($projection_pars) = $map_scale =~ /^$projection(\S+)/; #$projection_pars = sprintf ("$projection_pars$separator%1.5g", # $plot_width); $projection_pars="$plot_width"; } elsif ($use_scale) { # get landscape and portrait scales $plot_scale_landscape = $width_max_landscape/$dxx; if ($plot_scale_landscape*$dyy > $height_max_landscape) { $plot_scale_landscape = $height_max_landscape/$dyy; } $plot_scale_portrait = $width_max_portrait/$dxx; if ($plot_scale_portrait*$dyy > $height_max_portrait) { $plot_scale_portrait = $height_max_portrait/$dyy; } # decide which plot orientation to use if ($orientation == 1) { $portrait = 1; $plot_scale = $plot_scale_portrait; $width = $page_width_in{$pagesize}; $height = $page_height_in{$pagesize}; } elsif ($orientation == 2) { $landscape = 1; $plot_scale = $plot_scale_landscape; $width = $page_height_in{$pagesize}; $height = $page_width_in{$pagesize}; } elsif ($plot_scale_landscape > $plot_scale_portrait) { $landscape = 1; $plot_scale = $plot_scale_landscape; $width = $page_height_in{$pagesize}; $height = $page_width_in{$pagesize}; } else { $portrait = 1; $plot_scale = $plot_scale_portrait; $width = $page_width_in{$pagesize}; $height = $page_height_in{$pagesize}; } # set plot width $plot_width = $dxx * $plot_scale; $plot_height = $dyy * $plot_scale; # reset plot_scale if ratio required if ($use_ratio) { $top = int(1 / $plot_scale); $plot_scale = "1:" . "$top"; } # construct plot scale parameters ($projection_pars) = $map_scale =~ /^$projection(\S+)/; $projection_pars = sprintf ("$projection_pars$separator%1.5g", $plot_scale); # handle special case for linear projections if ($geographic) { $projection_pars = "$projection_pars" . "d"; } } elsif ($use_width) { # get landscape and portrait heights $plot_width_landscape = $height_max_landscape * $dxx / $dyy; if ($plot_width_landscape > $width_max_landscape) { $plot_width_landscape = $width_max_landscape; } $plot_width_portrait = $height_max_portrait * $dxx / $dyy; if ($plot_width_portrait > $width_max_portrait) { $plot_width_portrait = $width_max_portrait; } # decide which plot orientation to use if ($orientation == 1) { $portrait = 1; $plot_width = $plot_width_portrait; $plot_height = $plot_width * $dyy / $dxx; $width = $page_width_in{$pagesize}; $height = $page_height_in{$pagesize}; } elsif ($orientation == 2) { $landscape = 1; $plot_width = $plot_width_landscape; $plot_height = $plot_width * $dyy / $dxx; $width = $page_height_in{$pagesize}; $height = $page_width_in{$pagesize}; } elsif ($plot_width_landscape > $plot_width_portrait) { $landscape = 1; $plot_width = $plot_width_landscape; $plot_height = $plot_width * $dyy / $dxx; $width = $page_height_in{$pagesize}; $height = $page_width_in{$pagesize}; } else { $portrait = 1; $plot_width = $plot_width_portrait; $plot_height = $plot_width * $dyy / $dxx; $width = $page_width_in{$pagesize}; $height = $page_height_in{$pagesize}; } # construct plot scale parameters ($projection_pars) = $map_scale =~ /^$projection(\S+)/; $projection_pars = sprintf ("$projection_pars$separator%1.5g", $plot_width); # handle special case for linear projections if ($geographic) { $projection_pars = "$projection_pars" . "d"; } } # place the origin so plot is more or less centered $xoffset = ($width - $plot_width - $space_left - $space_right) / 2 + $space_left; $yoffset = ($height - $plot_height - $space_bottom - $space_top) / 2 + $space_bottom; # get plot degree annotation for geographic maps $degree_format = "ddd:mm"; if (!$gridprojected) { $xsize = $xmax - $xmin; $ysize = $ymax - $ymin; if ($xsize < $ysize) { $size = $xsize; } else { $size = $ysize; } if ($size > 1.0) { $degree_format = "ddd"; } elsif ($size > (1.0 / 60.0)) { $degree_format = "ddd:mm"; } else { $degree_format = "ddd:mm:ss"; } } # set pscoast control if ($coast_control && !$coast_wetfill && !$coast_dryfill && !$coast_pen && !$coast_boundary && !$coast_river) { $coast_dryfill = "200"; $coast_pen = "1p"; } if ($coast_control && !$coast_resolution) { $coast_resolution = "f"; } # come up with the filenames $cmdfile = "$root.cmd"; $psfile = "$root.ps"; $gmtfile = "gmtdefaults4\$\$"; # set some gmtisms $first_gmt = 1; $first = "-X\$X_OFFSET -Y\$Y_OFFSET -K -V >! \$PS_FILE"; $middle = "-K -O -V >> \$PS_FILE"; $end = "-O -V >> \$PS_FILE"; # set macro gmt default settings $gmt_def = "MEASURE_UNIT/inch"; push(@gmt_macro_defs, $gmt_def); $gmt_def = "PAPER_MEDIA/$page_gmt_name{$pagesize}+"; push(@gmt_macro_defs, $gmt_def); $gmt_def = "ANOT_FONT/Helvetica"; push(@gmt_macro_defs, $gmt_def); $gmt_def = "LABEL_FONT/Helvetica"; push(@gmt_macro_defs, $gmt_def); $gmt_def = "HEADER_FONT/Helvetica"; push(@gmt_macro_defs, $gmt_def); $gmt_def = "ANOT_FONT_SIZE/$page_anot_font{$pagesize}"; push(@gmt_macro_defs, $gmt_def); $gmt_def = "LABEL_FONT_SIZE/$page_anot_font{$pagesize}"; push(@gmt_macro_defs, $gmt_def); $gmt_def = "HEADER_FONT_SIZE/$page_header_font{$pagesize}"; push(@gmt_macro_defs, $gmt_def); $gmt_def = "FRAME_WIDTH/$frame_size"; push(@gmt_macro_defs, $gmt_def); $gmt_def = "TICK_LENGTH/$tick_size"; push(@gmt_macro_defs, $gmt_def); $gmt_def = "PAGE_ORIENTATION/LANDSCAPE"; push(@gmt_macro_defs, $gmt_def); $gmt_def = "COLOR_BACKGROUND/0/0/0"; push(@gmt_macro_defs, $gmt_def); $gmt_def = "COLOR_FOREGROUND/255/255/255"; push(@gmt_macro_defs, $gmt_def); $gmt_def = "COLOR_NAN/255/255/255"; push(@gmt_macro_defs, $gmt_def); $gmt_def = "PLOT_DEGREE_FORMAT/$degree_format"; push(@gmt_macro_defs, $gmt_def); # open the shellscript file if (!open(FCMD,">$cmdfile")) { print "\a"; die "Cannot open output file $cmdfile\nMacro $program_name aborted.\n"; } # write the shellscript header print FCMD "#! /bin/csh -f\n"; print FCMD "#\n# Shellscript to create Postscript plot of data in grd file\n"; print FCMD "# Created by macro $program_name\n"; print FCMD "#\n# This shellscript created by following command line:\n"; print FCMD "# $program_name $command_line\n"; # Define shell variables print FCMD "#\n# Define shell variables used in this script:\n"; print FCMD "set PS_FILE = $psfile\n"; print FCMD "set CPT_FILE = $cptfile\n"; print FCMD "set MAP_PROJECTION = $projection\n"; print FCMD "set MAP_SCALE = $projection_pars\n"; print FCMD "set MAP_REGION = $bounds_plot\n"; printf FCMD "set X_OFFSET = %1.5g\n", $xoffset; printf FCMD "set Y_OFFSET = %1.5g\n", $yoffset; # Reset GMT defaults, saving old defaults print FCMD "#\n# Save existing GMT defaults\n"; print FCMD "echo Saving GMT defaults...\n"; print FCMD "gmtdefaults -L >! $gmtfile\n"; print FCMD "#\n# Set new GMT defaults\n"; print FCMD "echo Setting new GMT defaults...\n"; foreach $gmt_def (@gmt_macro_defs) { ($gmt_par, $gmt_var) = $gmt_def =~ /^([^\/]+)\/(.+)/; print FCMD "gmtset $gmt_par $gmt_var\n"; } # Reset GMT defaults as per user commands if (@gmt_defs) { print FCMD "#\n# Set user defined GMT defaults\n"; print FCMD "echo Setting user defined GMT defaults...\n"; foreach $gmt_def (@gmt_defs) { ($gmt_par, $gmt_var) = $gmt_def =~ /^([^\/]+)\/(.+)/; print FCMD "gmtset $gmt_par $gmt_var\n"; } } # figure out labels $nlabels = 0; if ($labels) { @labels_split = split(/:/, $labels); $nlabels = @labels_split; if ($nlabels > 0) { $tlabel = shift(@labels_split); $tlabel =~ s/\\/\//g; } if ($nlabels > 1) { $xlabel = shift(@labels_split); $xlabel =~ s/\\/\//g; } if ($nlabels > 2) { $ylabel = shift(@labels_split); $ylabel =~ s/\\/\//g; } } if ($nlabels < 1) { if ($#xyfiles == 0) { $tlabel = "Data File $xyfiles[0]"; } else { $tlabel = "Data Files"; foreach $xyfile (@xyfiles) { $tlabel = $tlabel . " $xyfile"; } } } if ($nlabels < 2) { $xlabel = " "; } if ($nlabels < 3) { $ylabel = " "; } # set basemap axes annotation if ($tick_info) { $axes = $tick_info; if (!($tick_info =~ /.*:\..*/)) { $axes = "$axes:.\"$tlabel\":"; } } elsif ($projection =~ /^[Xx].*/ && !$geographic) { # figure out some reasonable tick intervals for the basemap &GetBaseTickLinear; $axes = "$base_tick_x:\"$xlabel\":/$base_tick_y:\"$ylabel\"::.\"$tlabel\":"; } else { # figure out some reasonable tick intervals for the basemap &GetBaseTick; $axes = "$base_tick:\"$xlabel\":/$base_tick:\"$ylabel\"::.\"$tlabel\":"; } # do xy plots for ($i = 0; $i < scalar(@xyfiles); $i++) { printf FCMD "#\n# Make xy data plot\n"; printf FCMD "echo Running psxy...\n"; printf FCMD "psxy $xyfiles[$i] \\\n\t"; printf FCMD "-J\$MAP_PROJECTION\$MAP_SCALE \\\n\t"; printf FCMD "-R\$MAP_REGION \\\n\t"; if ($xysymbols[$i] ne "N") { printf FCMD "-S$xysymbols[$i] \\\n\t"; } if ($xypens[$i] ne "N") { printf FCMD "-W$xypens[$i] \\\n\t"; } if ($xyfills[$i] ne "N") { printf FCMD "-G$xyfills[$i] \\\n\t"; } if ($xysegments[$i] ne "N") { printf FCMD "-M \\\n\t"; } if ($portrait) { printf FCMD "-P "; } if ($first_gmt == 1) { $first_gmt = 0; printf FCMD "$first\n"; } else { printf FCMD "$middle\n"; } } # do coastline plots if ($coast_control) { printf FCMD "#\n# Make coastline data plot\n"; printf FCMD "echo Running pscoast...\n"; printf FCMD "pscoast \\\n\t"; printf FCMD "-J\$MAP_PROJECTION\$MAP_SCALE \\\n\t"; printf FCMD "-R\$MAP_REGION \\\n\t"; if ($coast_lakefill) { printf FCMD "-C$coast_lakefill \\\n\t"; } if ($coast_resolution) { printf FCMD "-D$coast_resolution \\\n\t"; } if ($coast_dryfill) { printf FCMD "-G$coast_dryfill \\\n\t"; } if ($coast_river) { printf FCMD "-I$coast_river \\\n\t"; } for ($i = 0; $i < scalar(@coast_boundaries); $i++) { printf FCMD "-N$coast_boundaries[$i] \\\n\t"; } if ($coast_wetfill) { printf FCMD "-S$coast_wetfill \\\n\t"; } if ($coast_pen) { printf FCMD "-W$coast_pen \\\n\t"; } if ($portrait) { printf FCMD "-P "; } if ($first_gmt == 1) { $first_gmt = 0; printf FCMD "$first\n"; } else { printf FCMD "$middle\n"; } } # do pstext plot if (@text) { printf FCMD "#\n# Make text labels\n"; printf FCMD "echo Running pstext...\n"; printf FCMD "pstext -J\$MAP_PROJECTION\$MAP_SCALE \\\n\t"; printf FCMD "-R\$MAP_REGION \\\n\t"; printf FCMD "$middle < 1.0) { $pagescale = 1.0; } $view_pageflag = "--orientation=portrait --media=BBox -scale=$pagescale"; } elsif ($landscape) { $pagescale = 11.0 / $page_width_in{$pagesize}; if ($pagescale > 1.0) { $pagescale = 1.0; } $view_pageflag = "--orientation=landscape --media=BBox -scale=$pagescale"; } } elsif ($ps_viewer eq "ggv") { if ($portrait) { $view_pageflag = "--geometry=portrait"; } elsif ($landscape) { $view_pageflag = "--geometry=landscape"; } } if ($no_view_ps) { print FCMD "#echo Running $ps_viewer in background...\n"; print FCMD "#$ps_viewer $view_pageflag \$PS_FILE &\n"; } else { print FCMD "echo Running $ps_viewer in background...\n"; print FCMD "$ps_viewer $view_pageflag \$PS_FILE &\n"; } # claim it's all over print FCMD "#\n# All done!\n"; print FCMD "echo All done!\n"; # now close the shellscript and make it executable close FCMD; chmod 0775, $cmdfile; # tell program status if ($verbose) { print "\nProgram Status:\n"; print "--------------\n"; print "\n GMT Version:\n"; print " Version $gmt_version\n"; print "\n Plot Style:\n"; if (@xyfiles) { print " XY Plots of ", scalar(@xyfiles), " Datasets\n"; } if ($coast_control) { print " Coastline\n"; } if (@text) { print " ", scalar(@text), " Text labels\n"; } if ($length_scale) { print " Map distance scale\n"; } if ($unix_stamp_on && $unix_stamp) { print " Unix time stamp: $unix_stamp\n"; } elsif ($unix_stamp_on) { print " Unix time stamp\n"; } print "\n Input Files:\n"; foreach $xyfile (@xyfiles) { print " XY Data File: $xyfile\n"; } print "\n Output Files:\n"; print " Output plot name root: $root\n"; print " Plotting shellscript: $cmdfile\n"; print " Plot file: $psfile\n"; print "\n Plot Attributes:\n"; printf " Plot width: %.4f\n", $plot_width; printf " Plot height: %.4f\n", $plot_height; print " Page size: $pagesize\n"; print " Page width: $width\n"; print " Page height: $height\n"; print " Projection: -J$projection$projection_pars\n"; print " Axes annotation: $axes\n"; if ($portrait) { print " Orientation: portrait\n"; } else { print " Orientation: landscape\n"; } print "\n XY Data Attributes:\n"; printf " X min max: %9.4g %9.4g\n", $xmin, $xmax; printf " Y min max: %9.4g %9.4g\n", $ymin, $ymax; if (@xyfiles) { print "\n Primary XY Plotting Controls:\n"; printf " symbol pen fill segment file\n"; printf " ------ --- ---- ------- ----\n"; for ($i = 0; $i < scalar(@xyfiles); $i++) { printf " %-10s %-10s %-10s %-10s %s\n", $xysymbols[$i], $xypens[$i], $xyfills[$i], $xysegments[$i], $xyfiles[$i]; } } if ($coast_control) { print "\n Coastline Plotting Controls:\n"; } if ($coast_control && $coast_resolution) { printf " Coastline resolution: $coast_resolution\n"; } if ($coast_control && $coast_pen) { printf " Coastline pen: $coast_pen\n"; } if ($coast_control && $coast_wetfill) { printf " Ocean fill: $coast_wetfill\n"; } if ($coast_control && $coast_lakefill) { printf " Lake fill: $coast_lakefill\n"; } if ($coast_control && $coast_dryfill) { printf " Land fill: $coast_dryfill\n"; } if ($coast_control && $coast_river) { printf " Rivers: $coast_river\n"; } if ($coast_control && $coast_boundaries) { for ($i = 0; $i < scalar(@coast_boundaries); $i++) { printf " National Boundaries: $coast_boundaries[$i]\n"; } } if ($length_scale) { print "\n Miscellaneous Plotting Controls:\n"; } if ($length_scale) { printf " Length scale: $length_scale\n"; } print "\n GMT Default Values Reset in Script:\n"; foreach $gmt_def (@gmt_macro_defs) { ($gmt_par, $gmt_var) = $gmt_def =~ /^([^\/]+)\/(.+)/; printf " %-25s %s\n", $gmt_par, $gmt_var; } foreach $gmt_def (@gmt_defs) { ($gmt_par, $gmt_var) = $gmt_def =~ /^([^\/]+)\/(.+)/; printf " $gmt_par : $gmt_var\n"; } print "\n--------------\n"; } # print out final notes print "\nPlot generation shellscript <$cmdfile> created.\n"; print "\nInstructions:\n"; print " Execute <$cmdfile> to generate Postscript plot <$psfile>.\n"; if (!$no_view_ps) { print " Executing <$cmdfile> also invokes $ps_viewer "; print "to view the plot on the screen.\n"; } if ($verbose) { print "\n--------------\n\n"; } # execute shellscript if desired if ($execute) { if ($verbose) { print "Executing shellscript $cmdfile...\n"; } system "$cmdfile &"; } exit 0; #----------------------------------------------------------------------- sub min { # make local variables local ($min); # get the minimum of the arguments if ($_[0] < $_[1]) { $min = $_[0]; } else { $min = $_[1]; } $min; } #----------------------------------------------------------------------- sub max { # make local variables local ($max); # get the minimum of the arguments if ($_[0] > $_[1]) { $max = $_[0]; } else { $max = $_[1]; } $max; } #----------------------------------------------------------------------- sub GetDecimalDegrees { # make local variables local ($dec_degrees, $degrees, $minutes, $seconds); # deal with dd:mm:ss format if ($_[0] =~ /^\S+:\S+:\S+$/) { ($degrees, $minutes, $seconds) = $_[0] =~ /^(\S+):(\S+):(\S+)$/; if ($degrees =~ /^-\S+/) { $dec_degrees = $degrees - $minutes / 60.0 - $seconds / 3600.0; } else { $dec_degrees = $degrees + $minutes / 60.0 + $seconds / 3600.0; } } # deal with dd:mm format elsif ($_[0] =~ /^\S+:\S+$/) { ($degrees, $minutes) = $_[0] =~ /^(\S+):(\S+)$/; if ($degrees =~ /^-\S+/) { $dec_degrees = $degrees - $minutes / 60.0; } else { $dec_degrees = $degrees + $minutes / 60.0; } } # value already in decimal units of some sort else { $dec_degrees = $_[0]; } # return decimal degrees; $dec_degrees; } #----------------------------------------------------------------------- sub GetPageSize { # get space around edge of plot $space_top = 1.25 * $page_height_in{$pagesize} / $page_height_in{"a"}; $space_bottom = 1.50 * $page_height_in{$pagesize} / $page_height_in{"a"}; $space_left = 1.00 * $page_height_in{$pagesize} / $page_height_in{"a"}; $space_right = 1.00 * $page_height_in{$pagesize} / $page_height_in{"a"}; # set the relevent page width and height $width_max_landscape = $page_height_in{$pagesize} - $space_left - $space_right; $height_max_landscape = $page_width_in{$pagesize} - $space_bottom - $space_top; $width_max_portrait = $page_width_in{$pagesize} - $space_left - $space_right; $height_max_portrait = $page_height_in{$pagesize} - $space_bottom - $space_top; $frame_size = 0.075; if (($frame_size / $height_max_portrait) > 0.01) { $frame_size = 0.01 * $height_max_portrait; } $tick_size = 0.075; if (($tick_size / $height_max_portrait) > 0.01) { $tick_size = 0.01 * $height_max_portrait; } } #----------------------------------------------------------------------- sub GetProjection { # get the map projection flag ($projection) = $map_scale =~ /^(\w)/; ($projection_pars) = $map_scale =~ /^$projection(\S+)/; # see if plot scale or plot width defined $use_scale = 0; $use_width = 0; $separator = "/"; $trial_value = 1.0; # Cassini Projection if ($projection eq "c") { ($plot_scale) = $map_scale =~ /^c\S+\/\S+\/(\S+)$/; $use_scale = 1; } elsif ($projection eq "C") { ($plot_width) = $map_scale =~ /^C\S+\/\S+\/(\S+)$/; $use_width = 1; } # Mercator Projection elsif ($projection eq "m") { ($plot_scale) = $map_scale =~ /^m(\S+)$/; $use_scale = 1; $separator = ""; } elsif ($projection eq "M") { ($plot_width) = $map_scale =~ /^M(\S+)$/; $use_width = 1; $separator = ""; } # Oblique Mercator Projection elsif ($projection eq "o") { if ($map_scale =~ /^oa\S+/) { ($plot_scale) = $map_scale =~ /^oa\S+\/\S+\/\S+\/(\S+)$/; } elsif ($map_scale =~ /^ob\S+/) { ($plot_scale) = $map_scale =~ /^ob\S+\/\S+\/\S+\/\S+\/(\S+)$/; } elsif ($map_scale =~ /^oc\S+/) { ($plot_scale) = $map_scale =~ /^oc\S+\/\S+\/\S+\/\S+\/(\S+)$/; } $use_scale = 1; } elsif ($projection eq "O") { if ($map_scale =~ /^Oa\S+/) { ($plot_width) = $map_scale =~ /^Oa\S+\/\S+\/\S+\/(\S+)$/; } elsif ($map_scale =~ /^Ob\S+/) { ($plot_width) = $map_scale =~ /^Ob\S+\/\S+\/\S+\/\S+\/(\S+)$/; } elsif ($map_scale =~ /^Oc\S+/) { ($plot_width) = $map_scale =~ /^Oc\S+\/\S+\/\S+\/\S+\/(\S+)$/; } $use_width = 1; } # Equidistant Cylindrical Projection elsif ($projection eq "q") { ($plot_scale) = $map_scale =~ /^q\S+\/(\S+)$/; $use_scale = 1; } elsif ($projection eq "Q") { ($plot_width) = $map_scale =~ /^Q\S+\/(\S+)$/; $use_width = 1; } # Transverse Mercator Projection elsif ($projection eq "t") { ($plot_scale) = $map_scale =~ /^t\S+\/(\S+)$/; $use_scale = 1; } elsif ($projection eq "T") { ($plot_width) = $map_scale =~ /^T\S+\/(\S+)$/; $use_width = 1; } # Universal Transverse Mercator Projection elsif ($projection eq "u") { ($plot_scale) = $map_scale =~ /^u\S+\/(\S+)$/; $use_scale = 1; } elsif ($projection eq "U") { ($plot_width) = $map_scale =~ /^U\S+\/(\S+)$/; $use_width = 1; } # Basic Cylindrical Projection elsif ($projection eq "y") { ($plot_scale) = $map_scale =~ /^y\S+\/\S+\/(\S+)$/; $use_scale = 1; } elsif ($projection eq "Y") { ($plot_width) = $map_scale =~ /^Y\S+\/\S+\/(\S+)$/; $use_width = 1; } # Lambert Azimuthal Projection elsif ($projection eq "a") { ($plot_scale) = $map_scale =~ /^a\S+\/\S+\/(\S+)$/; $use_scale = 1; $trial_value = "1:1"; $use_ratio = 1; } elsif ($projection eq "A") { ($plot_width) = $map_scale =~ /^A\S+\/\S+\/(\S+)$/; $use_width = 1; } # Equidistant Projection elsif ($projection eq "e") { ($plot_scale) = $map_scale =~ /^e\S+\/\S+\/(\S+)$/; $use_scale = 1; $trial_value = "1:1"; $use_ratio = 1; } elsif ($projection eq "E") { ($plot_width) = $map_scale =~ /^E\S+\/\S+\/(\S+)$/; $use_width = 1; } # Orthographic Projection elsif ($projection eq "g") { ($plot_scale) = $map_scale =~ /^g\S+\/\S+\/(\S+)$/; $use_scale = 1; $trial_value = "1:1"; $use_ratio = 1; } elsif ($projection eq "G") { ($plot_width) = $map_scale =~ /^G\S+\/\S+\/(\S+)$/; $use_width = 1; } # General Sterographic Projection elsif ($projection eq "s") { ($plot_scale) = $map_scale =~ /^s\S+\/\S+\/(\S+)$/; $use_scale = 1; $trial_value = "1:1"; $use_ratio = 1; } elsif ($projection eq "S") { ($plot_width) = $map_scale =~ /^S\S+\/\S+\/(\S+)$/; $use_width = 1; } # Albers Projection elsif ($projection eq "b") { ($plot_scale) = $map_scale =~ /^b\S+\/\S+\/\S+\/\S+\/(\S+)$/; $use_scale = 1; } elsif ($projection eq "B") { ($plot_width) = $map_scale =~ /^B\S+\/\S+\/\S+\/\S+\/(\S+)$/; $use_width = 1; } # Lambert Conic Projection elsif ($projection eq "l") { ($plot_scale) = $map_scale =~ /^l\S+\/\S+\/\S+\/\S+\/(\S+)$/; $use_scale = 1; } elsif ($projection eq "L") { ($plot_width) = $map_scale =~ /^L\S+\/\S+\/\S+\/\S+\/(\S+)$/; $use_width = 1; } # Hammer Projection elsif ($projection eq "h") { ($plot_scale) = $map_scale =~ /^h\S+\/(\S+)$/; $use_scale = 1; } elsif ($projection eq "H") { ($plot_width) = $map_scale =~ /^H\S+\/(\S+)$/; $use_width = 1; } # Sinusoidal Projection elsif ($projection eq "i") { ($plot_scale) = $map_scale =~ /^i\S+\/(\S+)$/; $use_scale = 1; } elsif ($projection eq "I") { ($plot_width) = $map_scale =~ /^I\S+\/(\S+)$/; $use_width = 1; } # Eckert VI Projection elsif ($projection eq "k") { ($plot_scale) = $map_scale =~ /^k\S+\/(\S+)$/; $use_scale = 1; } elsif ($projection eq "K") { ($plot_width) = $map_scale =~ /^K\S+\/(\S+)$/; $use_width = 1; } # Robinson Projection elsif ($projection eq "n") { ($plot_scale) = $map_scale =~ /^n\S+\/(\S+)$/; $use_scale = 1; } elsif ($projection eq "N") { ($plot_width) = $map_scale =~ /^N\S+\/(\S+)$/; $use_width = 1; } # Winkel Tripel Projection elsif ($projection eq "r") { ($plot_scale) = $map_scale =~ /^r\S+\/(\S+)$/; $use_scale = 1; } elsif ($projection eq "R") { ($plot_width) = $map_scale =~ /^R\S+\/(\S+)$/; $use_width = 1; } # Mollweide Projection elsif ($projection eq "w") { ($plot_scale) = $map_scale =~ /^w\S+\/(\S+)$/; $use_scale = 1; } elsif ($projection eq "W") { ($plot_width) = $map_scale =~ /^W\S+\/(\S+)$/; $use_width = 1; } # Linear Polar Projection elsif ($projection eq "p") { ($plot_scale) = $map_scale =~ /^p(\S+)$/; $use_scale = 1; $separator = ""; } elsif ($projection eq "P") { ($plot_width) = $map_scale =~ /^P(\S+)$/; $use_width = 1; $separator = ""; } # Linear Projection elsif ($projection eq "x") { if ($map_scale =~ /^xd$/) { $geographic = 1; chop($map_scale); } else { ($plot_scale) = $map_scale =~ /^x(\S+)$/; } $use_scale = 1; $separator = ""; } elsif ($projection eq "X") { if ($map_scale =~ /^Xd$/) { $geographic = 1; chop($map_scale); } else { ($plot_width) = $map_scale =~ /^X(\S+)$/; } $use_width = 1; $separator = ""; } } #----------------------------------------------------------------------- sub GetBaseTickLinear { # figure out some reasonable tick intervals for the basemap $dxxx = $xmax - $xmin; if ($dxxx > 0) { $base = int((log($dxxx) / log(10.)) + 0.5); $base_tick_x = (10 ** $base) / 5; if ($dxxx / $base_tick_x < 5) { $base_tick_x = $base_tick_x / 4; } elsif ($dxxx / $base_tick_x < 10) { $base_tick_x = $base_tick_x / 2; } } $dyyy = $ymax - $ymin; if ($dyyy > 0) { $base = int((log($dyyy) / log(10.)) + 0.5); $base_tick_y = (10 ** $base) / 5; if ($dyyy / $base_tick_y < 5) { $base_tick_y = $base_tick_y / 4; } elsif ($dyyy / $base_tick_y < 10) { $base_tick_y = $base_tick_y / 2; } } } #----------------------------------------------------------------------- sub GetBaseTick { # figure out some reasonable tick intervals for the basemap $base_tick_x = ($xmax - $xmin) / 5; $base_tick_y = ($ymax - $ymin) / 5; $base_tick = &min($base_tick_x, $base_tick_y); if ($base_tick < 0.0002777777) { $base_tick = "1c"; } elsif ($base_tick < 0.0005555555) { $base_tick = "2c"; } elsif ($base_tick < 0.0013888889) { $base_tick = "5c"; } elsif ($base_tick < 0.0027777778) { $base_tick = "10c"; } elsif ($base_tick < 0.0041666667) { $base_tick = "15c"; } elsif ($base_tick < 0.0083333333) { $base_tick = "30c"; } elsif ($base_tick < 0.0166667) { $base_tick = "1m"; } elsif ($base_tick < 0.0333333) { $base_tick = "2m"; } elsif ($base_tick < 0.0833333) { $base_tick = "5m"; } elsif ($base_tick < 0.1666667) { $base_tick = "10m"; } elsif ($base_tick < 0.25) { $base_tick = "15m"; } elsif ($base_tick < 0.5) { $base_tick = "30m"; } elsif ($base_tick < 1.0) { $base_tick = "1"; } elsif ($base_tick < 2.0) { $base_tick = "2"; } elsif ($base_tick < 5.0) { $base_tick = "5"; } elsif ($base_tick < 10.0) { $base_tick = "10"; } elsif ($base_tick < 15.0) { $base_tick = "15"; } elsif ($base_tick < 30.0) { $base_tick = "30"; } elsif ($base_tick < 30.0) { $base_tick = "30"; } elsif ($base_tick < 360.0) { $base_tick = "60"; } } #----------------------------------------------------------------------- # This version of Getopts has been augmented to support multiple # calls to the same option. If an arg in argumentative is followed # by "+" rather than ":", then the corresponding scalar will # be concatenated rather than overwritten by multiple calls to # the same arg. # # Usage: # do Getopts('a:b+c'); # -a takes arg, -b concatenates args, # # -c does not take arg. Sets opt_* as a # # side effect. sub MBGetopts { local($argumentative) = @_; local(@args,$_,$first,$rest); local($errs) = 0; local($[) = 0; @args = split( / */, $argumentative ); while(@ARGV && ($_ = $ARGV[0]) =~ /^-(.)(.*)/) { ($first,$rest) = ($1,$2); $pos = index($argumentative,$first); if($pos >= $[) { if($args[$pos+1] eq ':') { shift(@ARGV); if($rest eq '') { ++$errs unless @ARGV; $rest = shift(@ARGV); } eval "\$opt_$first = \$rest;"; eval "\$flg_$first = 1;"; } elsif($args[$pos+1] eq '+') { shift(@ARGV); if($rest eq '') { ++$errs unless @ARGV; $rest = shift(@ARGV); } if (eval "\$opt_$first") { eval "\$opt_$first = \$opt_$first . \":::::::\" . \$rest;"; } else { eval "\$opt_$first = \$rest;"; } eval "\$flg_$first = 1;"; } elsif($args[$pos+1] eq '%') { shift(@ARGV); if($rest ne '') { eval "\$opt_$first = \$rest;"; } else { $rest = $ARGV[0]; ($one) = $rest =~ /^-(.).*/; $pos = index($argumentative,$one); if(!$one || $pos < $[) { eval "\$opt_$first = \$rest;"; shift(@ARGV); } } eval "\$flg_$first = 1;"; } else { eval "\$opt_$first = 1"; eval "\$flg_$first = 1;"; if($rest eq '') { shift(@ARGV); } else { $ARGV[0] = "-$rest"; } } } else { print STDERR "Unknown option: $first\n"; ++$errs; if($rest ne '') { $ARGV[0] = "-$rest"; } else { shift(@ARGV); } } } $errs == 0; } #-----------------------------------------------------------------------