Changeset 2413 for LMDZ5/trunk/libf/phylmd/albedo.F90

Timestamp:

Dec 18, 2015, 8:27:39 PM (9 years ago)

Author:

oboucher

Message:

limiting iflag_rrtm=0 to NSW=2 case (NSW=4 and 6 were possible with slightly different averaging procedures)
limiting iflag_rrtm choices to 0 and 1
limiting iflag_albedo choices to 0 and 1
cleaning up surface ocean albedo parametrisations
fmagic and pmagic corrections are now done in surf_ocean_albedo.F90
fmagic and pmagic corrections also possible for iflag_rrtm=1 now

File:

• LMDZ5/trunk/libf/phylmd/albedo.F90 (modified) (5 diffs)

LMDZ5/trunk/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