patch.v1.4.1_to_v1.4.2.tar

patch.cosp_modis_simulator.F90.txt0000644000446700022600000000111613211553676017367 0ustar  dswalesclimate--- a/cosp_modis_simulator.F90
+++ b/cosp_modis_simulator.F90
@@ -387,8 +387,6 @@ contains
     allocate(x%Optical_Thickness_vs_ReffICE(nPoints, numModisTauBins+1, numModisReffIceBins))
     allocate(x%Optical_Thickness_vs_ReffLIQ(nPoints, numModisTauBins+1, numModisReffLiqBins))
     x%Optical_Thickness_vs_Cloud_Top_Pressure(:, :, :) = R_UNDEF
-    x%Optical_Thickness_vs_ReffICE(:, :, :)            = R_UNDEF
-    x%Optical_Thickness_vs_ReffLIQ(:, :, :)            = R_UNDEF
   END SUBROUTINE CONSTRUCT_COSP_MODIS
 
   !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
patch.cosp_simulator.F90.txt0000644000446700022600000000176313211553722016174 0ustar  dswalesclimate--- a/cosp_simulator.F90
+++ b/cosp_simulator.F90
@@ -204,14 +204,6 @@ SUBROUTINE COSP_SIMULATOR(gbx,sgx,sghydro,cfg,vgrid,sgradar,sglidar,isccp,misr,m
     where(modis%Optical_Thickness_vs_Cloud_Top_Pressure /= R_UNDEF) modis%Optical_Thickness_vs_Cloud_Top_Pressure = &
                                                       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
   endif
patch.modis_simulator.F90.txt0000644000446700022600000000140213211553657016340 0ustar  dswalesclimate--- a/MODIS_simulator/modis_simulator.F90
+++ b/MODIS_simulator/modis_simulator.F90
@@ -862,10 +862,10 @@ contains
     g(:) = 0; w0(:) = 0. 
     tau(:) = ice_tau(:) + water_tau(:) 
     where (tau(:) > 0) 
-      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(:))
+      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(:))
     end where
     
     compute_nir_reflectance = compute_toa_reflectace(tau, g, w0)
README0000644000446700022600000000211613211553550011605 0ustar  dswalesclimate####################################################################################	 
COSPv1.4.2 patch for v1.4.1.
Purpose: This patch corrects a small error in the MODIS simulator. This error affects
         mixed-phase clouds visible at the top of the atmosphere. In addition, there
	 are two omissions from the previous version (v1.4.1) included. 1) The new
	 MODIS particle size diagnostics in v1.4.1 were output as fractional values,
	 not as percent. 2) These diagnsotics also needed to be set to undefined at
	 initialization.
Dustin Swales (dusitn.swales@noaa.gov)
####################################################################################	 
	 
There are 3 COSP source files impacted by this update:
1) MODIS_simulator/modis_simulator.F90
2) cosp_simulator.F90
3) cosp_modis_simulator.F90

The following commands (patch FileToPatch Patch) will apply the patches:
patch MODIS_simulator/modis_simulator.F90 patch.modis_simulator.F90.txt
patch cosp_simulator.F90                  patch.cosp_simulator.F90.txt
patch cosp_modis_simulator.F90            patch.cosp_modis_simulator.F90.txt