Index: LMDZ6/trunk/libf/phylmd/surf_landice_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/surf_landice_mod.F90	(revision 4537)
+++ LMDZ6/trunk/libf/phylmd/surf_landice_mod.F90	(revision 4538)
@@ -425,5 +425,5 @@
               rhod=p1lay(i)/RD/temp_air(i)
               fluxbs(i)=rhod*ws1(i)*cdragm(i)*zeta_bs*(AcoefQBS(i)-qsalt(i)) &
-                       / (1.-rhod*ws1(i)*zeta_bs*BcoefQBS(i)*dtime)
+                       / (1.-rhod*ws1(i)*zeta_bs*cdragm(i)*BcoefQBS(i)*dtime)
               !fluxbs(i)= zeta_bs*rhod*ws1(i)*cdragm(i)*(qbs1(i)-qsalt(i))
        enddo
Index: LMDZ6/trunk/libf/phylmdiso/surf_landice_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmdiso/surf_landice_mod.F90	(revision 4537)
+++ LMDZ6/trunk/libf/phylmdiso/surf_landice_mod.F90	(revision 4538)
@@ -514,5 +514,5 @@
               rhod=p1lay(i)/RD/temp_air(i)
               fluxbs(i)=rhod*ws1(i)*cdragm(i)*zeta_bs*(AcoefQBS(i)-qsalt(i)) &
-                       / (1.-rhod*ws1(i)*zeta_bs*BcoefQBS(i)*dtime)
+                       / (1.-rhod*ws1(i)*zeta_bs*cdragm(i)*BcoefQBS(i)*dtime)
               !fluxbs(i)= zeta_bs*rhod*ws1(i)*cdragm(i)*(qbs1(i)-qsalt(i))
        enddo