Hello,neks<br><blockquote>In my simulation, fluid density 'Rho', dynamic viscosity 'Visco', thermal diffusivity 'Conduct' are all the function of the temperature, so I define them in USERVP as below:<br><br>    subroutine uservp (ix,iy,iz,ieg) <br>      include 'SIZE'<br>      include 'TOTAL'<br>      include 'NEKUSE'<br>      <br>      common /nob_coef/Kn(1:3),Bn(1:3),Cn(1:3),Dn(1:3),<br>     $                    En(1:3),Fn(1:3)<br>      common /constant/Rhom,Cpm,Tm,Alpham,Kappam,Km,<br>     $                    Viscom,delta_T,rho0<br>      real*8 Kn,Bn,Cn,Dn,En,Fn,<br>     $       Rhom,Cpm,Tm,Alpham,Kappam,Km,<br>     $         Viscom,delta_T,rho0      <br>      integer i,j,k,e,n,iel<br><br>      <br>      Rhom = 0.9922*1e3<br>      Cpm = 4.1690*1e3<br>      Alpham = 3.8810*1e-4<br>      Km = 0.6297<br>      Kappam = 0.1528*1e-6<br>      Viscom = 0.669*1e-6<br>      Tm = 40 <br>      <br>      Kn=[ -175.9*1e-4,295.8*1e-6,-460*1e-8]<br>      Bn=[23.52*1e-4,-14.9*1e-6,0.]<br>      Cn= [-3.736*1e3,-3.98*1e3,0.]<br>      Dn= [0.084*1e3,4.6*1e3,0.]<br>      En = [195*1e-4,-159.8*1e-4,207*1e-4]<br>      Fn=[21.99,-17.8,0.]<br><br>      iel=gllel(ieg)<br>      if (IFIELD.eq. 1)then<br>      udiff  = Rhom<br>      utrans = 0.<br>      rho0 = 0.<br>      do i=1,3<br>      utrans = utrans + Kn(i)*((temp-0.5)*delta_T)**i<br>      rho0 = rho0 + Cn(i)*((temp-0.5)*delta_T)**i<br>      enddo<br>      <br>      utrans=utrans*Viscom+Viscom<br>      rho0=rho0*Rhom+Rhom<br>      utrans=utrans*rho0<br>      endif<br>      if (IFIELD.eq. 2)then<br>      udiff  = Rhom*Cpm<br>      utrans = 0.<br>      do i=1,3<br>      utrans = utrans + Fn(i)*((temp-0.5)*delta_T)**i<br>      enddo  <br>      utrans=utrans*Km+Km<br>      endif<br>      <br>      return<br>      end<br></blockquote><br>  then I set p30=1 in .rea file. But the program cannot run, it seems that some errors occured when reading .rea file. I know from the user guide that variable properties are supported but I am little confused about if the properties which are the functions of other flow field variables(like u, v, w or temperature) also can be set correctly in USERVP. Did anyone meet problems like this before? Any help is appreciated.<br><br><div align="right">regards,<br>Xu<br></div><span>--<br>Xu Bolun,<br>University of Science and Technology of China,<br>Hefei,Anhui,China<br></span>