Changeset 2435 for LMDZ5/branches/testing/libf/phylmd/albedo.F90

Timestamp:

Jan 28, 2016, 5:02:13 PM (9 years ago)

Author:

Laurent Fairhead

Message:

Merged trunk changes r2396:2434 into testing branch

Location:

LMDZ5/branches/testing

Files:

• . (modified) (1 prop)
• libf/phylmd/albedo.F90 (modified) (5 diffs)

LMDZ5/branches/testing/libf/phylmd/albedo.F90

@@ -22 +22 @@
     include "clesphys.h"
 
-    ! fmagic -> clesphys.h/.inc
-    ! REAL fmagic ! un facteur magique pour regler l'albedo
-    ! cc      PARAMETER (fmagic=0.7)
-    ! ccIM => a remplacer
-    ! PARAMETER (fmagic=1.32)
-    ! PARAMETER (fmagic=1.0)
-    ! PARAMETER (fmagic=0.7)
     INTEGER npts ! il controle la precision de l'integration
     PARAMETER (npts=120) ! 120 correspond a l'interval 6 minutes
@@ -74 +67 @@
           END DO
           IF (srmu/=0.0) THEN
-             albedo(i) = salb/srmu*fmagic + pmagic
+             albedo(i) = salb/srmu
           ELSE ! nuit polaire (on peut prendre une valeur quelconque)
-             albedo(i) = fmagic
+             albedo(i) = 1.0
           END IF
        END DO
@@ -119 +112 @@
           END DO
           IF (srmu/=0.0) THEN
-             albedo(i) = salb/srmu*fmagic + pmagic
+             albedo(i) = salb/srmu
           ELSE ! nuit polaire (on peut prendre une valeur quelconque)
-             albedo(i) = fmagic
+             albedo(i) = 1.0
           END IF
        END DO
@@ -146 +139 @@
     real, intent(out):: albedo(klon)
 
-    ! REAL fmagic ! un facteur magique pour regler l'albedo
-    ! cc      PARAMETER (fmagic=0.7)
-    ! ccIM => a remplacer
-    ! PARAMETER (fmagic=1.32)
-    ! PARAMETER (fmagic=1.0)
-    ! PARAMETER (fmagic=0.7)
-
     REAL fauxo
     INTEGER i
@@ -161 +147 @@
        DO i = 1, klon
           fauxo = (1.47-acos(max(rmu0(i), 0.0)))/0.15
-          albedo(i) = fmagic*(.03+.630/(1.+fauxo*fauxo)) + pmagic
+          albedo(i) = 0.03+.630/(1.+fauxo*fauxo)
           albedo(i) = max(min(albedo(i),0.60), 0.04)
        END DO
     ELSE
        DO i = 1, klon
-          albedo(i) = fmagic*0.058/(max(rmu0(i), 0.0)+0.30) + pmagic
+          albedo(i) = 0.058/(max(rmu0(i), 0.0)+0.30)
           albedo(i) = max(min(albedo(i),0.60), 0.04)
        END DO