/[lmdze]/trunk/phylmd/Interface_surf/alboc_cd.f
ViewVC logotype

Annotation of /trunk/phylmd/Interface_surf/alboc_cd.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 154 - (hide annotations)
Tue Jul 7 17:49:23 2015 UTC (8 years, 11 months ago) by guez
Original Path: trunk/Sources/phylmd/Interface_surf/alboc_cd.f
File size: 941 byte(s) 
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 guez 154 module alboc_cd_m
2    
3 guez 117 IMPLICIT NONE
4    
5 guez 154 contains
6 guez 117
7 guez 154 SUBROUTINE alboc_cd(rmu0, albedo)
8 guez 117
9 guez 154 ! From LMDZ4/libf/phylmd/albedo.F, version 1.2 2005/02/07 15:00:52
10 guez 117
11 guez 154 ! Author: Z. X. Li (LMD/CNRS)
12     ! Date: 1994/06/24
13 guez 117
14 guez 154 ! 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 guez 117
20 guez 154 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 guez 117
23 guez 154 ! Local:
24     REAL, PARAMETER:: fmagic = 1. ! facteur magique pour r\'egler l'alb\'edo
25 guez 117
26 guez 154 !----------------------------------------------------------
27 guez 117
28 guez 154 albedo = max(min(fmagic * 0.058 / (max(rmu0, 0.) + 0.3), 0.6), 0.04)
29 guez 117
30 guez 154 END SUBROUTINE alboc_cd
31 guez 117
32 guez 154 end module alboc_cd_m
  ViewVC Help
Powered by ViewVC 1.1.21