Changes between Version 12 and Version 13 of 1DPulsedJets


Ignore:
Timestamp:
01/25/12 13:31:48 (13 years ago)
Author:
ehansen
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • 1DPulsedJets

    v12 v13  
    1919== Zcooling.tab ==
    2020PHcooling.tab contains a lot of data that is not needed.  Reading in this data one line at a time is inefficient, especially when astrobear is being run on multiple processors.  This is why a separate script called ph2zformat.s^*^ is used to reformat PHcooling.tab into Zcooling.tab before running astrobear.  Remember that PHcooling.tab still requires the aforementioned properties in order to be reformatted correctly.  Zcooling.tab has the following format:
    21 
    22 nDensities [[BR]]
    23 nTemps [[BR]]
    24 nXs [[BR]]
    25 lognemin [[BR]]
    26 lognemax [[BR]]
    27 tempmin [[BR]]
    28 tempmax [[BR]]
    29 logxmin [[BR]]
    30 logxmax [[BR]]
    31 coolingrate(1,1,1) [[BR]]
    32 coolingrate(2,1,1) [[BR]]
    33 ... [[BR]]
     21{{{
     22nDensities
     23nTemps
     24nXs
     25lognemin
     26lognemax
     27tempmin
     28tempmax
     29logxmin
     30logxmax
     31coolingrate(1,1,1)
     32coolingrate(2,1,1)
     33...
    3434coolingrate(nDensities,nTemps,nxs)
    35 
     35}}}
    3636Where nDensities, nTemps, and nXs are the number of unique values for log(ne), T, and log(X) respectively.  So if these parameters are equal to 40, 30, and 21, then there are a total of 25,200 cooling rates (40*30*21 = 25200).  lognemin is the minimum value of the log(ne) column in PHcooling.tab, and the other minimum and maximum values are defined in the same way. Now the data does not need to be read in one line at a time for 25,200 or so values.  The entire array of cooling rates can be read in with a single READ statement.  No iteration or nested DO loops required.
    3737
     
    6565
    6666Where rate1,rate2,...,rate8 are assigned as aforementioned, and the "in" values represent the coordinates of our point of interest inside the unit cube.
     67
     68[[BR]]
     69
     70== 1D Jet Simulations ==