Changeset 3666 for LMDZ6/trunk/libf/phylmd/tracreprobus_mod.F90

Timestamp:

Apr 20, 2020, 12:13:34 PM (4 years ago)

Author:

lfalletti

Message:

Adding changes for Reprobus

File:

• LMDZ6/trunk/libf/phylmd/tracreprobus_mod.F90 (modified) (6 diffs)

LMDZ6/trunk/libf/phylmd/tracreprobus_mod.F90

@@ -16 +16 @@
     USE CHEM_REP, ONLY : pdt_rep, &  ! pas de temps reprobus
          daynum, iter, &             ! jourjulien, iteration chimie
-         pdel
+         pdel,&
+         d_q_rep,d_ql_rep,d_qi_rep
 #endif
     IMPLICIT NONE
@@ -46 +47 @@
 ! Local variables
 !----------------
-    INTEGER :: it, k
+    INTEGER :: it, k, niter
 
 #ifdef REPROBUS
     !   -- CHIMIE REPROBUS --
-    pdt_rep=pdtphys/2.
+!    pdt_rep=pdtphys/2.
+    niter=pdtphys/pdt_rep
+    write(*,*)'nb d appel de REPROBUS',niter
     
     DO k = 1, klev
@@ -60 +63 @@
        tr_seri(:,:,11)=tr_seri(:,:,8)
     END IF
+
+    d_q_rep(:,:)  =0.
+    d_ql_rep(:,:) =0.
+    d_qi_rep(:,:) =0.
     
-    DO  iter = 1,2
+    DO  iter = 1,niter
        daynum = FLOAT(julien) + gmtime + (iter-1)*pdt_rep/86400.
        
-       DO it=1, nbtr
+!       DO it=1, nbtr
 !     WRITE(lunout,*)it,' ',minval(tr_seri(:,:,it)),maxval(tr_seri(:,:,it))
 ! seulement pour les especes chimiques (pas l'age de l'air)
@@ -70 +77 @@
 ! correction: a 1.e-30 quand =0 ou negatif et
 ! call abort si >ou= 1.e10
-          WRITE(*,*)it,'nqtot',nqtot,'nbtr',nbtr
-          IF (it < nqtot) THEN
-             WRITE(*,*)'iciav',it,nqtot
-#ifdef REPROBUS
-             CALL minmaxqfi_chimie(it,tr_seri(1,1,it),0.,1.e10,'avant chimie ')
-#endif
-             WRITE(*,*)iter,'avpres'
-          ENDIF
-       ENDDO
+!          WRITE(*,*)it,'nqtot',nqtot,'nbtr',nbtr
+!          IF (it < nqtot) THEN
+!             WRITE(*,*)'iciav',it,nqtot
+!#ifdef REPROBUS
+!             CALL minmaxqfi_chimie(it,tr_seri(1,1,it),0.,1.e10,'avant chimie ')
+!#endif
+!             WRITE(*,*)iter,'avpres'
+!          ENDIF
+!       ENDDO
        
 #ifdef REPROBUS
@@ -95 +102 @@
        !                 et transporte par CHEM_REP
 
-       DO it=1, nbtr
+!       DO it=1, nbtr
 !     WRITE(lunout,*)it,' ',minval(tr_seri(:,:,it)),maxval(tr_seri(:,:,it))
 ! seulement pour les especes chimiques (pas l'age de l'air)
@@ -101 +108 @@
 ! correction: a 1.e-30 quand =0 ou negatif et
 ! call abort si >ou= 1.e10
-          WRITE(*,*)it,'nqtot',nqtot,'nbtr',nbtr
-          IF (it < nqtot) THEN
-             WRITE(*,*)'iciap',it,nqtot
-             CALL minmaxqfi_chimie(it,tr_seri(1,1,it),0.,1.e10,'apres chemmain')
-             WRITE(*,*)iter,'appres'
-          ENDIF
-       ENDDO
+!          WRITE(*,*)it,'nqtot',nqtot,'nbtr',nbtr
+!          IF (it < nqtot) THEN
+!             WRITE(*,*)'iciap',it,nqtot
+!             CALL minmaxqfi_chimie(it,tr_seri(1,1,it),0.,1.e10,'apres chemmain')
+!             WRITE(*,*)iter,'appres'
+!          ENDIF
+!       ENDDO
 
 #endif