tengo un programa que ya está listo pero lo que falta es que cuando quiero hacer un cálculo pues me manda a otra página ya que en la misma no puedo ya que es por php y lo que quiero es hacer el programa dentro de la misma página donde yo apriete calcular y me arroje de inmediato los resultados en la misma página...
este es el boton
<td <input type=submit name=calcular value='Calcular' onClick="window.open('02.php','UTMS','width=800, height=200, menubar=no, scrollbars=yes' )> </td>

estoy traduciendo un codigo matemático desde excel a php y me sale este símbolo que creo que es elebado y no se como hacerlo en php (^):
ejemplo.
$ni=(6421524,911/(1+0,006739497*(cos($rdlat))^2)^(1/2)*0.9996;

tengo un programa que se conecta a la base de datos para rescatar datos y en medio hay un boton calcular que me dirige a otra página donde se calcula el problema es que cuando lo hace pue me actualiza la página y me borra los datos anteriores
 si alguien me ayuda para ver que puedo hacer se lo agradesco montones
 los datos van adjuntos....

aqui si esta definido el form y las variables pero no me muestra resultado alguno que hago???

<?php
session_register('r');
include "verifica.php";
?>

<?php
$calcular=$_POST['calcula'];
if ($calcular) {
   
              }
       if ($orientalong=="1") { /* Si la orientacion de la Latitud esta al oeste*/
      $sd=(-1*((($segundoslong/60)/60)+($minutoslong/60)+($longitud)));/*Formula para hallar las sexas decimales*/
        $rd=$sd*3.1416/180;
        $aux=$sd/6+31;
        $huso=(int)$aux;
        $merhuso=($huso*6)-183;
        $delta=$rd-($merhuso*3.1416)/180;
        $a=(cos($rdlat))*(sin($delta));
        $xi=(1/2)*log((1+$a)/(1-$a));
        $eta=atan((tan($rdlat))/cos($delta))-$rdlat;
                              }

        else if ($orientalong=="2") {/* Si la orientacion de la Latitud esta al este del meridiano de Greenwich*/      
                 $sd=((($segundoslong/60)/60)+($minutoslong/60)+($longitud));
                $rd=$sd*3.1416/180;/*en radianes*/
                 $aux=$sd/6+31;/*Calculo huso*/
                 $huso=(int)$aux;;/*Calculo huso*/
                 $merhuso=($huso*6)-183;
                 $delta=$rd-($merhuso*3.1416)/180;
                 $a=cos($rdlat)*sin($delta);
                 $xi=(1/2)*log((1+$a)/(1-$a));
                 $eta=atan((tan($rdlat))/cos($delta))-$rdlat;
                                   }
   
        if ($orientalat=="1") { /*Calculando los valores de l Latitud*/
        $sd1=(-1*((($segundoslat/60)/60)+($minutoslat/60)+($latitud)));/*Formula para hallar las sexas decimales*/
        $rdlat=$sd1*3.1416/180;/*en radianes*/
        $eta=atan((tan($rdlat))/cos($delta))-$rdlat;
        $asemiejemayor=6378057.238;
        $bsemiejemenor=6334883.8;
        $a1=($asemiejemayor*$asemiejemayor);
        $b1=($bsemiejemenor*$bsemiejemenor);
        $c1=($a1-$b1)/$asemiejemayor;
        $excentricidad=sqrt((pow($asemiejemayor,2))-(pow($bsemiejemenor,2)))/$asemiejemayor;
        $segundaexcentricidad=sqrt((pow($asemiejemayor,2))-(pow($bsemiejemenor,2)))/$bsemiejemenor;
        $edos=pow($segundaexcentricidad,2);
        $rpc=pow($asemiejemayor,2)/$bsemiejemenor; /*radio polar de curvatura */
                /*$nuno=($rpc/(1+$edos(cos($rdlat)))*/
        $cero=$rpc/(1+$edos);
        $r=sqrt(cos(pow($rdlat,2)));
        $ni=($cero+$r)*0.9996;/*Ni*/
        $r1=($edos/2)*(pow($xi,2));
        $r2=(cos($rdlat));
        $r3=(pow($r2,2));
        $zeta=($r1*$r3);
        $res=2*$rdlat;/*Calculo de A1*/
        $auno=sin($res);
        $auxdos=cos($rdlat);/*Calculod e A2*/
        $auxtres=(pow($auxdos,2));
        $ados=$auno*$auxtres;
        $jdos=$rdlat+($auno/2);/*Calculod e J2*/
        $jcuatro=((3*$jdos)+$ados)/4;/*Calculod e J4*/
        $jseis1=5*$jcuatro+$ados;
        $jseis2=pow($rdlat,2);/*Calculod e J6 verificarlo*/
        $jseis3=cos($jseis2)/3;
        $jseis=$jseis1*$jseis3;
        $alfa=(3/4)*$edos;/*Calculod de Alfa/
        $beta=(5/3)*(pow($alfa,2));/*Calculod de Beta*/
        $gama=(35/27)*(pow($alfa,3));/*Calculod de Gama*/
        $bfi=(0.9996*$rpc)*($rdlat-($alfa*$jdos)+($beta*$jcuatro)-($gama*$jseis));/*Calculod de B(fi)*/
        $utme=($xi*$ni)*(1+($zeta/3))+500000;
        $utmnorte=($eta*$ni)*(1+$zeta)+$bfi+10000000;/*si la Latitues Norte*/
                                   }
   
        else if ($orientalat=="2") {
                 $sd1=((($segundoslat/60)/60)+($minutoslat/60)+($latitud));
                 $rdlat=$sd1*3.1416/180;/*en radianes*/
                 $eta=atan((tan($rdlat))/cos($delta))-$rdlat;
                 $asemiejemayor=6378057.238;
                 $bsemiejemenor=6334883.8;
                 $a1=($asemiejemayor*$asemiejemayor);
                 $b1=($bsemiejemenor*$bsemiejemenor);
                 $c1=($a1-$b1)/$asemiejemayor;
                 $excentricidad=sqrt((pow($asemiejemayor,2))-(pow($bsemiejemenor,2)))/$asemiejemayor;
                 $segundaexcentricidad=sqrt((pow($asemiejemayor,2))-(pow($bsemiejemenor,2)))/$bsemiejemenor;
                 $edos=pow($segundaexcentricidad,2);
                 $rpc=pow($asemiejemayor,2)/$bsemiejemenor; /*radio polar de curvatura */
                           /*$nuno=($rpc/(1+$edos(cos($rdlat)))*/
                 $cero=$rpc/(1+$edos);
                 $r=sqrt(cos(pow($rdlat,2)));
                 $ni=($cero+$r)*0.9996;/*Ni*/
                 $r1=($edos/2)*(pow($xi,2));
                 $r2=(cos($rdlat));
                 $r3=(pow($r2,2));
                 $zeta=($r1*$r3);
                 $res=2*$rdlat;/*Calculo de A1*/
                 $auno=sin($res);
                 $auxdos=cos($rdlat);/*Calculod e A2*/
                 $auxtres=(pow($auxdos,2));
                 $ados=$auno*$auxtres;
                 $jdos=$rdlat+($auno/2);/*Calculod e J2*/
                 $jcuatro=((3*$jdos)+$ados)/4;/*Calculod e J4*/
                 $jseis1=5*$jcuatro+$ados;
                 $jseis2=pow($rdlat,2);/*Calculod e J6 verificarlo*/
                 $jseis3=cos($jseis2)/3;
                 $jseis=$jseis1*$jseis3;
                 $alfa=(3/4)*$edos;/*Calculod de Alfa/
                 $beta=(5/3)*(pow($alfa,2));/*Calculod de Beta*/
                 $gama=(35/27)*(pow($alfa,3));/*Calculod de Gama*/
                 $bfi=(0.9996*$rpc)*($rdlat-($alfa*$jdos)+($beta*$jcuatro)-($gama*$jseis));/*Calculod de B(fi)*/
                 $utme=($xi*$ni)*(1+($zeta/3))+500000;
                 $utmnorte=(($eta*$ni)*(1+$zeta)+$bfi); /*Si la latitud es sur*/

                                        }
?>
<html>
<head>
<title>Ingreso  Incendio</title>
<LINK href="estilo.css" rel=STYLESHEET type="text/css">
<script type="text/javascript" src="funciones.js"></script>
<script type="text/javascript" src="funciones_antecedentes.js"></script>
<script type="text/javascript" src="select_dependientes.js"></script>
<style type="text/css">
<!--
.Estilo1 {color: #FF9900}
body,td,th {
   font-size: 10px;
}
body {
   background-image: url(../imagenes/fondo955x600.jpg);
}
-->
</style>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1"></head>
<?
if (!$_POST)
{
$sql="SELECT max(id_incendio) as id_incendio FROM antecedentes_generales";
$result=mysql_query($sql,$conexio);
if ($row= mysql_fetch_array($result))
{
$id_incendio=$row ["id_incendio"];
}
$sql="SELECT max(id_mensaje_incendio) as id_mensaje_incendio FROM mensaje_incendio";
$result=mysql_query($sql,$conexio);
if ($row= mysql_fetch_array($result))
{
$id_mensaje_incendio=$row ["id_mensaje_incendio"];
}
$sql="SELECT max(id_temporada) as id_temporada FROM temporada";
$result=mysql_query($sql,$conexio);
if ($row= mysql_fetch_array($result))
{
$id_temporada=($row ["id_temporada"]);
}
$sql="SELECT  temporada FROM temporada where id_temporada='$id_temporada'";
$result=mysql_query($sql,$conexio);
if ($row= mysql_fetch_array($result))
{
$temporada=($row ["temporada"]);
}
$sql="SELECT  max(situacion_nroreg) as situacion_nroreg FROM antecedentes_generales where id_temporada='$id_temporada' and id_region='$r' ";
$result=mysql_query($sql,$conexio);
if ($row= mysql_fetch_array($result))
{
$nroreg=($row ["situacion_nroreg"]);
}

$fecha_ac=date("20y-m-d");
$sql="SELECT  max(situacion_nrodias) as situacion_nrodias FROM antecedentes_generales where id_temporada='$id_temporada' and id_region='$r' and situacion_fecha='$fecha_ac'";
$result=mysql_query($sql,$conexio);
if ($row= mysql_fetch_array($result))
{
$nro_del_dia=($row ["situacion_nrodias"]);
}

?>
<body link="#000033" vlink="#000033" alink="#000033" onLoad="document.form.provincias.focus()">
<div id="Layer1" style="position:absolute; left:13px; top:5px; width:914px; height:24px; z-index:1"> <a href="/sidco/ingreso.php" title="Panel de control" target="_parent"><img src="../imagenes/panel%20de%20control.gif" width="120" height="22" border="0" style="cursor:pointer"></a>  <a href="/sidco/login.php" title="Cerrar Sesión" target="_parent"><img src="../imagenes/cerrar%20sesion.gif" width="120" height="22" border="0"></a>  <a href="/sidco/sub_menu_recursos.php" title="Recursos" target="_parent"><img src="../imagenes/recurso.gif" width="120" height="22" border="0"></a>    <a href="/sidco/foco_reg/ingreso.php" title="Foco" target="_parent"><img src="../imagenes/foco.gif" width="120" height="22" border="0"></a>    <a href="/sidco/mensaje_reg/ingreso.php" title="Mensaje" target="_parent"><img src="../imagenes/mensaje.gif" width="120" height="22" border="0"></a>       <a href="/sidco/sub_menu.php" title="Aperturas" target="_parent"><img src="../imagenes/aperturas.gif" width="120" height="22" border="0"></a>        <a href="buscar.php" title="Buscar" target="_parent"><img src="../imagenes/32538.ico" border="0" style="cursor:pointer"></a></div>
<br>
<br>
<br>
<FORM  method="post" ACTION="" name="form" onSubmit="return comprobar()">
 
   <table width="823" border="0"  align="left" cellpadding="0" >
    <tr>
      <td colspan="10"><strong>INGRESO  INCENDIO </strong></td>
    <tr>
      <td colspan="10"><hr></td>
    <tr>
      <td colspan="10"><strong>Antecedentes Generales </strong></td>
    <tr>
      <td width="124">Id</td>
      <td width="80"><input name="id_incendio" type="text" class="Estilo1" id="id_incendio2" value="<?echo $id_incendio + 1?>" size="10" readonly=""></td>
      <td width="61">Periodo</td>
      <td width="94"><input name="id_temporada" type="hidden" id="id_temporada" value="<?echo $id_temporada?>" size="10">
          <input name="temporada" type="text" id="temporada" value="<?echo $temporada?>" size="15" readonly="">          </td>
<td width="101">Provincia</td>
      <td width="158"><?php

   include 'conexion.php';
   conectar();
   $consulta=mysql_query("SELECT id_provincia, provincia_nom, id_region FROM provincia where id_region='$r' ");
   desconectar();

   // Voy imprimiendo el primer select compuesto por los paises
   echo "<select name='provincias' id='provincias'    onChange='cargaContenido(this.id)'>";
   echo "<option value=' '>Elige</option>";
   while($registro=mysql_fetch_row($consulta))
   {
      echo "<option value='".$registro[0]."'>".$registro[1]."</option>";
   }
   echo "</select>";
?>
      Comuna</td>
      <td colspan="4"><select  name="id_comuna"    id="select">
        <option value="0">Seleccion</option>
      </select></td>
     <tr>
      <td>N&ordm;
      Regional</td>
      <td><input name="situacion_nroreg" type="text" id="situacion_nroreg2" value="<?echo $nroreg + 1?>" size="10" readonly=""></td>
      <td>N&ordm; Local</td>
      <td><input name="situacion_nrolocal" type="text" id="situacion_nrolocal"  title="Elegir Provincia antes de hacer Clic "onClick="periodo()"  size="15"></td>
      <td>N&ordm; del Dia</td>
      <td colspan="5"><input name="situacion_nrodias" type="text" id="situacion_nrodias" value="<?echo $nro_del_dia+1?>" size="10" readonly=""></td>
     <tr>
      <td>Fecha</td>
      <td><input name="situacion_fecha"  readonly="" type="text" id="situacion_fecha2" value="<?echo date ( "j-m-20y" )?>" size="15"></td>
      <td>Registro</td>
      <td><input name="situacion_hora"  readonly="" type="text" id="situacion_hora2" value="<?echo date( "G:i:00" )?>" size="10"></td>
      <td>Estado</td>
      <td><select name="situacion_estado" size="1"   id="situacion_estado">
        <option value="0">Selecci&oacute;n</option>
        <option value="Bajo Observacion">Bajo Observacion</option>
        <option value="En Combate">En Combate</option>
        <option value="Controlado">Controlado</option>
        <option value="En Faena De Liquidacion">En Faena De Liquidacion</option>
        <option value="Extinguido">Extinguido</option>
        <option value="Rebrote">Rebrote</option>
                        </select></td>
      <td width="74">Ambito Conaf </td>
      <td colspan="3"><select name="situacion_ambito" id="situacion_ambito">
        <option value="0">Seleccion</option>
        <option value="Si">Si</option>
        <option value="No">No</option>
      </select></td>
    </tr>
    <tr>
      <td> Nombre  incendio </td>
      <td colspan="9"><input name="situacion_nombre" type="text" id="situacion_nombre" size="50" maxlength="50">
      </td>
    </tr>
        <tr>
          <td>Georef</td>
      <td colspan="9"><input name="id_georef" type="text" id="id_georef" value="<?echo $id_georef?>" "size=15 maxlength=7"></td>
     </tr>
        <tr>
      <td>Latitud</td>
          <td><input name="latitud" type="text"  size="10" maxlength="13" value="" ></td>
          <td> Minutos </td>
          <td><input name="minutoslat" type="text"  size="3" maxlength="3" value=""></td>
          <td>Segundos </td>
          <td><input name="segundoslat" type="text"  size="3" maxlength="3" value=""></td>
<td>
<select name="orientalat" >
<option value="1">N</option>
<option value="2">S</option>
</select></td>
</td>     </tr>
          <tr>
          <td>Longitud </td>
          <td><input name="longitud" type="text" value="" size="10" maxlength="13" ></td>
           <td> Minutos</td>
           <td><input name="minutoslong" type="text"  size="3" value=""  maxlength="3" ></td>
           <td>Segundos </td>
           <td><input name="segundoslong" type="text"  size="3" maxlength="3" value=""></td>
<td>
<select name="orientalong" >
<option value="1">E</option>
<option value="2">W</option>
</select>
</td> </tr>
         <td>
     <input type=submit name='calcular' value='calcular'  style="border:0px solid; color:#FFFFFF; background-color:#005B99" title="Calcular UTM">
        </td>
          <td>&nbsp;</td>
          <td colspan="3">&nbsp;</td>
        </tr>
        <tr>
          <td>Coord. UTM </td>
          <td>&nbsp;</td>
          <td>UTM Norte </td>
          <td><input name="utmnorte" type="text"  size="10" maxlength="15" value="<?php echo $utmnorte ?>" ></td>
          <td>UTM Este</td>
          <td><input name="utme" type="text"  size="10" maxlength="15" value="<?php echo $utme ?>" ></td>
          <td>&nbsp;</td>
          <td colspan="3">&nbsp;</td>
        </tr>
        <tr>
          <td>Ubic. Geografica </td>
          <td colspan="3"><input name="situacion_ugeografica" type="text" id="situacion_ugeografica"  size="50" maxlength="50"></td>
        <td >Cinta</td>
        <td colspan="6"><select name="situacion_cinta" id="situacion_cinta">
            <option value="0">Seleccion</option>
            <option value="Si">Si</option>
            <option value="No">No</option>
          </select></td>
        </tr>
        <tr>
          <td>P. Origen Incendio </td>
          <td><input name="situacion_poincendio" type="text" id="situacion_poincendio" size="20" maxlength="50"></td>
          <td>&nbsp;</td>
          <td>&nbsp;</td>
          <td>Propietario</td>
          <td><input name="situacion_propietario" type="text" id="situacion_propietario2" size="30"></td>
          <td colspan="4">&nbsp;</td>
        </tr>
        <tr>
          <td>Comb. Inicial </td>
          <td colspan="2"><select name="situacion_combinicial" id="situacion_combinicial">
              <option value="0">Seleccion</option>
              <option value="Arbolado natural">01-Arbolado natural</option>
              <option value="Plantacion pino">02-Plantacion pino</option>
           <option value="Plantacion eucaliptus">03-Plantacion eucaliptus</option>
           <option value="Otras plantaciones">04-Otras plantaciones</option>
           <option value="Matorral">05-Matorral</option>
           <option value="Pastizal">06-Pastizal</option>
           <option value="Cultivo agricola">07-Cultivo agricola</option>
           <option value="Desecho">08-Desecho</option>
           <option value="Basural">09-Basural</option>
           <option value="Otros">00-Otros</option>
            </select></td>
          <td>&nbsp;</td>
          <td>Inicio Cercano A </td>
          <td><select name="situacion_inccerca" id="situacion_inccerca">
              <option value="0">Seleccion</option>
              <option value="Camino principal">01-Camino principal</option>
              <option value="Camino secundario">02-Camino secundario</option>
           <option value="Via ferrea">03-Via ferrea</option>
           <option value="Inst. Industriales">04-Inst. Industriales</option>
           <option value="Casa o refugio">05-Casa o refugio</option>
           <option value="Lugares picnic">06-Lugares picnic</option>
           <option value="Senderos">07-Senderos</option>
           <option value="Faenas agricolas">08-Faenas agricolas</option>
           <option value="Faenas forestales">09-Faenas forestales</option>
           <option value="Otros">00-Otros</option>
            </select></td>
          <td>Investigada</td>
          <td colspan="3"><select name="situacion_investigada" id="select2">
            <option value="0">Seleccion</option>
            <option value="Si">Si</option>
            <option value="No">No</option>
          </select></td>
     </tr>
        <tr>
          <td>Rec. Detecci&oacute;n</td>
          <td colspan="3"><input name="id_deteccion" type="hidden" id="id_deteccion"  size="5" >
          <input name="deteccion_clave" type="text" id="deteccion_clave"   onChange="verificarecurso()" onClick="document.form.id_deteccion.value='',document.form.nombres.value='',document.form.deteccion_clave.value='',document.form.deteccion_nom.value=''" size="50"      >
          <input name="deteccion_nom" type="hidden" id="deteccion_nom"></td>
          <td>Tipo Detecci&oacute;n </td>
          <td colspan="5"><input name="nombres" type="text" id="nombres"></td>
        </tr>
        <tr>
          <td>Cod. Causa </td>
          <td colspan="3"><input name="id_causa" type="text" id="id_causa" onChange="verificarcausa2()" onClick="document.form.id_causa.value='',document.form.causa_desc.value='',document.form.tipocausa_desc.value='',document.form.grupocausa_desc.value=''" value="<?echo $id_causa?>" size=15 maxlength=15></td>
          <td rowspan="4">Observación
         
          </td>
          <td colspan="5" rowspan="4"><textarea name="situacion_obs" cols="40" rows="7" id="situacion_obs"><?echo $foco_obs?></textarea></td>
        </tr>
        <tr>
          <td>Causa Especifica </td>
          <td colspan="3"><textarea name="causa_desc" cols="30" rows="2" id="causa_desc" onChange="verificarcausa()" onClick="document.form.causa_desc.value='',document.form.id_causa.value='',document.form.tipocausa_desc.value='',document.form.grupocausa_desc.value=''"></textarea>          </td>
        </tr>
        <tr>
          <td>Causa general </td>
          <td colspan="3"><input name="tipocausa_desc" type="text" id="tipocausa_desc" size="50" maxlength="50"  readonly=""></td>
        </tr>
        <tr>
          <td>Grupo Causa </td>
          <td colspan="3"><input name="grupocausa_desc" type="text" id="grupocausa_desc" size="50" maxlength="50"  readonly=""></td>
        </tr>
        <tr>
          <td colspan="8"><hr></td>
          <td width="45">&nbsp;</td>
          <td width="36">&nbsp;</td>
        </tr>
      <tr><td colspan="10"><div align="center">
         <input name="submit" type="submit" class="botones"  onClick="this.form.action='ingreso.php?'"  value="Guardar y Salir.......">
         &nbsp;
         <input name="submit2" type="submit" class="botones"  onClick="this.form.action='guardar_y_siguiente.php?'"  value=" Guardar y Siguiente"  >
         &nbsp;         &nbsp;
         <input name="reset" type="reset" class="botones" value="Limpiar" >
        </div></td>
    </tr>
  </table>
</FORM>

tengo un programa en html y quiero hacer un scrip donde al ingresar la coordenada de latitud y longitud me transforme a su utm correspondiente si alguien me ayuda porfa desde ya gracias...
este es el programa aunque esta medio desordenado pero el al principio es el programa y lo de despues es el scrip que quiero hacer
         esta es una parte del programa donde hay que hacer el scrip mas adelante es lo que encontre para realizar la conversion...

<tr>
<td>&nbsp;</td>
<td><div align="left">
  </div></td>
<td>Latitud</td>
<td><input name="georef_latitud" type="TEXT" id="georef_latitud"  value="" size="10" maxlength="13"></td>
<td>Longitud</td>
<td colspan="8"><input name="georef_longitud" type="TEXT" id="georef_longitud"   value="" size="10" maxlength="13"> 
<div align="left"></div></td>
<tr>
<td>Coordenadas UTM </td>
  <td>&nbsp;</td>
  <td>UTM Norte</td>
  <td><input name="georef_utmy" type=TEXT value="" size="10" maxlength="15" ></td>
  <td> UTM Este</td>
  <td colspan="8"><input name="georef_utmx" type="TEXT" id="georef_utmx"  value=""  size="10" maxlength="15" ></td>
  <tr>

                  codigo a enlazar con scrip


 *                              DEFINES
 */

#define PI           3.14159265358979323e0    /* PI                        */
#define MIN_LAT      ( (-80.5 * PI) / 180.0 ) /* -80.5 degrees in radians    */
#define MAX_LAT      ( (84.5 * PI) / 180.0 )  /* 84.5 degrees in radians     */
#define MIN_EASTING  100000
#define MAX_EASTING  900000
#define MIN_NORTHING 0
#define MAX_NORTHING 10000000
GLOBAL DECLARATIONS
 */

static double UTM_a = 6378137.0;         /* Semi-major axis of ellipsoid in meters  */
static double UTM_f = 1 / 298.257223563; /* Flattening of ellipsoid                 */
static long   UTM_Override = 0;          /* Zone override flag                      */
FUNCTIONS
 *
 */

long Set_UTM_Parameters(double a,     
                        double f,
                        long   override)
{a                 : Semi-major axis of ellipsoid, in meters       (input)
 *    f                 : Flattening of ellipsoid                              (input)
 *    override          : UTM override zone, zero indicates no override (input)
 */

  double inv_f = 1 / f;
  long Error_Code = UTM_NO_ERROR;

  if (a <= 0.0)
  { /* Semi-major axis must be greater than zero */
    Error_Code |= UTM_A_ERROR;
  }
  if ((inv_f < 250) || (inv_f > 350))
  { /* Inverse flattening must be between 250 and 350 */
    Error_Code |= UTM_INV_F_ERROR;
  }
  if ((override < 0) || (override > 60))
  {
    Error_Code |= UTM_ZONE_OVERRIDE_ERROR;
  }
  if (!Error_Code)
  { /* no errors */
    UTM_a = a;
    UTM_f = f;
    UTM_Override = override;
  }
  return (Error_Code);
} /* END OF Set_UTM_Parameters */


void Get_UTM_Parameters(double *a,
                        double *f,
                        long   *override)
{*    a                 : Semi-major axis of ellipsoid, in meters       (output)
 *    f                 : Flattening of ellipsoid                              (output)
 *    override          : UTM override zone, zero indicates no override (output)
 */

  *a = UTM_a;
  *f = UTM_f;
  *override = UTM_Override;
} /* END OF Get_UTM_Parameters */


long Convert_Geodetic_To_UTM (double Latitude,
                              double Longitude,
                              long   *Zone,
                              char   *Hemisphere,
                              double *Easting,
                              double *Northing)
{    Latitude          : Latitude in radians                 (input)
 *    Longitude         : Longitude in radians                (input)
 *    Zone              : UTM zone                            (output)
 *    Hemisphere        : North or South hemisphere           (output)
 *    Easting           : Easting (X) in meters               (output)
 *    Northing          : Northing (Y) in meters              (output)
long Lat_Degrees;
  long Long_Degrees;
  long temp_zone;
  long Error_Code = UTM_NO_ERROR;
  double Origin_Latitude = 0;
  double Central_Meridian = 0;
  double False_Easting = 500000;
  double False_Northing = 0;
  double Scale = 0.9996;

  if ((Latitude < MIN_LAT) || (Latitude > MAX_LAT))
  { /* Latitude out of range */
    Error_Code |= UTM_LAT_ERROR;
  }
  if ((Longitude < -PI) || (Longitude > (2*PI)))
  { /* Longitude out of range */
    Error_Code |= UTM_LON_ERROR;
  }
  if (!Error_Code)
  { /* no errors */
    if((Latitude > -1.0e-9) && (Latitude < 0))
      Latitude = 0.0;
    if (Longitude < 0)
      Longitude += (2*PI) + 1.0e-10;

    Lat_Degrees = (long)(Latitude * 180.0 / PI);
    Long_Degrees = (long)(Longitude * 180.0 / PI);

    if (Longitude < PI)
      temp_zone = (long)(31 + ((Longitude * 180.0 / PI) / 6.0));
    else
      temp_zone = (long)(((Longitude * 180.0 / PI) / 6.0) - 29);

    if (temp_zone > 60)
      temp_zone = 1;
    /* UTM special cases */
    if ((Lat_Degrees > 55) && (Lat_Degrees < 64) && (Long_Degrees > -1)
        && (Long_Degrees < 3))
      temp_zone = 31;
    if ((Lat_Degrees > 55) && (Lat_Degrees < 64) && (Long_Degrees > 2)
        && (Long_Degrees < 12))
      temp_zone = 32;
    if ((Lat_Degrees > 71) && (Long_Degrees > -1) && (Long_Degrees < 9))
      temp_zone = 31;
    if ((Lat_Degrees > 71) && (Long_Degrees > && (Long_Degrees < 21))
      temp_zone = 33;
    if ((Lat_Degrees > 71) && (Long_Degrees > 20) && (Long_Degrees < 33))
      temp_zone = 35;
    if ((Lat_Degrees > 71) && (Long_Degrees > 32) && (Long_Degrees < 42))
      temp_zone = 37;

    if (UTM_Override)
    {
      if ((temp_zone == 1) && (UTM_Override == 60))
        temp_zone = UTM_Override;
      else if ((temp_zone == 60) && (UTM_Override == 1))
        temp_zone = UTM_Override;
      else if ((Lat_Degrees > 71) && (Long_Degrees > -1) && (Long_Degrees < 42))
      {
        if (((temp_zone-2) <= UTM_Override) && (UTM_Override <= (temp_zone+2)))
          temp_zone = UTM_Override;
        else
          Error_Code = UTM_ZONE_OVERRIDE_ERROR;
      }
      else if (((temp_zone-1) <= UTM_Override) && (UTM_Override <= (temp_zone+1)))
        temp_zone = UTM_Override;
      else
        Error_Code = UTM_ZONE_OVERRIDE_ERROR;
    }
    if (!Error_Code)
    {
      if (temp_zone >= 31)
        Central_Meridian = (6 * temp_zone - 183) * PI / 180.0;
      else
        Central_Meridian = (6 * temp_zone + 177) * PI / 180.0;
      *Zone = temp_zone;
      if (Latitude < 0)
      {
        False_Northing = 10000000;
        *Hemisphere = 'S';
      }
      else
        *Hemisphere = 'N';
      Set_Transverse_Mercator_Parameters(UTM_a, UTM_f, Origin_Latitude,
                                         Central_Meridian, False_Easting, False_Northing, Scale);
      Convert_Geodetic_To_Transverse_Mercator(Latitude, Longitude, Easting,
                                              Northing);
      if ((*Easting < MIN_EASTING) || (*Easting > MAX_EASTING))
        Error_Code = UTM_EASTING_ERROR;
      if ((*Northing < MIN_NORTHING) || (*Northing > MAX_NORTHING))
        Error_Code |= UTM_NORTHING_ERROR;
    }
  } /* END OF if (!Error_Code) */
  return (Error_Code);
} /* END OF Convert_Geodetic_To_UTM */