PARAMETER (NL=9,NLWF=9,NLCH=9,NVR=331,NIN=7,                              
     C NPT=100,NFR=50,NTZ=50,NCMX=40,NMC=421,MCW=600,MVW=60,                    
     C MPW=300,MAT=2,MSR=2,MSUB=2,MREC=13,M2D=60,MIN=63,MCSTR=130,              
     C MVSTR=12,MCST=30,MVPT=100,MPT=400,MFR=100,MMZ=100,MFJ=12,                
     C MFN=6000,MFV=600,MEN=50,                                                 
     C NL1=(NLCH*NLCH+NLCH)/2,MEL=(MPW*MPW+MPW)/2,MCN=MPW*NLWF*NLWF,            
     C MGW=3*MPW*MPW+2*MREC*MREC,MDS=2*MEN*(MSR+2))                             
      REAL KM(2,MPW),K(2)                                                       
      REAL FJ(MFJ)                                                            
      COMPLEX GLA(MCN),GLB(MCN),YLM(MCN),GW(MGW),ENS(MEN),DP                    
      DIMENSION INCST(2),MAPW(3),MCPW(3),         
     C CORD(3,MAT+2*MSUB),MESH(MSR),RMT(MSR+2*MSUB),                   
     C G1(2),G2(2),IG(3,MPW),IGR(MPW),         
     C IW(MPW),UMT(8*NL*MSR),CSEMB1(MPW),CSEMB2(MPW),XV(MIN),                   
     C WV(MIN),IC(3*MCW),IC2(2*MVW),PW(MCW),                                    
     C PS(2*NVR*NLWF*MSR),ICSTR(MCW),LCSTR(8*MCSTR),IVSTR(MVW),                 
     C LVSTR(8*MVSTR),NLMST(MSR),LATST(8*MSR),IST(2*MAT),                       
     C LST(NL1*MSR),MST(NL1*MSR),LCST(2*NL1*MSR),ST(32*NLCH),                   
     C CPSTR(MCSTR),ROLL(NVR*MCST),PGI(2*MVSTR),                                
     C VGI(2*MVSTR),PINT(MCSTR),VINT(MCSTR),                                    
     C DV(4*MVSTR),PMUF(NVR*MCST),                                              
     C VMUF(NVR*MCST),QINT(MCST),QMT(MCST),Q(MCSTR),VSUB1(4*MVPT),              
     C VSUB2(4*MVPT),P(MPT),PT(MFR),BV(NIN),YV(NIN),NCOR(MSR),                  
     C LNST(NLCH*MSR),ICV(MCSTR),ALA(MCN),BLA(MCN),H(MEL),S(MEL),               
     C CHW(MEL),SW(MEL*MSR),SI(MEL),HW(MEL),PL(2*NL),ZFUN(2*MMZ*M2D),           
     C DS(MDS),FN(MFN),FV(MFV),NL2J(NCMX*MSR)                                   
      CALL INTWRP(ITER,NAPW,NV1,IG,NCPW,IC,ICSTR,VINT,NVPW,IC2,IVSTR,           
     C VGI,NAT,NST,NSUB1,NSUB2,CORD,IST,RMT,G1,G2,G3,AR,D,ZEMB1,                
     C ZEMB2,MAPW,MCPW,NIN,XV,WV,BV,YV,FV,FN,HW,SI)                             
      CALL CHDEN(NEN,EBEG1,EBEG2,EIF1,EIF2,DP,S,SW,SI,H,NAT,NST,NAPW,           
     C NV1,NREC,IGR,CSEMB1,CSEMB2,D,DB,GW,IW,XV,WV,CHW,DS,WT,NK,NKK,PL,         
     C YLM,ALA,BLA,GLA,GLB,NL,NLWF,NLCH,NL1,AR*D,G1,G2,G3,IC,IG,CPW,            
     C NCPW,CORD,RMT,K,KM,ZFUN,NMZ,NOP,ST,NLMST,LATST,                          
     C LST,MST,LCST,NVR,MESH,PS,ROLL,NCSTR,CPSTR,LCSTR,INCST)                   
      CALL IFACDS(NEN,DP,S,SW,H,NST,NAPW,NREC,IGR,CSEMB1,CSEMB2,D,              
     C GW(1),GW(NAPW*NAPW+1),GW(2*NAPW*NAPW+1),GW(3*NAPW*NAPW+1),               
     C GW(3*NAPW*NAPW+NREC*NREC+1),IW,DS,ENS,WT,NK.EQ.NKK,K,IPLOT)              
      end