updating dustewrap_plot

Ilyes Choubani
1 parent 8bb63b77
Showing 1 changed file with 10 additions and 7 deletions Show diff stats
src/idl/dustemwrap_plot.pro
@@ -35,7 +35,7 @@ IF not keyword_set(st) THEN BEGIN
   st=dustem_run(p_dim)  
 ENDIF
  
-stop
+;stop
 ;If the interpolates arrays are not present (dustem_sed, etc...) run: 
 ;dustem_compute_sed etc to generate them. 
  
@@ -92,7 +92,10 @@ iswinext = !dustemcgwin_id.ext EQ la_undef()
 ;plotstrct = !p
 degtorad = !pi/180 ;factor that the arctan will be devided by in order to have an axis in degrees. Talk to JP maybe he wants angle in radians.
  
-dustem_psi_em = 0.5*atan(dustem_used,dustem_qsed)/degtorad
+;This should probably be moved to 'dustem_mpfit_run'. !!!!!!!
+;for now I am coding it like this. 
+
+if !run_pol then dustem_psi_em = 0.5*atan(dustem_used,dustem_qsed)/degtorad
  
 ;if iswinsed then begin ;I believe the iswin test should be below the case condition because I am not copying all the first 'HUGE LOOP' but rather replacing the positional arrays.            
 if test_m then begin
@@ -340,11 +343,11 @@ endif
  
  
 ;Putting everything required for the plotting of the legends (column density,model chi2 and rchi2) and the two windows for the parameters and the plugins respectively
-
-;     xxpos = 
-;     yypos =
-;     ltit  =  'Model: '+
-;     xyouts,xxpos,yypos,ltit+!dustem_model,color=0,/normal,charsize=legend_charsize
+    ;p_sed = [0.05,0.67,0.475,0.90]
+;     xxpos = 0.03
+;     yypos = 0.95
+;     ltit  =  'Model: '
+;     cgwindow,'xyouts',xxpos,/addcmd,windid=winid_m,yypos,ltit+!dustem_model,color=0,/normal,charsize=legend_charsize
 ;     xxpos =
 ;     yypos =
 ;     ltit = 'N_{H}:'+strmid(string(*!dustem_HCD,format='(1E10.2)'),2)+'(H/cm^2)'