Change code in dynamics/nonhydro/solve.F to read as follows (thanks: Wes Jones)


C                                                                                24SEP99.358
C     INITIALIZE DIFFUSION COEFFS                                                24SEP99.359
C                                                                                24SEP99.360
cmic$ do all autoscope
cmic$2 private(i,j,k)
c$doacross
c$& local(i,j,k)
c$omp parallel do default(shared)
c$omp&private(i,j,k)
      DO J=1,MJX
        DO K=1,MKX
          DO I=1,MIX                                                             24SEP99.363
            CRSXK(I,J,K)=0.                                                      24SEP99.364
            DOTXK(I,J,K)=0.                                                      24SEP99.365
          ENDDO                                                                  24SEP99.366
        ENDDO                                                                    24SEP99.367
      ENDDO                                                                      24SEP99.368


Effect of this change is to parallelize this set of loops and to avoid 
redistributing data over caches.

-Rotang, Sept 27, 1999