/*** ^^A -*-C++-*- **********************************************/
/*	ltog 		29ix09					*/
/****************************************************************/
/*	Short Description :					*/
/*	Set LB and GB for Lorentz-Gauss conversion		*/
/****************************************************************/
/*
$Id: ltog,v 1.10.16.1 2006/03/16 16:27:06 svcu Exp $
*/

float tlb,tgb,taq;

GETCURDATA
tlb=0; tgb=0;
if (strlen(cmd) == 0) 
{
    GETFLOAT("Enter Lorentzian width to remove : ",tlb);
    GETFLOAT("Enter Gaussian width to add : ",tgb);
}
else
{
    if (2 != sscanf(cmd, "%f  %f", &tlb, &tgb))
    {
 
	STOPMSG("illegal input, two real numbers needed, Lorentzian width and Gaussian width");
    }
}

FETCHPAR("AQ",&taq);
d1=tgb;
tgb=0.4413*tlb/(taq*d1*d1);
tlb=-tlb;
STOREPAR("LB",tlb);
if (tlb>0.0)
{
	STOPMSG("illegal input, Lorentzian width must be positive");
}
if (tgb>1.0)
{
	STOPMSG("illegal input, attempted resolution enhancement too extreme, please increase Gaussian width");
}
STOREPAR("GB",tgb);
GM
FT
PK
QUITMSG("--- ltog finished ---");