#include "dims.h"
      SUBROUTINE ADHEEL
      implicit none
C     ****
C     ****     MODIFY STENCILS AND RHS SO AS TO INSERT HEELIS
C     ****       POTENTIAL FOR GEOMAGNETIC LATITUDE .GE. 60.0 DEG.
C     ****
C     ****
#include "params.h"
#include "ceee.h"
#include "cterp.h"
#include "dynphi.h"
#include "consdyn.h"
#include "ioncr.h"
!
! Locals:
      integer :: ncc,imx,jmx,n
      
      crit(1) = 0.261799387
      crit(2) = 0.523598775
C     ****
C     ****     FILL ARRAY COLATC WITH NORTHERN AURORAL COLATITUDES
C     ****       CORRESPONDING TO EACH NORTHERN GEOMAGNETIC GRID
C     ****       POINT.  VALUES IN COLATC(IMX0,JMX0).
C     ****       ALSO CALCULATE FRACTIONAL PRESENCE OF DYNAMO
C     ****       EQUATION GIVEN CRITICAL COLATITUDES CRIT(2).
C     ****       VALUES IN P(IMX0,JMX0).
C     ****
      CALL COLATH
C     ****
C     ****     CALCULATE PHIHM, THE HEELIS POTENTIAL IN
C     ****       GEOMAGNETIC COORDINATES.
C     ****
      istar = 0
      CALL POTM
C     ****
C     ****     MODIFY STENCILS
C     ****
      NCC = 1
      IMX = IMX0
      JMX = JMX0
      DO 5 N = 1,5
        CALL STENMOD(IMX,JMX,CEE(NCC),PHIHM(1,JMX0),P,COLATC)
        NCC = NCC+9*IMX*JMX
        IF(N.EQ.1)NCC = NCC+IMX*JMX
        IMX = (IMX+1)/2
        JMX = (JMX+1)/2
    5 CONTINUE
      RETURN
      END
C
C--------------------------------------------------------------------------
#include "dims.h"
      SUBROUTINE ADHEELMD
!      
! am_02/02:   adheelmd: adheel for modified mudpack solver modmud
!             (isolve = 2 in dyncal)
!             subroutine stenmd: calculate modified and unmodified
!             coefficient stencil
!
      implicit none
C     ****
C     ****     MODIFY STENCILS AND RHS SO AS TO INSERT HEELIS
C     ****       POTENTIAL FOR GEOMAGNETIC LATITUDE .GE. 60.0 DEG.
C     ****
C     ****
#include "params.h"
#include "ceee.h"
#include "cterp.h"
#include "dynphi.h"
#include "consdyn.h"
#include "ioncr.h"
!
! Locals:
      integer :: ncc,imx,jmx,n
!
      crit(1) = 0.261799387
      crit(2) = 0.523598775
!
C     ****
C     ****     FILL ARRAY COLATC WITH NORTHERN AURORAL COLATITUDES
C     ****       CORRESPONDING TO EACH NORTHERN GEOMAGNETIC GRID
C     ****       POINT.  VALUES IN COLATC(IMX0,JMX0).
C     ****       ALSO CALCULATE FRACTIONAL PRESENCE OF DYNAMO
C     ****       EQUATION GIVEN CRITICAL COLATITUDES CRIT(2).
C     ****       VALUES IN P(IMX0,JMX0).
C     ****
      CALL COLATH
C     ****
C     ****     CALCULATE PHIHM, THE HEELIS POTENTIAL IN
C     ****       GEOMAGNETIC COORDINATES.
C     ****
      istar = 0
      CALL POTM
C     ****
C     ****     MODIFY STENCILS
C     ****
      NCC = 1
      IMX = IMX0
      JMX = JMX0
      DO 5 N = 1,5
       CALL STENMD(IMX,JMX,CEE(NCC),PHIHM(1,JMX0),P,COLATC)
	NCC = NCC+9*IMX*JMX
	IF(N.EQ.1)NCC = NCC+IMX*JMX
	IMX = (IMX+1)/2
	JMX = (JMX+1)/2
    5 CONTINUE
      RETURN
      END
C