do_eqstate - equation of state

INTERFACE:

    subroutine do_eqstate()

DESCRIPTION:

Here, the equation of state is calculated for every 3D grid point. USES:

    use domain, only: imin,imax,jmin,jmax,kmax,az
    use variables_3d, only: kmin,T,S,rho,buoy,hn,alpha,beta
    use getm_timers, only: tic, toc, TIM_EQSTATE
  $ use omp_lib
    IMPLICIT NONE

LOCAL VARIABLES:

    integer                   :: i,j,k
    REALTYPE                  :: x
    REALTYPE                  :: p1,s1,t1
    REALTYPE                  :: th,densp