[nils] r435 committed - Added the script for the MIE results

[ni...@googlecode.com]

Revision: 435
Author: Jen...@gmail.com
Date: Thu Aug 4 10:41:11 2011
Log: Added the script for the MIE results
http://code.google.com/p/nils/source/detail?r=435

Added:
/trunk/students/jenstarj/tex/figs/mie_rs_log_plot3.pro

=======================================
--- /dev/null
+++ /trunk/students/jenstarj/tex/figs/mie_rs_log_plot3.pro Thu Aug 4
10:41:11 2011
@@ -0,0 +1,106 @@
+Size=2
+if !d.name eq 'PS' then begin
+ device,xsize=30,ysize=40,yoffset=3
+ !p.charthick=4 & !p.thick=4 & !x.thick=4 & !y.thick=4
+ !p.charsize=size
+endif else begin
+ !p.charsize=size
+ device,decompose=0
+ device,retain=2
+ loadct,5
+ Window, XSize=850, YSize=1000
+ !p.charsize=size
+end
+
+!p.multi=[0,1,1]
+
+rs1d= [0.01, 0.015, 0.02, 0.03, 0.04, 0.05, 0.1]
+rs3d = [0.015, 0.02, 0.03, 0.04, 0.05, 0.1]
+rsmaas2D=[0.003, 0.005, 0.01, 0.02, 0.05, 0.1]
+rsmaas3D=[0.005, 0.01, 0.02, 0.05, 0.1]
+rshong=[0.01,0.05,0.1]
+
+;#####The pencil code###############
+mie1d=[9.791, 13.23, 16.51, 22.90, 29.12, 35.25, 65.55]
+mie2d=[0.2985, 0.558, 0.889, 1.753, 2.880, 4.266, 15.03]
+mie3d=[0.0229, 0.0457, 0.127, 0.271, 0.493, 3.301]
+
+;######Maas#########################
+
+mie2dmaas=[0.05, 0.08, 0.205, 0.7, 3.3, 12]
+mie3dmaas=[0.001, 0.004, 0.02, 0.24, 2]
+
+;######Hong#########################
+
+mie3dHongAir=[0.01, 0.55, 7]
+
+;######The plots###############
+
+plot,rs1d,mie1d,ytitle='!8E!d!6min!n!6',xtitle='!8r!d!6s0!n
[cm]!6',ps=-2,xrange=[0.009,0.11],xstyle=1,yrange=[0.008,100], ystyle=1 ,
/xlog, /ylog, Li=0
+
+oplot, rs1d, mie2d, Li=0, ps=-4
+oplot, rs3d, mie3d, Li=0, ps=-5
+
+oplot, rsmaas3D, mie3dmaas, Li=0, ps=-5, col = 50
+oplot, rsmaas2D, mie2dmaas, Li=0, ps=-4, col = 50
+
+oplot, rshong, mie3dHongAir, Li=0, ps=-5, col = 100
+
+;########## REG ############### commeted out
+
+;nx=1000
+;v = double(linspace(0.01, 0.1, nx))
+
+;reg1d = 422.422*v^0.825
+;reg2d = 713.197*v^1.702
+;reg3d = 1310.448*v^2.623
+
+;oplot, v, reg1d, Li=0, col=180
+;oplot, v, reg2d, Li=2, col=180
+;oplot, v, reg3d, Li=3, col=180
+
+;###########LEGEND############
+
+
+x=0.0035
+dx=0.002
+y= 10
+
+xx=0.03
+dxx=0.013
+yy= 3
+
+xxx=0.0035
+dxxx=0.002
+yyy=0.05
+
+xxxx = 0.0059
+dxxxx = 0
+
+legend,xx,dxx,yy^(-3.5),0,'!6Pencil Code (H!d2!n/air)!6'
+legend,xx,dxx,yy^(-3.8),0,'!6Maas et al. (H!d2!n/O!d2!n)!6', col = 50
+legend,xx,dxx,yy^(-4.1),0,'!6Hong et al. (H!d2!n/air)!6', col = 100
+
+;legend,xxx,dxxx,yyy^1.0,0,'!6Curve fitting for', col=180
+;legend,xxxx,dxxxx,yyy^1.1,0,'!6the Pencil Code', col=180
+
+;xyouts,0.04,0.02,'!8E!6!dmin!n(Glassmann)!6'
+
+xyouts,0.01,17,'!61D: [mJ/cm!u2!n]!6'
+xyouts,0.01,0.8,'!62D: [mJ/cm]!6'
+xyouts,0.01,0.03,'!63D: [mJ]!6'
+
+xyouts,0.024,23,'!8Ar!u!60.83'
+xyouts,0.024,1.7,'!8Br!u!61.70'
+xyouts,0.024,0.15,'!8Cr!u!62.62'
+
+;######(Glassmann)############
+xmin=[0.01, 0.10]
+min=[0.018, 0.018]
+;oplot, xmin, min
+
+
+
+print,'mv idl.ps
~/numerics/pencil-code/nils/students/jenstarj/tex/figs/mie_rs_log_plot3.ps'
+
+END