- Timestamp:
- Jul 19, 2024, 5:41:58 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_sno_albedo.F
r3980 r5082 218 218 219 219 220 IF (iflag_albcalc .GE.1) THEN ! Albedo calculation according to Kokhanovsky and Zege 2004220 IF (iflag_albcalc >= 1) THEN ! Albedo calculation according to Kokhanovsky and Zege 2004 221 221 222 222 dalbed = 0.0 … … 752 752 alb30 = exp(-(alpha**0.5)* 3.0 / 7.0 * (1.0 + 20 * cos30)) 753 753 754 IF ((wavelengths(i) .GE.lambdamin).AND.755 . (wavelengths(i) .LT.lambdamax)) THEN754 IF ((wavelengths(i)>=lambdamin).AND. 755 . (wavelengths(i)<lambdamax)) THEN 756 756 albint=albint+R*Pas ! rectangle integration -> can be improved with trapezintegration 757 757 norm=norm+Pas
Note: See TracChangeset
for help on using the changeset viewer.