anchor.diff:        REAL, PARAMETER, PRIVATE:: D0s = 200.E-6
         REAL, PARAMETER, PRIVATE:: D0g = 250.E-6
         REAL, PRIVATE:: D0i, xm0s, xm0g
  +      !$omp threadprivate(D0i, xm0s, xm0g)
   
   !..Lookup table dimensions
         INTEGER, PARAMETER, PRIVATE:: nbins = 100
  
positive.code: !..Schmidt number
        REAL, PARAMETER, PRIVATE:: Sc = 0.632
        REAL, PRIVATE:: Sc3
  
  !..Homogeneous freezing temperature
        REAL, PARAMETER, PRIVATE:: HGFR = 235.16
  
negative.code: 
        enddo i_loop
        enddo j_loop
        !$omp end do
        !$omp end parallel
  
  ! DEBUG - GT
      !   write(mp_debug,'(a,7(a,e13.6,1x,a,i3,a,i3,a,i3,a,1x))') 'MP-GT:', &
  