C **** C **** /CONSTS/ FOR CONSTANTS USED DURING THE COURSE OF THE C **** DYNAMO CALCULATION C **** COMMON/ECONS/RE,R0,R00,PI,DTR,RTD,DLATM,DLONM,YLATM(JMAXM), 1 YLONM(IMAXMP),RCOS0S(JMAXM),DT0DTS(JMAXM),TABLE(181,2), 2 DT1DTS(JMAXM) C **** C **** WHERE: C **** IMAXG=72, JMAXG=36, IMAXM=641, JMAXM=385, 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 **** R0 = RE+H0 C **** R00 = RE+H00 C **** PI = 3.142------- C **** DTR = DEGREES TO RADIANS = PI/180. C **** RTD = RADIANS TO DEGREES = 180./PI C **** DLATM = PI/(JMAXM-1) C **** DLONM = 2.*PI/IMAXM 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 **** TABLE = TABLE GIVING VALUES THETA0 FOR EQUALLY SPACED C **** VALUES OF THETAS C **** DT1DTS = DT0DTS/abs(sinIm) (is non-zero at equator) C **** C