#!/bin/sh
#
# This file is part of Rheolef.
#
# Copyright (C) 2000-2009 Pierre Saramito 
#
# Rheolef 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.
#
# Rheolef 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 Rheolef; if not, write to the Free Software
# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
#
# -------------------------------------------------------------------------

#Prog:qmg2geo
#NAME: qmg2geo - convert qmg mesh in geo format
#@cindex mesh
#@pindex qmg2geo
#@pindex geo
#@fiindex @file{.qmg} qmg mesh
#@fiindex @file{.dmn} domain names
#@fiindex @file{.geo} mesh
#@toindex @code{qmg}
#DESCRIPTION:
#Convert a qmg @file{.qmg} into @file{.geo} one.
#The output goes to standart output.
#
#SYNOPSIS:
#@example
#  qmg2geo [-qmg] @var{input}.qmg @var{input}.dmn
#@end example
#
#DOMAIN NAME FILE:
#This auxilliary @file{.dmn} file defines the boundary domain names
#as used by Rheolef, since @code{qmg} uses numeric
#labels for domains.
#@example
#  EdgeDomainNames
#    4
#    bottom
#    right
#    top
#    left
#@end example
#
#End:

#============================================================
# Usage
#============================================================

PROG=`expr $0 : '.*/\(.*\)' \| $0`

USAGE="usage: $PROG input[.qmg] input.dmn"

#============================================================
# Arguments
#============================================================

if test $# = 0; then
    echo $USAGE 1>&2; exit 0
fi

tmpdir=${TMPDIR-"/tmp"}
qmg_file=""
geobin="geo"
while test $# -ne 0; do
    case $1 in
	-help)        echo $USAGE 1>&2; exit 0;;
      	*.qmg)        qmg_file=$1
		      ;;
      	-qmg)        if test x"$2" = x""; then
			echo "$PROG: no argument to -qmg option found." 1>&2
		        echo $USAGE 1>&2; exit 1
                      fi
      		      qmg_file=$2
		      shift
		      ;;
        *.dmn)        dmn_file=$1;;
        -dmn)        if test x"$2" = x""; then
                        echo "$PROG: no argument to -dmn option found." 1>&2
                        echo $USAGE 1>&2; exit 1
                      fi
                      dmn_file=$2
                      shift
                      ;;
	-geobin)      if test x"$2" = x""; then
			echo "$PROG: no argument to -geobin option found." 1>&2
			echo $USAGE 1>&2; exit 1
		      fi
                      geobin=$2
		      shift
		      ;;
	*) 	      
		      echo "$PROG: invalid $1 option." 1>&2
		      echo $USAGE 1>&2; exit 1
		      ;;
    esac
    shift
done

if test x"${qmg_file}" = x""; then
   echo $USAGE 1>&2
   exit 1
fi

#============================================================
# buid the `.geo' mesh file
#============================================================

cat ${qmg_file} ${dmn_file} | col | ${geobin} -input-qmg -upgrade -geo -noverbose -
status=$?
if test $status -ne 0; then
    echo "$PROG(fatal): ${geobin} exit with $status status." 1>&2
    exit $status
fi

