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

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/watercloud_mod.F (modified) (3 diffs)

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

@@ -27 +27 @@
      &                      qperemin
       use dimradmars_mod, only: naerkind
+      use write_output_mod, only: write_output
       IMPLICIT NONE
 
@@ -268 +269 @@
      &                           cloudfrac(:,:)
       !! CLFvarying outputs
-      CALL WRITEDIAGFI(ngrid,'pqeffice','pqeffice',
-     &             'kg/kg',3,pqeff(:,:,igcm_h2o_ice))
-      CALL WRITEDIAGFI(ngrid,'pteff','pteff',
-     &             'K',3,pteff(:,:))
-      CALL WRITEDIAGFI(ngrid,'tcond','tcond',
-     &             'K',3,tcond(:,:))
-      CALL WRITEDIAGFI(ngrid,'cloudfrac','cloudfrac',
-     &             'K',3,cloudfrac(:,:))
+      CALL write_output('pqeffice','pqeffice',
+     &             'kg/kg',pqeff(:,:,igcm_h2o_ice))
+      CALL write_output('pteff','pteff',
+     &             'K',pteff(:,:))
+      CALL write_output('tcond','tcond',
+     &             'K',tcond(:,:))
+      CALL write_output('cloudfrac','cloudfrac',
+     &             'K',cloudfrac(:,:))
       END IF ! end if (CLFvarying)
 c------------------------------------------------------------------
@@ -713 +714 @@
 #ifndef MESOSCALE
 c=======================================================================
-      call WRITEDIAGFI(ngrid,"pdqice2","pdqcloudice apres microphysique"
-     &      ,"kg/kg.s-1",3,pdqcloud(1:ngrid,1:nlay,igcm_h2o_ice))
-      call WRITEDIAGFI(ngrid,"pdqvap2","pdqcloudvap apres microphysique"
-     &      ,"kg/kg.s-1",3,pdqcloud(1:ngrid,1:nlay,
+      call write_output("pdqice2","pdqcloudice apres microphysique"
+     &      ,"kg/kg.s-1",pdqcloud(:,:,igcm_h2o_ice))
+      call write_output("pdqvap2","pdqcloudvap apres microphysique"
+     &      ,"kg/kg.s-1",pdqcloud(:,:,
      &      igcm_h2o_vap))
       if (hdo) then
-      call WRITEDIAGFI(ngrid,"pdqiceD","pdqiceD apres microphysique"
-     &      ,"kg/kg.s-1",3,pdqcloud(1:ngrid,1:nlay,igcm_hdo_ice))
-      call WRITEDIAGFI(ngrid,"pdqvapD","pdqvapD apres microphysique"
-     &      ,"kg/kg.s-1",3,pdqcloud(1:ngrid,1:nlay,
+      call write_output("pdqiceD","pdqiceD apres microphysique"
+     &      ,"kg/kg.s-1",pdqcloud(:,:,igcm_hdo_ice))
+      call write_output("pdqvapD","pdqvapD apres microphysique"
+     &      ,"kg/kg.s-1",pdqcloud(:,:,
      &      igcm_hdo_vap))
       endif
-      call WRITEDIAGFI(ngrid,"pdqccn2","pdqcloudccn apres microphysique"
-     &      ,"kg/kg.s-1",3,pdqcloud(1:ngrid,1:nlay,
+      call write_output("pdqccn2","pdqcloudccn apres microphysique"
+     &      ,"kg/kg.s-1",pdqcloud(:,:,
      &      igcm_ccn_mass)) 
-      call WRITEDIAGFI(ngrid,"pdqccnN2","pdqcloudccnN apres "//
-     &      "microphysique","nb/kg.s-1",3,pdqcloud(1:ngrid,1:nlay,
+      call write_output("pdqccnN2","pdqcloudccnN apres "//
+     &      "microphysique","nb/kg.s-1",pdqcloud(:,:,
      &      igcm_ccn_number))
-      call WRITEDIAGFI(ngrid,"pdqdust2", "pdqclouddust apres "//
-     &      "microphysique","kg/kg.s-1",3,pdqcloud(1:ngrid,1:nlay,
+      call write_output("pdqdust2", "pdqclouddust apres "//
+     &      "microphysique","kg/kg.s-1",pdqcloud(:,:,
      &      igcm_dust_mass))
-      call WRITEDIAGFI(ngrid,"pdqdustN2", "pdqclouddustN apres "//
-     &      "microphysique","nb/kg.s-1",3,pdqcloud(1:ngrid,1:nlay,
+      call write_output("pdqdustN2", "pdqclouddustN apres "//
+     &      "microphysique","nb/kg.s-1",pdqcloud(:,:,
      &      igcm_dust_number)) 
 c=======================================================================