Changeset 2087

Show
Ignore:
Timestamp:
05/19/08 16:41:03 (2 months ago)
Author:
gary
Message:

merge bugfix r2075 - do not include turbulent width when defining fine continuum resolution

Files:

Legend:

Unmodified
Added
Removed
Modified
Copied
Moved

  • branches/c08_branch/source/cont_createmesh.cpp

    r2034 r2087  
    457457         *  
    458458         * TeLowestFineOpacity of 1e4 K is temperature were line width is  
    459          * evaluationed.  Tests were done using the stop temperature in its place 
     459         * evaluated.  Tests were done using the stop temperature in its place 
    460460         * Te below 1e4 K made fine opacity grid huge  
    461461         * do not let temp get higher than 1e4 either - code run with stop temp 10 set 
    462462         * stop temp of 1e10K and assert thrown at line 204 of cont_createpointers.c  
    463463         * simply use 1e4 K as a characteristic temperature */ 
     464        /** \todo       1       set temp of 1e4K will be too coarse a line for PDRs where 
     465         * H2 line overlap is very important */ 
    464466        double TeLowestFineOpacity = 1e4; 
    465467        rfield.fine_opac_velocity_width =  
    466468                (realnum)sqrt(2.*BOLTZMANN/ATOMIC_MASS_UNIT*TeLowestFineOpacity/ 
    467                 /* TurbVel/5 will divide line into 5 cells when turbulence-dominated */ 
    468                 dense.AtomicWeight[rfield.fine_opac_nelem] + POW2(DoppVel.TurbVel/5.)) /  
     469                dense.AtomicWeight[rfield.fine_opac_nelem] ) /  
    469470                /* we want fine_opac_nresolv continuum elements across this line 
    470471                 * default is 1, changed with SET FINE CONTINUUM command */