INTEGER mpiTypeXFaceThread_xy_r4(MAX_NO_THREADS)

      INTEGER mpiTypeXFaceThread_xy_r8(MAX_NO_THREADS)

      INTEGER mpiTypeYFaceThread_xy_r4(MAX_NO_THREADS)

      INTEGER mpiTypeYFaceThread_xy_r8(MAX_NO_THREADS)

      INTEGER mpiTypeXFaceThread_xyz_r4(MAX_NO_THREADS)

      INTEGER mpiTypeXFaceThread_xyz_r8(MAX_NO_THREADS)

      INTEGER mpiTypeYFaceThread_xyz_r4(MAX_NO_THREADS)

      INTEGER mpiTypeYFaceThread_xyz_r8(MAX_NO_THREADS)

      LOGICAL thError(MAX_NO_THREADS)

      LOGICAL threadIsRunning(MAX_NO_THREADS)

      LOGICAL threadIsComplete(MAX_NO_THREADS)

      LOGICAL allMyEdgesAreSharedMemory(MAX_NO_THREADS)

      Real*8  phiGMR8(lShare8, 0:MAX_NO_THREADS )

      Real*4  phiGMR4(lShare4, 0:MAX_NO_THREADS )

      INTEGER phiGMI (lShare4, 0:MAX_NO_THREADS )

      Real*8  phiGSR8 (lShare8, 0:MAX_NO_THREADS )

      Real*4  phiGSR4 (lShare4, 0:MAX_NO_THREADS )

      INTEGER phiGSI  (lShare4, 0:MAX_NO_THREADS )

      _RL hh(4*imax+1,4*imax,MAX_NO_THREADS)

      _RL c(4*imax,MAX_NO_THREADS),s(4*imax,MAX_NO_THREADS)

      _RL rs(4*imax+1,MAX_NO_THREADS),t,ro

      integer i1Thid(MAX_NO_THREADS)

      _RL r0, gam, eps1(MAX_NO_THREADS)

      INTEGER BAR2_level(cacheLineSize/4,MAX_NO_THREADS)

      INTEGER BAR2_barrierCount(cacheLineSize/4,MAX_NO_THREADS)

      INTEGER BAR2_spinsCount(cacheLineSize/4,MAX_NO_THREADS)

      INTEGER BAR2_spinsMax(cacheLineSize/4,MAX_NO_THREADS)

      INTEGER BAR2_spinsMin(cacheLineSize/4,MAX_NO_THREADS)

      DO I=1, MAX_NO_THREADS

      INTEGER key1(lShare4,MAX_NO_THREADS)

      INTEGER key2(lShare4,MAX_NO_THREADS)

      INTEGER key3(lShare4,MAX_NO_THREADS)

      INTEGER bCount(MAX_NO_THREADS)

      INTEGER masterSet(MAX_NO_THREADS)

      _RL fOrogr(NGP,MAX_NO_THREADS)

      _RL WVSurf(NGP,MAX_NO_THREADS)

      _RL snLat (NGP,MAX_NO_THREADS)

      _RL csLat (NGP,MAX_NO_THREADS)

      _RL field(sNx,sNy,nNr,MAX_NO_THREADS)

      _RL fmask1(NGP,3,MAX_NO_THREADS)

      _RL sst1   (NGP,MAX_NO_THREADS)

      _RL sti1   (NGP,MAX_NO_THREADS)

      _RL stl1   (NGP,MAX_NO_THREADS)

      _RL soilw1 (NGP,MAX_NO_THREADS)

      _RL alb1   (NGP,0:3,MAX_NO_THREADS)

      _RL dTsurf (NGP,3,MAX_NO_THREADS)

      _RL PRECNV (NGP,MAX_NO_THREADS)

      _RL PRECLS (NGP,MAX_NO_THREADS)

      _RL EnPrec (NGP,MAX_NO_THREADS)

      _RL CBMF   (NGP,MAX_NO_THREADS)

      _RL TSR    (NGP,MAX_NO_THREADS)

      _RL SSR    (NGP,0:3,MAX_NO_THREADS)

      _RL SLR    (NGP,0:3,MAX_NO_THREADS)

      _RL OLR    (NGP,MAX_NO_THREADS)

      _RL SHF    (NGP,0:3,MAX_NO_THREADS)

      _RL EVAP   (NGP,0:3,MAX_NO_THREADS)

      _RL SPEED0 (NGP,MAX_NO_THREADS)

      _RL DRAG   (NGP,0:3,MAX_NO_THREADS)

      _RL TT_SWclr(NGP,NLEV,MAX_NO_THREADS)

      _RL TSWclr  (NGP,MAX_NO_THREADS)

      _RL SSWclr  (NGP,MAX_NO_THREADS)

      _RL TT_LWclr(NGP,NLEV,MAX_NO_THREADS)

      _RL OLWclr  (NGP,MAX_NO_THREADS)

      _RL SLWclr  (NGP,MAX_NO_THREADS)

      _RL RH   (NGP,NLEV,MAX_NO_THREADS)

      _RL TS    (NGP,MAX_NO_THREADS)

      _RL TSKIN (NGP,MAX_NO_THREADS)

      _RL T0    (NGP,MAX_NO_THREADS)

      _RL Q0    (NGP,MAX_NO_THREADS)

      _RL CLTOP (NGP,MAX_NO_THREADS)

      _RL CLOUDC(NGP,MAX_NO_THREADS)

      _RL TT_CNV (NGP,NLEV,MAX_NO_THREADS)

      _RL QT_CNV (NGP,NLEV,MAX_NO_THREADS)

      _RL TT_LSC (NGP,NLEV,MAX_NO_THREADS)

      _RL QT_LSC (NGP,NLEV,MAX_NO_THREADS)

      _RL TT_RSW (NGP,NLEV,MAX_NO_THREADS)

      _RL TT_RLW (NGP,NLEV,MAX_NO_THREADS)

      _RL TT_PBL (NGP,NLEV,MAX_NO_THREADS)

      _RL QT_PBL (NGP,NLEV,MAX_NO_THREADS)

      Real*8 GSR8_value(lShare8,MAX_NO_THREADS)

      Real*8 GSR8_value(lShare8,MAX_NO_THREADS)

      DO I=1, MAX_NO_THREADS

      DO I = 2, MAX_NO_THREADS

      IF   ( nThreads .GT. MAX_NO_THREADS ) THEN

     &   nTx*nTy, MAX_NO_THREADS

      PARAMETER (size = nmaxevents*maxTimers*MAX_NO_THREADS)

      INTEGER*8 i_result(nmaxevents, maxTimers, MAX_NO_THREADS)

      INTEGER*8 i_result1(nmaxevents, maxTimers, MAX_NO_THREADS)

      INTEGER*8 i_result2(nmaxevents, maxTimers, MAX_NO_THREADS)

      REAL*8 fp_result(nmaxevents, maxTimers, MAX_NO_THREADS)

      INTEGER timerStarts( maxTimers , MAX_NO_THREADS)

      INTEGER timerStops ( maxTimers , MAX_NO_THREADS)

      Real*8 timerUser  ( maxTimers , MAX_NO_THREADS)

      Real*8 timerWall  ( maxTimers , MAX_NO_THREADS)

      Real*8 timerSys   ( maxTimers , MAX_NO_THREADS)

      Real*8 timerT0User( maxTimers , MAX_NO_THREADS)

      Real*8 timerT0Wall( maxTimers , MAX_NO_THREADS)

      Real*8 timerT0Sys ( maxTimers , MAX_NO_THREADS)

      INTEGER timerStatus( maxTimers , MAX_NO_THREADS)

      INTEGER timerNameLen( maxTimers , MAX_NO_THREADS)

      CHARACTER*(maxString) timerNames( maxTimers , MAX_NO_THREADS)

      INTEGER nTimers(MAX_NO_THREADS)

      INTEGER*8 values(nmaxevents, maxTimers , MAX_NO_THREADS),

     $     values1(nmaxevents, maxTimers, MAX_NO_THREADS),

     $     values2(nmaxevents, maxTimers, MAX_NO_THREADS)

      INTEGER*8 i_result(nmaxevents, maxTimers, MAX_NO_THREADS)

      INTEGER*8 i_result1(nmaxevents, maxTimers, MAX_NO_THREADS)

      INTEGER*8 i_result2(nmaxevents, maxTimers, MAX_NO_THREADS)

      REAL*8 fp_result(nmaxevents, maxTimers, MAX_NO_THREADS)

      DATA nTimers  /MAX_NO_THREADS*0/

      PARAMETER (size = 3*nmaxevents*maxTimers*MAX_NO_THREADS)

      INTEGER*8 values(nmaxevents, maxTimers , MAX_NO_THREADS),

     $     values1(nmaxevents, maxTimers, MAX_NO_THREADS),

     $     values2(nmaxevents, maxTimers, MAX_NO_THREADS)

      INTEGER ioErrorCount(MAX_NO_THREADS)

      INTEGER myBxLo(MAX_NO_THREADS)

      INTEGER myBxHi(MAX_NO_THREADS)

      INTEGER myByLo(MAX_NO_THREADS)

      INTEGER myByHi(MAX_NO_THREADS)

      INTEGER MAX_NO_THREADS

      PARAMETER ( MAX_NO_THREADS =  4 )