Changeset 1717 for branches/newmole/source

Timestamp:

12/23/07 06:39:42 (11 months ago)

Author:

rjrw

Message:

Code consensus view on the implementation of photodesorption rates.

Location:

branches/newmole/source

Files:

• mole.h (modified) (1 diff)
• mole_reactions.cpp (modified) (1 diff)
• mole_species.cpp (modified) (2 diffs)

branches/newmole/source/mole.h

@@ -107 +107 @@
         int num_total, num_calc, num_compacted;
         
-        realnum eden_f, grain_area, grain_avg_area;
+        realnum eden_f, grain_area;
 
         /** total charge in molecules */

branches/newmole/source/mole_reactions.cpp

@@ -997 +997 @@
          * (cf the d'Hendercourt reference in Bergin et al) */
 
-        fixit(); // Should this be 4* (i.e. the total not the cs area)?  Depends on definition of radiation field.
-        return mole.grain_avg_area * hmi.UV_Cont_rel2_Draine_DB96_depth *(1.232e7f * 1.71f);
+        fixit(); 
+        // 2e-15 is mean adsorbed species cross section, from
+  // >>refer Greenberg, L. T., 1973, IAUS, 52, 413
+        // 1.232e7f * 1.71f is from DB96 referred to below 
+        // will be multiplied by yield factor (~1e-4) to get overall photodesorption rate
+        return 2e-15 * hmi.UV_Cont_rel2_Draine_DB96_depth *(1.232e7f * 1.71f);
 }
 

branches/newmole/source/mole_species.cpp

@@ -596 +596 @@
         int i;
         int nd;
-        double den_times_area, dengrains;
+        double den_times_area;
 
         DEBUG_ENTRY("mole_update_species_cache()");
@@ -613 +613 @@
         if(gv.lgDustOn)
         {
-                den_times_area = dengrains = 0.0;
+                den_times_area = 0.0;
                 for( nd=0; nd < gv.nBin; nd++ )
                 {
                         /* >>chng 06 mar 04, update expression for projected grain surface area, PvH */
                         den_times_area += gv.bin[nd]->IntArea/4.*gv.bin[nd]->cnv_H_pCM3;
-                        dengrains += gv.bin[nd]->cnv_GR_pCM3;
                 }
                 
                 mole.grain_area = den_times_area;
-                mole.grain_avg_area = den_times_area/dengrains;
         }
         else
         {
-                mole.grain_area = mole.grain_avg_area = 0.0;
+                mole.grain_area = 0.0;
         }