Changeset 2270 for LMDZ5/trunk/libf/dyn3d/advtrac.F90

Timestamp:

May 7, 2015, 5:45:04 PM (9 years ago)

Author:

crisi

Message:

Adding isotopes in the dynamics and more generally tracers of tracers.
CRisi

File:

• LMDZ5/trunk/libf/dyn3d/advtrac.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/dyn3d/advtrac.F90

@@ -9 +9 @@
   !            M.A Filiberti (04/2002)
   !
-  USE infotrac, ONLY: nqtot, iadv
+  USE infotrac, ONLY: nqtot, iadv,nqperes,ok_iso_verif
   USE control_mod, ONLY: iapp_tracvl, day_step
 
@@ -223 +223 @@
      !     Appel des sous programmes d'advection
      !-----------------------------------------------------------
-     do iq=1,nqtot
+
+     if (ok_iso_verif) then
+           call check_isotopes_seq(q,1,ip1jmp1,'advtrac 162')
+     endif !if (ok_iso_verif) then
+
+     do iq=1,nqperes
         !        call clock(t_initial)
         if(iadv(iq) == 0) cycle 
@@ -230 +235 @@
         !   ----------------------------------------------------------------
         if(iadv(iq).eq.10) THEN
-           call vlsplt(q(1,1,iq),2.,massem,wg,pbarug,pbarvg,dtvr)
-
+           ! CRisi: on fait passer tout q pour avoir acces aux fils
+           
+           write(*,*) 'advtrac 239: iq,q(1721,19,:)=',iq,q(1721,19,:)     
+           call vlsplt(q,2.,massem,wg,pbarug,pbarvg,dtvr,iq)
+           
            !   ----------------------------------------------------------------
            !   Schema "pseudo amont" + test sur humidite specifique
@@ -238 +246 @@
         else if(iadv(iq).eq.14) then
            !
-           CALL vlspltqs( q(1,1,1), 2., massem, wg , &
-                pbarug,pbarvg,dtvr,p,pk,teta )
+           write(*,*) 'advtrac 248: iq,q(1721,19,:)=',iq,q(1721,19,:)
+           CALL vlspltqs( q, 2., massem, wg , &
+                pbarug,pbarvg,dtvr,p,pk,teta,iq)
+           
            !   ----------------------------------------------------------------
            !   Schema de Frederic Hourdin
@@ -388 +398 @@
      end DO
 
+     if (ok_iso_verif) then
+           call check_isotopes_seq(q,ip1jmp1,'advtrac 397')
+     endif !if (ok_iso_verif) then
 
      !------------------------------------------------------------------