Changeset 2932 for trunk/LMDZ.MARS/libf/phymars/compute_dtau_mod.F90

Timestamp:

Apr 3, 2023, 6:17:05 PM (23 months ago)

Author:

romain.vande

Message:

Mars PCM:
Add a new routine to write output called write_output.
It needs to be imported (for example : use write_output_mod, only: write_output)
Then, it requires only 4 arguments : the name of the variable, its title, its units and the variable itself.
It detects the dimension of the variable and decide to output either in diagfi or diagsoil.
It is also compatible with XIOS (xml file needs to be adapted)
Writediagfi and writediagsoil routines are still available in case.
RV

File:

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

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

@@ -25 +25 @@
         USE dust_param_mod, only: odpref, t_scenario_sol
         USE read_dust_scenario_mod, only: read_dust_scenario
+        use write_output_mod, only: write_output
         
         IMPLICIT NONE
@@ -80 +81 @@
         endif
        ! for diagnostics
-        call WRITEDIAGFI(ngrid,"tau_pref_target", &
+        call write_output("tau_pref_target", &
                           "target visible dust opacity column at 610Pa", &
-                          "",2,tau_pref_target)
+                          "",tau_pref_target(:))
 
         ! 2. Compute dtau() and dustliftday()
@@ -102 +103 @@
 
        ! for diagnostics
-        call WRITEDIAGFI(ngrid,"dtau","opacity difference wrt scenario",&
-                          "",2,dtau)
-        call WRITEDIAGFI(ngrid,"dustliftday","dust injection rate",     &
-                          "s-1",2,dustliftday)
+        call write_output("dtau","opacity difference wrt scenario",&
+                          "",dtau(:))
+        call write_output("dustliftday","dust injection rate",     &
+                          "s-1",dustliftday(:))
          
         ! 4. Save local time