Changeset 2514 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

May 4, 2021, 10:23:43 AM (4 years ago)

Author:

romain.vande

Message:

Mars Dynamico : _ update of r2507: correct ztime_fin for dynamico

_ Test for no writting Diagfi for Dynamico
_ Correction for Xios output : Change of calendar : start_date=0 as a diagfi, timestep for output physic=1s, day_length=ndtphys, it is needed to have output in days in the xml to have the correct time

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• physiq_mod.F (modified) (1 diff)
• writediagfi.F (modified) (2 diffs)
• xios_output_mod.F90 (modified) (3 diffs)

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

@@ -2322 +2322 @@
      .                  - day_ini - time_phys
              ENDIF
+             if (ecritstart==0) then
+                ztime_fin = ztime_fin-(day_end-day_ini)
+             endif
 
           ELSE ! IF LMDZ

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

@@ -45 +45 @@
      &                               is_master, gather
       USE mod_grid_phy_lmdz, only : klon_glo, Grid1Dto2D_glo,
-     &                              nbp_lon, nbp_lat, nbp_lev
+     &                              nbp_lon, nbp_lat, nbp_lev,
+     &                              grid_type, unstructured
       implicit none
 
@@ -114 +115 @@
       real areafi_glo(ngrid) ! mesh area on global physics grid
 #endif
+
+      if (grid_type==unstructured) then
+           return
+      endif
 
 !***************************************************************

trunk/LMDZ.MARS/libf/phymars/xios_output_mod.F90

@@ -104 +104 @@
      write(lunout,*) "initialize_xios_output: build calendar"
     endif
-    timestep%second=dtphys
+    timestep%second=1
     if (nint(dtphys).ne.dtphys) then
       write(*,*) "initialize_xios_output: warning physics timestep is not an integer!"
-      timestep%second=nint(dtphys)
     endif
     if (nint(daysec).ne.daysec) then
@@ -116 +115 @@
     CALL xios_define_calendar(type="user_defined", &
                               timestep=timestep, &
-                              day_length=nint(daysec), &
+                              day_length=nint(daysec/dtphys), &
                               month_lengths=[61,66,66,65,60,54,50,46,47,47,51,56])
     !NB: it would make more sense to define months and their length in the
@@ -130 +129 @@
     ! Now define the start time of this simulation
     ! NB: we substract dtphys because we want to set the origin of the time axis
-    start_date=time_origin+xios_duration(0,0,day,0,0,timeofday*daysec-dtphys)
+    start_date=time_origin
     call xios_set_start_date(start_date=start_date)
     if (prt_level>=10) then