Changeset 1549 for trunk/LMDZ.VENUS/libf/phyvenus

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

Location:

trunk/LMDZ.VENUS/libf/phyvenus

Files:

• dyn1d/rcm1d.F (modified) (3 diffs)
• gr_fi_ecrit.F90 (added)
• physiq_mod.F (moved) (moved from trunk/LMDZ.VENUS/libf/phyvenus/physiq.F) (2 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/dyn1d/rcm1d.F

@@ -14 +14 @@
       USE iniphysiq_mod, ONLY: iniphysiq
       USE mod_const_mpi, ONLY: comm_lmdz
+      USE physiq_mod, ONLY: physiq
       IMPLICIT NONE
 
@@ -68 +69 @@
       REAL du(llm),dv(llm),dtemp(llm)
       REAL dudyn(llm),dvdyn(llm),dtempdyn(llm)
-      REAL dpsurf    
+      REAL dpsurf(1)    
       REAL,allocatable :: dq(:,:)
 
@@ -449 +450 @@
 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)

trunk/LMDZ.VENUS/libf/phyvenus/physiq_mod.F

@@ -1 +1 @@
 !
-! $Header: /home/cvsroot/LMDZ4/libf/phylmd/physiq.F,v 1.8 2005/02/24 09:58:18 fairhead Exp $
+! $Id: $
 !
-c
+      MODULE physiq_mod
+
+      IMPLICIT NONE
+
+      CONTAINS
+
       SUBROUTINE physiq (nlon,nlev,nqmax,
      .            debut,lafin,rjourvrai,gmtime,pdtphys,
@@ -1767 +1772 @@
       ENDIF
       
-      RETURN
-      END
-
-
-
-***********************************************************************
-***********************************************************************
-***********************************************************************
-***********************************************************************
-***********************************************************************
-***********************************************************************
-***********************************************************************
-***********************************************************************
-
-      SUBROUTINE gr_fi_ecrit(nfield,nlon,iim,jjmp1,fi,ecrit)
-      IMPLICIT none
-c
-c Tranformer une variable de la grille physique a
-c la grille d'ecriture
-c
-      INTEGER nfield,nlon,iim,jjmp1, jjm
-      REAL fi(nlon,nfield), ecrit(iim*jjmp1,nfield)
-c
-      INTEGER i, n, ig
-c
-      jjm = jjmp1 - 1
-      DO n = 1, nfield
-         DO i=1,iim
-            ecrit(i,n) = fi(1,n)
-            ecrit(i+jjm*iim,n) = fi(nlon,n)
-         ENDDO
-         DO ig = 1, nlon - 2
-           ecrit(iim+ig,n) = fi(1+ig,n)
-         ENDDO
-      ENDDO
-      RETURN
-      END
+      END SUBROUTINE physiq
+
+      END MODULE physiq_mod
+