Index: LMDZ6/trunk/libf/phylmd/cosp/cosp_simulator.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/cosp/cosp_simulator.F90	(revision 3182)
+++ LMDZ6/trunk/libf/phylmd/cosp/cosp_simulator.F90	(revision 3184)
@@ -205,4 +205,13 @@
                                                       modis%Optical_Thickness_vs_Cloud_Top_Pressure*100.0
   endif
+  if (cfg%Lcrimodis) then
+     where(modis%Optical_Thickness_vs_ReffICE /= R_UNDEF) modis%Optical_Thickness_vs_ReffICE = &
+                                                      modis%Optical_Thickness_vs_ReffICE*100.0
+  endif
+  if (cfg%Lcrlmodis) then
+     where(modis%Optical_Thickness_vs_ReffLIQ /= R_UNDEF) modis%Optical_Thickness_vs_ReffLIQ = &
+                                                      modis%Optical_Thickness_vs_ReffLIQ*100.0
+  endif
+
   if (cfg%Lcltmodis) then
     where(modis%Cloud_Fraction_Total_Mean /= R_UNDEF) modis%Cloud_Fraction_Total_Mean = modis%Cloud_Fraction_Total_Mean*100.0
Index: LMDZ6/trunk/libf/phylmd/cosp/modis_simulator.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/cosp/modis_simulator.F90	(revision 3182)
+++ LMDZ6/trunk/libf/phylmd/cosp/modis_simulator.F90	(revision 3184)
@@ -863,8 +863,8 @@
     tau(:) = ice_tau(:) + water_tau(:) 
     where (tau(:) > 0) 
-      g(:)  = (water_tau(:) * water_g(:)                + ice_tau(:) * ice_g(:)            ) / & 
-              tau(:) 
-      w0(:) = (water_tau(:) * water_g(:) * water_w0(:)  + ice_tau(:) * ice_g(:) * ice_w0(:)) / &
-              (g(:) * tau(:))
+      w0(:) = (water_tau(:) * water_w0(:)  + ice_tau(:) * ice_w0(:)) / &
+              tau(:)
+      g(:) = (water_tau(:) * water_g(:) * water_w0(:)  + ice_tau(:) * ice_g(:) * ice_w0(:)) / &
+             (w0(:) * tau(:))
     end where