Changeset 2932 for trunk/LMDZ.MARS/libf/phymars/co2condens_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/co2condens_mod.F (modified) (4 diffs)

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

@@ -30 +30 @@
        use comcstfi_h, only: cpp, g, r, pi 
        use dust_param_mod, only: freedust
-
+       use write_output_mod, only: write_output
 #ifndef MESOSCALE
        USE vertical_layers_mod, ONLY: ap, bp
@@ -332 +332 @@
      &                       rdust,zcondicea,zfallice,zdq_scav,pdqsc)
       endif
-             call WRITEdiagfi(ngrid,"zfallice",
+             call write_output("zfallice",
      &         "",
-     &         " ",2,zfallice(ngrid,1))
+     &         " ",zfallice(:,1))
       ELSE ! if co2 clouds
         condens_layer(:,:) = 0.
@@ -349 +349 @@
           piceco2(ig) = piceco2(ig) + zdqssed_co2(ig)*ptimestep
         ENDDO
-       call WRITEdiagfi(ngrid,"zfallice",
+       call write_output("zfallice",
      &         "",
-     &         " ",2,zdqssed_co2) ! otherwise we have not 1 day(1proc) = 1 day (24procs) test
+     &         " ",zdqssed_co2(:)) ! otherwise we have not 1 day(1proc) = 1 day (24procs) test
       ENDIF ! end of if co2clouds
 
-      call WRITEdiagfi(ngrid,"pdtc_atm",
+      call write_output("pdtc_atm",
      &         "temperature tendency due to CO2 condensation",
-     &         " ",3,pdtc) 
+     &         " ",pdtc(:,:)) 
      
-       call WRITEdiagfi(ngrid,"condens_layer",
+       call write_output("condens_layer",
      &         "",
-     &         " ",3,condens_layer)
+     &         " ",condens_layer(:,:))
 
 c     *************************
@@ -716 +716 @@
 c        END DO
 c     END DO
-c     call WRITEDIAGFI(ngrid,'tconda1',
+c     call write_output('tconda1',
 c    &'Taux de condensation CO2 atmospherique /Pa',
-c    & 'kg.m-2.Pa-1.s-1',3,tconda1)
-c     call WRITEDIAGFI(ngrid,'tconda2',
+c    & 'kg.m-2.Pa-1.s-1',tconda1)
+c     call write_output('tconda2',
 c    &'Taux de condensation CO2 atmospherique /m',
-c    & 'kg.m-3.s-1',3,tconda2)
+c    & 'kg.m-3.s-1',tconda2)
 
 ! output falling co2 ice in 1st layer:
-!      call WRITEDIAGFI(ngrid,'fallice',
+!      call write_output('fallice',
 !     &'Precipitation of co2 ice',
-!     & 'kg.m-2.s-1',2,zfallice(1,1))
+!     & 'kg.m-2.s-1',zfallice(1,1))
 
 #ifndef MESOSCALE