Changeset 2697 for LMDZ5/trunk/libf/phylmd/ocean_albedo.F90

Timestamp:

Nov 2, 2016, 10:45:49 PM (8 years ago)

Author:

oboucher

Message:

Correcting declaration of ZUE and ZUE2 and one vector calculation

File:

• LMDZ5/trunk/libf/phylmd/ocean_albedo.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/ocean_albedo.F90

@@ -70 +70 @@
 INTEGER                         :: JWL, INU    ! indexes
 REAL                            :: ZWL         ! input parameter: wavelength and diffuse/direct fraction of light
-REAL:: ZCHLABS, ZAW, ZBW, ZREFM, ZYLMD         ! scalar computation variables
+REAL:: ZCHLABS, ZAW, ZBW, ZREFM, ZYLMD, ZUE, ZUE2 ! scalar computation variables
 !
 REAL, DIMENSION(klon) :: ZAP, ZXX2, ZR00, ZRR0, ZRRR               ! computation variables
-REAL, DIMENSION(klon) :: ZR22, ZUE, ZUE2, ZR11DF, ZALBT            ! computation variables
+REAL, DIMENSION(klon) :: ZR22, ZR11DF                              ! computation variables
 REAL, DIMENSION(klon) :: ZBBP, ZNU, ZHB                            ! computation variables
 REAL, DIMENSION(klon) :: ZR11, ZRW, ZRWDF, ZRDF                    ! 4 components of the OSA
@@ -156 +156 @@
   WHERE ( ZCHL(1:knon) > 0.02 )
     ZNU(:)=MIN(0.0,0.5*(ZWORK3(:)-0.3))
-    ZBBP(:)=(0.002+0.01*(0.5-0.25*ZWORK3(:))*(ZWL/550.)**ZNU)*ZWORK2(:)
+    ZBBP(:)=(0.002+0.01*(0.5-0.25*ZWORK3(:))*(ZWL/550.)**ZNU(:))*ZWORK2(:)
   ELSEWHERE
     ZBBP(:)=0.019*(550./ZWL)*ZWORK2(:)       !ZBBPf=0.0113 at chl<=0.02