Changeset 1549 for trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F

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.GENERIC/libf/phystd/dyn1d/rcm1d.F (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F

@@ -32 +32 @@
       use mod_interface_dyn_phys, only: init_interface_dyn_phys
       use inifis_mod, only: inifis
+      use physiq_mod, only: physiq
       implicit none
 
@@ -96 +97 @@
       REAL du(llm),dv(llm),dtemp(llm)
       REAL dudyn(llm),dvdyn(llm),dtempdyn(llm)
-      REAL dpsurf    
+      REAL dpsurf(1)    
       REAL,ALLOCATABLE :: dq(:,:)
       REAL,ALLOCATABLE :: dqdyn(:,:)
@@ -949 +950 @@
 c       ----------------------------------------------------------
 
-           psurf=psurf+dtphys*dpsurf   ! evolution de la pression de surface
+           psurf=psurf+dtphys*dpsurf(1)   ! evolution de la pression de surface
            DO ilevel=1,nlevel
               plev(ilevel)=ap(ilevel)+psurf*bp(ilevel)