Code to provide micro physics subroutines for radiation calculation
fms_mod
constants_mod
rad_utilities_mod
longwave_params_mod
esfsw_parameters_mod
call microphys_rad_init (cldhm_abs_in, cldml_abs_in)
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] |
call microphys_sw_driver (is, ie, js, je, Cloud_microphysics, & Cloud_rad_props, donner_flag )
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] |
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] |
call microphys_lw_driver (is, ie, js, je, Cloud_microphysics, & Cloud_rad_props, donner_flag )
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] |
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] |
call isccp_microphys_sw_driver (is, js, iswband, Cloud_microphysics,cldext )
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] |
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] |
call isccp_microphys_lw_driver (is, js, ilwband, & Cloud_microphysics,abscoeff )
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] |
abscoeff | abscoeff absorption coefficient for
clouds in each of the longwave frequency bands [ km**(-1) ] [real] |
call lwemiss_calc ( deltaz, abscoeff, cldemiss)
deltaz | Pressure layer thickness [real] |
abscoeff | Absorption coefficient [real] |
cldemiss | Emissivity calculated from absorption coefficient [real] |
call comb_cldprops_calc ( deltaz, & cldext, cldsct, cldasymm, abscoeff, & Lsc_microphys, Cell_microphys, & Meso_microphys, Lscrad_props, & Cellrad_props, Mesorad_props)
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] |
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] |
call microphys_rad_end
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)
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] |
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] |
call slingo & (conc_drop , size_drop , & cldextivlliq, cldssalbivlliq, cldasymmivlliq, & starting_band, ending_band )
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 |
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] |
call savijarvi & (conc_rain , size_rain , & cldextivlrain, cldssalbivlrain, cldasymmivlrain)
conc_rain | the rain drop water concentration in grams / meter**3 [real] |
size_rain | the rain drop effective diameter in microns [real] |
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] |
call fu & (conc_ice , size_ice , & cldextivlice, cldssalbivlice, cldasymmivlice)
conc_ice | the ice water concentation in grams / meter**3 [real] |
size_ice | the ice crystal effective size in microns [real] |
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] |
call icesolar & (conc_ice , size_ice , & cldextivlice, cldssalbivlice, cldasymmivlice)
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] |
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] |
call snowsw & (conc_snow, & cldextivlsnow, cldssalbivlsnow, cldasymmivlsnow)
conc_snow | the snow concentration in grams / meter**3 [real] |
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] |
subroutine cloud_lwpar (nonly, nbmax, nnn, & size_drop, size_ice, size_rain, & conc_drop, conc_ice, conc_rain, conc_snow, & do_dge_lw, abscoeff)
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] |
abscoeff | cloud absorption coefficient [real] |
subroutine cloud_lwem_oneband & (conc_drop, conc_ice, size_drop, size_ice, & abscoeff)
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] |
abscoeff | cloud absorption coefficient [real] |
call el (conc_ice , size_ice , & cldextbndicelw, cldssalbbndicelw, cldasymmbndicelw)
conc_ice | the ice crystal concentation in grams / meter**3 [real] |
size_ice | the ice crystal effective size in microns [real] |
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] |
call el_dge ( conc_ice , size_ice , & cldextbndicelw, cldssalbbndicelw, cldasymmbndicelw)
conc_ice | the ice crystal concentation in grams / meter**3 [real] |
size_ice | the ice crystal effective size in microns [real] |
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] |
call cliqlw (conc_drop, cldextbnddroplw)
conc_drop | the cloud drop concentration in grams / meter**3 [real] |
cldextbnddroplw | the specified values of the extinction
coefficient for cloud drops in kilometers**(-1)
over wavenumber bands used by the radiation code [real] |
subroutine furainlw & (conc_rain , & cldextbndrainlw, cldssalbbndrainlw, cldasymmbndrainlw)
conc_rain | the rain drop water concentration in grams / meter**3 [real] |
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] |
subroutine fusnowlw & (conc_snow , & cldextbndsnowlw, cldssalbbndsnowlw, cldasymmbndsnowlw)
conc_snow | the snow drop water concentration in grams / meter**3 [real] |
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] |