Changeset 1621 for trunk/LMDZ.MARS/libf/phymars/iniwrite.F

Timestamp:

Oct 25, 2016, 9:23:21 AM (8 years ago)

Author:

emillour

Message:

Further work on full dynamics/physics separation.

LMDZ.COMMON:

• added phy_common/vertical_layers_mod.F90 to store information on vertical grid. This is where routines in the physics should get the information.
• The contents of vertical_layers_mod intialized via dynphy_lonlat/inigeomphy_mod.F90.

LMDZ.MARS:

• physics now completely decoupled from dynamics; the physics package may now be compiled as a library (-libphy option of makelmdz_fcm).
• created an "ini_tracer_mod" routine in module "tracer_mod" for a cleaner initialization of the later.
• removed some purely dynamics-related outputs (etot0, zoom parameters, etc.) from diagfi.nc and stats.nc outputs as these informations are not available in the physics.

LMDZ.GENERIC:

• physics now completely decoupled from dynamics; the physics package may now be compiled as a library (-libphy option of makelmdz_fcm).
• added nqtot to tracer_h.F90.
• removed some purely dynamics-related outputs (etot0, zoom parameters, etc.) from diagfi.nc and stats.nc outputs as these informations are not available in the physics.

LMDZ.VENUS:

• physics now completely decoupled from dynamics; the physics package may now be compiled as a library (-libphy option of makelmdz_fcm).
• added infotrac_phy.F90 to store information on tracers in the physics. Initialized via iniphysiq.
• added cpdet_phy_mod.F90 to store t2tpot etc. functions to be used in the physics. Initialized via iniphysiq. IMPORTANT: there are some hard-coded constants! These should match what is in cpdet_mod.F90 in the dynamics.
• got rid of references to moyzon_mod module within the physics. The required variables (tmoy, plevmoy) are passed to the physics as arguments to physiq.

LMDZ.TITAN:

• added infotrac_phy.F90 to store information on tracers in the physics. Initialized via iniphysiq.
• added cpdet_phy_mod.F90 to store t2tpot etc. functions to be used in the physics.
• Extra work required to completely decouple physics and dynamics: moyzon_mod should be cleaned up and information passed from dynamics to physics as as arguments. Likewise moyzon_ch and moyzon_mu should not be queried from logic_mod (which is in the dynamics).

EM

File:

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

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

@@ -3 +3 @@
       use comsoil_h, only: mlayer, nsoilmx
       USE comcstfi_h, only: g, mugaz, omeg, rad, rcp, pi 
-      USE comvert_mod, ONLY: ap,bp,aps,bps,pseudoalt
-      USE logic_mod, ONLY: fxyhypb,ysinus
-      USE serre_mod, ONLY: clon,clat,grossismx,grossismy,dzoomx,dzoomy
+      USE vertical_layers_mod, ONLY: ap,bp,aps,bps,pseudoalt
+!      USE logic_mod, ONLY: fxyhypb,ysinus
+!      USE serre_mod, ONLY: clon,clat,grossismx,grossismy,dzoomx,dzoomy
       USE time_phylmdz_mod, ONLY: hour_ini, daysec, dtphys
-      USE ener_mod, ONLY: etot0,ptot0,ztot0,stot0,ang0
+!      USE ener_mod, ONLY: etot0,ptot0,ztot0,stot0,ang0
       USE regular_lonlat_mod, ONLY: lon_reg, lat_reg
       USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, nbp_lev
@@ -75 +75 @@
       tab_cntrl(10) = daysec
       tab_cntrl(11) = dtphys
-      tab_cntrl(12) = etot0
-      tab_cntrl(13) = ptot0
-      tab_cntrl(14) = ztot0
-      tab_cntrl(15) = stot0
-      tab_cntrl(16) = ang0
+!      tab_cntrl(12) = etot0
+!      tab_cntrl(13) = ptot0
+!      tab_cntrl(14) = ztot0
+!      tab_cntrl(15) = stot0
+!      tab_cntrl(16) = ang0
 
       tab_cntrl(27) = hour_ini
@@ -85 +85 @@
 c    ..........    P.Le Van  ( ajout le 8/04/96 )    .........
 c         .....        parametres  pour le zoom          ......   
-      tab_cntrl(17)  = clon
-      tab_cntrl(18)  = clat
-      tab_cntrl(19)  = grossismx
-      tab_cntrl(20)  = grossismy
+!      tab_cntrl(17)  = clon
+!      tab_cntrl(18)  = clat
+!      tab_cntrl(19)  = grossismx
+!      tab_cntrl(20)  = grossismy
 c
 c     .....   ajout  le 6/05/97 et le 15/10/97  .......
 c
-      IF ( fxyhypb )   THEN
-        tab_cntrl(21) = 1.
-        tab_cntrl(22) = dzoomx
-        tab_cntrl(23) = dzoomy
-      ELSE
-        tab_cntrl(21) = 0.
-        tab_cntrl(22) = dzoomx
-        tab_cntrl(23) = dzoomy
-        tab_cntrl(24) = 0.
-        IF( ysinus )  tab_cntrl(24) = 1.
-      ENDIF
+!      IF ( fxyhypb )   THEN
+!        tab_cntrl(21) = 1.
+!        tab_cntrl(22) = dzoomx
+!        tab_cntrl(23) = dzoomy
+!      ELSE
+!        tab_cntrl(21) = 0.
+!        tab_cntrl(22) = dzoomx
+!        tab_cntrl(23) = dzoomy
+!        tab_cntrl(24) = 0.
+!        IF( ysinus )  tab_cntrl(24) = 1.
+!      ENDIF
 
 c    .........................................................