Changeset 2641 for LMDZ5/branches/testing/libf/dyn3dpar/leapfrog_p.F

Timestamp:

Sep 29, 2016, 11:26:46 PM (8 years ago)

Author:

Laurent Fairhead

Message:

Merged trunk changes r2593:2640 into testing branch

Location:

LMDZ5/branches/testing

Files:

• . (modified) (1 prop)
• libf/dyn3dpar/leapfrog_p.F (modified) (7 diffs)

LMDZ5/branches/testing/libf/dyn3dpar/leapfrog_p.F

@@ -28 +28 @@
      &                       periodav, ok_dyn_ave, output_grads_dyn,
      &                       iapp_tracvl
+       USE comvert_mod, ONLY: ap,bp,pressure_exner,presnivs
+       USE comconst_mod, ONLY: cpp, dtvr, ihf, dtphys, pi, jmp1
+       USE logic_mod, ONLY: iflag_phys,ok_guide,forward,leapf,apphys,
+     &                      statcl,conser,apdiss,purmats,ok_strato
+       USE temps_mod, ONLY: itaufin,jD_ref,jH_ref,day_ini,
+     &                        day_ref,start_time,dt
+
       IMPLICIT NONE
 
@@ -63 +70 @@
 #include "dimensions.h"
 #include "paramet.h"
-#include "comconst.h"
 #include "comdissnew.h"
-#include "comvert.h"
 #include "comgeom.h"
-#include "logic.h"
-#include "temps.h"
-#include "ener.h"
 #include "description.h"
-#include "serre.h"
 !#include "com_io_dyn.h"
 #include "iniprint.h"
@@ -1455 +1456 @@
                 call Gather_Field(teta,ip1jmp1,llm,0)
                 call Gather_Field(pk,ip1jmp1,llm,0)
-                call Gather_Field(phi,ip1jmp1,llm,0)
+                call Gather_Field(phi,ip1jmp1,llm,0)
                 do iq=1,nqtot
                   call Gather_Field(q(1,1,iq),ip1jmp1,llm,0)
                 enddo
-                call Gather_Field(masse,ip1jmp1,llm,0)
+                call Gather_Field(masse,ip1jmp1,llm,0)
                 call Gather_Field(ps,ip1jmp1,1,0)
-                call Gather_Field(phis,ip1jmp1,1,0)
+                call Gather_Field(phis,ip1jmp1,1,0)
                 if (mpi_rank==0) then
                  CALL writedynav(itau,vcov,
      &                 ucov,teta,pk,phi,q,masse,ps,phis)
-                endif
+                endif
 #endif
 !$OMP END MASTER
@@ -1477 +1478 @@
             IF( MOD(itau,iecri).EQ.0) THEN
              ! Ehouarn: output only during LF or Backward Matsuno
-             if (leapf.or.(.not.leapf.and.(.not.forward))) then
+             if (leapf.or.(.not.leapf.and.(.not.forward))) then
 c$OMP BARRIER
 c$OMP MASTER
@@ -1515 +1516 @@
                 call Gather_Field(ucov,ip1jmp1,llm,0)
                 call Gather_Field(teta,ip1jmp1,llm,0)
-                call Gather_Field(phi,ip1jmp1,llm,0)
+                call Gather_Field(phi,ip1jmp1,llm,0)
                 do iq=1,nqtot
                   call Gather_Field(q(1,1,iq),ip1jmp1,llm,0)
                 enddo
-                call Gather_Field(masse,ip1jmp1,llm,0)
+                call Gather_Field(masse,ip1jmp1,llm,0)
                 call Gather_Field(ps,ip1jmp1,1,0)
-                call Gather_Field(phis,ip1jmp1,1,0)
+                call Gather_Field(phis,ip1jmp1,1,0)
                 if (mpi_rank==0) then
-                 CALL writehist(itau,vcov,ucov,teta,phi,q,masse,ps,phis)
-                endif
+                 CALL writehist(itau,vcov,ucov,teta,phi,q,masse,ps,phis)
+                endif
 !              CALL writehist_p(histid,histvid, itau,vcov, 
 !     &                         ucov,teta,phi,q,masse,ps,phis)
@@ -1650 +1651 @@
                 call Gather_Field(teta,ip1jmp1,llm,0)
                 call Gather_Field(pk,ip1jmp1,llm,0)
-                call Gather_Field(phi,ip1jmp1,llm,0)
+                call Gather_Field(phi,ip1jmp1,llm,0)
                 do iq=1,nqtot
                   call Gather_Field(q(1,1,iq),ip1jmp1,llm,0)
                 enddo
-                call Gather_Field(masse,ip1jmp1,llm,0)
+                call Gather_Field(masse,ip1jmp1,llm,0)
                 call Gather_Field(ps,ip1jmp1,1,0)
-                call Gather_Field(phis,ip1jmp1,1,0)
+                call Gather_Field(phis,ip1jmp1,1,0)
                 if (mpi_rank==0) then
                  CALL writedynav(itau,vcov,
      &                 ucov,teta,pk,phi,q,masse,ps,phis)
-                endif
+                endif
 #endif
 !$OMP END MASTER
@@ -1706 +1707 @@
                 call Gather_Field(ucov,ip1jmp1,llm,0)
                 call Gather_Field(teta,ip1jmp1,llm,0)
-                call Gather_Field(phi,ip1jmp1,llm,0)
+                call Gather_Field(phi,ip1jmp1,llm,0)
                 do iq=1,nqtot
                   call Gather_Field(q(1,1,iq),ip1jmp1,llm,0)
                 enddo
-                call Gather_Field(masse,ip1jmp1,llm,0)
+                call Gather_Field(masse,ip1jmp1,llm,0)
                 call Gather_Field(ps,ip1jmp1,1,0)
-                call Gather_Field(phis,ip1jmp1,1,0)
+                call Gather_Field(phis,ip1jmp1,1,0)
                 if (mpi_rank==0) then
                  CALL writehist(itau,vcov,ucov,teta,phi,q,masse,ps,phis)
-                endif
+                endif
 !                CALL writehist_p(histid, histvid, itau,vcov , 
 !     &                           ucov,teta,phi,q,masse,ps,phis)