Changeset 24

Timestamp:

06/24/10 10:23:43 (2 years ago)

Author:

tabel

Message:

Some cleaning up. Comments in example script to go with documentation.

Location:

trunk

Files:

• TOOLS/construct_histogram.pro (modified) (1 diff)
• TOOLS/construct_interpolated_slice_data.pro (modified) (2 diffs)
• scripts/time_averaged_density_pdf.pro (modified) (4 diffs)

trunk/TOOLS/construct_histogram.pro

@@ -184 +184 @@
                  xstyle=1, ystyle=1, $
                  charsize=2, thick=2
+           meanx = obin1
+           meanh = mh
         endif else begin
            if ((dims eq 1) or (hi.cumoption eq 3)) then begin ; requested plot instead of histogram

trunk/TOOLS/construct_interpolated_slice_data.pro

@@ -190 +190 @@
               < (index_range[const_sub,*]-1)
 
-  ind = where((ii_s eq (index_range[const_sub,*]-1)) eq 1, count)
-  if count gt 0 then  ii_s[ind] -= 1
+  read_two = 0
+  read_two = (total(["div V", "density slope", "DM density slope"] eq list_str[var_index]) gt 0) < 1
+  print, "read_two", read_two
+  if read_two then begin
+     ind = where((ii_s eq (index_range[const_sub,*]-1)) eq 1, count)
+     if count gt 0 then  ii_s[ind] -= 1
+  endif
 
   ii_sp1 = (ii_s+1) ; plus one to have two indices to interpolate from
@@ -208 +213 @@
 start[*,other_subs[0]] = ii_l1
 start[*,other_subs[1]] = ii_l2
-read_two = 0
-read_two = (total(["div V", "density slope", "DM density slope"] eq list_str[var_index]) gt 0) < 1
-print, "read_two", read_two
 if read_two then righti[*,const_sub]   = ii_sp1[*] $
 else righti[*,const_sub]     = ii_s[*]

trunk/scripts/time_averaged_density_pdf.pro

@@ -2 +2 @@
 ; use the loaded snapshots to make a time average density probability
 ; distributions function (pdf)
+; at the end the variable mh holds the function and obin1 the
+; corresponding density values.
+
 pro time_averaged_density_pdf
 @TOOLS/common_blocks.inc
@@ -20 +23 @@
 ;  if time_index lt ntimes-1 then starttime = time_index
 
+; tells jacques to switch to the histogram tab
   J_event, 'HISTOGRAMS'
   hi.field1 = this_field
@@ -25 +29 @@
 
   for time_index=starttime,ntimes-1,astep DO BEGIN ; loop over all loaded timesteps
-     set_list_widget_values
-     wait, 1.Ma
-     select_current_grids 
-     check_cancel
+     set_list_widget_values                        ; update gui with the new time
+     wait, 1.
+     select_current_grids       ;
+     check_cancel               ; allow user to type  c in command line to stop 
      if stopit then begin
         stopit = 0
         break
      endif
+; tells jacques to draw the histogram
      J_event, 'DRAW_BUTTON'
 
@@ -40 +45 @@
      if verbose then print, 'memory use in Mb', memory()/1024./1e3
   endfor
-  time_index = ntimes-1
-  tapdf /= totaltime
+  time_index = ntimes-1         ; keep time_index in loaded range
+
+  tapdf /= totaltime            ; compute the average
+
   mh = tapdf                    ;  store in this common block variable
                                 ; to make it accesible from
-                                ; ccommandline
+                                ; commandline
 
+; show the plot of the final result
   titlestring1=this_field+' ['+get_unit(this_field,/ulabel)+']'
   plot, obin1, (mh) , xlog=shouldbe_logged(this_field), $