Changeset 2709 for LMDZ5/trunk/libf/phylmd

Timestamp:

Nov 22, 2016, 10:05:03 AM (8 years ago)

Author:

Laurent Fairhead

Message:

knon/klon does not seem to agree with the WHERE construct and crashes MPI
parallelisation.
LF

File:

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

LMDZ5/trunk/libf/phylmd/ocean_albedo.F90

@@ -69 +69 @@
 !
 INTEGER                         :: JWL, INU    ! indexes
+INTEGER                         :: JI
 REAL                            :: ZWL         ! input parameter: wavelength and diffuse/direct fraction of light
 REAL:: ZCHLABS, ZAW, ZBW, ZREFM, ZYLMD, ZUE, ZUE2 ! scalar computation variables
@@ -154 +155 @@
   ZAP(1:knon) = 0.06*ZCHLABS*ZWORK1(1:knon) +0.2*(XAW440+0.06*ZWORK1(1:knon))*ZYLMD
    
-  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(:)
-  ELSEWHERE
-    ZBBP(:)=0.019*(550./ZWL)*ZWORK2(:)       !ZBBPf=0.0113 at chl<=0.02
-  ENDWHERE
-   
+!!  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(:)
+!!  ELSEWHERE
+!!    ZBBP(:)=0.019*(550./ZWL)*ZWORK2(:)       !ZBBPf=0.0113 at chl<=0.02
+!!  ENDWHERE
+
+    do JI = 1, knon
+      IF (ZCHL(JI) > 0.02) THEN
+        ZNU(JI)=MIN(0.0,0.5*(ZWORK3(JI)-0.3))
+        ZBBP(JI)=(0.002+0.01*(0.5-0.25*ZWORK3(JI))*(ZWL/550.)**ZNU(JI)) &
+                  *ZWORK2(JI)
+      ELSE
+        ZBBP(JI)=0.019*(550./ZWL)*ZWORK2(JI)       !ZBBPf=0.0113 at chl<=0.02 
+      ENDIF
+    ENDDO
+
   ! Morel-Gentili(1991), Eq (12)
   ! ZHB=h/(h+2*ZBBPf*(1.-h))