# 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':
3.141593     0.600000
# number of Fourier moments for surface 'nm_srf':
1
# number of relative azimuths, values in degrees and radians:
4
  0.00   0.000000000000e+00
180.00   3.141592653590e+00
 30.00   5.235987755983e-01
150.00   2.617993877991e+00
# number of solar zenith angles, values in degrees and mu0=cos(sza):
2
60.00   0.500000000000000
84.26   0.100000000000000
# number of view zenith angles, values in degrees and mu=-cos(vza):
3
 0.00   -1.000000000000000
60.00   -0.500000000000000
84.26   -0.100000000000000
# number of Gauss nodes in hemisphere, zeros, and weights:
24
0.002406390001489   0.006170614899993
0.012635722014345   0.014265694314467
0.030862723998634   0.022138719408710
0.056792236497799   0.029649292457718
0.089999007013049   0.036673240705536
0.129937904210723   0.043095080765976
0.175953174031512   0.048809326052057
0.227289264305580   0.053722135057983
0.283103246186977   0.057752834026863
0.342478660151918   0.060835236463902
0.404440566263192   0.062918728173414
0.467971553568697   0.063969097673376
0.532028446431303   0.063969097673376
0.595559433736808   0.062918728173414
0.657521339848082   0.060835236463902
0.716896753813023   0.057752834026863
0.772710735694420   0.053722135057983
0.824046825968488   0.048809326052057
0.870062095789277   0.043095080765976
0.910000992986951   0.036673240705536
0.943207763502200   0.029649292457718
0.969137276001366   0.022138719408710
0.987364277985655   0.014265694314467
0.997593609998511   0.006170614899993
# EOF