Changeset 1549 for trunk/LMDZ.MARS/libf/phymars/dyn1d

Timestamp:

May 6, 2016, 12:30:29 PM (9 years ago)

Author:

emillour

Message:

All GCMs:
Further adaptations to keep up with changes in LMDZ5 concerning
physics/dynamics separation (up to rev r2420 of LMDZ5)

• all physics packages:
• added module callphysiq_mod.F90 in dynphy_lonlat/phy* which contains the routine "call_physiq" which is called by calfis* and calls the physics. This way different "physiq" routine from different physics packages may be called: The calfis* routines now exposes all available fields that might be transmitted to physiq but which is actually send (ie: expected/needed by physiq) is decided in call_physiq.
• turned "physiq.F[90]" into module "physiq_mod.F[90]" for better control of "physiq" arguments. for phyvenus/phytitan, extracted gr_fi_ecrit from physiq.F as gr_fi_ecrit.F90 (note that it can only work in serial).

• misc:
• updated wxios.F90 to keep up with LMDZ5 modifications.

• dyn3d_common:
• infotrac.F90 keep up with LMDZ5 modifications (cosmetics)

• dyn3d:
• gcm.F90: cosmetic cleanup.
• leapfrog.F90: fix computation of date as function of itau.

• dyn3dpar:
• gcm.F: cosmetic cleanup.
• leapfrog_p.F90: fix computation of date as function of itau.

NB: physics are given the date corresponding to the end of the
physics step.

• dynphy_lonlat:
• calfis.F : added computation of relative wind vorticity.
• calfis_p.F: added computation of relative wind vorticity (input required by Earth physics)

EM

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F (modified) (3 diffs)

trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F

@@ -27 +27 @@
       use mod_interface_dyn_phys, only: init_interface_dyn_phys
       USE phys_state_var_init_mod, ONLY: phys_state_var_init
+      USE physiq_mod, ONLY: physiq
       IMPLICIT NONE
 
@@ -95 +96 @@
       REAL du(nlayer),dv(nlayer),dtemp(nlayer)
       REAL dudyn(nlayer),dvdyn(nlayer),dtempdyn(nlayer)
-      REAL dpsurf    
+      REAL dpsurf(1)    
       REAL,ALLOCATABLE :: dq(:,:)
       REAL,ALLOCATABLE :: dqdyn(:,:)
@@ -792 +793 @@
 c       ----------------------------------------------------------
 
-           psurf=psurf+dtphys*dpsurf   ! surface pressure change
+           psurf=psurf+dtphys*dpsurf(1)   ! surface pressure change
            DO ilevel=1,nlevel
              plev(ilevel)=ap(ilevel)+psurf*bp(ilevel)