Changeset 3870 for trunk/LMDZ.MARS/libf/phymars/writediagmicrofi.F

Timestamp:

Jul 29, 2025, 3:51:19 PM (3 months ago)

Author:

jmauxion

Message:

Mars PCM:
In writediag* routines, the files are now opened and closed only at dump time instead of every physical step (sensible speed up in 1D mainly).
JM

File:

• trunk/LMDZ.MARS/libf/phymars/writediagmicrofi.F (modified) (5 diffs)

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

@@ -81 +81 @@
       REAL area((nbp_lon+1),nbp_lat)
 
-      integer isample
+      integer, save :: isample
       integer ierr,ierr2
       integer i,j,l, ig0
@@ -124 +124 @@
 #endif
 
-!***************************************************************
-! Compute the output rate
-
-      isample = steps_per_sol/outputs_per_sol
-
-!***************************************************************
-
 ! At very first call, check if there is a "diagfi.def" to use and read it
 ! -----------------------------------------------------------------------
       IF (firstcall) THEN
          firstcall=.false.
+         
+         ! Compute the output rate
+         isample = steps_per_sol/outputs_per_sol
 
 !$OMP MASTER
@@ -264 +260 @@
          ENDIF
 
+         ierr= NF_CLOSE(nid) ! Close the NETCDF file once initialized
+
          endif ! of if (is_master)
-
-      else
-
-         if (is_master) then
-           ! only the master is required to do this
-
-           ! Open the NetCDF file
-           ierr = NF_OPEN(fichnom,NF_WRITE,nid)
-         endif ! of if (is_master)
-
       endif ! if (firstnom.eq.'1234567890')
 
@@ -297 +285 @@
         if (is_master) then
            ! only the master is required to do this
+
+           ! Time to write data, open NETCDF file
+           ierr=NF_OPEN(fichnom,NF_WRITE,nid)
+
         if ((nom.eq.firstnom) .and. (microstep.eq.1)) then
         ! We have identified a "first call" (at given date)
@@ -669 +661 @@
         endif ! of if (dim.eq.3) elseif(dim.eq.2)...
 
+        if (is_master) then
+          ! Close until next variable to dump or next dump iteration
+          ierr= NF_CLOSE(nid)
+        endif
+
       endif ! of if ( MOD(zitau+1,isample) .eq.0.)
-
-      if (is_master) then
-        ierr= NF_CLOSE(nid)
-      endif
-
+      
       end