Changeset 1047 for trunk/LMDZ.MARS/libf/aeronomars/deposition.F

Timestamp:

Sep 23, 2013, 9:56:47 AM (11 years ago)

Author:

emillour

Message:

Mars GCM:

• IMPORTANT CHANGE: Removed all reference/use of ngridmx (dimphys.h) in routines (necessary prerequisite to using parallel dynamics); in most cases this just means adding 'ngrid' as routine argument, and making local saved variables allocatable (and allocated at first call). In the process, had to convert many *.h files to equivalent modules: yomaer.h => yomaer_h.F90 , surfdat.h => surfdat_h.F90 , comsaison.h => comsaison_h.F90 , yomlw.h => yomlw_h.F90 , comdiurn.h => comdiurn_h.F90 , dimradmars.h => dimradmars_mod.F90 , comgeomfi.h => comgeomfi_h.F90, comsoil.h => comsoil_h.F90 , slope.h => slope_mod.F90
• Also updated EOF routines, everything is now in eofdump_mod.F90
• Removed unused routine lectfux.F (in dyn3d)

EM

File:

• trunk/LMDZ.MARS/libf/aeronomars/deposition.F (modified) (2 diffs)

trunk/LMDZ.MARS/libf/aeronomars/deposition.F

@@ -1 @@
-      subroutine deposition(ig, ig_vl1, pplay, pplev, zzlay, zzlev,
+      subroutine deposition(ngrid, nlayer, nq,
+     &                      ig, ig_vl1, pplay, pplev, zzlay, zzlev,
      $                      zu, zv, zt, zycol, ptimestep, co2ice)
 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
@@ -7 +8 @@
 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
 c
-      use tracer_mod, only: nqmx
+      use surfdat_h, only: z0 ! surface roughness
+      use conc_mod, only: rnew ! specific gas constant
       implicit none
 c
-#include "dimensions.h"
-#include "dimphys.h"
-#include "planete.h"
-#include "chimiedata.h"
-#include "conc.h"
-#include "surfdat.h"
+!#include "dimensions.h"
+!#include "dimphys.h"
+!#include "planete.h"
+!#include "chimiedata.h"
+!#include "conc.h"
+!#include "surfdat.h"
 c
 c     input
 c 
-      integer ig                         ! grid point index
-      integer ig_vl1                     ! viking 1 grid point
-      real    pplay(ngridmx,nlayermx)    ! pressure at the middle of the layers (pa)
-      real    pplev(ngridmx,nlayermx+1)  ! pressure at layer boundaries (pa)
-      real    zzlay(ngridmx,nlayermx)    ! altitude at the middle of the layers (m)
-      real    zzlev(ngridmx,nlayermx+1)  ! altitude at layer boundaries (m)
-      real    zu(ngridmx,nlayermx)       ! u component of the wind (m.s-1)
-      real    zv(ngridmx,nlayermx)       ! v component of the wind (m.s-1)
-      real    zt(ngridmx,nlayermx)       ! temperature (k)
-      real    zycol(nlayermx,nqmx)       ! composition (volume mixing ratio)
-      real    ptimestep                  ! physical timestep (s)
-      real    co2ice(ngridmx)            ! co2 ice surface layer (kg.m-2)
+      integer,intent(in) :: ngrid ! number of atmospheric columns
+      integer,intent(in) :: nlayer ! number of atmospheric layers
+      integer,intent(in) :: nq ! number of tracers
+      integer ig                     ! grid point index
+      integer ig_vl1                 ! viking 1 grid point
+      real    pplay(ngrid,nlayer)    ! pressure at the middle of the layers (pa)
+      real    pplev(ngrid,nlayer+1)  ! pressure at layer boundaries (pa)
+      real    zzlay(ngrid,nlayer)    ! altitude at the middle of the layers (m)
+      real    zzlev(ngrid,nlayer+1)  ! altitude at layer boundaries (m)
+      real    zu(ngrid,nlayer)       ! u component of the wind (m.s-1)
+      real    zv(ngrid,nlayer)       ! v component of the wind (m.s-1)
+      real    zt(ngrid,nlayer)       ! temperature (k)
+      real    zycol(nlayer,nq)       ! composition (volume mixing ratio)
+      real    ptimestep              ! physical timestep (s)
+      real    co2ice(ngrid)          ! co2 ice surface layer (kg.m-2)
 c
 c     local