Changeset 2344 for LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

Timestamp:

Aug 21, 2015, 9:23:13 AM (9 years ago)

Author:

Ehouarn Millour

Message:

Physics/dynamics separation: get rid of all the 'include "temps.h"' in the physics; variables in module time_phylmdz_mod must be used instead. Also added JD_cur, JH_cur and JD_ref in module phys_cal_mod, in preparation for having physics handle its calendar internally.
EM

File:

• LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90 (modified) (6 diffs)

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -25 +25 @@
 
     USE dimphy, only: klon, klev, klevp1, nslay
-    USE control_mod, only: day_step, iphysiq
+    USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
+    USE time_phylmdz_mod, only: day_step_phy, start_time, itau_phy
     USE phys_output_ctrlout_mod, only: o_phis, o_aire, is_ter, is_lic, is_oce, &
          is_ave, is_sic, o_contfracATM, o_contfracOR, &
@@ -268 +269 @@
 
 
-    INCLUDE "temps.h"
     INCLUDE "clesphys.h"
     INCLUDE "thermcell.h"
     INCLUDE "compbl.h"
     INCLUDE "YOMCST.h"
-    INCLUDE "dimensions.h"
 
     ! Input
@@ -288 +287 @@
     REAL, DIMENSION(klon, klev+1) :: paprs
     REAL, DIMENSION(klon,klev,nqtot) :: qx, d_qx
-    REAL, DIMENSION(klon, llm) :: zmasse
+    REAL, DIMENSION(klon, klev) :: zmasse
     LOGICAL :: flag_aerosol_strat
     INTEGER :: flag_aerosol 
@@ -295 +294 @@
 
     ! Local
-    INTEGER, PARAMETER :: jjmp1=jjm+1-1/jjm
     INTEGER :: itau_w
     INTEGER :: i, iinit, iinitend=1, iff, iq, nsrf, k, ll, naero
@@ -302 +300 @@
     REAL, DIMENSION (klon,klev+1) :: zx_tmp_fi3d1
     CHARACTER (LEN=4)              :: bb2
-    INTEGER, DIMENSION(iim*jjmp1)  :: ndex2d
-    INTEGER, DIMENSION(iim*jjmp1*klev) :: ndex3d
+    INTEGER, DIMENSION(nbp_lon*nbp_lat)  :: ndex2d
+    INTEGER, DIMENSION(nbp_lon*nbp_lat*klev) :: ndex3d
     REAL, PARAMETER :: dobson_u = 2.1415e-05 ! Dobson unit, in kg m-2
 !   REAL, PARAMETER :: missing_val=nf90_fill_real
@@ -312 +310 @@
 
     ! On calcul le nouveau tau:
-    itau_w = itau_phy + itap + start_time * day_step / iphysiq
+    itau_w = itau_phy + itap + start_time * day_step_phy
     ! On le donne à iophy pour que les histwrite y aient accès:
     CALL set_itau_iophy(itau_w)