PUBLIC INTERFACE ~ PUBLIC DATA ~ PUBLIC ROUTINES ~ NAMELIST ~ DIAGNOSTIC FIELDS ~ ERROR MESSAGES ~ REFERENCES ~ NOTES

Module microphys_rad_mod

Contact:  fil
Reviewers:  smf
Change History: WebCVS Log


OVERVIEW

Code to provide micro physics subroutines for radiation calculation



OTHER MODULES USED

             fms_mod
constants_mod
rad_utilities_mod
longwave_params_mod
esfsw_parameters_mod

PUBLIC INTERFACE

microphys_rad_init:
The microphys_rad module constructor
microphys_sw_driver:
Subroutine to deploy micro physics radiation calculation, particularly for cloud parameterizations in the shortwave
microphys_lw_driver:
Subroutine to deploy micro physics radiation calculation, particularly for cloud parameterizations in the longwave
isccp_microphys_sw_driver:
Subroutine to deploy micro physics radiation calculation, particularly for cloud parameterizations in the shortwave
isccp_microphys_lw_driver:
Subroutine to deploy micro physics radiation calculation, particularly for cloud parameterizations in the longwave
lwemiss_calc:
Subroutine to compute infrared emissivity from the absorption coefficient
comb_cldprops_calc:
Subroutine to define the total-cloud radiative properties
microphys_rad_end:
microphys_rad_end is the destructor for microphys_rad_mod.
cloudpar:
Subroutine to determine cloud single scattering parameters
slingo:
Subroutine to determine single scattering parameters for clouds
savijarvi:
Subroutine to define the single scattering parameters for rain drop
fu:
Subroutine to define the single scattering parameters for ice crystals
icesolar:
Subroutine to define the single scattering parameters for ice crystals
snowsw:
Subroutine to define the single scattering parameters for snow
cloud_lwpar:
Subroutine to determine cloud infrared emissivity
cloud_lwem_oneband:
Subroutine to determine a single broadband cloud infrared emissivity
el:
Subroutine to calculates total optical depth and scattering optical depth
el_dge:
Subroutine to calculate total optical depth and scattering optical depth in infrared.
cliqlw:
Subroutine to calculate longwave absorption optical depth for liquid
furainlw:
Subroutine to calculate total optical depth and scattering optical depth in infrared for cloud rain water
fusnowlw:
Subroutine to calculate total optical depth and scattering optical depth in infrared for cloud rain water


PUBLIC DATA

None.


