Changeset 4150 for trunk/LMDZ.MARS/libf/phymars/iostart.F90

Timestamp:

Mar 24, 2026, 5:27:57 PM (7 weeks ago)

Author:

yluo

Message:

Mars PCM:
Refactored the tracer mass-conservation fixer for dynamics introduced in r4142:
Added a runtime flag to enable/disable the correction of tracer mass non-conservation in dynamics.
Reorganized the implementation into a modular structure.
Restricted the correction to a subset of tracers: CO, O2, H2, HO2, H2O2, N2, Ar, and He.
Enabled correction for dynamics timesteps preceding the first physics timestep of a simulation.
YCL

File:

• trunk/LMDZ.MARS/libf/phymars/iostart.F90 (modified) (3 diffs)

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

@@ -8 +8 @@
     
     ! restartfi.nc file dimension identifiers: (see open_restartphy())
+    INTEGER,SAVE :: idim0  ! "scalar" dimension
     INTEGER,SAVE :: idim1  ! "index" dimension
     INTEGER,SAVE :: idim2  ! "physical_points" dimension
@@ -510 +511 @@
         write(*,*)'open_restartphy: problem writing title '
         write(*,*)trim(nf90_strerror(ierr))
+      ENDIF
+
+      ierr=NF90_DEF_DIM(nid_restart,"scalar",1,idim0)
+      IF (ierr/=NF90_NOERR) THEN
+        write(*,*)'open_restartphy: problem defining scalar dimension '
+        write(*,*)trim(nf90_strerror(ierr))
+        CALL abort_physic("open_restartphy","Failed defining scalar",1)
       ENDIF
 
@@ -1061 +1069 @@
         ENDIF
         return ! nothing left to do
+      ELSEIF (var_size==1) THEN
+        ! We know it is a scalar
+        idim1d=idim0
       ELSEIF (var_size==length) THEN
         ! We know it is a "controle" kind of 1D array