anchor.diff:                  ! $omp end parallel
               endif
           else
  -                domain%nsquared = options%lt_options%N_squared
  +                domain%nsquared = log(options%lt_options%N_squared)
           endif
   
           ! smooth array has it's own parallelization, so this probably can't go in a critical section
  
positive.code:         nz  = size(domain%u,2)
          nxu = size(domain%u,1)
          nyv = size(domain%v,3)
          
          ! Here we just set externalN based on two options.
          ! if NfromForcing is True and nvar is defined then N is to be
          ! read from the forcing data set instead of being calculated by ICAR,
  
negative.code:                         ! look up table is computed in log space
                          domain%nsquared(:,j,k) = log(domain%nsquared(:,j,k))
                      end do
  
                      if (smooth_nsq) then
                          do j=1,nz
                              ! compute window as above.
                              top = min(j+vsmooth,nz)
                              bottom = max(1, j - (vsmooth - (top-j)) )
  
                              do smoothz = bottom, j-1
                                  domain%nsquared(:,j,k) = domain%nsquared(:,j,k) + domain%nsquared(:,smoothz,k)
                              end do
  