Changeset 406 for LMDZ.3.3/branches/rel-LF/libf

Timestamp:

Oct 15, 2002, 4:38:42 PM (22 years ago)

Author:

lmdzadmin

Message:

Albedo ECHAM3 (SB)
IM/LF

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/albedo.F (modified) (4 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/albedo.F

@@ -18 +18 @@
       REAL fmagic ! un facteur magique pour regler l'albedo
 ccc      PARAMETER (fmagic=0.7)
-      PARAMETER (fmagic=1.1)
+cccIM => a remplacer  
+cccIM   PARAMETER (fmagic=1.1)
+c       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
@@ -26 +29 @@
       REAL rmu,alb, srmu, salb, fauxo, aa, bb
       INTEGER i, k
+cccIM
+      LOGICAL ancien_albedo
+      PARAMETER(ancien_albedo=.FALSE.) 
+c     SAVE albedo
+c
+      IF ( ancien_albedo ) THEN
 c
       zpi = 4. * ATAN(1.)
@@ -63 +72 @@
       ENDIF
   999 CONTINUE
+c
+c nouvel albedo 
+c
+      ELSE
+c
+      zpi = 4. * ATAN(1.)
+c
+c Calculer la longitude vraie de l'orbite terrestre:
+      CALL orbite(rjour,zlonsun,zdist)
+c
+c Calculer la declinaison du soleil (qui varie entre + et - R_incl):
+      zdeclin = ASIN(SIN(zlonsun*zpi/180.0)*SIN(R_incl*zpi/180.0))
+c
+      DO 1999 i=1,klon
+      aa = SIN(rlat(i)*zpi/180.0) * SIN(zdeclin)
+      bb = COS(rlat(i)*zpi/180.0) * COS(zdeclin)
+c
+c Midi local (angle du temps = 0.0):
+      rmu = aa + bb * COS(0.0)
+      rmu = MAX(0.0, rmu)
+      alb = 0.058/(rmu + 0.30)
+      srmu = rmu
+      salb = alb * rmu
+c
+c Faire l'integration numerique de midi a minuit (le facteur 2
+c prend en compte l'autre moitie de la journee):
+      DO k = 1, npts
+         rmu = aa + bb * COS(FLOAT(k)/FLOAT(npts)*zpi)
+         rmu = MAX(0.0, rmu)
+         alb = 0.058/(rmu + 0.30)
+         srmu = srmu + rmu * 2.0
+         salb = salb + alb*rmu * 2.0
+      ENDDO
+      IF (srmu .NE. 0.0) THEN
+         albedo(i) = salb / srmu * fmagic
+      ELSE ! nuit polaire (on peut prendre une valeur quelconque)
+         albedo(i) = fmagic
+      ENDIF
+1999  CONTINUE
+      ENDIF
       RETURN
       END
@@ -86 +135 @@
       REAL fmagic ! un facteur magique pour regler l'albedo
 ccc      PARAMETER (fmagic=0.7)
-      PARAMETER (fmagic=1.1)
+cccIM => a remplacer  
+cccIM   PARAMETER (fmagic=1.1)
+c       PARAMETER (fmagic=1.0)
+        PARAMETER (fmagic=0.7) 
 c
       REAL fauxo
       INTEGER i
+cccIM
+      LOGICAL ancien_albedo
+      PARAMETER(ancien_albedo=.FALSE.) 
+c     SAVE albedo
+c
+      IF ( ancien_albedo ) THEN
 c
       DO i = 1, klon
+c
+         rmu0(i) = MAX(rmu0(i),0.0)
+c
          fauxo = ( 1.47 - ACOS( rmu0(i) ) )/0.15
          albedo(i) = fmagic*( .03 + .630/( 1. + fauxo*fauxo))
          albedo(i) = MAX(MIN(albedo(i),0.60),0.04)
       ENDDO
+c
+c nouvel albedo 
+c
+      ELSE
+c
+      DO i = 1, klon
+         rmu0(i) = MAX(rmu0(i),0.0)
+         albedo(i) = 0.058/(rmu0(i) + 0.30)
+         albedo(i) = MAX(MIN(albedo(i),0.60),0.04)
+      ENDDO
+c
+      ENDIF
 c
       RETURN