#include <math.h>
#include "dimensiony.h"
#include "/data/25/gung/,eigfun/premy_st.h"
#include "/data/25/gung/,src/,header/const.h"


void corfac(nod,w,q,pm,xac,xf,xln)
int nod;
float w,*xac,*xf,*xln;
char q;
prem_st *pm;
{
float fac,xmu,xka,rt;

/*corrected by XDL 3/14/1996*/         
  fac=2.*log((double)(0.5*w)/PI)/PI;
/*  fac=2.*log10((double)(0.5*w)/PI)/PI; */
  xmu=pm->qmu[nod]*fac;
  *xln=1.+xmu;
  
  if(q=='S')
  { xka=fac*pm->qkappa[nod];
    rt=(4./3.)* (pm->a[nod]+pm->c[nod]-2.*pm->f[nod]+5.*pm->n[nod]+
        6.*pm->l[nod])/
       (8.*pm->a[nod]+3.*pm->c[nod]+4.*pm->f[nod]+8.*pm->l[nod]);
    *xf=1.+((1.-rt)*xka-.5*rt*xmu)/(1.-1.5*rt);
    *xac=1.+(1.-rt)*xka+rt*xmu;
  }
  else *xac=*xf=0.0;
}