!
! am_02/02: be3 field added
!
! This is included in transf.f and intgrls.f for the dynamo:
! Params.h must also be included for dimension parameters.
!
      integer,parameter :: kmax=zkmx, kmaxp=kmax+1
      real :: ssigma1,ssigma2,zz,ww,adotv,axv,vxb,bxa,adota,adotam,
     |  a1dta2,sini,be3,sigma1m,sigma2m,zm,wm,adotvm,axvm,vxbm,bxam,
     |  a1a2m,siniam,bmodm,azm,pm
C     ****
C     ****     FIELDS TO BE TRANSFORMED TO GEOMAGNETIC SPACE
C     ****
! Global common:
      COMMON /transmag/ SSIGMA1(IMAXGP,0:JMAXGP,-2:KMAX),
     |  SSIGMA2(IMAXGP,0:JMAXGP,-2:KMAX),ZZ(IMAXGP,0:JMAXGP,-2:KMAXP),
     |  WW(IMAXGP,0:JMAXGP,-2:KMAX),ADOTV(IMAXGP,0:JMAXGP,-2:KMAX,2),
     |  AXV(IMAXGP,0:JMAXGP,-2:KMAX,2),VXB(IMAXGP,0:JMAXGP,-2:KMAX),
     |  BXA(IMAXGP,0:JMAXGP,2),ADOTA(IMAXGP,0:JMAXGP,2),
     |  A1DTA2(IMAXGP,0:JMAXGP),SINI(IMAXGP,0:JMAXGP),
! cga am_1200
     |  BE3(IMAXGP,0:JMAXGP)
! cge am_1200
C     ****
C     ****     THESE SAME QUANTITIES IN GEOMAGNETIC SPACE
C     ****       (ONE LATITUDE AT A TIME)
!
! Thread-private:
      common /transmag_priv/
     |  SIGMA1M(IMAXMP,-2:KMAX),SIGMA2M(IMAXMP,-2:KMAX),
     |  ZM(IMAXMP,-2:KMAXP),WM(IMAXMP,-2:KMAX),
     |  ADOTVM(IMAXMP,-2:KMAX,2),AXVM(IMAXMP,-2:KMAX,2),
     |  VXBM(IMAXMP,-2:KMAX),BXAM(IMAXMP,2),ADOTAM(IMAXMP,2),
     |  A1A2M(IMAXMP),SINIAM(IMAXMP),BMODM(IMAXMP),AZM(IMAXMP,2),
     |  PM(IMAXMP)
! cga am_1200
!$OMP THREADPRIVATE (/transmag_priv/)
!DIR$ TASKCOMMON transmag_priv
! cge am_1200