Changeset 4046 for LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

Timestamp:

Dec 15, 2021, 11:18:49 PM (2 years ago)

Author:

dcugnet

Message:

First commit for new tracers.

• parser routines readTracFiles, strings_mod and tracer_types added in misc using revision 4 of ​https://svn.lmd.jussieu.fr/tracers-parser
• tested in sequential and parallel mode using ioipsl.
• for now, only two fields of "tracers(:)" derived type vector are used: "name" and "longName".

File:

• LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90 (modified) (8 diffs)

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -25 +25 @@
 
     USE dimphy, ONLY: klon, klev, klevp1
-    USE infotrac_phy, ONLY: nbtr
+    USE infotrac_phy, ONLY: nbtr, nqtot, nqo, type_trac, tracers, niadv
+    USE strings_mod,  ONLY: maxlen
     USE mod_phys_lmdz_para, ONLY: is_north_pole_phy,is_south_pole_phy
     USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
@@ -366 +367 @@
     USE pbl_surface_mod, ONLY: snow
     USE indice_sol_mod, ONLY: nbsrf
-    USE infotrac_phy, ONLY: nqtot, nqo, type_trac, tname, niadv
     USE geometry_mod, ONLY: cell_area, latitude_deg, longitude_deg
     USE surface_data, ONLY: type_ocean, version_ocean, ok_veget, landice_opt
@@ -441 +441 @@
     INTEGER ISW
     CHARACTER*1 ch1
-    CHARACTER*20 varname
+    CHARACTER(LEN=maxlen) :: varname, dn
 
 #ifdef CPP_XIOS
@@ -454 +454 @@
 
     INTEGER :: nqup
+
 
     ! On calcul le nouveau tau:
@@ -512 +513 @@
           DO iq=nqo+1, nqtot
             iiq=niadv(iq)
-            varname=trim(tname(iiq))
-            WRITE (lunout,*) 'XIOS var=', nqo, iq, nqtot, varname
+            dn = 'd'//TRIM(tracers(iiq)%name)//'_'
+            WRITE (lunout,*) 'XIOS var=', nqo, iq, nqtot, tracers(iiq)%name
+
+            unt = "kg kg-1"
+            varname=trim(tracers(iiq)%name)
             CALL xios_add_child(group_handle, child, varname)
-            CALL xios_set_attr(child, name=varname, unit="kg kg-1")
-            varname='d'//trim(tname(iiq))//'_vdf'
+            CALL xios_set_attr(child, name=varname, unit=unt)
+
+            unt = "kg kg-1 s-1"
+            varname=TRIM(dn)//'vdf'
             CALL xios_add_child(group_handle, child, varname)
-            CALL xios_set_attr(child, name=varname, unit="kg kg-1 s-1")
-            varname='d'//trim(tname(iiq))//'_the'
+            CALL xios_set_attr(child, name=varname, unit=unt)
+            varname=TRIM(dn)//'the'
             CALL xios_add_child(group_handle, child, varname)
-            CALL xios_set_attr(child, name=varname, unit="kg kg-1 s-1")
-            varname='d'//trim(tname(iiq))//'_con'
+            CALL xios_set_attr(child, name=varname, unit=unt)
+            varname=TRIM(dn)//'con'
             CALL xios_add_child(group_handle, child, varname)
-            CALL xios_set_attr(child, name=varname, unit="kg kg-1 s-1")
-            varname='d'//trim(tname(iiq))//'_lessi_impa'
+            CALL xios_set_attr(child, name=varname, unit=unt)
+            varname=TRIM(dn)//'lessi_impa'
             CALL xios_add_child(group_handle, child, varname)
-            CALL xios_set_attr(child, name=varname, unit="kg kg-1 s-1")
-            varname='d'//trim(tname(iiq))//'_lessi_nucl'
+            CALL xios_set_attr(child, name=varname, unit=unt)
+            varname=TRIM(dn)//'lessi_nucl'
             CALL xios_add_child(group_handle, child, varname)
