SUBROUTINE FLUX_CALC

 -----------------------------------------------------------------
 What it does
 ============

Calculate radiative fluxes.

1. Properties independent of the spectral bands are set.
2. A loop over bands is then entered.
   1. Grey optical properties are set
   2. An appropriate subroutine is called to treat the gaseous overlaps.
   3. The final fluxes are assigned.

 -----------------------------------------------------------------
 Structure of Routine
 ====================

1. Many variables are passed into subroutine.
2. Many variables defined
3. DATA statements
4. Initial determination of flags and switches
5. INITIAL CALCULATIONS FOR SURFACE PROPERTIES
   • CALL COLLECT_SURFACE
6. INITIAL CALCULATIONS FOR AEROSOLS
   • CALL SET_MOIST_AEROSOL_PROPERTIES
7. INITIAL CALCULATIONS FOR CLOUDS
   • CALL SET_CLOUD_POINTER
   • CALL SET_CLOUD_GEOMETRY
   • CALL AGGREGATE_CLOUD
8. CALCULATE ENERGY TRANSFER COEFFICIENTS IN A MIXED COLUMN, OR SPLIT THE ATMOSPHERE INTO COLUMNS WITH A COLUMN MODEL
   • CALL OVERLAP_MIX_MAXIMUM or
   • CALL OVERLAP_MIX_RANDOM or
   • CALL OVERLAP_TRIPLE or
   • CALL SPLIT_MAXIMUM
9. CALCULATE THE ATMOSPHERIC DENSITIES
   • CALL CALCULATE_DENSITY
10. INITIALIZE THE TOTAL FLUXES.
    • CALL INITIALIZE_FLUX
11. INITIALIZATION OF DIAGNOSTICS AND COUPLING ARRAYS FOR THE UNIFIED MODEL.
    • CALL R2_INIT_COUPLE_DIAG
12. SOLVE THE EQUATION OF TRANSFER IN EACH BAND AND INCREMENT THE FLUXES. - DO I_BAND=I_FIRST_BAND, I_LAST_BAND
    1. CALL SET_SCATTERING
    2. CALL SCALE_ABSORB
    3. CALL RESCALE_CONTINUUM
    4. CALCULATE THE GREY EXTINCTION WITHIN THE BAND.
       • CALL GREY_EXTINCTION
    5. CALL RESCALE_ASYMMETRY
    6. CALL DIFF_PLANCK_SOURCE
    7. SET THE SURFACE PROPERTIES
       • CALL SET_SURFACE_PROPERTIES
    8. CALL A SOLVER APPROPRIATE TO THE PRESENCE OF GASES AND THE OVERLAP ASSUMED
       • CALL SOLVE_BAND_WITHOUT_GAS or
       • CALL INITIALIZE_FLUX and
       • ( CALL SOLVE_BAND_ONE_GAS or
       • CALL SOLVE_BAND_RANDOM_OVERLAP or
       • CALL SOLVE_BAND_FESFT or
       • CALL SOLVE_BAND_CLR_FESFT or
       • CALL SOLVE_BAND_K_EQV )
    9. INCREMENT THE TOTAL FLUXES
       • CALL AUGMENT_TOTAL_FLUX
    10. INCREMENT THE BAND-DEPENDENT DIAGNOSTICS FOR THE UNIFIED MODEL
        • CALL R2_COUPLE_DIAG
    11. PASS THE CALCULATED FLUXES INTO THE OUTPUT ARRAYS.
        • CALL ASSIGN_FLUX

Called by

Calls

No radiation routines radintlw (core lem code) radintsw (core lem code)

aggregate_cloud (summ) augment_total_flux (summ) assign_flux (summ) calculate_density (summ) collect_surface (summ) diff_planck_source grey_extinction (summ) initialize_flux overlap_mix_maximum overlap_mix_random overlap_triple r2_couple_diag r2_init_couple_diag rescale_asymmetry rescale_contiuum scale_absorb set_cloud_geometry set_cloud_pointer set_moist_aerosol_properties set_scattering set_surface_properties split_maximum solve_band_one_gas solve_band_clr_fesft solve_band_fesft solve_band_k_eqv solve_band_random_overlap solve_band_without_gas

 -----------------------------------------------------------------

Return to radiation home page