/*** ^^A -*-C++-*- **********************************************/
/*	setgtc2d 		26xi09					*/
/****************************************************************/
/*	Short Description :					*/
/*	Set LB and GB for simple Gaussian multiplication in 2D	*/
/****************************************************************/
/*
$Id: setgtc2d,v 1.10.16.1 2006/03/16 16:27:06 svcu Exp $
*/

float tlb,tgb,tgtc,taq,tdeff,td;

GETCURDATA
tlb=-0.01; tgtc=0;
FETCHPARS("AQ",&taq);
FETCHPARS("TDEFF",&tdeff);
FETCHPARS("TD",&td);
if (&tdeff == 0){}
else{taq=taq*tdeff/td;}

if (strlen(cmd) == 0) 
{
  tgtc=taq*0.5;
}
else
{
    if (2 != sscanf(cmd, "%f", &tgtc))
    {
	STOPMSG("illegal input: real number needed, Gaussian time constant");
    }
}
if (tgtc<0)
{
	STOPMSG("illegal input: Gaussian time constant must be positive");
}
tgb=-3.141593*tlb*tgtc*tgtc/(taq*2);
STOREPAR("LB",tlb);
STOREPAR("GB",tgb);
STOREPAR("WDW",2);
QUIT;