extern void pgsah_c(fint *,float *,float *);