void sjacg(long*,long,long,float[],float[],float*,long,float[],float[],long[],
      float[],long,long[],long);