Scheduled Downtime
On Tuesday 24 October 2023 @ 5pm MT the forums will be in read only mode in preparation for the downtime. On Wednesday 25 October 2023 @ 5am MT, this website will be down for maintenance and expected to return online later in the morning.
Normal Operations
The forums are back online with normal operations. If you notice any issues or errors related to the forums, please reach out to help@ucar.edu

To modify the spectrum, do I need to change frcsol in radconstants.F90?

Hi,I am using WACCM (compset: F_1850_WACCM, CESM version: 1.2.2) to simulate the atmospheric chemistry on exoplanets around M dwarf stars. To modify the spectrum, there are several files in the namelist of the atmospheric component (solar_data_file, solar_parms_file, etc) to be rewritten. But not long a ago I found another file in the CESM source code, namely, $CESM_ROOT/models/atm/cam/src/physics/cam/radconstants.F90, which contains a variable named frcsol:   ! minimum wavelength of band in micrometers  real(r8), parameter :: wavmin(nswbands) = &       (/  0.200_r8, 0.245_r8,   0.265_r8,   0.275_r8,   0.285_r8, &           0.295_r8, 0.305_r8,   0.350_r8,   0.640_r8,   0.700_r8,   0.700_r8, &           0.700_r8, 0.700_r8,   0.700_r8,   0.700_r8,   0.700_r8, &           2.630_r8, 4.160_r8,   4.160_r8/)    ! maximum wavelength of band in micrometers  real(r8), parameter :: wavmax(nswbands) = &       (/  0.245_r8, 0.265_r8,   0.275_r8,   0.285_r8,   0.295_r8, &           0.305_r8, 0.350_r8,   0.640_r8,   0.700_r8,   5.000_r8,   5.000_r8, &           5.000_r8, 5.000_r8,   5.000_r8,   5.000_r8,   5.000_r8, &           2.860_r8, 4.550_r8,   4.550_r8/)    ! Fraction of solar flux in each stated spectral interval  real(r8), public, parameter :: frcsol(nswbands) = &                             (
 

marsh

Member
I think this is more of a CAM question than a WACCM question. My guess is that frcsol is used when only TSI is specified and some assumption is necessary on how that solar irradiance is apportioned across the spectrum. In WACCM we provide the spectral flux (SSI) and do not use the TSI. In WACCM4 the SSI is used in CAMRT to scale heating rates below 65 km. It does that by ratioing the speficied flux within a CAMRT spectral band to a reference spectrum and applying that to the heating rate in that band. 
 
Thank you for your reply.I did some tests about this and reviewed the source code. In WACCM frcsol acts like a "baseline" for the irradiance from the Sun. Modifying the tsi and ssi in solar_data_file is definitely enough to change the spectrum.However, in other compsets this may not be true, because there isn't always a "solar_data_file" specifying the energy flux in every band. Modifying "frcsol" may still be the only way to change the spectrum.Yangcheng
 

mmills

CSEG and Liaisons
Staff member
There is a namelist variable that is set in WACCM compsets to use spectral heating rates: solar_htng_spctrl_scl          = .true.When this flag is set, TSI is never used.
 
Top