PUBLIC ROUTINES

  1. microphys_rad_init

    call microphys_rad_init (cldhm_abs_in, cldml_abs_in)
    DESCRIPTION
    This subroutine initializes micro physics module data, determines micro physics parameterization scheme based on initialization input data file.


    INPUT
    cldhm_abs_in    boundaries in sigma pressure level between high and middle clouds
       [cldhm_abs_in]
    cldml_abs_in    boundaries in sigma pressure level between middle and low clouds
       [cldml_abs_in]

  2. microphys_sw_driver

    call microphys_sw_driver (is, ie, js, je, Cloud_microphysics, & Cloud_rad_props, donner_flag )
    DESCRIPTION
    This subroutine takes cloud micro physics parameters and calculate broad band cloud radiation parameters. For example the input parameters are cloud droplet concentration, size, and composition; the output parameters are cloud extinction coefficient, scattering coefficient, and assymmetry parameters.


    INPUT
    is    starting indice of x dimension in current physics window
       [integer]
    ie    ending indice of x dimension in current physics window
       [integer]
    js    starting indice of y dimension in current physics window
       [integer]
    je    ending indice of y dimension in current physics window
       [integer]
    Cloud_microphysics    derived type variable containing cloud microphysical specification information
       [microphysics_type]
    donner_flag    OPTIONAL: logical flag which if present indicates that clouds from donner_deep_mod are being processed, and that an ice parameterization associated with that scheme (which differs from that used by strat_cloud_mod) is to be used.
       [logical]

    INPUT/OUTPUT
    Cloud_rad_props    derived type variable containing the micro- physically-based sw cloud radiative proper- ties [ microrad_properties_type ] the components defined in this routine: %cldext parameterization band values of the cloud extinction coefficient [ km**(-1) ] %cldsct parameterization band values of the cloud scattering coefficient [ km**(-1) ] %cldasymm parameterization band values of the asymmetry factor [ dimensionless ]
       [microrad_properties_type]

  3. microphys_lw_driver

    call microphys_lw_driver (is, ie, js, je, Cloud_microphysics, & Cloud_rad_props, donner_flag )
    DESCRIPTION
    This subroutine takes cloud micro physics parameters and calculate broad band cloud radiation parameters. For example the input parameters are cloud droplet concentration, size, and composition; the output parameters are cloud extinction coefficient, scattering coefficient, and assymmetry parameters.


    INPUT
    is    starting indice of x dimension in current physics window
       [integer]
    ie    ending indice of x dimension in current physics window
       [integer]
    js    starting indice of y dimension in current physics window
       [integer]
    je    ending indice of y dimension in current physics window
       [integer]
    Cloud_microphysics    derived type variable containing cloud microphysical specification information
       [microphysics_type]
    donner_flag    OPTIONAL: logical flag which if present indicates that clouds from donner_deep_mod are being processed, and that an ice parameterization associated with that scheme (which differs from that used by strat_cloud_mod) is to be used.
       [logical]

    INPUT/OUTPUT
    Cloud_rad_props    derived type variable containing the micro- physically-based sw cloud radiative proper- ties [ microrad_properties_type ] the components defined in this routine: %cldext parameterization band values of the cloud extinction coefficient [ km**(-1) ] %cldsct parameterization band values of the cloud scattering coefficient [ km**(-1) ] %cldasymm parameterization band values of the asymmetry factor [ dimensionless ]
       [microrad_properties_type]

  4. isccp_microphys_sw_driver

    call isccp_microphys_sw_driver (is, js, iswband, Cloud_microphysics,cldext )
    DESCRIPTION
    This subroutine takes cloud micro physics parameters and calculate broad band cloud radiation parameters. For example the input parameters are cloud droplet concentration, size, and composition; the output parameters are cloud extinction coefficient, scattering coefficient, and assymmetry parameters.


    INPUT
    is    starting indice of x dimension in current physics window
       [integer]
    js    starting indice of y dimension in current physics window
       [integer]
    iswband    swband whose extinction we desire
       [integer]
    Cloud_microphysics    derived type variable containing cloud microphysical specification information
       [microphysics_type]

    OUTPUT
    cldext     derived type variable containing the micro- physically-based sw cloud radiative proper- ties [ microrad_properties_type ] the components defined in this routine: %cldext parameterization band values of the cloud extinction coefficient [ km**(-1) ]
       [real]

  5. isccp_microphys_lw_driver

    call isccp_microphys_lw_driver (is, js, ilwband, & Cloud_microphysics,abscoeff )
    DESCRIPTION
    This subroutine takes cloud micro physics parameters and calculate broad band cloud radiation parameters. For example the input parameters are cloud droplet concentration, size, and composition; the output parameters are cloud extinction coefficient, scattering coefficient, and assymmetry parameters.


    INPUT
    is    starting indice of x dimension in current physics window
       [integer]
    js    starting indice of y dimension in current physics window
       [integer]
    ilwband    lwband whose absorption we desire
       [integer]
    Cloud_microphysics    derived type variable containing cloud microphysical specification information
       [microphysics_type]

    OUTPUT
    abscoeff    abscoeff absorption coefficient for clouds in each of the longwave frequency bands [ km**(-1) ]
       [real]

  6. lwemiss_calc

    call lwemiss_calc ( deltaz, abscoeff, cldemiss)
    DESCRIPTION
    Subroutine to compute infrared emissivity from the absorption coefficient


    INPUT
    deltaz    Pressure layer thickness
       [real]
    abscoeff    Absorption coefficient
       [real]

    OUTPUT
    cldemiss    Emissivity calculated from absorption coefficient
       [real]

  7. comb_cldprops_calc

    call comb_cldprops_calc ( deltaz, & cldext, cldsct, cldasymm, abscoeff, & Lsc_microphys, Cell_microphys, & Meso_microphys, Lscrad_props, & Cellrad_props, Mesorad_props)
    DESCRIPTION
    Subroutine to define the total-cloud radiative properties to be seen by the radiation package, obtained by the appropriate combination of the large-scale, mesoscale and cell-scale clouds present in a grid box.


    INPUT
    Lsc_microphys    microphysical specification for large-scale clouds
       [microphysics_type]
    Cell_microphys    microphysical specification for convective cell clouds associated with donner convection
       [microphysics_type]
    Meso_microphys    microphysical specification for meso-scale clouds assciated with donner convection
       [microphysics_type]
    Lscrad_props    cloud radiative properties for the large-scale clouds
       [microrad_properties_type]
    Cellrad_props    cloud radiative properties for the convective cell clouds associated with donner convection
       [microrad_properties_type]
    Mesorad_props    cloud radiative properties for the meso-scale clouds associated with donner convection
       [microrad_properties_type]

    INPUT/OUTPUT
    cldext    parameterization band values of the cloud extinction coefficient [ km**(-1) ]
       [real]
    cldsct    parameterization band values of the cloud scattering coefficient [ km**(-1) ]
       [real]
    cldasymm    parameterization band values of the asymmetry factor [ dimensionless ]
       [real]
    abscoeff    combined absorption coefficient for clouds in each of the longwave frequency bands [ km**(-1) ]
       [real]

  8. microphys_rad_end

    call microphys_rad_end 
    
    DESCRIPTION
    microphys_rad_end is the destructor for microphys_rad_mod.


  9. cloudpar

    call cloudpar ( nonly, nbmax, nnn, & size_drop, size_ice, size_rain, conc_drop, & conc_ice, conc_rain, conc_snow, do_dge_sw, & cldext, cldsct, cldasymm)
    DESCRIPTION
    determine the parameterization band values of the single scattering parameters (extinction coefficient, scattering coefficient and asymmetry factor) for clouds from the size and/or concentration of each constituent (cloud drops, rain drops, ice crystals and snow) present.


    INPUT
    nonly    The single band for calculations. Note that this is used only in the case of a call from cloudrad_diagnostics to do isccp simulator work. For all other calls, nonly should be 0 and will have no effect on the calculations below
       [integer]
    nbmax    The number of individual bands to do calculations over. Note that for normal GCM calls this will be 1. For calls using stochastic clouds with or without the isccp simulator this will be equal to the number of shortwave bands
       [integer]
    nnn    This integer controls which cloud state to access for radiation calculations. For normal GCM applications this will be 1. For Full Independent Column Approximation calculations with stochast- ic clouds this will be the profile number being accessed.
       [integer]
    conc_drop    total cloud droplet concentration
       [real]
    conc_ice    ice crystal concentration
       [real]
    conc_rain    rain droplet concetration
       [real]
    conc_snow    snow concentration
       [real]
    size_drop    cloud droplet size distribution
       [real]
    size_ice    ice crystal size distribution
       [real]
    size_rain    rain drop size distribution
       [real]
    do_dge_sw    use sw parameterizations using generalized effective size developed by Fu et al (1998) (if true). otherwise use parameterizations by Fu et al using effective size.
       [logical]

    OUTPUT
    cldext    the parameterization band values of the cloud extinction coefficient in kilometer**(-1)
       [real]
    cldsct    the parameterization band values of the cloud scattering coefficient in kilometer**(-1)
       [real]
    cldasymm    the parameterization band values of the asymmetry factor
       [real]

  10. slingo

    call slingo & (conc_drop , size_drop , & cldextivlliq, cldssalbivlliq, cldasymmivlliq, & starting_band, ending_band )
    DESCRIPTION
    define the single scattering parameters for cloud drops using the Slingo parameterization for his spectral intervals. slingo, a., a gcm parameterization of the shortwave properties of water clouds., j. atmos. sci.,46, 1419-1427, 1989.


    INPUT
    conc_drop    the cloud drop liquid water concentration in grams meter**3
       [real]
    size_drop    the cloud drop effective diameter in microns
       [real]
    starting_band   
    ending_band   

    OUTPUT
    cldextivlliq    The specified spectral values of the extinction coefficient in kilometer**(-1) for drops
       [real]
    cldssalbivlliq    the specified spectral values of the single-scattering albedo for drops
       [real]
    cldasymmivlliq    the specified spectral values of the asymmetry factor for drops
       [real]

  11. savijarvi

    call savijarvi & (conc_rain , size_rain , & cldextivlrain, cldssalbivlrain, cldasymmivlrain)
    DESCRIPTION
    define the single scattering parameters for rain drops using the Savijarvi parameterization for his spectral intervals.


    INPUT
    conc_rain    the rain drop water concentration in grams / meter**3
       [real]
    size_rain    the rain drop effective diameter in microns
       [real]

    OUTPUT
    cldextivlrain    the specified spectral values of the extinction coefficient for rain in kilometers**(-1)
       [real]
    cldssalbivlrain    the specified spectral values of the single- scattering albedo for rain
       [real]
    cldasymmivlrain    the specified spectral values of the asymmetry factor for rain
       [real]

  12. fu

    call fu & (conc_ice , size_ice , & cldextivlice, cldssalbivlice, cldasymmivlice)
    DESCRIPTION
    define the single scattering parameters for ice crystals using the Fu parameterization for his spectral intervals.

    references:

    fu, q., an accurate parameterization of the solar radiative properties of cirrus clouds for climate models., j. climate, 9, 2058-2082, 1996.

    notes: the ice crystal effective size (D^sub^ge in his paper) can only be 18.6 <= D^sub^ge <= 130.2 microns.

    the single scattering properties for wavenumbers < 2000 cm-1 are assigned the values in the first interval, since the formulation is not defined for those wavenumbers.

    the extinction coefficient is converted to kilometer**(-1) the unit utilized by the shortwave routine Swresf.

    a value of 1.0E-100 is added to the size so that no division by zero occurs when the size is zero, in defining the extinction coefficient.


    INPUT
    conc_ice    the ice water concentation in grams / meter**3
       [real]
    size_ice    the ice crystal effective size in microns
       [real]

    OUTPUT
    cldextivlice    the specified spectral values of the extinction coefficient for ice particles in kilometers**(-1)
       [real]
    cldssalbivlice    the specified spectral values of the single- scattering albedo for ice particles
       [real]
    cldasymmivlice    the specified spectral values of the asymmetry factor for ice particles
       [real]

  13. icesolar

    call icesolar & (conc_ice , size_ice , & cldextivlice, cldssalbivlice, cldasymmivlice)
    DESCRIPTION
    define the single scattering parameters for ice crystals using the Fu parameterization for his spectral intervals.

    references:

    fu and Liou (1993, JAS)

    notes: the ice crystal effective size (D^sub^e in paper) can only be 18.6 <= D^sub^e <= 130.2 microns.

    the single scattering properties for wavenumbers < 2000 cm-1 are assigned the values in the first interval, since the formulation is not defined for those wavenumbers.

    the extinction coefficient is converted to kilometer**(-1) the unit utilized by the shortwave routine Swresf.

    a value of 1.0E-100 is added to the size so that no division by zero occurs when the size is zero, in defining the extinction coefficient.


    INPUT
    conc_ice    the ice water concentation in grams / meter**3
       [real]
    size_ice    the ice crystal effective size in microns Corresponds to minimum dimension of hexagonal crystal.
       [real]

    OUTPUT
    cldextivlice    the specified spectral values of the extinction coefficient for ice particles in kilometers**(-1)
       [real]
    cldssalbivlice    the specified spectral values of the single- scattering albedo for ice particles
       [real]
    cldasymmivlice    the specified spectral values of the asymmetry factor for ice particles
       [real]

  14. snowsw

    call snowsw & (conc_snow, & cldextivlsnow, cldssalbivlsnow, cldasymmivlsnow)
    DESCRIPTION
    define the single scattering parameters for snow using the Fu parameterization for his spectral intervals. author: leo donner, gfdl, 11 Sept 98

    references:

    fu, q., et al., (See notes from Kuo-Nan Liou, 1 Sept 98). (SNOW)

    notes: the single scattering properties for wavenumbers < 2500 cm-1 are assigned the values in the first interval, since the formulation is not defined for those wavenumbers.

    the extinction coefficient is in units of kilometer**(-1)


    INPUT
    conc_snow    the snow concentration in grams / meter**3
       [real]

    OUTPUT
    cldextivlsnow    the specified spectral values of the extinction coefficient for snow in kilometers**(-1)
       [real]
    cldssalbivlsnow    the specified spectral values of the single- scattering albedo for snow
       [real]
    cldasymmivlsnow    the specified spectral values of the asymmetry factor for snow
       [real]

  15. cloud_lwpar

    subroutine cloud_lwpar (nonly, nbmax, nnn, & size_drop, size_ice, size_rain, & conc_drop, conc_ice, conc_rain, conc_snow, & do_dge_lw, abscoeff)
    DESCRIPTION
    determine the infrared cloud emissivities for specified wavenumber bands from parameterizations for absorption coefficients due to cloud drops, cloud ice crystals, rain and snow. conceptually one could have separate concentrations and sizes for "thin" or randomly overlapped and for maximally overlapped clouds. for now, there is one concentration and size, therefore the two emissivities are set equal.


    INPUT
    nonly    The single band for calculations. Note that this is used only in the case of a call from cloudrad_diagnostics to do isccp simulator work. For all other calls, nonly should be 0 and will have no effect on the calculations below
       [integer]
    nbmax    The number of individual bands to do calculations over. Note that for normal GCM calls this will be 1. For calls using stochastic clouds with or without the isccp simulator this will be equal to the number of longwave bands
       [integer]
    nnn    This integer controls which cloud state to access for radiation calculations. For normal GCM applications this will be 1. For Full Independent Column Approximation calculations with stochast- ic clouds this will be the profile number being accessed.
       [integer]
    conc_drop    total cloud droplet concentration
       [real]
    conc_ice    ice cloud droplet concentration
       [real]
    conc_rain    rain droplet concetration
       [real]
    conc_snow    snow concentration
       [real]
    size_drop    cloud droplet size distribution
       [real]
    size_ice    ice droplet size distribution
       [real]
    size_rain    rain droplet size distribution
       [real]
    do_dge_lw    flag for using dge longwave parameterization
       [logical]

    OUTPUT
    abscoeff    cloud absorption coefficient
       [real]

  16. cloud_lwem_oneband

    subroutine cloud_lwem_oneband & (conc_drop, conc_ice, size_drop, size_ice, & abscoeff)
    DESCRIPTION
    determine the infrared cloud emissivities for a single broadband from parameterizations for absorption coefficients due to cloud drops and cloud ice crystals. the parameterization comes from from the 5-band formulation given by Ebert and Curry (1992, J. Geophys. Res., vol. 97, pp. 3831-3836). S. Klein derived the coefficients for the 1-band version used here.


    INPUT
    conc_drop    total cloud droplet concentration
       [real]
    conc_ice    ice cloud droplet concentration
       [real]
    size_drop    cloud droplet size distribution
       [real]
    size_ice    ice droplet size distribution
       [real]

    OUTPUT
    abscoeff    cloud absorption coefficient
       [real]

  17. el

    call el (conc_ice , size_ice , & cldextbndicelw, cldssalbbndicelw, cldasymmbndicelw)
    DESCRIPTION
    This subroutine calculates total optical depth and scattering optical depth for infrared radiation using Fu and Liou (1993, JAS). To be used for crystal effective sizes from 20 to 130 um. limits changed to 18.6 to 130.2 um on 2 April 1999 to match shortwave limits.


    INPUT
    conc_ice    the ice crystal concentation in grams / meter**3
       [real]
    size_ice    the ice crystal effective size in microns
       [real]

    OUTPUT
    cldextbndicelw    the specified values of the extinction coefficient for ice particles in kilometers**(-1) over wavenumber bands used by the radiation code
       [real]
    cldssalbbndicelw    the specified values of the single- scattering albedo for ice particles over wavenumber bands used by the radiation code
       [real]
    cldasymmbndicelw    the specified values of the asymmetry factor for ice particles over wavenumber bands used by the radiation code
       [real]

  18. el_dge

    call el_dge ( conc_ice , size_ice , & cldextbndicelw, cldssalbbndicelw, cldasymmbndicelw)
    DESCRIPTION
    calculates total optical depth and scattering optical depth for infrared radiation using Fu et al (J. Clim., 11,2223 (1998)). To be used for crystal generalized effective diameters from 18.6 to 130.2 um to match shortwave limits.


    INPUT
    conc_ice    the ice crystal concentation in grams / meter**3
       [real]
    size_ice    the ice crystal effective size in microns
       [real]

    OUTPUT
    cldextbndicelw    the specified values of the extinction coefficient for ice particles in kilometers**(-1) over wavenumber bands used by the radiation code
       [real]
    cldssalbbndicelw    the specified values of the single- scattering albedo for ice particles over wavenumber bands used by the radiation code
       [real]
    cldasymmbndicelw    the specified values of the asymmetry factor for ice particles over wavenumber bands used by the radiation code
       [real]

  19. cliqlw

    call cliqlw (conc_drop, cldextbnddroplw)
    DESCRIPTION
    Calculates longwave absorption optical depth for liquid. Follows Held et al. (J. Atmos. Sci., 1993).

    Leo Donner, GFDL, 1 Feb 1999


    INPUT
    conc_drop    the cloud drop concentration in grams / meter**3
       [real]

    OUTPUT
    cldextbnddroplw    the specified values of the extinction coefficient for cloud drops in kilometers**(-1) over wavenumber bands used by the radiation code
       [real]

  20. furainlw

    subroutine furainlw & (conc_rain , & cldextbndrainlw, cldssalbbndrainlw, cldasymmbndrainlw)
    DESCRIPTION
    Calculates absorption coefficient for cloud rain water for longwave radiation (Fu et al., 1995, J. Atmos. Sci.) To be used for rain water with radii between 60 um and 1.8 mm. See also notes from Q. Fu (4 Sept 98) note: the size of the rain water from the microphysical model (if existing) does not enter into the calculations.

    Leo Donner, GFDL, 20 Mar 99


    INPUT
    conc_rain    the rain drop water concentration in grams / meter**3
       [real]

    OUTPUT
    cldextbndrainlw    the specified values of the extinction coefficient for rain water in kilometers**(-1) over wavenumber bands used by the radiation code
       [real]
    cldssalbbndrainlw    the specified values of the single- scattering albedo for rain water over wavenumber bands used by the radiation code
       [real]
    cldasymmbndrainlw    the specified values of the asymmetry factor for rain water over wavenumber bands used by the radiation code
       [real]

  21. fusnowlw

    subroutine fusnowlw & (conc_snow , & cldextbndsnowlw, cldssalbbndsnowlw, cldasymmbndsnowlw)
    DESCRIPTION
    Calculates absorption coefficient for cloud rain water for longwave radiation (Fu et al., 1995, J. Atmos. Sci.) To be used for rain water with radii between 60 um and 1.8 mm. See also notes from Q. Fu (4 Sept 98) note: the size of the rain water from the microphysical model (if existing) does not enter into the calculations.

    Leo Donner, GFDL, 20 Mar 99


    INPUT
    conc_snow    the snow drop water concentration in grams / meter**3
       [real]

    OUTPUT
    cldextbndsnowlw    the specified values of the extinction coefficient for snow drop water in kilometers**(-1) over wavenumber bands used by the radiation code
       [real]
    cldssalbbndsnowlw    the specified values of the single- scattering albedo for snow drop water over wavenumber bands used by the radiation code
       [real]
    cldasymmbndsnowlw    the specified values of the asymmetry factor for snow drop water over wavenumber bands used by the radiation code
       [real]


DATA SETS

None.


ERROR MESSAGES

None.


top