

#ifdef ALLREAL
#ifdef ALLDOUBLE
      implicit real*8 (a-h,o-z)
#else
      implicit real (a-h,o-z)
#endif
#endif
      integer iorder, lz, mz, mlz, nz, kz,
     &     ntz, npz, jz,iz, jzntz, lHz, jHz,
     &        lz4, lz9, kz4, kz8, kz9, ireal
      parameter( iorder= 4 )
#ifdef ALLREAL
      parameter( lz=4000,mz=01,mlz=0,nz=32,kz=48, ntz=200,npz=48 )
      parameter( lHz= 200 )
      parameter( jz=093,iz= lHz+3,jzntz=jz+ntz+3,jHz=1 )
#else
#ifdef STRUC
      parameter( lz=103,mz=034,mlz=-10,nz=017,kz=022, ntz=130,npz=074 )
      parameter( jz=113,iz= lz+3,jzntz=jz+ntz+3)
      parameter( lHz= lz, jHz = 1 )
#else
      parameter( lz=10000,mz=01,mlz=0,nz=09,kz=2, ntz=001,npz=048 )
      parameter( lHz= 50 )
      parameter( jz=093,iz= lHz+3,jzntz=jz+ntz+3,jHz=1 )
#endif
#endif
      parameter( lz4= lz,kz4= kz )
      parameter( lz9= lz,kz9= kz )
#ifdef HALF
      parameter( kz8= 2*kz )
#else
      parameter( kz8= kz )
#endif
     
c
c     dif( iorder-4 ) 1,1,1
c     parameter( lz4= 1,kz4=1 )
c*      mlz must be < = 0
c*      nz-1, mz-1, mlz are the extrimum values allowed for n and m.

#ifdef ALLREAL
#ifdef ALLDOUBLE
#define DATATYPE2 real*8
#define DATATYPE real*8
#define conjg dble
#define aimag dble
#define cabs  abs
#define cexp  exp
#define atan  datan
#define atan2  datan2
#define amax1  dmax1
#define amin1  dmin1
#else
#define DATATYPE2 real
#define DATATYPE real
#define conjg real
#define aimag real
#define RCON real
#define cabs abs
#define cexp exp
#endif
      parameter( ireal = 1 )
#else
#define DATATYPE2 real
#define DATATYPE complex
      parameter( ireal = 0 )
#endif

