Changes between Version 41 and Version 42 of ModulesOnAstroBear
- Timestamp:
- 01/24/13 10:16:26 (12 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
ModulesOnAstroBear
v41 v42 287 287 Info%q(i,j,k,irho) = 1d0 288 288 Info%q(i,j,k,ivx)=10d0 289 IF (iE /= 0) Info%q(i,j,k, iE)=gamma7+50d0 289 290 IF (ivy /= 0) Info%q(1:Info%mX(1),1:Info%mX(2),1:Info%mX(3), ivy)=0d0 290 291 IF (ivz /= 0) Info%q(1:Info%mX(1),1:Info%mX(2),1:Info%mX(3), ivz)=0d0 291 IF (iE /= 0) Info%q(1:Info%mX(1),1:Info%mX(2),1:Info%mX(3), iE)=gamma7+50d0292 292 END IF 293 293 END DO … … 367 367 Info%q(i,j,k,irho) = 1d0 368 368 Info%q(i,j,k,ivx)=velocity 369 IF (iE /= 0) Info%q(i,j,k, iE)=gamma7+half*velocity**2 369 370 IF (ivy /= 0) Info%q(1:Info%mX(1),1:Info%mX(2),1:Info%mX(3), ivy)=0d0 370 371 IF (ivz /= 0) Info%q(1:Info%mX(1),1:Info%mX(2),1:Info%mX(3), ivz)=0d0 371 IF (iE /= 0) Info%q(1:Info%mX(1),1:Info%mX(2),1:Info%mX(3), iE)=gamma7+half*velocity**2372 372 END IF 373 373 END DO … … 466 466 If we want to get more complicated - we can modify other clump/wind/ambient attributes. All of the available (and default) options for the various objects should be documented on the AstroBearObjects page. 467 467 468 To try the problem module we just built, we can copy the data files from {{{$(ASTROBEAR)/modules/Template}}} directory to the current problem directory {{{$(ASTROBEAR)/modules/Problem}}} and add these three lines to problem.data 469 {{{ 470 &ProblemData 471 rho = 10.0 472 radius = 1.0 473 velocity = 10.0 474 / 475 }}} 476 477 If we follow the procedure [BuildingAstroBear#SettingupandCompilingaProblem "Setting up and Compiling A Problem"] to try to compile and run our new problem module. If you set the dimensions and final time in global.data as 478 {{{ 479 nDim = 2 ! number of dimensions for this problem (1-3) 480 GmX = 30,30,1 ! Base grid resolution [x,y,z] 481 MaxLevel = 0 ! Maximum level for this simulation (0 is fixed grid) 482 }}} 483 484 and 485 {{{ 486 final_time = 15d-1 ! The final time in computational units. 487 final_frame = 20 ! The final frame [10] 488 }}} 489 490 and follow [BuildingAstroBear#RunningAstroBear "First Run"] to run it. Following Chapter 2 to use VisIT to analyze the results, we can get a movie like this [http://www.pas.rochester.edu/~bliu/ProblemModule/ClumpOBJ.mpeg ClumpMovie][[BR]] 468 491 469 492