C -----------------------------------------------------------------
#if defined (DOC)
C
C* COMDECK COMLEG: Description of Gaussian grid related parameters
C
C RMU : mu ----| GAUS(NJBEG:NJEND,1)
C RWT : weights of the gaussian | GAUS(NJBEG:NJEND,2)
C quadrature |
C RWOCS : weight/(1 - mu**2) | GAUS(NJBEG:NJEND,3)
C R1MU2 : 1.- MU**2 | GAUS(NJBEG:NJEND,4)
C RSQM2 : SQRT(1 - MU**2) | GAUS(NJBEG:NJEND,5)
C RCOLAT: colatitude | GAUS(NJBEG:NJEND,6)
C R1QM2 : 1./SQRT(R1MU2) | GAUS(NJBEG:NJEND,7)
C R1MUI : 1./(1 - MU**2) | GAUS(NJBEG:NJEND,8)
C R1MUA : 1./(A*(1 - MU**2)) | GAUS(NJBEG:NJEND,9)
C RLATI : latitude : arcsin(mu) ----| GAUS(NJBEG:NJEND,10)
C --------------------------------------------------------
C GAUS(NJBEG:NJEND,10): Global field containing all of the parameters
C --------------------------------------------------------
C
C Positioning within spectral arrays (triangular truncation only)
C ---------------------------------
C a) When dimensioned as SX(NLA)
C NIND(0:NTRUNC) : position of spectral component (m,m)
C NCLM(0:NTRUNC) : number of components for a given m
C
C b) When dimensioned as SX(NLARH)
C NINDRH(0:NTRUNC) : position of spectral component (m,m)
#endif
C
REAL*8 RMU(NJBEG:NJEND), RWT(NJBEG:NJEND), RWOCS(NJBEG:NJEND)
S ,R1MU2(NJBEG:NJEND),RSQM2(NJBEG:NJEND),RCOLAT(NJBEG:NJEND)
S ,R1QM2(NJBEG:NJEND),R1MUI(NJBEG:NJEND),R1MUA(NJBEG:NJEND)
S ,RLATI(NJBEG:NJEND)
C
REAL*8 GAUS(NJBEG:NJEND,10)
INTEGER NIND(0:NTRUNC), NINDRH(0:NTRUNC), NCLM(0:NTRUNC)
C
C POINTER STATEMENTS
POINTER (PTMU, RMU), (PTWT, RWT), (PTWOCS,RWOCS)
S , (PT1MU2,R1MU2), (PTSQM2,RSQM2), (PTCOLA,RCOLAT)
S , (PT1QM2,R1QM2), (PT1MUI,R1MUI), (PT1MUA, R1MUA)
S , (PTLATI,RLATI)
C
POINTER (PTGAUS,GAUS)
POINTER (PTIND,NIND),(PTINDRH,NINDRH), (PTCLM,NCLM)
C
C Common containing the pointers
C ------------------------------
COMMON /POMLEG/PTMU, PTWT, PTWOCS, PT1MU2, PTSQM2, PTCOLA
S ,PT1QM2, PT1MUI, PT1MUA, PTLATI, PTGAUS
S ,PTIND, PTINDRH, PTCLM