Bug correction in SNOW3L

src/SURFEX/snow3L_isba.F90

@@ -352,10 +352,9 @@ IF (SIZE(DMK%XSNOWDEND)>0) THEN
   DMK%XIMPUR_CONC (:,:,:) = XUNDEF 
 ENDIF
 !
-
-IF (SIZE(DMK%XDIFF_RATIO)>1) THEN
-  DMK%XSPEC_ALB(:,:) = XUNDEF
-  DMK%XDIFF_RATIO(:,:) = XUNDEF
+IF (PEK%TSNOW%SCHEME=='CRO') THEN
+   DMK%XSPEC_ALB(:,:) = XUNDEF
+   DMK%XDIFF_RATIO(:,:) = XUNDEF
 ENDIF
 !
 DEK%XSNDRIFT(:)    = 0.0
@@ -367,7 +366,10 @@ PFLSN_COR(:)   = 0.0
 PTHRUFAL(:)    = 0.0
 PEVAPCOR(:)    = 0.0
 PQS(:)         = XUNDEF
-IF (SIZE(DMK%XSYTMASS)>1) DMK%XSYTMASS(:)    = 0.0
+!
+IF (PEK%TSNOW%SCHEME=='CRO' .AND. IO%LSNOWSYTRON) THEN
+   DMK%XSYTMASS(:) = 0.0
+ENDIF
 !
 ZSNOW(:)       = 0.0
 ZSNOWD(:)      = 0.0
@@ -1582,15 +1584,15 @@ IF (PEK%TSNOW%SCHEME=='CRO') THEN
     ENDDO
   ENDDO      
 !
-  IF (SIZE(DMK%XDIFF_RATIO)>1) THEN
-     DO JWRK=1,SIZE(P_DIR_SW,2)
+IF (PEK%TSNOW%SCHEME=='CRO' .AND. GCOMPUTECRODIAG) THEN
+   DO JWRK=1,SIZE(P_DIR_SW,2)
       DO JJ=1,KSIZE1
-        JI = KMASK(JJ)
-        DMK%XDIFF_RATIO(JI,JWRK)=ZP_DIFF_RATIO(JJ,JWRK)
-        DMK%XSPEC_ALB(JI,JWRK)=ZP_SPEC_ALB(JJ,JWRK) 
+         JI = KMASK(JJ)
+         DMK%XDIFF_RATIO(JI,JWRK)=ZP_DIFF_RATIO(JJ,JWRK)
+         DMK%XSPEC_ALB(JI,JWRK)=ZP_SPEC_ALB(JJ,JWRK) 
       ENDDO
-    ENDDO
-  ENDIF
+   ENDDO
+ENDIF
 !
   IF (GCOMPUTECRODIAG)THEN
   ! This is equivalent to test the value of DGMI%LPROSNOW which does not enter in ISBA