Changeset 4058

Timestamp:

Jan 13, 2022, 11:47:02 PM (2 years ago)

Author:

dcugnet

Message:

Fix the call to vlspltgen_loc in advtrac_loc (one wrong additional calling argument).

File:

• LMDZ6/trunk/libf/dyn3dmem/advtrac_loc.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/dyn3dmem/advtrac_loc.F90

@@ -10 +10 @@
    !            M.A Filiberti (04/2002)
    !
-   USE infotrac,     ONLY: nqtot, tracers,ok_iso_verif
+   USE infotrac,     ONLY: nqtot, tracers, ok_iso_verif
    USE control_mod,  ONLY: iapp_tracvl, day_step, planet_type
    USE comconst_mod, ONLY: dtvr
@@ -106 +106 @@
    !---------------------------------------------------------------------------
 !ym   CALL massbar_p(massem,massebx,masseby)
-!ym   ----> Normalement, inutile pour les schémas classiques
-!ym   ----> Revérifier lors de la parallélisation des autres schemas
+!ym   ----> Normalement, inutile pour les schemas classiques
+!ym   ----> Reverifier lors de la parallelisation des autres schemas
 
 #ifdef DEBUG_IO    
@@ -128 +128 @@
 !  CALL WaitRequest(TestRequest)
 !$OMP BARRIER
-                 
+
 !  WRITE(*,*) 'advtrac 157: appel de vlspltgen_loc'
-   CALL vlspltgen_loc(q, tracers(:)%iadv, 2., massem, wg, pbarug, pbarvg, dtvr, p, pk, teta )
+   CALL vlspltgen_loc(q, 2., massem, wg, pbarug, pbarvg, dtvr, p, pk, teta )
 
 #ifdef DEBUG_IO     
@@ -167 +167 @@
          CASE(12)  !--- Schema de Frederic Hourdin
          !--------------------------------------------------------------------
+            STOP 'advtrac : appel a vlspltqs :schema non parallelise'
             CALL adaptdt(iadv,dtbon,n,pbarug,massem)   ! pas de temps adaptatif
             IF(n > 1) WRITE(*,*) 'WARNING horizontal dt=',dtbon,'dtvr=',dtvr,'n=',n
@@ -259 +260 @@
 
       !----------------------------------------------------------------------
-      ! On impose une seule valeur du traceur au pôle Sud j=jjm+1=jjp1 et Nord j=1
+      ! On impose une seule valeur du traceur au pole Sud j=jjm+1=jjp1 et Nord j=1
       !----------------------------------------------------------------------
       !  CALL traceurpole(q(1,1,iq),massem)