anchor.diff:       radequiltemp(:,:,k) = (prs(:,:,k)/p0)**(2./7.)*&
            (315.-deltmp*sin(lats)**2-delthz*log(prs(:,:,k)/p0)*cos(lats)**2)
        enddo
  -     !$omp end parallel do 
  +     !$omp end parallel do
        blprof = (blprof-sigbot)/(1.-sigbot)
  -     where (blprof < 0) 
  +     where (blprof < 0)
            blprof = 0
        endwhere
        where (radequiltemp < tempstrat)
  
positive.code:      call grdtospec(forcingg(:,:,k), spectral_div_tend(:,k))
       !ddivspecdt(:,k) = forcingspec(:,k)
       enddo
       !$omp end parallel do 
       spectral_lnps_tend = 0.
  
       !print *, 'Maximum tendency in divergence: '
  
negative.code:          dtracerspecdt(:,k,trac) = spectral_tracer_tend(:,k,trac)
       enddo
       enddo
       !$omp end parallel do
       dlnpsspecdt(:) = spectral_lnps_tend(:)
  
  
  