Index: /trunk/LMDZ.VENUS/libf/phyvenus/moldiff_red.F90
===================================================================
--- /trunk/LMDZ.VENUS/libf/phyvenus/moldiff_red.F90	(revision 3217)
+++ /trunk/LMDZ.VENUS/libf/phyvenus/moldiff_red.F90	(revision 3218)
@@ -106,63 +106,64 @@
 
 ! Initializations at first call
+
       if (firstcall) then
 
-! Liste des especes qui sont diffuses si elles sont presentes
-
-	ListeDiff(1)='co2'
-	ListeDiff(2)='o'
-	ListeDiff(3)='n2'
-	ListeDiff(4)='ar'
-	ListeDiff(5)='co'
-	ListeDiff(6)='h2'
-	ListeDiff(7)='h'
-	ListeDiff(8)='d2'
-	ListeDiff(9)='hd'
-	ListeDiff(10)='d'
-	ListeDiff(11)='o2'
-	ListeDiff(12)='h2o'
-	ListeDiff(13)='o3'
-	ListeDiff(14)='n'
-	ListeDiff(15)='he'
-        ListeDiff(16)='n2d'
-	i_esc_h=1000
-	i_esc_h2=1000
-
-! On regarde quelles especes diffusables sont presentes
-
-	ncompdiff=0
-	indic_diff(1:nq)=0
-	
-	do nn=1,nq
-
-	do n=1,nb_esp_diff
-	if (ListeDiff(n) .eq. tname(nn)) then
-	indic_diff(nn)=1
-	if (tname(nn) .eq. 'h') i_esc_h=n
-	if (tname(nn) .eq. 'h2') i_esc_h2=n
-	endif
-	enddo
-
-	if (indic_diff(nn) .eq. 1) then
-	print*,'specie ', tname(nn), 'diffused in moldiff_red'
-	ncompdiff=ncompdiff+1
-	endif
-	enddo
-	print*,'number of diffused species:',ncompdiff
-	allocate(gcmind(ncompdiff))
-	allocate(mol_tr(ncompdiff))
-
-!   Store gcm indexes in gcmind and molar masses in mol_tr
-	n=0
-	do nn=1,nq
-	if (indic_diff(nn) .eq. 1) then
-	n=n+1
-	gcmind(n)=nn
-        mol_tr(n) = M_tr(nn)
-	endif
-	enddo
-
-!	print*,'gcmind' ,gcmind
-!       STOP
+!        list of diffused species (if present)
+
+         ListeDiff(1)  = 'co2'
+         ListeDiff(2)  = 'o'
+         ListeDiff(3)  = 'n2'
+         ListeDiff(4)  = 'ar'
+         ListeDiff(5)  = 'co'
+         ListeDiff(6)  = 'h2'
+         ListeDiff(7)  = 'h'
+         ListeDiff(8)  = 'd2'
+         ListeDiff(9)  = 'hd'
+         ListeDiff(10) = 'd'
+         ListeDiff(11) = 'o2'
+         ListeDiff(12) = 'h2o'
+         ListeDiff(13) = 'o3'
+         ListeDiff(14) = 'n'
+         ListeDiff(15) = 'he'
+         ListeDiff(16) = 'n2d'
+
+         i_esc_h  = 1000
+         i_esc_h2 = 1000
+
+!        search for species that can be diffused
+
+         ncompdiff = 0
+         indic_diff(1:nq) = 0
+
+         do nn = 1,nq
+            do n = 1,nb_esp_diff
+               if (ListeDiff(n) == tname(nn)) then
+                  indic_diff(nn) = 1
+               end if
+            end do
+
+            if (indic_diff(nn) == 1) then
+               print*,'species ', tname(nn), 'diffused in moldiff_red'
+               ncompdiff = ncompdiff + 1
+            end if
+         end do
+
+         print*,'number of diffused species:',ncompdiff
+         allocate(gcmind(ncompdiff))
+         allocate(mol_tr(ncompdiff))
+
+!        store gcm indexes in gcmind and molar masses in mol_tr
+
+         n = 0
+         do nn = 1,nq
+            if (indic_diff(nn) == 1) then
+               n = n + 1
+               gcmind(n) = nn
+               mol_tr(n) = M_tr(nn)
+               if (tname(nn) == 'h') i_esc_h = n
+               if (tname(nn) == 'h2') i_esc_h2 = n
+            end if
+         end do
+
 ! find vertical index above which diffusion is computed