Ignore:
Timestamp:
Mar 24, 2004, 2:10:56 PM (21 years ago)
Author:
lmdzadmin
Message:

Remplace fmagic*1.2 par fmagic dans alboc_cd
IM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/branches/rel-LF/libf/phylmd/albedo.F

    r493 r512  
    2222ccc      PARAMETER (fmagic=0.7)
    2323cccIM => a remplacer 
    24         PARAMETER (fmagic=1.2)
    25 c       PARAMETER (fmagic=1.1)
     24        PARAMETER (fmagic=1.32)
    2625c       PARAMETER (fmagic=1.0)
    2726c       PARAMETER (fmagic=0.7)
     
    9998c     alb = 0.058/(rmu + 0.30) * 1.5
    10099      alb = 0.058/(rmu + 0.30) * 1.2
     100c     alb = 0.058/(rmu + 0.30) * 1.3
    101101      srmu = rmu
    102102      salb = alb * rmu
     
    110110c        alb = 0.058/(rmu + 0.30) * 1.5
    111111         alb = 0.058/(rmu + 0.30) * 1.2
     112c        alb = 0.058/(rmu + 0.30) * 1.3
    112113         srmu = srmu + rmu * 2.0
    113114         salb = salb + alb*rmu * 2.0
     
    144145ccc      PARAMETER (fmagic=0.7)
    145146cccIM => a remplacer 
    146         PARAMETER (fmagic=1.2)
    147 c       PARAMETER (fmagic=1.1)
     147        PARAMETER (fmagic=1.32)
    148148c       PARAMETER (fmagic=1.0)
    149149c       PARAMETER (fmagic=0.7)
     
    173173      DO i = 1, klon
    174174         rmu0(i) = MAX(rmu0(i),0.0)
    175 cIM cf. PB      albedo(i) = 0.058/(rmu0(i) + 0.30)
    176 c        albedo(i) = 0.058/(rmu0(i) + 0.30) * 1.5
    177          albedo(i) = fmagic*0.058/(rmu0(i) + 0.30) * 1.2
     175cIM:orig albedo(i) = 0.058/(rmu0(i) + 0.30)
     176         albedo(i) = fmagic * 0.058/(rmu0(i) + 0.30)
    178177         albedo(i) = MAX(MIN(albedo(i),0.60),0.04)
    179178      ENDDO
Note: See TracChangeset for help on using the changeset viewer.