Changeset 4046 for LMDZ6/trunk/libf/phylmdiso/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/phylmdiso/phys_output_write_mod.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmdiso/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, &
+        nqtottr,itr_indice
+    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
@@ -383 +385 @@
     USE pbl_surface_mod, ONLY: snow
     USE indice_sol_mod, ONLY: nbsrf
-    USE infotrac_phy, ONLY: nqtot, nqo, type_trac, tname, niadv, &
-        nqtottr,itr_indice ! C Risi: ajout nqtottr
 #ifdef ISO
     USE infotrac_phy, only: ntraciso,niso,nqtottr
@@ -467 +467 @@
     INTEGER ISW
     CHARACTER*1 ch1
-    CHARACTER*20 varname
+    CHARACTER(LEN=maxlen) :: varname, dn
 
 #ifdef CPP_XIOS
@@ -538 +538 @@
           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 
@@ -587 +592 @@
           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