Changeset 2344 for LMDZ5/trunk/libf/phylmd/phys_output_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_mod.F90 (modified) (8 diffs)

LMDZ5/trunk/libf/phylmd/phys_output_mod.F90

@@ -44 +44 @@
     USE surface_data, ONLY : ok_snow
     USE phys_output_ctrlout_mod
-    USE mod_grid_phy_lmdz, only: klon_glo
+    USE mod_grid_phy_lmdz, only: klon_glo,nbp_lon,nbp_lat
     USE print_control_mod, ONLY: prt_level,lunout
     USE vertical_layers_mod, ONLY: ap,bp,preff,presnivs
-
+    USE time_phylmdz_mod, ONLY: day_ini, itau_phy, start_time, annee_ref, day_ref
 #ifdef CPP_XIOS
     ! ug Pour les sorties XIOS
@@ -54 +54 @@
 
     IMPLICIT NONE
-    include "dimensions.h"
-    include "temps.h"
     include "clesphys.h"
     include "thermcell.h"
@@ -69 +67 @@
     REAL, DIMENSION(klon, klev+1), INTENT(IN)   :: paprs
     REAL, DIMENSION(klon,klev,nqtot), INTENT(IN):: qx, d_qx
-    REAL, DIMENSION(klon, llm), INTENT(IN)      :: zmasse
+    REAL, DIMENSION(klon, klev), INTENT(IN)      :: zmasse
 
 
@@ -106 +104 @@
     CHARACTER(LEN=2)                      :: bb3
     CHARACTER(LEN=6)                      :: type_ocean
-    INTEGER, DIMENSION(iim*jjmp1)         ::  ndex2d
-    INTEGER, DIMENSION(iim*jjmp1*klev)    :: ndex3d
+    INTEGER, DIMENSION(nbp_lon*jjmp1)         ::  ndex2d
+    INTEGER, DIMENSION(nbp_lon*jjmp1*klev)    :: ndex3d
     INTEGER                               :: imin_ins, imax_ins
     INTEGER                               :: jmin_ins, jmax_ins
@@ -343 +341 @@
           IF (phys_out_regfkey(iff)) then
              imin_ins=1
-             imax_ins=iim
+             imax_ins=nbp_lon
              jmin_ins=1
              jmax_ins=jjmp1
 
              ! correction abderr        
-             do i=1,iim
+             do i=1,nbp_lon
                 WRITE(lunout,*)'io_lon(i)=',io_lon(i)
                 IF (io_lon(i).le.phys_out_lonmin(iff)) imin_ins=i
@@ -367 +365 @@
                   io_lat(jmax_ins),io_lat(jmin_ins)
 
-             CALL histbeg(phys_out_filenames(iff),iim,io_lon,jjmp1,io_lat, &
+             CALL histbeg(phys_out_filenames(iff),nbp_lon,io_lon,jjmp1,io_lat, &
                   imin_ins,imax_ins-imin_ins+1, &
                   jmin_ins,jmax_ins-jmin_ins+1, &
@@ -523 +521 @@
     use ioipsl
     USE phys_cal_mod
+    USE time_phylmdz_mod, ONLY: day_ref, annee_ref
     USE print_control_mod, ONLY: lunout
 
@@ -532 +531 @@
     real                :: ttt,xxx,timestep,dayseconde,dtime
     parameter (dayseconde=86400.)
-    include "temps.h"
     include "comconst.h"