User: rhaas Date: 2012/07/05 12:36 AM
Modified: /trunk/src/ GRHydro_Prim2Con.F90
Log: * add limit on values of Y_e also in Prim2Con
File Changes:
Directory: /trunk/src/ ======================
File [modified]: GRHydro_Prim2Con.F90 Delta lines: +2 -0 =================================================================== --- trunk/src/GRHydro_Prim2Con.F90 2012-06-26 23:48:01 UTC (rev 374) +++ trunk/src/GRHydro_Prim2Con.F90 2012-07-05 05:36:29 UTC (rev 375) @@ -270,6 +270,8 @@ *dvelz*dvelz + 2*gxy*dvelx*dvely + 2*gxz*dvelx *dvelz + 2*gyz& *dvely*dvelz))
+ ye = max(min(ye,GRHydro_Y_e_max),GRHydro_Y_e_min) + call EOS_Omni_press(handle,keytemp,GRHydro_eos_rf_prec,n,& drho,deps,temp,ye,dpress,keyerr,anyerr) ! error handling