! am_02/02: adding DT1DTS = DT0DTS/abs(sinIm)
C     ****
C     ****     /consdyn/ FOR CONSTANTS USED DURING THE COURSE OF THE
C     ****       DYNAMO CALCULATION
C     ****
      real :: re,h0,h00,hs,r0,r00,rs,r1,dtr,rtd,dlatg,dlong,dlatm,
     |  dlonm,ylatg,ylong,ylatm,ylonm,rcos0s,dt0dts,unit,table,dt1dts
      COMMON/consdyn/RE,H0,H00,HS,R0,R00,RS,R1,DTR,RTD,DLATG,DLONG,
     1  DLATM,DLONM,YLATG(0:JMAXGP),YLONG(IMAXGP),YLATM(JMAXM),
     2  YLONM(IMAXMP),RCOS0S(JMAXM),DT0DTS(JMAXM),UNIT(IMAXM),
     3  TABLE(91,2),DT1DTS(JMAXM)
C     ****
C     ****     WHERE:
C     ****       IMAXG=72, JMAXG=36, IMAXM=160, JMAXM=193,
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     ****       TABLE = TABLE GIVING VALUES THETA0 FOR EQUALLY SPACED
C     ****               VALUES OF THETASC     
C     ****       DT1DTS= DT0DTS/abs(sinIm) (is non-zero at equator)
C     ****
C