sample.mco

蒙特卡罗算法。用盟特卡罗算法模拟光在组织钟的传播

A1 	# Version number of the file format.

####
# Data categories include: 
# InParm, RAT, 
# A_l, A_z, Rd_r, Rd_a, Tt_r, Tt_a, 
# A_rz, Rd_ra, Tt_ra 
####

# User time:  1488.13 sec =     0.41 hr.   Simulation time of this run.

InParm 			# Input parameters. cm is used.
sample.mco 	A		# output file name, ASCII.
1000000 			# No. of photons
0.01	0.01		# dz, dr [cm]
40	50	1	# No. of dz, dr, da.

3					# Number of layers
#n	mua	mus	g	d	# One line for each layer
1					# n for medium above
1.37	1	100	0.9	0.1	# layer 1
1.37	1	10	0	0.1	# layer 2
1.37	2	10	0.7	0.2	# layer 3
1					# n for medium below

RAT #Reflectance, absorption, transmission. 
0.0243729      	#Specular reflectance [-]
0.237544       	#Diffuse reflectance [-]
0.641594       	#Absorbed fraction [-]
0.0964896      	#Transmittance [-]

A_l #Absorption as a function of layer. [-]
      0.2612
      0.1488
      0.2316

A_z #A[0], [1],..A[nz-1]. [1/cm]
  2.7990E+00
  2.8099E+00
  2.8125E+00
  2.7845E+00
  2.7265E+00
  2.6497E+00
  2.5552E+00
  2.4446E+00
  2.3323E+00
  2.2037E+00
  2.0549E+00
  1.9209E+00
  1.7779E+00
  1.6612E+00
  1.5438E+00
  1.4102E+00
  1.2974E+00
  1.1851E+00
  1.0694E+00
  9.5971E-01
  1.7314E+00
  1.6297E+00
  1.5446E+00
  1.4724E+00
  1.4007E+00
  1.3238E+00
  1.2726E+00
  1.2104E+00
  1.1565E+00
  1.1089E+00
  1.0698E+00
  1.0307E+00
  9.9160E-01
  9.5890E-01
  9.3477E-01
  9.0365E-01
  8.7446E-01
  8.6162E-01
  8.4770E-01
  8.3724E-01

Rd_r #Rd[0], [1],..Rd[nr-1]. [1/cm2]
  3.0924E+01
  1.0071E+01
  6.1187E+00
  4.4015E+00
  3.5138E+00
  2.8528E+00
  2.4344E+00
  2.0692E+00
  1.7715E+00
  1.5409E+00
  1.3425E+00
  1.2006E+00
  1.0319E+00
  8.9424E-01
  7.9300E-01
  7.0183E-01
  6.1493E-01
  5.4537E-01
  4.8518E-01
  4.2576E-01
  3.7759E-01
  3.2900E-01
  2.9997E-01
  2.6421E-01
  2.3595E-01
  2.0775E-01
  1.8657E-01
  1.6498E-01
  1.4860E-01
  1.4064E-01
  1.2364E-01
  1.0955E-01
  9.9937E-02
  8.9281E-02
  8.3917E-02
  7.1952E-02
  6.5441E-02
  5.8384E-02
  5.3626E-02
  4.8753E-02
  4.5503E-02
  3.8753E-02
  3.6688E-02
  3.1858E-02
  3.1348E-02
  2.6778E-02
  2.5402E-02
  2.2791E-02
  2.0812E-02
  3.1424E-01

Rd_a #Rd[0], [1],..Rd[na-1]. [sr-1]
  3.4038E-02

Tt_r #Tt[0], [1],..Tt[nr-1]. [1/cm2]
  5.8942E-01
  6.0085E-01
  5.6453E-01
  5.5132E-01
  5.3655E-01
  5.3603E-01
  5.2027E-01
  4.9193E-01
  4.6582E-01
  4.3591E-01
  4.1109E-01
  3.9854E-01
  3.6636E-01
  3.3945E-01
  3.2122E-01
  2.9694E-01
  2.8041E-01
  2.6182E-01
  2.3866E-01
  2.2141E-01
  2.1004E-01
  1.9047E-01
  1.7161E-01
  1.6091E-01
  1.5083E-01
  1.3239E-01
  1.2549E-01
  1.1543E-01
  9.9789E-02
  9.4186E-02
  8.4806E-02
  7.8703E-02
  7.4215E-02
  6.6476E-02
  6.0723E-02
  5.6045E-02
  5.1798E-02
  4.6792E-02
  4.3551E-02
  4.0498E-02
  3.6332E-02
  3.4432E-02
  3.0875E-02
  2.9627E-02
  2.6649E-02
  2.3876E-02
  2.2584E-02
  2.1318E-02
  1.8962E-02
  3.3392E-01

Tt_a #Tt[0], [1],..Tt[na-1]. [sr-1]
  1.3826E-02

# A[r][z]. [1/cm3]
# A[0][0], [0][1],..[0][nz-1]
# A[1][0], [1][1],..[1][nz-1]
# ...
# A[nr-1][0], [nr-1][1],..[nr-1][nz-1]
A_rz
  3.2948E+03   3.0386E+03   2.5033E+03   1.8514E+03   1.2758E+03 
  8.5363E+02   5.5268E+02   3.6460E+02   2.4357E+02   1.7068E+02 
  1.1574E+02   8.3167E+01   6.3684E+01   4.9714E+01   3.7366E+01 
  3.1753E+01   2.3568E+01   2.0270E+01   1.5639E+01   1.3584E+01 
  1.9635E+01   1.7311E+01   1.4253E+01   1.2972E+01   9.3743E+00 
  9.7259E+00   8.8083E+00   7.1041E+00   6.9629E+00   5.8110E+00 
  4.5151E+00   5.3302E+00   3.9514E+00   3.7088E+00   3.5495E+00 
  3.2786E+00   3.3136E+00   3.2377E+00   2.6409E+00   2.3877E+00 
  7.7127E+01   1.3435E+02   2.1857E+02   2.8337E+02   3.0338E+02 
  2.8525E+02   2.4567E+02   1.9975E+02   1.5713E+02   1.2392E+02 
  9.4105E+01   7.2111E+01   5.4411E+01   4.3708E+01   3.5286E+01 
  2.7410E+01   2.2960E+01   1.8669E+01   1.5932E+01   1.1772E+01 
  1.9931E+01   1.6372E+01   1.3586E+01   1.2065E+01   1.1066E+01 
  9.0454E+00   8.0976E+00   7.3499E+00   6.8339E+00   5.2720E+00 
  5.5565E+00   4.3873E+00   4.1344E+00   3.8387E+00   3.8148E+00 
  3.2499E+00   3.0570E+00   2.8812E+00   2.6316E+00   2.3044E+00