Changeset 3221

Timestamp:

Feb 16, 2024, 1:13:16 PM (10 months ago)

Author:

flefevre

Message:

Venus : a bit more cosmetics in molecular diffusion

File:

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

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

@@ -166 +166 @@
          end do
 
-! find vertical index above which diffusion is computed
-
-        do l=1,nlayer
-        if (pplay(1,l) .gt. Pdiff) then
-        il0=l
-        endif
-        enddo
-        il0=il0+1
-        print*,'vertical index for diffusion',il0,pplay(1,il0)
-
-        allocate(dij(ncompdiff,ncompdiff))
-
-        call moldiffcoeff_red(dij,indic_diff,gcmind,ncompdiff)
-        print*,'MOLDIFF  EXO'
-
-! allocatation des tableaux dependants du nombre d especes diffusees
-        allocate(qq(nlayer,ncompdiff))
-        allocate(qnew(nlayer,ncompdiff))
-        allocate(qint(nlayer,ncompdiff))
-        allocate(FacMass(nlayer,ncompdiff))
-        allocate(rhok(nlayer,ncompdiff))
-        allocate(rhokinit(nlayer,ncompdiff))
-
-        allocate(wi(ncompdiff))
-        allocate(wad(ncompdiff))
-        allocate(uthermal(ncompdiff))
-        allocate(lambdaexo(ncompdiff))
-        allocate(Hspecie(ncompdiff))
-        allocate(Mtot1(ncompdiff))
-        allocate(Mtot2(ncompdiff))
-        allocate(Mraf1(ncompdiff))
-        allocate(Mraf2(ncompdiff))
-
-        firstcall= .false.
-        step=1
-      endif ! of if (firstcall)
-
-!
+!        find vertical index above which diffusion is computed
+
+         do l = 1,nlayer
+            if (pplay(1,l) > Pdiff) then
+               il0 = l
+            end if
+         end do
+
+         il0 = il0 + 1
+
+         print*,'vertical index for diffusion',il0,pplay(1,il0)
+
+         allocate(dij(ncompdiff,ncompdiff))
+
+         call moldiffcoeff_red(dij,indic_diff,gcmind,ncompdiff)
+
+         print*,'MOLDIFF  EXO'
+
+!       allocation of arrays depending on the number of diffused species
+
+         allocate(qq(nlayer,ncompdiff))
+         allocate(qnew(nlayer,ncompdiff))
+         allocate(qint(nlayer,ncompdiff))
+         allocate(FacMass(nlayer,ncompdiff))
+         allocate(rhok(nlayer,ncompdiff))
+         allocate(rhokinit(nlayer,ncompdiff))
+
+         allocate(wi(ncompdiff))
+         allocate(wad(ncompdiff))
+         allocate(uthermal(ncompdiff))
+         allocate(lambdaexo(ncompdiff))
+         allocate(Hspecie(ncompdiff))
+         allocate(Mtot1(ncompdiff))
+         allocate(Mtot2(ncompdiff))
+         allocate(Mraf1(ncompdiff))
+         allocate(Mraf2(ncompdiff))
+
+         firstcall = .false.
+         step = 1
+
+      end if ! of if (firstcall)
+
 !ccccccccccccccccccccccccccccccccccccccccccccccccccccccc