Changeset 1541 for trunk/LMDZ.MARS/libf/phymars/writediagfi.F

Timestamp:

Apr 17, 2016, 10:23:21 PM (9 years ago)

Author:

emillour

Message:

Mars GCM:

• fix for 1D in writediagfi to enable writing at "ecritphy" rate.
• removed iniprint.h from phymars/dyn1d since it is in "misc"
• Some code cleanup in anticipation of future updates:
  • changed variable names in comgeomphy.F90: give them more explicit names: rlond => longitude , rlatd => latitude, airephy => cell_area, cuphy => dx , cvphy => dy
  • removed long(), lati() and area() from comgeomfi_h.F90, use longitude(), latitude() and cell_are() from comgeomphy.F90 instead

EM

File:

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

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

@@ -40 +40 @@
 !=================================================================
       use surfdat_h, only: phisfi
-      use comgeomphy, only: airephy
+      use comgeomphy, only: cell_area
       use time_phylmdz_mod, only: ecritphy, day_step, iphysiq, day_ini
       USE mod_phys_lmdz_para, only : is_parallel, is_mpi_root,
@@ -118 +118 @@
 !Sortie des variables au rythme voulu
 
-      irythme = int(ecritphy) ! sortie au rythme de ecritphy
-!     irythme = iconser  ! sortie au rythme des variables de controle
-!     irythme = iphysiq  ! sortie a tous les pas physique
-!     irythme = iecri*day_step ! sortie au rythme des fichiers histoires
-!     irythme = periodav*day_step ! sortie au rythme des fichiers histmoy
+      irythme = int(ecritphy) ! output rate set by ecritphy
 
 !***************************************************************
@@ -189 +185 @@
           ! Gather phisfi() geopotential on physics grid
           call Gather(phisfi,phisfi_glo)
-          ! Gather airephy() mesh area on physics grid
-          call Gather(airephy,areafi_glo)
+          ! Gather cell_area() mesh area on physics grid
+          call Gather(cell_area,areafi_glo)
 #else
          phisfi_glo(:)=phisfi(:)
-         areafi_glo(:)=airephy(:)
+         areafi_glo(:)=cell_area(:)
 #endif
 
@@ -260 +256 @@
 
       endif ! if (firstnom.eq.'1234567890')
-
-      if (klon_glo.eq.1) then
-        ! in testphys1d, for the 1d version of the GCM, iphysiq and irythme
-        ! are undefined; so set them to 1
-        iphysiq=1 
-        irythme=1
-      endif
 
 ! Increment time index 'zitau' if it is the "fist call" (at given time level)