Title of the Web Page

#!/usr/bin/env bash # # These functions can be used from any sh/bash script by specifying # . gmt_shell_functions.sh # in your script. Placing it in .bashrc makes the functions available # on the command line as well. See documentation for usage. # # Note: Several of the original functions are now obsolete due to # new functionality in GMT. These are still available for use # but no longer documented and are placed at the end of this file # under the OBSOLETE tag. #----GMT SHELL FUNCTIONS-------------------- # Creates a unique temp directory and points GMT_TMPDIR to it gmt_init_tmpdir () { export GMT_TMPDIR=$(mktemp -d ${TMPDIR:-/tmp}/gmt.XXXXXX) } # Remove the temp directory created by gmt_init_tmpdir gmt_remove_tmpdir () { rm -rf $GMT_TMPDIR unset GMT_TMPDIR } # Remove all files and directories in which the current process number is part of the file name gmt_cleanup() { rm -rf *$$* if [ $# -eq 1 ]; then rm -rf ${1}* fi } # Send a message to stderr gmt_message() { echo "$*" >&2 } # Print a message to stderr and exit gmt_abort() { echo "$*" >&2 exit } # Returns the number of fields or arguments gmt_get_nfields() { echo $* | awk '{print NF}' } # Same with backwards compatible name... gmt_nfields() { echo $* | awk '{print NF}' } # Returns the given field (arg 1) in current record (arg 2) # Must pass arg 2 inside double quotes to preserve it as one item gmt_get_field() { echo $2 | cut -f$1 -d ' ' } # Return w/e/s/n from given table file(s) # May also add -Idx/dy to round off answer gmt_get_region() { printf "%s/%s/%s/%s\n" $(gmt info -C $* | cut -d' ' -f1-4) } # Return the w/e/s/n from the header in grd file gmt_get_gridregion() { printf "%s/%s/%s/%s\n" $(gmt grdinfo -Cn $* | cut -d' ' -f1-4) } # Make output PostScript file name based on script base name gmt_set_psfile() { echo $(basename $1 '.sh').ps } # Make output PDF file name based on script base name gmt_set_pdffile() { echo $(basename $1 '.sh').pdf } # For animations: Create a lexically increasing file namestem (no extension) based on prefix and frame number # i.e., prefix_###### gmt_set_framename() { echo $1 $2 | awk '{printf "%s_%06d\n", $1, $2}' } # For animations: Increment frame counter by one gmt_set_framenext() { echo $(($1 + 1)) } # For KMZ: Package a bunch of *.kml [and *.png] into a single kmz gmt_build_kmz() { if [ $# -eq 0 ]; then cat << EOF >&2 gmt_build_kmz - Create a single KMZ file from many KML files usage: gmt_build_kmz -p [-r] *.kml [*.png] *.kml: The KML files you want to include in the KML file. If these link to local PNG files then list those too. -p Specify a prefix for naming the .kmz file -r Remove all provided files after KMZ is built [leave alone] EOF return fi if [ ! "X$1" = "X-p" ]; then echo "gmt_build_kmz: Requires -p as first argument" >&2 return fi remove=0 name=$2 ; shift ; shift; if [ "X$1" = "X-r" ]; then remove=1 shift fi mkdir -p kml; mv -f $* kml cat <<- EOF > doc.kml EOF ls kml/*.kml > /tmp/$$.lis while read file; do name=$(basename $file .kml) cat << EOF >> doc.kml $name $file EOF done < /tmp/$$.lis cat << EOF >> doc.kml EOF zip -rq9 $name.kmz doc.kml kml if [ $remove -eq 0 ]; then mv -f kml/* .. fi rm -rf kml doc.kml /tmp/$$.lis } # For animations: Build animated gif from stills gmt_build_gif() { if [ $# -eq 0 ]; then cat << EOF >&2 gmt_build_gif - Process stills to animated gif with convert usage: gmt_build_gif [-d ] [-l ] [-r ] is the prefix of the still images and the resulting movie -d Specify path to of directory with the stills [current directory] -l Specify number of times to loop, 0 for continuous [0] -r Set frame rate in images per second [24] -n Dry-run. Do not run convert command, just print it EOF return fi command -v ${GRAPHICSMAGICK-gm} >/dev/null 2>&1 || { echo "gmt_build_gif: Cannot find graphicsmagick in your path - exiting." >&2; return; } rate=24; dir=.; loop=0; dryrun=0 while [ $# -ne 1 ]; do case "$1" in "-d") dir=$2 ; shift ;; "-l") loop=$2 ; shift ;; "-n") dryrun=1 ;; "-r") rate=$2 ; shift ;; *) echo "gmt_build_gif: No such option ($1)" >&2 ;; esac shift done delay=$(gmt math -Q 100 $rate DIV RINT =) # Delay to nearest 1/100 s if [ $dryrun -eq 1 ]; then cat <<- EOF ${GRAPHICSMAGICK-gm} convert -delay $delay -loop $loop +dither "$dir/${1}_*.*" ${1}.gif EOF else ${GRAPHICSMAGICK-gm} convert -delay $delay -loop $loop +dither "$dir/${1}_*.*" ${1}.gif echo "gmt_build_gif: GIF file is called ${1}.gif" >&2 fi } # For animations: Build a m4v movie from stills gmt_build_movie() { if [ $# -eq 0 ]; then cat << EOF >&2 gmt_build_movie - Process stills to MP4 movie with ffmpeg Note: M4V Requires images to have an even pixel width usage: gmt_build_movie [-d ] [-n] [-r ] [-v] is the prefix of the still images and the resulting movie -d Specify path to of directory with the stills [current directory] -r Set frame rate in images per second [24] -n Dry-run. Do not run ffmpeg command, just print it -v Verbose. Give progress messages EOF return fi command -v ffmpeg >/dev/null 2>&1 || { echo "gmt_build_movie: Cannot find ffmpeg in your path - exiting." >&2; return; } rate=24; dir=.; dryrun=0; blabber=quiet while [ $# -ne 1 ]; do case "$1" in "-d") dir=$2 ; shift ;; "-n") dryrun=1 ;; "-r") rate=$2 ; shift ;; "-v") blabber=verbose ;; *) echo "gmt_build_movie: No such option ($1)" >&2 ;; esac shift done if [ $dryrun -eq 1 ]; then cat <<- EOF ffmpeg -loglevel $blabber -f image2 -pattern_type glob -framerate $rate -y -i "$dir/$1_*.*" -pix_fmt yuv420p ${1}.m4v EOF else ffmpeg -loglevel $blabber -f image2 -pattern_type glob -framerate $rate -y -i "$dir/$1_*.*" -pix_fmt yuv420p ${1}.m4v echo "gmt_build_movie: MP4 file is called ${1}.m4v" >&2 fi } # For animations: Build an animation script template gmt_movie_script() { if [ $# -eq 0 ]; then cat << EOF >&2 gmt_movie_script - Create template script for animation usage: gmt_movie_script [-c