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 |