Changeset 2932 for trunk/LMDZ.MARS/libf/phymars/improvedclouds_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/improvedclouds_mod.F (modified) (2 diffs)

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

@@ -19 +19 @@
       use conc_mod, only: mmean
       use comcstfi_h, only: pi, cpp
+      use write_output_mod, only: write_output
       implicit none
       
@@ -540 +541 @@
 #ifndef MESOSCALE
 !      IF (ngrid.ne.1) THEN ! 3D
-!         call WRITEDIAGFI(ngrid,"satu","ratio saturation","",3,
+!         call write_output("satu","ratio saturation","",
 !     &                    satu_out)
-!         call WRITEDIAGFI(ngrid,"dM","ccn variation","kg/kg",3,
+!         call write_output("dM","ccn variation","kg/kg",
 !     &                    dM_out)
-!         call WRITEDIAGFI(ngrid,"dN","ccn variation","#",3,
+!         call write_output("dN","ccn variation","#",
 !     &                    dN_out)
-!         call WRITEDIAGFI(ngrid,"error","dichotomy max error","%",2,
+!         call write_output("error","dichotomy max error","%",
 !     &                    error2d)
-!         call WRITEDIAGFI(ngrid,"zqsat","zqsat","kg",3,
+!         call write_output("zqsat","zqsat","kg",
 !     &                    zqsat)
 !      ENDIF
 
 !      IF (ngrid.eq.1) THEN ! 1D
-!         call WRITEDIAGFI(ngrid,"error","incertitude sur glace","%",1,
+!         call write_output("error","incertitude sur glace","%",
 !     &                    error_out)
-         call WRITEdiagfi(ngrid,"resist","resistance","s/m2",1,
-     &                    res_out)
-         call WRITEdiagfi(ngrid,"satu_bf","satu before","kg/kg",1,
-     &                    satubf)
-         call WRITEdiagfi(ngrid,"satu_af","satu after","kg/kg",1,
-     &                    satuaf)
-         call WRITEdiagfi(ngrid,"vapbf","h2ovap before","kg/kg",1,
-     &                    zq0(1,1,igcm_h2o_vap))
-         call WRITEdiagfi(ngrid,"vapaf","h2ovap after","kg/kg",1,
-     &                    zq(1,1,igcm_h2o_vap))
-         call WRITEdiagfi(ngrid,"icebf","h2oice before","kg/kg",1,
-     &                    zq0(1,1,igcm_h2o_ice))
-         call WRITEdiagfi(ngrid,"iceaf","h2oice after","kg/kg",1,
-     &                    zq(1,1,igcm_h2o_ice))
-         call WRITEdiagfi(ngrid,"ccnbf","ccn before","/kg",1,
-     &                    zq0(1,1,igcm_ccn_number))
-         call WRITEdiagfi(ngrid,"ccnaf","ccn after","/kg",1,
-     &                    zq(1,1,igcm_ccn_number))
-c         call WRITEDIAGFI(ngrid,"growthrate","growth rate","m^2/s",1,
+         call write_output("resist","resistance","s/m2",
+     &                    res_out(:,:))
+         call write_output("satu_bf","satu before","kg/kg",
+     &                    satubf(:,:))
+         call write_output("satu_af","satu after","kg/kg",
+     &                    satuaf(:,:))
+         call write_output("vapbf","h2ovap before","kg/kg",
+     &                    zq0(:,:,igcm_h2o_vap))
+         call write_output("vapaf","h2ovap after","kg/kg",
+     &                    zq(:,:,igcm_h2o_vap))
+         call write_output("icebf","h2oice before","kg/kg",
+     &                    zq0(:,:,igcm_h2o_ice))
+         call write_output("iceaf","h2oice after","kg/kg",
+     &                    zq(:,:,igcm_h2o_ice))
+         call write_output("ccnbf","ccn before","/kg",
+     &                    zq0(:,:,igcm_ccn_number))
+         call write_output("ccnaf","ccn after","/kg",
+     &                    zq(:,:,igcm_ccn_number))
+c         call write_output("growthrate","growth rate","m^2/s",
 c     &                    gr_out)
-c         call WRITEDIAGFI(ngrid,"nuclearate","nucleation rate","",1,
+c         call write_output("nuclearate","nucleation rate","",
 c     &                    rate_out)
-c         call WRITEDIAGFI(ngrid,"dM","ccn variation","kg",1,
+c         call write_output("dM","ccn variation","kg",
 c     &                    dM_out)
-c         call WRITEDIAGFI(ngrid,"dN","ccn variation","#",1,
+c         call write_output("dN","ccn variation","#",
 c     &                    dN_out)
-         call WRITEdiagfi(ngrid,"zqsat","p vap sat","kg/kg",1,
-     &                    zqsat)
-!         call WRITEDIAGFI(ngrid,"satu","ratio saturation","",1,
-!     &                    satu_out)
-         call WRITEdiagfi(ngrid,"rice","ice radius","m",1,
-     &                    rice)
-!         call WRITEDIAGFI(ngrid,"rdust_sca","rdust","m",1,
+         call write_output("zqsat","p vap sat","kg/kg",
+     &                    zqsat(:,:))
+!         call write_output("satu","ratio saturation","",
+!     &                    satu_out(:,:))
+         call write_output("rice","ice radius","m",
+     &                    rice(:,:))
+!         call write_output("rdust_sca","rdust","m",
 !     &                    rdust)
-!         call WRITEDIAGFI(ngrid,"rsedcloud","rsedcloud","m",1,
+!         call write_output("rsedcloud","rsedcloud","m",
 !     &                    rsedcloud)
-!         call WRITEDIAGFI(ngrid,"rhocloud","rhocloud","kg.m-3",1,
+!         call write_output("rhocloud","rhocloud","kg.m-3",
 !     &                    rhocloud)
 !      ENDIF