fixed typo in polyk2 for ieos=7

src/main/eos.F90

@@ -217,6 +217,7 @@ subroutine equationofstate(eos_type,ponrhoi,spsoundi,rhoi,xi,yi,zi,eni,tempi,gam
  case(7)
 !
 !-- z-dependent locally isothermal eos
+!
     r2 = xi**2 + yi**2
     cs2mid = polyk * r2**(-qfacdisc)
     cs2atm = polyk2 * r2**(-qfacdisc2)

src/setup/setup_disc.f90

@@ -627,7 +627,7 @@ subroutine equation_of_state(gamma)
                print "(/,a)",' setting ieos=7 for locally isothermal disc with stratification'
                call temp_to_HR(temp_mid0,H_R(onlydisc),R_ref(onlydisc),mcentral,cs)
                call temp_to_HR(temp_atm0,H_R_atm,R_ref(onlydisc),mcentral,cs)
-               polyk2 = (cs*(1./R_ref(onlydisc))**qfacdisc2)**2
+               polyk2 = (cs/R_ref(onlydisc)**(-qfacdisc2))**2
                z0 = z0_ref/R_ref(onlydisc)**beta_z
              else
                ieos = 3