Changeset 1144 for trunk/source/radius_next.cpp
- Timestamp:
- 05/21/07 20:25:16 (20 months ago)
- Files:
-
- 1 modified
-
trunk/source/radius_next.cpp (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/radius_next.cpp
r1111 r1144 567 567 568 568 /*>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>*/ 569 /* check whether change in wind velocity constrains DRAD */ 570 /*>>chng 04 may 25, 571 * WJH 22 May 2004: disable when we are near the sonic point since 569 /* check whether change in wind velocity constrains DRAD 570 * WJH 22 May 2004: disable when we are near the sonic point since 572 571 * the velocity may be jumping all over the place but we just want 573 572 * to push through it as quickly as we can */ … … 575 574 { 576 575 double v = fabs(wind.windv); 577 /* this is fractional accel over length */576 /* this is fractional acceleration over length */ 578 577 dVeldRad = fabs(wind.windv-OldWindVelocity)/ 579 578 MAX2(v,0.1*timesc.sound_speed_isothermal)/radius.drad; … … 1615 1614 { 1616 1615 fprintf( punch.ipDRout, 1617 "Wind, dVeldRad=%10.3e\n", 1618 dVeldRad ); 1616 "Wind, dVeldRad=%.3e windv=%.3e\n", 1617 dVeldRad , 1618 wind.windv); 1619 1619 } 1620 1620
