Changeset 3218

Timestamp:

Feb 16, 2024, 9:15:46 AM (10 months ago)

Author:

flefevre

Message:

VENUS: bug fix in Jeans escape (wrong index for H and H2) + some cosmetics

File:

• trunk/LMDZ.VENUS/libf/phyvenus/moldiff_red.F90 (modified) (1 diff)

trunk/LMDZ.VENUS/libf/phyvenus/moldiff_red.F90

@@ -106 +106 @@
 
 ! 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