more work on lunar illumination

modules/util/util.py

@@ -664,6 +664,7 @@ rho_ice = 0.917
 
 
 def compute_lwc_iwc(iphase, reff, opd):
+    iphase = iphase.flatten()
     lwp_dcomp = np.zeros(reff.shape[0])
     iwp_dcomp = np.zeros(reff.shape[0])
     lwp_dcomp[:] = np.nan