Changeset 1130 for trunk/LMDZ.MARS/libf/phymars/eofdump_mod.F90

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/eofdump_mod.F90 (modified) (4 diffs)

trunk/LMDZ.MARS/libf/phymars/eofdump_mod.F90

@@ -27 +27 @@
       integer,intent(in) :: ngrid ! total number of physics grid points
       integer,intent(in) :: nlayer ! number of atmospheric layers
-      real*4 u(ngrid,nlayer)
-      real*4 v(ngrid,nlayer)
-      real*4 t(ngrid,nlayer)
-      real*4 rho(ngrid,nlayer)
-      real*4 ps(ngrid)
+      real,intent(in) :: u(ngrid,nlayer)
+      real,intent(in) :: v(ngrid,nlayer)
+      real,intent(in) :: t(ngrid,nlayer)
+      real,intent(in) :: rho(ngrid,nlayer)
+      real,intent(in) :: ps(ngrid)
       integer,save :: count=0
       integer i,j,l, ig
@@ -56 +56 @@
           do j=1+eofskip/2,jjm+1,eofskip
             ig = 1+ (j-2)*iim +i
+#ifdef NC_DOUBLE
+            write(uedata) (real(u(ig,l)),l=1,nlayer)
+            write(uedata) (real(v(ig,l)),l=1,nlayer)
+            write(uedata) (real(t(ig,l)),l=1,nlayer)
+            write(uedata) (real(rho(ig,l)),l=1,nlayer)
+            write(uedata) real(ps(ig))
+#else
             write(uedata) (u(ig,l),l=1,nlayer)
             write(uedata) (v(ig,l),l=1,nlayer)
@@ -61 +68 @@
             write(uedata) (rho(ig,l),l=1,nlayer)
             write(uedata) ps(ig)
+#endif
           enddo
         enddo
@@ -112 +120 @@
           if(j.eq.1) stop 'Problem in ineofdump.F'
           if(j.eq.jjm+1) stop 'Problem in ineofdump.F'
+#ifdef NC_DOUBLE
+          write(uehead,*) real(long(ig)*180./pi),real(lati(ig)*180./pi)
+#else
           write(uehead,*) long(ig)*180./pi, lati(ig)*180./pi
+#endif
 !         write(*,*) 'eof grid j=',j,' lat= ', lati(ig)*180./pi
         enddo
       enddo
 
+#ifdef NC_DOUBLE
+      write(uehead,*) real(aps)
+      write(uehead,*) real(bps)
+#else
       write(uehead,*) aps
       write(uehead,*) bps
+#endif
       close(uehead)
 !