Show
Ignore:
Timestamp:
05/21/07 20:25:16 (20 months ago)
Author:
gary
Message:

most files, remove dead code, correct comments

init_coreload.cpp - start of reorganizing punch commands

Files:
1 modified

Legend:

Unmodified
Added
Removed

  • trunk/source/pressure_total.cpp

    r1113 r1144  
    172172        dense.xMassDensity = (float)(dense.wmole*ATOMIC_MASS_UNIT*dense.pden); 
    173173         
    174         /*>>chng 04 may 25, introduce folling test on xMassDensity0  
     174        /*>>chng 04 may 25, introduce following test on xMassDensity0  
    175175         * WJH 21 may 04, this is the mass density that corresponds to the hden  
    176176         * specified in the init file. It is used to calculate the mass flux in the  
    … …  
    282282                                                TotalInsanity(); 
    283283                                        /* >>chng 01 jan 08, do not include highest level since maser can occur, 
    284                                         * and pops are set to small number in this case */ 
     284                                         * and pops are set to small number in this case */ 
    285285                                        /* >>chng 06 aug 17, should go to numLevels_local instead of _max. */ 
    286286                                        for( ipLo=0; ipLo < (iso.numLevels_local[ipISO][nelem] - 2); ipLo++ )  
    … …  
    330330                                                        /* >>chng 04 may 13, add Aul > 10 to prevent forbidden lines with 
    331331                                                         * tiny optical depths - taumin doubled - occur, but still opt thin  
    332                                                          * sympton was very large radiation pressure from Paschen beta */ 
     332                                                         * symptom was very large radiation pressure from Paschen beta */ 
    333333                                                        /* >>chng 04 aug 28, add test on ipCont being positive to skip bogus lines*/ 
    334                                                         /*NB - this code must be kep parrallel with code in prt_linepres */ 
     334                                                        /*NB - this code must be keep parallel with code in prt_linepres */ 
    335335                                                        if( Transitions[ipISO][nelem][ipHi][ipLo].Emis.ipCont > 0 && 
    336336                                                                Transitions[ipISO][nelem][ipHi][ipLo].Emis.Aul > 10. && 
    337337                                                            Transitions[ipISO][nelem][ipHi][ipLo].Hi.Pop > smallfloat && 
    338338                                                                /* >>chng 04 may 21, change PopOpc to poplo as defined above, 
    339                                                                  * to not overcount 23P of He */ 
     339                                                                 * to not over count 23P of He */ 
    340340                                                                /*Transitions[ipISO][nelem][ipHi][ipLo].PopOpc > smallfloat &&*/ 
    341341                                                                poplo > smallfloat && 
    … …  
    351351                                                                        StatesElem[ipISO][nelem][ipHi].g)/ 
    352352                                                                        /* >>chng 04 may 21, change PopOpc to poplo as defined above, 
    353                                                                          * to not overcount 23P of He */ 
     353                                                                         * to not over count 23P of He */ 
    354354                                                                        /*(Transitions[ipISO][nelem][ipHi][ipLo].PopOpc/ */ 
    355355                                                                        (poplo/ 
    … …  
    709709        { 
    710710                /* this branch is time dependent AND constant density */ 
    711                 /*fprintf(ioQQQ,"DEBUG enthalph HIT1\n");*/ 
     711                /*fprintf(ioQQQ,"DEBUG enthalpy HIT1\n");*/ 
    712712                /* this is the time-varying case where density is constant */ 
    713                 /* \todo        1       this has 3/2 on the PresGasCurr while trule dynamics case below 
     713                /* \todo        1       this has 3/2 on the PresGasCurr while true dynamics case below 
    714714                 * has 5/2 - so this is not really the enthalpy density - better 
    715715                 * would be to always use this term and add the extra PresGasCurr 
    … …  
    723723        { 
    724724                /* this branch either advective or constant pressure */ 
    725                 /*fprintf(ioQQQ,"DEBUG enthalph HIT2\n");*/ 
     725                /*fprintf(ioQQQ,"DEBUG enthalpy HIT2\n");*/ 
    726726                /* this is usual dynamics case, or time-varying case where pressure 
    727727                 * is kept constant and PdV work is added to the cell */ 
    … …  
    740740                /* >>chng 05 feb 03, do not update RadBetaMax in this block, wait until 
    741741                 * following, which also protects against the evaluation before zone 1, 
    742                  * when soln may not be stable 
     742                 * when solution may not be stable 
    743743                pressure.RadBetaMax = (float)MAX2(pressure.RadBetaMax,pressure.pbeta); */ 
    744744                /* stop RadPres from exceeding the gas pressure on first iteration */ 
    … …  
    767767                                /*@+redef@*/ 
    768768                                /*lint -e644 Piso_seq not initialized */ 
    769                                 /* trace verious constituents in radiation pressure */ 
     769                                /* trace serious constituents in radiation pressure */ 
    770770                                if( DEBUG_LOC  ) 
    771771                                { 
    … …  
    823823                else if( ipSet == 8 ) 
    824824                { 
    825                         /* hyperfile struct lines */ 
     825                        /* hyperfine struct lines */ 
    826826                        strcpy( pressure.chLineRadPres, "hyperf "); 
    827827                        ASSERT( ip1>=0  ); 
    … …  
    875875 
    876876        { 
    877                 /*@-redef@*/ 
    878877                enum{DEBUG_LOC=false}; 
    879                 /*@+redef@*/ 
    880878                if( DEBUG_LOC && pressure.PresTotlCurr > SMALLFLOAT /*&& iteration > 1*/ ) 
    881879                {