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: 
  end subroutine calculate_heldsuarez
  
  subroutine set_tendencies(dvrtgdt,ddivgdt,dvirtempgdt, dqdt, dlnpsgdt,dtracergdt) bind(c, name='gfs_set_tendencies')
  
      real(c_double), intent(in), dimension(nlons,nlats,nlevs) :: &
          dvrtgdt,ddivgdt,dvirtempgdt, dqdt
      real(c_double), intent(in), dimension(nlons,nlats,nlevs,ntrac) :: &
          dtracergdt
      real(c_double), intent(in), dimension(nlons,nlats) :: dlnpsgdt
  
negative.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: '
  