Mesoscale Parallel Experiment Log

Experiment Name
Eta-32 with modified cloud microphysics
Parallel Slot
Control Slot
Start date of parallel experiment
12Z 2002/03/19
End date of parallel experiment
00Z 2002/07/09
Environmental Modeling Center scientists
Brad Ferrier, Yi Jin, Ying Lin, Eric Rogers
Abstract (including Motivation, Hypothesis and Method)
The Eta-32 is being used to test the impact modifications to the cloud microphysics

  • I. Model changes

    - 1. New grid-scale cloud properties are written and read in the restart files, allowing for their cycling in the forecast system. Old arrays from the Zhao et al. (1997) microphysics were removed from all code.
    - 2. Arrays defining convective cloud-base and cloud-top properties used in radiation are now properly reinitialized after being written to the model restart files. Code was restructured to have convective cloud arrays used in radiation, then written to restart files, and then reinitialized to no-cloud values using a single set of arrays. Additional arrays that separate shallow from deep convection have also been added to the model restart files.
    - 3. Code is restructured so that cycling will produce the same forecast as a longer forecast. For example, the forecasts will be the same if the model is run for two 3-h forecasts (e.g., two EDAS cycles with no assimilation of observations) compared to one 6-h forecast. In other words, the model forecast does not depend on the reading of arrays that are read in at the end of the restart files when restarting the model from the middle of the forecast (i.e., a nonzero start time).
    - 4. Whole array operations were found to be inefficient and were replaced with explicit do-loop indexing.
    - 5. Lower limits were placed to turbulent kinetic energy (Q2) throughout the code, set to values of 0.2. This change affects the writing of Q2 to restart files but does not affect the model integration, since lower limits of Q2 are internally enforced in TURBL.
    - 6. Longwave radiation is updated every hour along with shortwave radiation.
    - 7. Inline compiler directives have been added at the beginning of numerous subroutines, which allow for easy debugging of code. These directives override compiler directives specified in the make files, and they are included in those subroutines that call MPI routines. Now the code can be run with full trapping, checks on array bounds, and checks on argument lists by just manipulating the makefile.
    - 8. The lower limit for all hydrometeor mixing ratios is set to 1.e-12 for all model routines (parameter EPSQ).
    - 9. Convective precipitation rates used to define optical properties of convective clouds are now modified in the precipitation assimilation code.
    - 10. Some of the values in array ZEFFIJ along the eastern boundary were uninitialized. A small patch of code was added in GOSSIP to assign values at these locations to values at the next grid point west (in from) the eastern boundary.
    - 11. Tunable parameters in the microphysics defined based on grid resolution in subroutine GSMCONST are restructured more clearly. These quantities are assumed to vary linearly as a function grid resolution, but this functional dependence can be easily changed in the future.
    - 12. Probablistic freezing of rain as a function of temperature, as parameterized by Lin et al. (1983) following Biggs (1953), has been correctly coded in the grid-scale microphysics in GSMCONST and in GSMCOLUMN. What is currently coded in operations is producing underflow errors and is effectively not activated, however, this process is an extremely slow process and is not expected to have much impact on precipitation forecasts.
    - 13. A small bug was fixed in SFCDIF, in which on rare occasions subscript ML for array ZEFF was out of range (i.e., larger than the upper limit of 4). No discernible impact on the forecast is expected.
    - 14. No dependence of cloud type (convective or grid scale) was seen in the GFDL radiation driver, so the array ITYP was removed from RADFS.
    - 15. The following changes were made in the representation of clouds in radiation.
    a. Convective cloud fraction as a function of precipitation rate of Slingo (1987) was increased by 20%, allowing for shallow, nonprecipitating convective clouds to have an assumed cloud fraction of 10%.
    b. Grid-scale cloud fraction is now parameterized following Xu and Randall (1996), replacing the Randall (1994) formulation. Also, relative humidity is now calculated using the same approach as in the EGCP01 grid-scale microphysics. One change was made with respect to Xu and Randall, however, in that a relative humidity of 100% or more will not produce a cloud fraction of 1.0. Grid-scale cloud must be present in order for nonzero cloud fractions to be represented. Grid-scale cloud is defined as the sum of the mixing ratios of cloud water and total ice (i.e., ignoring rain).
    c. The radiative effects from low-level clouds in the lowest 100 mb and upper-level clouds above the tropopause are no longer ignored.
    d. Absorption coefficients for convective clouds is assumed to be 0.16 at temperatures <= -10C assuming ice is the dominant hydrometeor category as in the grid scheme, 0.08 at warmer temperatures assuming water is dominant at warmer temperatures. The values for these coefficients is based on Harshvardhan et al. (1989). Absorption coefficients for grid-scale cloud water is 0.08*min(1., Qw/Q0), where Qw is the cloud water mixing ratio and Q0=0.1e-3 kg/kg. Absorption coefficients for grid-scale ice is equal to 500*Qi, where Qi is the mixing ratio for all ice particles (kg/kg). This value was derived based on the optical properties for snow provided by Q. Fu (U. Washington, personal communication) for ice particles. Operational code currently assumes the absorption coefficients vary as functions of temperature for convective and grid-scale clouds.
    e. Many of the algorithms for calculating cloud fraction were simplified and numerous comments were added to allow future modifications to be made much more easily and effectively.

  • II. Post changes

    -1. Read in new, expanded restart files.
    -2. Post new cloud fields (cloud water, cloud ice, rain, snow, and total condensate) on pressure and Eta surfaces,.as well as original model arrays (F_rain, F_ice, F_RimeF) on Eta surfaces for debugging.
    -3. Post total cloud fraction from grid-scale and convective clouds on Eta surfaces using the same algorithms as in the Eta radiation driver (this is not true in the operational Eta post). Also post cloud-base and cloud-top pressures from shallow, nonprecipitating convection, deep convection, and grid-scale convection. The cloud efficiency parameter used in the convective parameterization can now also be posted.
    -4. Visibility calculations have been changed to use the new cloud fields, responding to mixing ratios of cloud water, cloud ice, rain, and snow.
    -5. The same method is used to calculate relative humidity (RH) throughout the post. This currently is defined as RH with respect to water at all levels, however, it will be very easy to calculate RH with respect to ice(if desired) in the future. There were MANY subroutines that used the ice/water flag from Zhao and Carr (1997) for calculating RH, and the code was poorly structured and not always uniform throughout the post. These blocks of code have been substituted with a simple subroutine call.
    -6. The same lower limit for the condensate mixing ratio of EPSQ=1.e-12 is used throughout the post as it is used throughout the model. Small changes were made to the include parameter file "params" that removed the need to use the parameter file "cuparm" in WETBULB. The "cuparm" file is no longer used in the post.
    -7. Inline compiler directives have been added to the beginning of numerous subroutines, as in the model, to allow for easy compiling of code with or without traps, checks on array bounds, and checks on argument lists by simply changing lines in the make file.
    -8. Bug fixes identified by Tuccillo have been made to the J index in FIXED, and an argument to MPI_ISEND in GRIBIT was changed that solved the problem of the post occasionally failing due to segmentation faults.

  • III. Brief summary of changes

    Changes were made allows the new cloud arrays to be written to the model restart files and posted to grib files. Numerous changes were made in the optical properties of grid-scale and convective clouds that should increase the areal coverage of clouds, including modifications in the calculation of cloud fraction and the absorption coefficient of clouds. These changes will result in the radiative properties of the clouds being more consistent with the other physical parameterizations. These changes are also accurately reflected in the post, allowing the end user to review many more aspects of the convective and grid-scale cloud properties simulated in the Eta.
  • Experiment changes log
    Implementation changes list
    Product changes list
    Background links
    Evaluation of parallel results
    Daily forecast maps
    Daily forecast stats
    Verification of precipitation and against rawinsondes / surface data

    Turned off to conduct clean test of new 3DVAR w/updated radiance processing (including NOAA-16 HIRS channel)

    Mesoscale Parallel Home Page


    Page Last Modified: March 19, 2002