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

Timestamp:

Apr 2, 2016, 4:09:43 PM (9 years ago)

Author:

emillour

Message:

Mars GCM:

• Got rid of references to "dimensions.h" from physics packages: use nbp_lon (=iim), nbp_lat (==jjp1) and nbp_lev from module mod_grid_phy_lmdz (in phy_common) instead.
• Added "ioipsl_getin_p_mod.F90" (getin_p routine) in phy_common to correctly read in parameters from *.def files in a parallel environment.

EM

File:

• trunk/LMDZ.MARS/libf/phymars/eofdump_mod.F90 (modified) (7 diffs)

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

@@ -1 +1 @@
 module eofdump_mod
 ! this module controls the production of data for EOFs
+! it won't work if run in parallel (but it's OK, we don't use it anymore...)
+! Mainly kept for reference.
 implicit none
 ! Dump profiles for EOFs every ieofs physics timesteps,
@@ -15 +17 @@
       subroutine eofdump(ngrid,nlayer,u,v,t,rho,ps)
 
+      use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat
       implicit none
 !
@@ -21 +24 @@
 !     Corrected small bug in sampling rate/count, EM 11/2007
 !
-#include "dimensions.h"
 !
 
@@ -52 +54 @@
       if (mod(count+1,ieofs).eq.0) then
 !        write(*,*)'eofdump: dump --> ps(1)=',ps(1)
-        do i=1,iim,eofskip
-          do j=1+eofskip/2,jjm+1,eofskip
-            ig = 1+ (j-2)*iim +i
+        do i=1,nbp_lon,eofskip
+          do j=1+eofskip/2,nbp_lat,eofskip
+            ig = 1+ (j-2)*nbp_lon +i
 #ifdef NC_DOUBLE
             write(uedata) (real(u(ig,l)),l=1,nlayer)
@@ -82 +84 @@
       use time_phylmdz_mod, only: daysec, dtphys
       USE comvert_mod, ONLY: aps,bps
+      use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat
       implicit none
 !
 !     Initialise dumping of profiles for EOF calculations
 !
-#include "dimensions.h"
 
       integer,intent(in) :: ngrid ! total number of physics grid points
@@ -96 +98 @@
 
       if (firstcall) then
-         npgrid=ngrid+2*(iim-1)
+         npgrid=ngrid+2*(nbp_lon-1)
          firstcall=.false.
       endif
@@ -111 +113 @@
       open(uehead,file='profiles.hdr',form='formatted')
       write(uehead,*) 0.E+0,0,0,ieofs,1,0
-      write(uehead,*) iim,npgrid/iim,npgrid,nlayer
+      write(uehead,*) nbp_lon,npgrid/nbp_lon,npgrid,nlayer
 
-      do i=1,iim,eofskip
-        do j=1+eofskip/2,jjm+1,eofskip    
-          ig = 1+ (j-2)*iim +i
+      do i=1,nbp_lon,eofskip
+        do j=1+eofskip/2,nbp_lat,eofskip    
+          ig = 1+ (j-2)*nbp_lon +i
           if(j.eq.1) stop 'Problem in ineofdump.F'
-          if(j.eq.jjm+1) stop 'Problem in ineofdump.F'
+          if(j.eq.nbp_lat) stop 'Problem in ineofdump.F'
 #ifdef NC_DOUBLE
           write(uehead,*) real(long(ig)*180./pi),real(lati(ig)*180./pi)