additional changes

Author: mlee03

full/atm_land_ice_flux_exchange.F90

@@ -1998,12 +1998,14 @@ subroutine sfc_boundary_layer ( dt, Time, Atm, Land, Ice, Land_Ice_Atmos_Boundar
 
 
     !> ON THE EXCHANGE GRID, COMPUTE T_SURF**4
-    !! NOTE, TO COMPUTE FLUXES, T_SURF**4 IS SENT TO THE EXCHANGE GRID AND NOT T_SURF DUE TO NONLINEARITY IN
-    !! THE STEFAN-BOLTZMANN LAW WHERE LOGWAVE_FLUX = STEFAN_BOLTZMANN _CONSTANT * T**4.
+    !! NOTE, TO COMPUTE LONGWAVE RADIATION FLUXES MORE ACCURATELY, T_SURF**4 IS
+    !! PUT ON THE EXCHANGE GRID RATHER THAN T_SURF IN ORDER TO ACCOUNT FOR THE NONLINEARITY IN
+    !! THE STEFAN-BOLTZMANN LAW WHERE LONGWAVE_FLUX = STEFAN_BOLTZMANN _CONSTANT * T**4.
     !! ON THE EXCHANGE GRID, AS QUANTITIES ARE REMAPPED, FIELDS ARE AREA-WEIGHTED (AVERAGED)
     !! SUCH THAT OUTPUT_TEMPERATURE = SUM(INPUT_TEMPERATURE * (XGRID_AREA)/(INPUT_GRID_AREA))
     !! WHERE THE SUM IS OVER ALL XGRID CELLS THAT OVERLAP WITH THE OUTPUT CELL. 
     !! BECAUSE OF THIS WEIGHTING, THE COMPUTED FLUX WOULD DIFFER FROM USING <T**4> VS <T>**4
+    !! (WHERE <> IS USED TO DENOTE AVERAGING)
     !{
     !$OMP parallel do default(none) shared(my_nblocks,block_start,block_end,ex_avail,  &
     !$OMP                                  ex_drag_q,ex_wind,ex_cd_q,ex_t_surf4,ex_t_surf ) &