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

How to find the code to calculate top of atmosphere solar radiation at each grid cell

Eric

Eric
Member
Hello,

According to the CAM6.0 scientific support document, the top of atmosphere solar radiation at each grid cell is calculated as:
1619837577914.png
Could you help me find the location of this calculation? Which .F90 file should I look into? Thank you!

Best
Eric
 

andrew

Member
Try looking for "solin" in ./src/physics/rrtmg/radsw.F90


solin(i) = sum(sfac(:)*solar_band_irrad(:)) * eccf * coszrs(i)

I think is going to be the solar broadband flux across bands and solar zenith angle, etc.
 

Eric

Eric
Member
Hi Andrew,

Thanks for the help! I noticed that in the radsw.F90, under "! Initialize output fields:", each output field is expressed as an array. For example, fsds(1:ncol) = 0.0_r8. I wonder what does ncol mean here? What is the physical meaning of ncol?

In my understanding, at each time step, there should be a spatial loop allowing the radiation calculation for each grid cell. I guess the radiation calculation in the radsw.F90 and radiation.F90 is inside only one grid cell. Could you help me find the do loop for the spatial dimension?

Thank you!

Best
Eric
 
Top