User: rhaas Date: 2012/06/20 11:19 AM

Modified: /trunk/src/ GRHydro_ParamCheck.F90

Log: GRHydro: do not access whichpsidcspeed in ParamCheck unless evolve_MHD

File Changes:

Directory: /trunk/src/ ======================

File [modified]: GRHydro_ParamCheck.F90 Delta lines: +5 -1 =================================================================== --- trunk/src/GRHydro_ParamCheck.F90 2012-06-20 16:18:04 UTC (rev 359) +++ trunk/src/GRHydro_ParamCheck.F90 2012-06-20 16:19:10 UTC (rev 360) @@ -126,7 +126,7 @@ end if

- if(clean_divergence.ne.0) then + if(evolve_MHD.ne.0 .and. clean_divergence.ne.0) then if (CCTK_EQUALS(psidcspeed,"char speed")) then whichpsidcspeed = 0 else if (CCTK_EQUALS(psidcspeed,"light speed")) then @@ -138,6 +138,10 @@ end if end if

+ if(clean_divergence.eq.0 .and. evolve_MHD.ne.0) then + call CCTK_WARN(1, "You activated divergence cleaning but do not evolve a magnetic field. At best, nothing should happen, be preapered for a segfault otherwise.") + end if + if (Tmunu_damping_radius_min .gt. Tmunu_damping_radius_max) then call CCTK_WARN(0, "Minimum damping radius is greater than maximum damping radius!") end if