anchor.diff:        INTEGER, PARAMETER, PRIVATE:: ntb_ark = 4
         INTEGER, PARAMETER, PRIVATE:: ntb_IN = 55
         INTEGER, PRIVATE:: niIN2
  -      !$omp threadprivate(niIN2)
  +      ! $omp threadprivate(niIN2)
   
         DOUBLE PRECISION, DIMENSION(nbins+1):: xDx
         DOUBLE PRECISION, DIMENSION(nbc):: Dc, dtc
  
positive.code:       INTEGER, PARAMETER, PRIVATE:: ntb_i1 = 55
        INTEGER, PARAMETER, PRIVATE:: ntb_t = 9
        INTEGER, PRIVATE:: nic1, nic2, nii2, nii3, nir2, nir3, nis2, nig2, nig3
        !$omp threadprivate(nic1, nic2, nii2, nii3, nir2, nir3, nis2, nig2, nig3)
        INTEGER, PARAMETER, PRIVATE:: ntb_arc = 7
        INTEGER, PARAMETER, PRIVATE:: ntb_arw = 9
        INTEGER, PARAMETER, PRIVATE:: ntb_art = 7
  
negative.code: !..Schmidt number
        REAL, PARAMETER, PRIVATE:: Sc = 0.632
        REAL, PRIVATE:: Sc3
        ! $omp threadprivate(Sc3)
  
  !..Homogeneous freezing temperature
        REAL, PARAMETER, PRIVATE:: HGFR = 235.16
  