-
Notifications
You must be signed in to change notification settings - Fork 41
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
add simplified albedo calculation model (#1283)
* add simplified albedo calculation model * fix the whitespace check error * add cam5 data to the repo * fix the double albedo definition error
- Loading branch information
Showing
6 changed files
with
36 additions
and
67 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,22 @@ | ||
/*------------------------------------------------------------ | ||
Computes surface albedos | ||
For more details , see Briegleb, Bruce P., 1992: Delta-Eddington | ||
Approximation for Solar Radiation in the NCAR Community Climate Model, | ||
Journal of Geophysical Research, Vol 97, D7, pp7603-7612). | ||
NOTE: this is the simplest formula for albedo calculation, we have to | ||
update to the latest implementation from microphysics output. | ||
*/ | ||
#ifndef ERF_ALBEDO_H_ | ||
#define ERF_ALBEDO_H_ | ||
void set_albedo(const real1d& coszrs, real2d& albedo_dir, real2d& albedo_dif) | ||
{ | ||
// Albedos for land type I (Briegleb) | ||
auto nswbands = albedo_dir.extent(0); | ||
auto ncol = albedo_dif.extent(1); | ||
yakl::c::parallel_for(yakl::c::Bounds<2>(nswbands,ncol), YAKL_LAMBDA (int ibnd, int icol) { | ||
albedo_dir(ibnd, icol) = 1.4 * 0.24 / ( 1. + 0.8 * coszrs(icol)); | ||
albedo_dif(ibnd, icol) = 1.2 * 0.24; | ||
}); | ||
} | ||
#endif |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Binary file not shown.
Binary file not shown.