# shape of arrays: Jm[naz*nmu0*nmu3],     Jp[nm_srf*nmu0*ng3],
#                  Tm[nm_srf*nmu3*ng3],   Rp[nm_srf*ng3*ng3]
# WARNING: Jp[ ] does *not* contain Kronecker delta, (2-d0m)
# py: data = np.fromfile('*.bin', dtype=np.float64)
# py: arr[0:ny, 0:nx] = np.reshape(data, (ny, nx))
# py: dtype(*.bin) = float_64bit

# To flux, 'Fo', and total optical thickness, 'tau0':
1.000000     0.630310
# number of Fourier moments for surface 'nm_srf':
3
# number of relative azimuths, values in degrees and radians:
4
  0.00   0.000000000000e+00
180.00   3.141592653590e+00
 60.00   1.047197551197e+00
240.00   4.188790204786e+00
# number of solar zenith angles, values in degrees and mu0=cos(sza):
2
30.00   0.866025403784439
60.00   0.500000000000000
# number of view zenith angles, values in degrees and mu=-cos(vza):
13
 0.00   -1.000000000000000
 5.00   -0.996194698091746
10.00   -0.984807753012208
15.00   -0.965925826289068
20.00   -0.939692620785908
25.00   -0.906307787036650
30.00   -0.866025403784439
35.00   -0.819152044288992
40.00   -0.766044443118978
45.00   -0.707106781186547
50.00   -0.642787609686539
55.00   -0.573576436351046
60.00   -0.500000000000000
# number of Gauss nodes in hemisphere, zeros, and weights:
8
0.019855071751232   0.050614268145185
0.101666761293187   0.111190517226687
0.237233795041836   0.156853322938944
0.408282678752175   0.181341891689181
0.591717321247825   0.181341891689181
0.762766204958164   0.156853322938944
0.898333238706813   0.111190517226687
0.980144928248768   0.050614268145185
# EOF