source: LMDZ4/trunk/libf/phylmd/albedo.F @ 862

Last change on this file since 862 was 830, checked in by lmdzadmin, 17 years ago

Le parametre "fmagic" reglage albedo est lu dans physiq.def (conf_phys.F90)
IM

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 5.5 KB
Line 
1!
2! $Header$
3!
4c
5c
6      SUBROUTINE alboc(rjour,rlat,albedo)
7      USE dimphy
8      IMPLICIT none
9c======================================================================
10c Auteur(s): Z.X. Li (LMD/CNRS) (adaptation du GCM du LMD)
11c Date: le 16 mars 1995
12c Objet: Calculer l'albedo sur l'ocean
13c Methode: Integrer numeriquement l'albedo pendant une journee
14c
15c Arguments;
16c rjour (in,R)  : jour dans l'annee (a compter du 1 janvier)
17c rlat (in,R)   : latitude en degre
18c albedo (out,R): albedo obtenu (de 0 a 1)
19c======================================================================
20cym#include "dimensions.h"
21cym#include "dimphy.h"
22#include "YOMCST.h"
23#include "clesphys.h"
24c
25c fmagic -> clesphys.h/.inc
26c     REAL fmagic ! un facteur magique pour regler l'albedo
27ccc      PARAMETER (fmagic=0.7)
28cccIM => a remplacer 
29c       PARAMETER (fmagic=1.32)
30c       PARAMETER (fmagic=1.0)
31c       PARAMETER (fmagic=0.7)
32      INTEGER npts ! il controle la precision de l'integration
33      PARAMETER (npts=120) ! 120 correspond a l'interval 6 minutes
34c
35      REAL rlat(klon), rjour, albedo(klon)
36      REAL zdist, zlonsun, zpi, zdeclin
37      REAL rmu,alb, srmu, salb, fauxo, aa, bb
38      INTEGER i, k
39cccIM
40      LOGICAL ancien_albedo
41      PARAMETER(ancien_albedo=.FALSE.)
42c     SAVE albedo
43c
44      IF ( ancien_albedo ) THEN
45c
46      zpi = 4. * ATAN(1.)
47c
48c Calculer la longitude vraie de l'orbite terrestre:
49      CALL orbite(rjour,zlonsun,zdist)
50c
51c Calculer la declinaison du soleil (qui varie entre + et - R_incl):
52      zdeclin = ASIN(SIN(zlonsun*zpi/180.0)*SIN(R_incl*zpi/180.0))
53c
54      DO 999 i=1,klon
55      aa = SIN(rlat(i)*zpi/180.0) * SIN(zdeclin)
56      bb = COS(rlat(i)*zpi/180.0) * COS(zdeclin)
57c
58c Midi local (angle du temps = 0.0):
59      rmu = aa + bb * COS(0.0)
60      rmu = MAX(0.0, rmu)
61      fauxo = (1.47-ACOS(rmu))/.15
62      alb = 0.03+0.630/(1.+fauxo*fauxo)
63      srmu = rmu
64      salb = alb * rmu
65c
66c Faire l'integration numerique de midi a minuit (le facteur 2
67c prend en compte l'autre moitie de la journee):
68      DO k = 1, npts
69         rmu = aa + bb * COS(FLOAT(k)/FLOAT(npts)*zpi)
70         rmu = MAX(0.0, rmu)
71         fauxo = (1.47-ACOS(rmu))/.15
72         alb = 0.03+0.630/(1.+fauxo*fauxo)
73         srmu = srmu + rmu * 2.0
74         salb = salb + alb*rmu * 2.0
75      ENDDO
76      IF (srmu .NE. 0.0) THEN
77         albedo(i) = salb / srmu * fmagic
78      ELSE ! nuit polaire (on peut prendre une valeur quelconque)
79         albedo(i) = fmagic
80      ENDIF
81  999 CONTINUE
82c
83c nouvel albedo
84c
85      ELSE
86c
87      zpi = 4. * ATAN(1.)
88c
89c Calculer la longitude vraie de l'orbite terrestre:
90      CALL orbite(rjour,zlonsun,zdist)
91c
92c Calculer la declinaison du soleil (qui varie entre + et - R_incl):
93      zdeclin = ASIN(SIN(zlonsun*zpi/180.0)*SIN(R_incl*zpi/180.0))
94c
95      DO 1999 i=1,klon
96      aa = SIN(rlat(i)*zpi/180.0) * SIN(zdeclin)
97      bb = COS(rlat(i)*zpi/180.0) * COS(zdeclin)
98c
99c Midi local (angle du temps = 0.0):
100      rmu = aa + bb * COS(0.0)
101      rmu = MAX(0.0, rmu)
102cIM cf. PB  alb = 0.058/(rmu + 0.30)
103c     alb = 0.058/(rmu + 0.30) * 1.5
104      alb = 0.058/(rmu + 0.30) * 1.2
105c     alb = 0.058/(rmu + 0.30) * 1.3
106      srmu = rmu
107      salb = alb * rmu
108c
109c Faire l'integration numerique de midi a minuit (le facteur 2
110c prend en compte l'autre moitie de la journee):
111      DO k = 1, npts
112         rmu = aa + bb * COS(FLOAT(k)/FLOAT(npts)*zpi)
113         rmu = MAX(0.0, rmu)
114cIM cf. PB      alb = 0.058/(rmu + 0.30)
115c        alb = 0.058/(rmu + 0.30) * 1.5
116         alb = 0.058/(rmu + 0.30) * 1.2
117c        alb = 0.058/(rmu + 0.30) * 1.3
118         srmu = srmu + rmu * 2.0
119         salb = salb + alb*rmu * 2.0
120      ENDDO
121      IF (srmu .NE. 0.0) THEN
122         albedo(i) = salb / srmu * fmagic
123      ELSE ! nuit polaire (on peut prendre une valeur quelconque)
124         albedo(i) = fmagic
125      ENDIF
1261999  CONTINUE
127      ENDIF
128      RETURN
129      END
130c=====================================================================
131      SUBROUTINE alboc_cd(rmu0,albedo)
132      USE dimphy
133      IMPLICIT none
134c======================================================================
135c Auteur(s): Z.X. Li (LMD/CNRS)
136c date: 19940624
137c Calculer l'albedo sur l'ocean en fonction de l'angle zenithal moyen
138c Formule due a Larson and Barkstrom (1977) Proc. of the symposium
139C on radiation in the atmosphere, 19-28 August 1976, science Press,
140C 1977 pp 451-453, ou These de 3eme cycle de Sylvie Joussaume.
141c
142c Arguments
143c rmu0    (in): cosinus de l'angle solaire zenithal
144c albedo (out): albedo de surface de l'ocean
145c======================================================================
146cym#include "dimensions.h"
147cym#include "dimphy.h"
148#include "clesphys.h"
149      REAL rmu0(klon), albedo(klon)
150c
151c     REAL fmagic ! un facteur magique pour regler l'albedo
152ccc      PARAMETER (fmagic=0.7)
153cccIM => a remplacer 
154c       PARAMETER (fmagic=1.32)
155c       PARAMETER (fmagic=1.0)
156c       PARAMETER (fmagic=0.7)
157c
158      REAL fauxo
159      INTEGER i
160cccIM
161      LOGICAL ancien_albedo
162      PARAMETER(ancien_albedo=.FALSE.)
163c     SAVE albedo
164c
165      IF ( ancien_albedo ) THEN
166c
167      DO i = 1, klon
168c
169         rmu0(i) = MAX(rmu0(i),0.0)
170c
171         fauxo = ( 1.47 - ACOS( rmu0(i) ) )/0.15
172         albedo(i) = fmagic*( .03 + .630/( 1. + fauxo*fauxo))
173         albedo(i) = MAX(MIN(albedo(i),0.60),0.04)
174      ENDDO
175c
176c nouvel albedo
177c
178      ELSE
179c
180      DO i = 1, klon
181         rmu0(i) = MAX(rmu0(i),0.0)
182cIM:orig albedo(i) = 0.058/(rmu0(i) + 0.30)
183         albedo(i) = fmagic * 0.058/(rmu0(i) + 0.30)
184         albedo(i) = MAX(MIN(albedo(i),0.60),0.04)
185      ENDDO
186c
187      ENDIF
188c
189      RETURN
190      END
191c========================================================================
Note: See TracBrowser for help on using the repository browser.