Index: trunk/LMDZ.MARS/libf/aeronomars/moldiff_red.F90
===================================================================
--- trunk/LMDZ.MARS/libf/aeronomars/moldiff_red.F90	(revision 2156)
+++ trunk/LMDZ.MARS/libf/aeronomars/moldiff_red.F90	(revision 2157)
@@ -1173,5 +1173,6 @@
 	endif
 	enddo
-	D(l)=1d0/interm
+!	D(l)=1d0/interm
+        D(l)=(1D0-Nk(l,nn)/N(l))/interm
 	enddo
 	END