Changeset 2932 for trunk/LMDZ.MARS/libf/phymars/nonoro_gwd_ran_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/nonoro_gwd_ran_mod.F90 (modified) (7 diffs)

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

@@ -47 +47 @@
       use vertical_layers_mod, only : presnivs
       use geometry_mod, only: cell_area
+      use write_output_mod, only: write_output
 #ifdef CPP_XIOS
      use xios_output_mod, only: send_xios_field
@@ -222 +223 @@
     PH(:, 1) = 2. * PP(:, 1) - PH(:, 2)
 
-    call writediagfi(ngrid,'nonoro_pp','nonoro_pp', 'm',3,PP(:,1:nlayer))
-    call writediagfi(ngrid,'nonoro_ph','nonoro_ph', 'm',3,PH(:,1:nlayer))
+    call write_output('nonoro_pp','nonoro_pp', 'm',PP(:,:))
+    call write_output('nonoro_ph','nonoro_ph', 'm',PH(:,:))
 
     ! Launching level for reproductible case
@@ -250 +251 @@
         ZH(:, 1) = H0 * LOG(PR / (PH(:, 1) + PSEC))
 
-    call writediagfi(ngrid,'nonoro_zh','nonoro_zh', 'm',3,ZH(:,2:nlayer+1))
+    call write_output('nonoro_zh','nonoro_zh', 'm',ZH(:,2:nlayer+1))
 
     ! Winds and Brunt Vaisala frequency
@@ -275 +276 @@
        MAX_K(II)=MAX(kmin,kstar)
     ENDDO
-    call writediagfi(ngrid,'nonoro_bv','nonoro_bv', 'm',3,BV(:,2:nlayer+1))
+    call write_output('nonoro_bv','nonoro_bv', 'm',BV(:,2:nlayer+1))
 
 !-----------------------------------------------------------------------------------------------------------------------
@@ -444 +445 @@
     ! (upward) that should be compensated by circulation 
     ! and induce additional friction at the surface 
-    call writediagfi(ngrid,'nonoro_u_epflux_tot','nonoro_u_epflux_tot', '',3,u_epflux_tot(:,2:nlayer+1))
-    call writediagfi(ngrid,'nonoro_v_epflux_tot','nonoro_v_epflux_tot', '',3,v_epflux_tot(:,2:nlayer+1))
+    call write_output('nonoro_u_epflux_tot','nonoro_u_epflux_tot', '',u_epflux_tot(:,2:nlayer+1))
+    call write_output('nonoro_v_epflux_tot','nonoro_v_epflux_tot', '',v_epflux_tot(:,2:nlayer+1))
 
     ! 5.2 AR-1 RECURSIVE FORMULA (13) IN VERSION 4
@@ -457 +458 @@
     ENDDO
     d_t(:,:) = 0.
-    call writediagfi(ngrid,'nonoro_d_u','nonoro_d_u', '',3,d_u)
-    call writediagfi(ngrid,'nonoro_d_v','nonoro_d_v', '',3,d_v)
+    call write_output('nonoro_d_u','nonoro_d_u', '',d_u(:,:))
+    call write_output('nonoro_d_v','nonoro_d_v', '',d_v(:,:))
 
     ! 5.3 Update tendency of wind with the previous (and saved) values
@@ -469 +470 @@
     dv_nonoro_gwd(:,:) = d_v(:,:)
 
-    call writediagfi(ngrid,'du_nonoro_gwd','du_nonoro_gwd', '',3,du_nonoro_gwd)
-    call writediagfi(ngrid,'dv_nonoro_gwd','dv_nonoro_gwd', '',3,dv_nonoro_gwd)
+    call write_output('du_nonoro_gwd','du_nonoro_gwd', '',du_nonoro_gwd(:,:))
+    call write_output('dv_nonoro_gwd','dv_nonoro_gwd', '',dv_nonoro_gwd(:,:))
     
     ! Cosmetic: evaluation of the cumulated stress