Changeset 512 for LMDZ.3.3/branches
- Timestamp:
- Mar 24, 2004, 2:10:56 PM (21 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/branches/rel-LF/libf/phylmd/albedo.F
r493 r512 22 22 ccc PARAMETER (fmagic=0.7) 23 23 cccIM => a remplacer 24 PARAMETER (fmagic=1.2) 25 c PARAMETER (fmagic=1.1) 24 PARAMETER (fmagic=1.32) 26 25 c PARAMETER (fmagic=1.0) 27 26 c PARAMETER (fmagic=0.7) … … 99 98 c alb = 0.058/(rmu + 0.30) * 1.5 100 99 alb = 0.058/(rmu + 0.30) * 1.2 100 c alb = 0.058/(rmu + 0.30) * 1.3 101 101 srmu = rmu 102 102 salb = alb * rmu … … 110 110 c alb = 0.058/(rmu + 0.30) * 1.5 111 111 alb = 0.058/(rmu + 0.30) * 1.2 112 c alb = 0.058/(rmu + 0.30) * 1.3 112 113 srmu = srmu + rmu * 2.0 113 114 salb = salb + alb*rmu * 2.0 … … 144 145 ccc PARAMETER (fmagic=0.7) 145 146 cccIM => a remplacer 146 PARAMETER (fmagic=1.2) 147 c PARAMETER (fmagic=1.1) 147 PARAMETER (fmagic=1.32) 148 148 c PARAMETER (fmagic=1.0) 149 149 c PARAMETER (fmagic=0.7) … … 173 173 DO i = 1, klon 174 174 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 175 cIM:orig albedo(i) = 0.058/(rmu0(i) + 0.30) 176 albedo(i) = fmagic * 0.058/(rmu0(i) + 0.30) 178 177 albedo(i) = MAX(MIN(albedo(i),0.60),0.04) 179 178 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.