Changeset 900 for LMDZ4/trunk/libf/phylmd/albedo.F

Timestamp:

Feb 7, 2008, 3:26:01 PM (17 years ago)

Author:

Laurent Fairhead

Message:

Inclusion de la constante pmagic dans le calcul de l'albedo MAFO/JLD
LF

File:

• LMDZ4/trunk/libf/phylmd/albedo.F (modified) (4 diffs)

LMDZ4/trunk/libf/phylmd/albedo.F

@@ -75 +75 @@
       ENDDO
       IF (srmu .NE. 0.0) THEN
-         albedo(i) = salb / srmu * fmagic
+         albedo(i) = salb / srmu * fmagic+pmagic
       ELSE ! nuit polaire (on peut prendre une valeur quelconque)
          albedo(i) = fmagic
@@ -120 +120 @@
       ENDDO
       IF (srmu .NE. 0.0) THEN
-         albedo(i) = salb / srmu * fmagic
+         albedo(i) = salb / srmu * fmagic+pmagic
       ELSE ! nuit polaire (on peut prendre une valeur quelconque)
          albedo(i) = fmagic
@@ -170 +170 @@
 c
          fauxo = ( 1.47 - ACOS( rmu0(i) ) )/0.15
-         albedo(i) = fmagic*( .03 + .630/( 1. + fauxo*fauxo))
+         albedo(i) = fmagic*( .03 + .630/( 1. + fauxo*fauxo))+pmagic
          albedo(i) = MAX(MIN(albedo(i),0.60),0.04)
       ENDDO
@@ -181 +181 @@
          rmu0(i) = MAX(rmu0(i),0.0)
 cIM:orig albedo(i) = 0.058/(rmu0(i) + 0.30)
-         albedo(i) = fmagic * 0.058/(rmu0(i) + 0.30)
+         albedo(i) = fmagic * 0.058/(rmu0(i) + 0.30)+pmagic
          albedo(i) = MAX(MIN(albedo(i),0.60),0.04)
       ENDDO