C     ****
C     ****     /CONSTS/ FOR CONSTANTS USED DURING THE COURSE OF THE
C     ****       DYNAMO CALCULATION
C     ****
      COMMON/CONSTS/RE,H0,H00,HS,R0,R00,RS,R1,PI,DTR,RTD,DLATG,DLONG,
     1  DLATM,DLONM,YLATG(0:JMAXGP),YLONG(IMAXGP),YLATM(JMAXM),
     2  YLONM(IMAXMP),RCOS0S(JMAXM),DT0DTS(JMAXM),UNIT(IMAXM),DZ,
     3  TABLE(91,2)
C     ****
C     ****     WHERE:
C     ****       IMAXG=72, JMAXG=36, IMAXM=80, JMAXM=97,
C     ****         IMAXGP=IMAXG+1, JMAXGP=JMAXG+1, IMAXMP=IMAXM+1
C     ****       RE = EARTHS EQUATORIAL RADIUS = 6.378165E8 CM
C     ****       H0 = LOWER BOUND FOR QUADRATURES = 9.0E6
C     ****       H00 = LOWER BOUNDARY OF MODEL = 9.7E6
C     ****       HS = REFERENCE LEVEL FOR MAGNETIC FIELD =1.3E7 CM
C     ****       R0 = RE+H0
C     ****       R00 = RE+H00
C     ****       RS = RE+HS
C     ****       R1 = CONST. USED IN TRANSFORMATION FROM THATA0 TO
C     ****            THEATAS = 1.06E7
C     ****       PI = 3.142-------
C     ****       DTR = DEGREES TO RADIANS = PI/180.
C     ****       RTD = RADIANS TO DEGREES = 180./PI
C     ****       DLATG = PI/JMAXG
C     ****       DLONG = 2.*PI/IMAXG
C     ****       DLATM = PI/(JMAXM-1)
C     ****       DLONM = 2.*PI/IMAXM
C     ****       YLATG = ARRAY OF GEOGRAPHIC LATITUDE VALUES (RADIANS)
C     ****       YLONG = ARRAY OF GEOGRAPHIC LONGITUDE VALUES (RADIANS)
C     ****       YLATM = ARRAY OF GEOMAGNETIC LATITUDE VALUES EQUALLY
C     ****                SPACED BY 2 DEGREES IN THETA0 BUT EXPRESSED
C     ****                INTERMS OF THETAS.     (RADIANS)
C     ****       YLONM = ARRAY OF GEOMAGNETIC LONGITUDE VALUES (RADIANS)
C     ****       RCOS0S = COS(THETA0)/COS(THETAS)
C     ****       DT0DTS = D(THETA0)/D(THETAS)
C     ****       UNIT = VECTOR OF WHICH ALL ELEMENTS ARE UNITY
C     ****       DZ = SEPARATION OF MODEL PRRESSURE LEVELS
C     ****       TABLE = TABLE GIVING VALUES THETA0 FOR EQUALLY SPACED
C     ****               VALUES OF THETAS
C     ****
C