Changeset 4058 for trunk/LMDZ.MARS/libf/aeronomars

Timestamp:

Feb 10, 2026, 9:46:10 AM (8 days ago)

Author:

emillour

Message:

Mars PCM:
Add a "output_diagfi" (default .true.) flag so that the user can decide if
a diagfi.nc file should be outputted.
Also do some cleaning around the few remaining calls to writediagfi
or reference to it througout the code; write_output() should be used.
EM

Location:

trunk/LMDZ.MARS/libf/aeronomars

Files:

• calchim_mod.F90 (modified) (3 diffs)
• surfacearea.F (modified) (2 diffs)

trunk/LMDZ.MARS/libf/aeronomars/calchim_mod.F90

@@ -40 +40 @@
       use iono_h, only: temp_elect
       use wstats_mod, only: wstats
+      use write_output_mod, only: write_output
       use callkeys_mod, only: photochem
 
@@ -223 +224 @@
 !     for output:
 
-      logical,save :: output        ! to issue calls to writediagfi and stats
+      logical,save :: output        ! to issue calls to write_output and stats
       real :: jo3_3d(ngrid,nlayer)  ! Photodissociation rate O3->O1D (s-1)
       real :: jh2o_3d(ngrid,nlayer)  ! Photodissociation rate H2O->H+OH (s-1)
@@ -956 +957 @@
       if (photochem .and. output) then
          if (ngrid > 1) then
-            call writediagfi(ngrid,'jo3','j o3->o1d',    &
-                             's-1',3,jo3_3d(1,1))
-            call writediagfi(ngrid,'jh2o','jh2o',    &
-                             's-1',3,jh2o_3d(1,1))
-            call writediagfi(ngrid,'iter','iterations',  &
-                             ' ',3,iter_3d(1,1))
-
-!            if (.not. unichim) then
-               call writediagfi(ngrid,'emission_no',        &
-                    'NO nightglow emission rate','cm-3 s-1',3,emission_no)
-               call writediagfi(ngrid,'emission_o2',        &
-                    'O2 nightglow emission rate','cm-3 s-1',3,emission_o2)
-!            endif
+           call write_output('jo3','j o3->o1d','s-1',jo3_3d)
+           call write_output('jh2o','jh2o','s-1',jh2o_3d)
+           call write_output('iter_chemistry','iterations',' ',iter_3d)
+
+           call write_output('emission_no', &
+                             'NO nightglow emission rate','cm-3 s-1', &
+                             emission_no)
+           call write_output('emission_o2', &
+                             'O2 nightglow emission rate','cm-3 s-1', &
+                             emission_o2)
             
-               call wstats(ngrid,'jo3','j o3->o1d',       &
+           call wstats(ngrid,'jo3','j o3->o1d',       &
                            's-1',3,jo3_3d(1,1))
-               call wstats(ngrid,'emission_no',           &
+           call wstats(ngrid,'emission_no',           &
                    'NO nightglow emission rate','cm-3 s-1',3,emission_no)
-               call wstats(ngrid,'emission_o2',           &
+           call wstats(ngrid,'emission_o2',           &
                    'O2 nightglow emission rate','cm-3 s-1',3,emission_o2)
-               call wstats(ngrid,'mmean','mean molecular mass',       &
+           call wstats(ngrid,'mmean','mean molecular mass',       &
                            'g.mole-1',3,mmean(1,1))
          end if ! of if (ngrid.gt.1)

trunk/LMDZ.MARS/libf/aeronomars/surfacearea.F

@@ -16 +16 @@
       use comcstfi_h, only: pi
       use wstats_mod, only: wstats
+      use write_output_mod, only: write_output
       use callkeys_mod, only: microphys
       implicit none
@@ -126 +127 @@
      $            "micron2 cm-3",3,surfice*1.e6)
 
-      call writediagfi(ngrid,"surfdust", "Dust surface area",
-     $            "micron2 cm-3",3,surfdust*1.e6)
-      call writediagfi(ngrid,"surfice", "Ice cloud surface area",
-     $            "micron2 cm-3",3,surfice*1.e6)
+      call write_output("surfdust","Dust surface area",
+     &                  "micron2 cm-3",surfdust*1.e6)
+      call write_output("surfice","Ice cloud surface area",
+     &                  "micron2 cm-3",surfice*1.e6)
 
       end subroutine surfacearea