Changeset 2932 for trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F

Timestamp:

Apr 3, 2023, 6:17:05 PM (20 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/vdifc_mod.F (modified) (3 diffs)

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

@@ -28 +28 @@
 c      use geometry_mod, only: longitude_deg,latitude_deg !  Joseph
       use dust_param_mod, only: doubleq, submicron, lifting
+      use write_output_mod, only: write_output
 
       IMPLICIT NONE
@@ -395 +396 @@
 ! Some usefull diagnostics for the new surface layer parametrization :
 
-!        call WRITEDIAGFI(ngrid,'vdifc_zcdv_true',
+!        call write_output('vdifc_zcdv_true',
 !     &              'momentum drag','no units',
-!     &                         2,zcdv_true)
-!        call WRITEDIAGFI(ngrid,'vdifc_zcdh_true',
+!     &                         zcdv_true(:))
+!        call write_output('vdifc_zcdh_true',
 !     &              'heat drag','no units',
-!     &                         2,zcdh_true)
-!        call WRITEDIAGFI(ngrid,'vdifc_ust',
-!     &              'friction velocity','m/s',2,ust)
-!       call WRITEDIAGFI(ngrid,'vdifc_tst',
-!     &              'friction temperature','K',2,tst)
-!        call WRITEDIAGFI(ngrid,'vdifc_zcdv',
+!     &                         zcdh_true(:))
+!        call write_output('vdifc_ust',
+!     &              'friction velocity','m/s',ust(:))
+!       call write_output('vdifc_tst',
+!     &              'friction temperature','K',tst(:))
+!        call write_output('vdifc_zcdv',
 !     &              'aerodyn momentum conductance','m/s',
-!     &                         2,zcdv)
-!        call WRITEDIAGFI(ngrid,'vdifc_zcdh',
+!     &                         zcdv(:))
+!        call write_output('vdifc_zcdh',
 !     &              'aerodyn heat conductance','m/s',
-!     &                         2,zcdh)
+!     &                         zcdh(:))
 
 c    ** schema de diffusion turbulente dans la couche limite
@@ -1057 +1058 @@
 c --------- end of tracers  ----------------------------
 
-         call WRITEDIAGFI(ngrid,"surf_h2o_lh",
+         call write_output("surf_h2o_lh",
      &                          "Ground ice latent heat flux",
-     &                               "W.m-2",2,surf_h2o_lh(:))
+     &                               "W.m-2",surf_h2o_lh(:))
 
 C       Diagnostic output for HDO
         if (hdo) then
-            CALL WRITEDIAGFI(ngrid,'hdoflux',
+            CALL write_output('hdoflux',
      &                       'hdoflux',
-     &                       ' ',2,hdoflux)  
-            CALL WRITEDIAGFI(ngrid,'h2oflux',
+     &                       ' ',hdoflux(:))  
+            CALL write_output('h2oflux',
      &                       'h2oflux',
-     &                       ' ',2,h2oflux)
+     &                       ' ',h2oflux(:))
         endif