Parent Directory
| 
Revision Log
Removed argument dtphys of physiq. Use it directly from comconst in physiq instead. Donwgraded variables eignfnu, eignfnv of module inifgn_m to dummy arguments of SUBROUTINE inifgn. They were not used elsewhere than in the calling procedure inifilr. Renamed argument dv of inifgn to eignval_v. Made alboc and alboc_cd independent of the size of arguments. Now we can call them only at indices knindex in interfsurf_hq, where we need them. Fixed a bug in alboc_cd: rmu0 was modified, and the corresponding actual argument in interfsurf_hq is an intent(in) argument of interfsurf_hq. Variables of size knon instead of klon in interfsur_lim and interfsurf_hq. Removed argument alb_new of interfsurf_hq because it was the same than alblw. Simplified test on cycle_diurne, following LMDZ. Moved tests on nbapp_rad from physiq to read_clesphys2. No need for separate counter itaprad, we can use itap. Define lmt_pas and radpas from integer input parameters instead of real-type computed values.
| 1 | module alboc_cd_m |
| 2 | |
| 3 | IMPLICIT NONE |
| 4 | |
| 5 | contains |
| 6 | |
| 7 | SUBROUTINE alboc_cd(rmu0, albedo) |
| 8 | |
| 9 | ! From LMDZ4/libf/phylmd/albedo.F, version 1.2 2005/02/07 15:00:52 |
| 10 | |
| 11 | ! Author: Z. X. Li (LMD/CNRS) |
| 12 | ! Date: 1994/06/24 |
| 13 | |
| 14 | ! Calculer l'alb\'edo sur l'oc\'ean en fonction de l'angle |
| 15 | ! z\'enithal moyen. Formule due \`a Larson and Barkstrom, |
| 16 | ! Proceedings of the symposium on radiation in the atmosphere, |
| 17 | ! 19-28 August 1976, science Press, 1977, pages 451-453, ou |
| 18 | ! th\`ese de 3\`eme cycle de Sylvie Joussaume. |
| 19 | |
| 20 | REAL, intent(in):: rmu0(:) ! cosinus de l'angle solaire z\'enithal |
| 21 | real, intent(out):: albedo(:) ! alb\'edo de surface de l'oc\'ean |
| 22 | |
| 23 | ! Local: |
| 24 | REAL, PARAMETER:: fmagic = 1. ! facteur magique pour r\'egler l'alb\'edo |
| 25 | |
| 26 | !---------------------------------------------------------- |
| 27 | |
| 28 | albedo = max(min(fmagic * 0.058 / (max(rmu0, 0.) + 0.3), 0.6), 0.04) |
| 29 | |
| 30 | END SUBROUTINE alboc_cd |
| 31 | |
| 32 | end module alboc_cd_m |
| ViewVC Help | |
| Powered by ViewVC 1.1.21 |