Index | Rundeck: P1SoM40 | Created: Wed May 8 03:30:18 EDT 2024 |
RAD_COM | File: RAD_COM.f |
Author : Original Development Team
Version:
Subroutines:
Global Variables:
- im : used from resolution
- jm : used from resolution
- lm : used from resolution
- lm_req : used from atm_com
- s0 : used from radpar
- solar 'constant' needs to be saved between calls to radiation
- itrmax : used from radpar
- s0 : used from radpar
- solar 'constant' needs to be saved between calls to radiation
- abstractorbit : used from abstractorbit_mod
- nrad :
integer
- : DT_Rad = NRad*DTsrc
Initial Value= 5
Used by: | :atm_phase1 | :init_rad | :radia |
- modrd :
integer
- : if MODRD=0 do radiation, else skip
Used by: | :atm_phase1 | :init_rad | :radia |
- variable_orb_par :
integer
- 1 if orbital parameters are time dependent
Initial Value= -2
Used by: | :daily_orbit |
- orb_par_year_bp :
integer
- = offset from model_year or 1950 (fixed case)
Initial Value= 0
Used by: | :daily_orbit |
- dimrad_sv :
integer, parameter
- dimension sum of input fields saved for radia_only runs
Initial Value= IM*JM*(7*LM+3*LM_REQ+24)
- rqt :
real*8, allocatable, dimension(:, :, :)
- Radiative equilibrium temperatures above model top
Used by: | :alloc_rad_com | :diaga | :init_rad | :radia |
- tchg :
real*8, allocatable, dimension(:, :, :)
- Total temperature change in adjusted forcing runs
Used by: | :alloc_rad_com | :radia |
- srhr :
real*8, allocatable, dimension(:, :, :)
- (1->LM) Solar raditive heating rate (W/m^2) (short wave)
Used by: | :alloc_rad_com | :radia | SUBDAILY:get_subdd |
- trhr :
real*8, allocatable, dimension(:, :, :)
- (1->LM) Thermal raditive heating rate (W/m^2) (long wave)
Used by: | :alloc_rad_com | :atm_phase1_exports | :radia | :surface | :surface_diag1 | :surface_diag1a | SUBDAILY:get_subdd |
- srhr :
real*8, allocatable, dimension(:, :, :)
- (1->LM) Solar raditive heating rate (W/m^2) (short wave)
Used by: | :alloc_rad_com | :radia | SUBDAILY:get_subdd |
- trhr :
real*8, allocatable, dimension(:, :, :)
- (1->LM) Thermal raditive heating rate (W/m^2) (long wave)
Used by: | :alloc_rad_com | :atm_phase1_exports | :radia | :surface | :surface_diag1 | :surface_diag1a | SUBDAILY:get_subdd |
- trsurf :
real*8, allocatable, dimension(:, :, :)
- upward thermal radiation at the surface from rad step W/m2
Used by: | :alloc_rad_com | :atm_phase1_exports | :radia | :reset_surf_fluxes | :surface | :surface_diag1a |
- fsf :
real*8, allocatable, dimension(:, :, :)
- Solar Forcing over each type (W/m^2)
Used by: | :alloc_rad_com | :atm_phase1_exports | :radia | :reset_surf_fluxes | :surface |
- fsrdir :
real*8, allocatable, dimension(:, :)
- Solar incident at surface, direct fraction (1)
Used by: | :alloc_rad_com | :radia | SOIL_DRV:earth |
- dirvis :
real*8, allocatable, dimension(:, :)
- Direct beam solar incident at surface (W/m^2)
Used by: | :alloc_rad_com | :radia |
- srvissurf :
real*8, allocatable, dimension(:, :)
- Incident solar direct+diffuse visible at surface (W/m^2)
Used by: | :alloc_rad_com | :radia | SOIL_DRV:earth |
- srdn :
real*8, allocatable, dimension(:, :)
- Total incident solar at surface (W/m^2)
Used by: | :alloc_rad_com | :radia | SUBDAILY:get_subdd |
- fsrdif :
real*8, allocatable, dimension(:, :)
- diffuse visible incident solar at surface
Used by: | :alloc_rad_com | :radia |
- dirnir :
real*8, allocatable, dimension(:, :)
- direct nir incident solar at surface
Used by: | :alloc_rad_com | :radia |
- difnir :
real*8, allocatable, dimension(:, :)
- diffuse nir incident solar at surface
Used by: | :alloc_rad_com | :radia |
- srnflb_save :
real*8, allocatable, dimension(:, :, :)
- Net solar radiation (W/m^2)
Used by: | :ahourly | :alloc_rad_com |
- trnflb_save :
real*8, allocatable, dimension(:, :, :)
- Net thermal radiation (W/m^2)
Used by: | :ahourly | :alloc_rad_com |
- tausumi :
real*8, dimension(:, :), allocatable
- column-sum water,ice cloud opt. depths (for diags)
Used by: | :alloc_rad_com | :radia | SUBDAILY:get_subdd |
- tausumw :
real*8, dimension(:, :), allocatable
- column-sum water,ice cloud opt. depths (for diags)
Used by: | :alloc_rad_com | :radia | SUBDAILY:get_subdd |
- ttausv_save :
real*8, allocatable, dimension(:, :, :, :)
- Tracer optical thickness
Used by: | :ahourly | :alloc_rad_com |
- ttausv_cs_save :
real*8, allocatable, dimension(:, :, :, :)
- Tracer optical thickness clear sky
Used by: | :ahourly | :alloc_rad_com |
- cfrac :
real*8, allocatable, dimension(:, :)
- Total cloud fraction as seen be radiation
Used by: | :alloc_rad_com | :radia | SUBDAILY:get_subdd |
- rcld :
real*8, allocatable, dimension(:, :, :)
- Total cloud optical depth as seen be radiation
Used by: | :alloc_rad_com | :radia |
- chem_tracer_save :
real*8, allocatable, dimension(:, :, :, :)
- 3D O3, CH4 saved elsewhere for use in radiation
Used by: | :alloc_rad_com | :radia |
- rad_to_chem :
real*8, allocatable, dimension(:, :, :, :)
- save 3D quantities from radiation code for use in
Used by: | :alloc_rad_com | :radia |
- kliq :
integer, allocatable, dimension(:, :, :, :)
- Flag indicating dry(0)/wet(1) atmosphere (memory feature)
Used by: | :alloc_rad_com | :radia |
- rad_to_file :
real*8, allocatable, dimension(:, :, :, :)
Used by: | :alloc_rad_com |
- ikliq :
integer
- 0,1,-1 initialize kliq as dry,equil,current model state
Initial Value= -1
- rhfix :
real*8
- const.rel.humidity passed to radiation for aeros. tests
Initial Value= -1.
Used by: | :init_rad | :radia |
- dalbsnx :
real*8
- global coeff for snow alb change by black carbon depos
Initial Value= 0.
Used by: | :daily_rad | :init_rad | :radia |
- albsn_yr :
integer
- year of blk carb depos used for snow alb. reduction
Initial Value= 1951
Used by: | :daily_rad | :init_rad |
- cc_cdncx :
real*8
- scaling factor relating cld cvr change and CDNC change
Initial Value= .0000d0
Used by: | :init_rad | :radia |
- oc_cdncx :
- scaling factor relating cld opt depth and CDNC change
- vcdnc :
real*8, parameter, dimension(7)
- pressure,vertical profile for cld.cvr change
Initial Value= (/ .35d0, .20d0, .10d0, .17d0, .10d0, .08d0, 0.d0/)
Used by: | :init_rad |
- pcdnc :
real*8, parameter, dimension(7)
- pressure,vertical profile for cld.cvr change
Initial Value= (/984.d0, 964.d0, 934.d0, 884.d0, 810.d0, 710.d0, 550.d0/)
Used by: | :init_rad |
- od_cdncx :
real*8
- Initial Value
= .0000d0
Used by: | :init_rad | :radia |
- cdncl :
real*8, dimension(LM)
- = vcdnc interpolated to current vertical resolution
Used by: | :init_rad | :radia |
- cosz1 :
real*8, allocatable, dimension(:, :)
- Mean Solar Zenith angle for curr. physics(not rad) time step
Used by: | :ahourly | :alloc_rad_com | :atm_phase1 | :atm_phase1_exports | :condse | :init_rad | :radia | :run_radonly | :surface | SUBDAILY:get_subdd |
- cosz_day :
real*8, allocatable, dimension(:, :)
- Mean Solar Zenith angle for current day
Used by: | :alloc_rad_com | :daily_orbit |
- sunset :
real*8, allocatable, dimension(:, :)
- Time of sunset for current day (radians from local noon)
Used by: | :alloc_rad_com | :daily_orbit |
- s0x :
real*8
- solar constant multiplication factor
Initial Value= 1.
Used by: | :init_rad | :radia |
- s0_day :
integer
- obs.date of solar constant (if 0: time var)
Initial Value= 182
Used by: | :init_rad |
- s0_yr :
integer
- obs.date of solar constant (if 0: time var)
Initial Value= 1951
Used by: | :init_rad |
- co2x :
real*8
- scaling factor for CO2 with respect to value in "GHG" file
Initial Value= 1.
Used by: | :daily_rad | :init_rad | :radia | SOIL_DRV:earth |
- n2ox :
real*8
- scaling factor for N2O with respect to value in "GHG" file
Initial Value= 1.
Used by: | :daily_rad | :init_rad | :radia |
- ch4x :
real*8
- scaling factor for CH4 with respect to value in "GHG" file
Initial Value= 1.
Used by: | :daily_rad | :init_rad | :radia |
- cfc11x :
real*8
- scaling factor for CFC11 with respect to value in "GHG" file
Initial Value= 1.
Used by: | :daily_rad | :init_rad | :radia |
- cfc12x :
real*8
- scaling factor for CFC12 with respect to value in "GHG" file
Initial Value= 1.
Used by: | :daily_rad | :init_rad | :radia |
- xghgx :
real*8
- scaling factor for XGHG with respect to value in "GHG" file
Initial Value= 1.
Used by: | :daily_rad | :init_rad | :radia |
- o2x :
real*8
- scaling factor for O2 with respect to value in "GHG" file
Initial Value= 1.
Used by: | :daily_rad | :init_rad |
- no2x :
real*8
- scaling factor for NO2 with respect to value in "GHG" file
Initial Value= 1.
Used by: | :daily_rad | :init_rad |
- yghgx :
real*8
- scaling factor for YGHG with respect to value in "GHG" file
Initial Value= 2.
Used by: | :daily_rad | :init_rad |
- so2x :
real*8
- scaling factor for SO2 with respect to value in "GHG" file
Initial Value= 0.
Used by: | :daily_rad | :init_rad |
- ch4x_radoverchem :
real*8
- scaling factor for chemistry
Initial Value= 1.d0
Used by: | :init_rad | :radia |
- n2cx :
real*8
- Initial Value
= 1.
Used by: | :daily_rad | :init_rad |
- ghg_day :
integer
- obs.date of well-mixed GHgases (if 0: time var)
Initial Value= 182
Used by: | :init_rad |
- ghg_yr :
integer
- obs.date of well-mixed GHgases (if 0: time var)
Initial Value= 1951
Used by: | :daily_ch4ox | :daily_rad | :ghghst | :init_rad | :radia |
- ref_mult :
real*8
- Initial Value
= 1.
Used by: | :init_rad |
- volc_day :
integer
- obs.date of Volc.Aerosols (if 0: time var)
Initial Value= 182
Used by: | :init_rad |
- volc_yr :
integer
- obs.date of Volc.Aerosols (if 0: time var)
Initial Value= 1951
Used by: | :daily_rad | :init_rad |
- aero_yr :
integer
- obs.year of troposph.Aerosols (if 0: use current yr)
Initial Value= 1951
Used by: | :init_rad |
- o3_yr :
integer
- obs.year of Ozone (if 0: use current year)
Initial Value= 1951
Used by: | :daily_rad | :init_rad |
- crops_yr :
integer
- obs.year of crops (if 0: time var, -1: default)
Initial Value= -1
- h2ostratx :
real*8
- strat_water_vapor, cloud, Ozone scaling factor
Initial Value= 1.
Used by: | :daily_rad | :init_rad |
- h2obych4 :
real*8
- if not 0: add CH4 produced H2O into layers 1->LM
Initial Value= 1.
Used by: | :daily_ch4ox | :init_rad |
- cldx :
real*8
- Initial Value
= 1.
Used by: | :init_rad | :radia |
- o3x :
real*8
- Initial Value
= 1.
Used by: | :daily_rad | :init_rad |
- dh2o :
real*8, allocatable, dimension(:, :, :)
- zonal H2O-prod.rate in kg/m^2/ppm_CH4/second in layer L
Used by: | :alloc_rad_com | :daily_ch4ox | :init_rad |
- cosd :
real*8
- orbit related variables computed once a day
Used by: | :daily_orbit |
- sind :
real*8
- orbit related variables computed once a day
Used by: | :daily_orbit |
- rsdist :
real*8
- orbit related variables computed once a day
Used by: | :daily_orbit | :radia |
- alb :
real*8, allocatable, dimension(:, :, :), target
- is SRNFLB(1)/(SRDFLB(1)+1.D-20),PLAVIS,PLANIR,ALBVIS,ALBNIR,
Used by: | :alloc_rad_com | :radia |
- salb :
real*8, pointer, dimension(:, :)
- (1.-broadband surface albedo) - saved in rsf
Used by: | :alloc_rad_com | SUBDAILY:get_subdd |
- rad_interact_aer :
integer
- =1 for radiatively active non-chem tracers
Initial Value= 0
Used by: | :init_rad | :radia |
- clim_interact_chem :
integer
- =1 for radiatively active chem tracers
Initial Value= 0
Used by: | :init_rad | :radia |
- nradfrc :
integer
- sets frequency of inst. rad. forcing calculations
Initial Value= 1
Used by: | :ij_defs | :init_rad | :radia |
- rad_forc_lev :
integer
- = 0 for TOA, 1 for LTROPO (default=0)
Initial Value= 0
Used by: | :init_rad | :radia |
- cloud_rad_forc :
integer
- = 1 for calculation of cloud + aerosol + O3
Initial Value= 0
Used by: | :init_rad | :radia | DIAG_SERIAL:diagij |
- cloud_rad_forc :
integer
- = 1 for calculation of cloud + aerosol + O3
Initial Value= 0
Used by: | :init_rad | :radia | DIAG_SERIAL:diagij |
- aer_rad_forc :
integer
- = 1 for calculation of aerosol radiative forcing
Initial Value= 0
Used by: | :init_rad | :radia | DIAG_SERIAL:diagij |
- cloud_aer_o3_rad_forc :
integer
- Initial Value
= 0
Used by: | :init_rad | :radia |
- co2ppm :
real*8
- Current CO2 level as seen by radiation
Initial Value= 280.
Used by: | :daily_rad | SOIL_DRV:earth |
- ql0 :
- global parts of local arrays (to avoid OMP-copyin)
- plb0 :
real*8, dimension(lm_req)
- global parts of local arrays (to avoid OMP-copyin)
Used by: | :init_rad | :radia |
- ntrix :
integer, dimension(itrmax)
- Indexing array for optional aerosol interaction
Initial Value= 0
Used by: | :init_rad | :radia |
- shl0 :
real*8, dimension(lm_req)
Used by: | :init_rad | :radia |
- wttr :
real*8, dimension(itrmax)
- weighting array for optional aerosol interaction
Initial Value= 1.
Used by: | :init_rad | :radia |
- nrad_clay :
integer
- index of clay in arrays for optional aerosol interaction
Used by: | :init_rad | :radia |
- snoage_def :
integer
- determines how snowage is calculated:
Initial Value= 0
Used by: | :daily_rad | :init_rad |
- snoage :
real*8, allocatable, dimension(:, :, :)
Used by: | :alloc_rad_com | :condse | :daily_rad | :radia | :read_rad_ic | SOIL_DRV:init_land_surface |