Changeset 1130 for trunk/LMDZ.MARS/libf/phymars/vdif_kc.F

Timestamp:

Dec 20, 2013, 4:04:56 PM (11 years ago)

Author:

emillour

Message:

Mars GCM:
Series of changes to enable running in parallel (using LMDZ.COMMON dynamics);
Current LMDZ.MARS can still notheless be compiled and run in serial mode
"as previously".
Summary of main changes:

• Main programs (newstart, start2archive, xvik) that used to be in dyn3d have been moved to phymars.
• dyn3d/control.h is now module control_mod.F90
• rearanged input/outputs routines everywhere to handle serial/MPI cases. physdem.F => phyredem.F90 , phyetat0.F => phyetat0.F90 ; all read/write routines for startfi files are gathered in module iostart.F90
• added parallelism related routines init_phys_lmdz.F90, comgeomphy.F90, dimphy.F90, iniphysiq.F90, mod_grid_phy_lmdz.F90, mod_phys_lmdz_mpi_data.F90, mod_phys_lmdz_mpi_transfert.F90, mod_phys_lmdz_omp_data.F90, mod_phys_lmdz_omp_transfert.F90, mod_phys_lmdz_para.F90, mod_phys_lmdz_transfert_para.F90 in phymars and mod_const_mpi.F90 in dyn3d (for compliance with parallel case)
• created generic routines 'planetwide_maxval' and 'planetwide_minval', in module "planetwide_mod", that enable obtaining the min and max of a field over the whole planet.

EM

File:

• trunk/LMDZ.MARS/libf/phymars/vdif_kc.F (modified) (4 diffs)

trunk/LMDZ.MARS/libf/phymars/vdif_kc.F

@@ -32 +32 @@
 c 
 c.......................................................................
-      REAL dt,g
-      REAL zlev(ngrid,nlay+1)
-      REAL zlay(ngrid,nlay)
-      REAL u(ngrid,nlay)
-      REAL v(ngrid,nlay)
-      REAL teta(ngrid,nlay)
-      REAL cd(ngrid)
-      REAL q2(ngrid,nlay+1)
-      REAL km(ngrid,nlay+1)
-      REAL kn(ngrid,nlay+1)
-      REAL zq(ngrid,nlay,nq)
+      REAL,INTENT(IN) :: dt,g
+      REAL,INTENT(IN) :: zlev(ngrid,nlay+1)
+      REAL,INTENT(IN) :: zlay(ngrid,nlay)
+      REAL,INTENT(IN) :: u(ngrid,nlay)
+      REAL,INTENT(IN) :: v(ngrid,nlay)
+      REAL,INTENT(IN) :: teta(ngrid,nlay)
+      REAL,INTENT(IN) :: cd(ngrid)
+      REAL,INTENT(INOUT) :: q2(ngrid,nlay+1)
+      REAL,INTENT(OUT) :: km(ngrid,nlay+1)
+      REAL,INTENT(OUT) :: kn(ngrid,nlay+1)
+      REAL,INTENT(IN) :: zq(ngrid,nlay,nq)
 c.......................................................................
 c
@@ -247 +247 @@
       nlev=nlay+1
       
+      ! Other initializations of local variables (to be clean)
+      long(:,:)=0.
+      n2(:,:)=0.
+      sn(:,:)=0.
+      snq2(:,:)=0.
+      sm(:,:)=0.
+      smq2(:,:)=0.
 c.......................................................................
 c  Special treatment for co2
@@ -381 +388 @@
         gninf=.false.
         gnsup=.false.
-        long(igrid,ilev)=long(igrid,ilev)
-        long(igrid,ilev)=long(igrid,ilev)
+        long(igrid,ilev)=long(igrid,ilev) ! not very useful...
+        long(igrid,ilev)=long(igrid,ilev) ! not very useful...
 c
         IF (gn.lt.gnmin) THEN
@@ -609 +616 @@
 c.......................................................................
 c
+
+!        call writediagfi(ngrid,'vdif_kc_q2','','',3,q2(:,1:nlay))
+!        call writediagfi(ngrid,'vdif_kc_km','','',3,km(:,1:nlay))
+!        call writediagfi(ngrid,'vdif_kc_kn','','',3,kn(:,1:nlay))
+!        call writediagfi(ngrid,'vdif_kc_unsdz','','',3,unsdz(:,1:nlay))
+!        call writediagfi(ngrid,'vdif_kc_unsddecz','','',3,
+!     &                   unsdzdec(:,1:nlay))
+!        call writediagfi(ngrid,'vdif_kc_q','','',3,q(:,1:nlay))
+!        call writediagfi(ngrid,'vdif_kc_kmpre','','',3,
+!     &                   kmpre(:,1:nlay))
+!        call writediagfi(ngrid,'vdif_kc_long','','',3,long(:,1:nlay))
+!        call writediagfi(ngrid,'vdif_kc_sn','','',3,sn(:,1:nlay))
+!        call writediagfi(ngrid,'vdif_kc_sm','','',3,sm(:,1:nlay))
+
       RETURN
       END