-            CALL xios_set_attr(child, name=varname, unit="kg kg-1 s-1")
-            varname='d'//trim(tname(iiq))//'_insc'
+            CALL xios_set_attr(child, name=varname, unit=unt)
+            varname=TRIM(dn)//'insc'
             CALL xios_add_child(group_handle, child, varname)
-            CALL xios_set_attr(child, name=varname, unit="kg kg-1 s-1")
-            varname='d'//trim(tname(iiq))//'_bcscav'
+            CALL xios_set_attr(child, name=varname, unit=unt)
+            varname=TRIM(dn)//'bcscav'
             CALL xios_add_child(group_handle, child, varname)
-            CALL xios_set_attr(child, name=varname, unit="kg kg-1 s-1")
-            varname='d'//trim(tname(iiq))//'_evapls'
+            CALL xios_set_attr(child, name=varname, unit=unt)
+            varname=TRIM(dn)//'evapls'
             CALL xios_add_child(group_handle, child, varname)
-            CALL xios_set_attr(child, name=varname, unit="kg kg-1 s-1")
-            varname='d'//trim(tname(iiq))//'_ls'
+            CALL xios_set_attr(child, name=varname, unit=unt)
+            varname=TRIM(dn)//'ls'
             CALL xios_add_child(group_handle, child, varname)
-            CALL xios_set_attr(child, name=varname, unit="kg kg-1 s-1")
-            varname='d'//trim(tname(iiq))//'_trsp'
+            CALL xios_set_attr(child, name=varname, unit=unt)
+            varname=TRIM(dn)//'trsp'
             CALL xios_add_child(group_handle, child, varname)
-            CALL xios_set_attr(child, name=varname, unit="kg kg-1 s-1")
-            varname='d'//trim(tname(iiq))//'_sscav'
+            CALL xios_set_attr(child, name=varname, unit=unt)
+            varname=TRIM(dn)//'sscav'
             CALL xios_add_child(group_handle, child, varname)
-            CALL xios_set_attr(child, name=varname, unit="kg kg-1 s-1")
-            varname='d'//trim(tname(iiq))//'_sat'
+            CALL xios_set_attr(child, name=varname, unit=unt)
+            varname=TRIM(dn)//'sat'
             CALL xios_add_child(group_handle, child, varname)
-            CALL xios_set_attr(child, name=varname, unit="kg kg-1 s-1")
-            varname='d'//trim(tname(iiq))//'_uscav'
+            CALL xios_set_attr(child, name=varname, unit=unt)
+            varname=TRIM(dn)//'uscav'
             CALL xios_add_child(group_handle, child, varname)
-            CALL xios_set_attr(child, name=varname, unit="kg kg-1 s-1")
+            CALL xios_set_attr(child, name=varname, unit=unt)
           ENDDO
           !On ajoute les variables 2D traceurs par l interface fortran 
@@ -561 +567 @@
           DO iq=nqo+1, nqtot
             iiq=niadv(iq)
-            varname='cum'//trim(tname(iiq))
+
+            unt = "kg m-2"
+            varname='cum'//trim(tracers(iiq)%name)
             WRITE (lunout,*) 'XIOS var=', iq, nqtot, varname
             CALL xios_add_child(group_handle, child, varname)
-            CALL xios_set_attr(child, name=varname, unit="kg m-2")
-            varname='cumd'//trim(tname(iiq))//'_dry'
+            CALL xios_set_attr(child, name=varname, unit=unt)
+
+            unt = "kg m-2 s-1"
+            varname='cumd'//trim(tracers(iiq)%name)//'_dry'
             CALL xios_add_child(group_handle, child, varname)
-            CALL xios_set_attr(child, name=varname, unit="kg m-2 s-1")
+            CALL xios_set_attr(child, name=varname, unit=unt)
           ENDDO
     ENDIF
@@ -573 +583 @@
 #endif
 #endif
+
     ! ug la boucle qui suit ne sert qu'une fois, pour l'initialisation, sinon il n'y a toujours qu'un seul passage:
     DO iinit=1, iinitend
@@ -600 +611 @@
        ENDIf
        CALL histwrite_phy(o_aire, zx_tmp_fi2d)
-
        IF (vars_defined) THEN
           DO i=1, klon