Skip to content
Snippets Groups Projects
Commit 609b6c04 authored by WAUTELET Philippe's avatar WAUTELET Philippe
Browse files

Philippe + Pierre T 06/12/2023: write_lfin: improve metadata for aerosol fluxes

parent 461e264e
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 65 additions and 66 deletions
src/MNH/write_lfin.f90
+ 65
− 66
View file @ 609b6c04
...@@ -223,7 +223,6 @@ USE MODD_EOL_SHARED_IO ...@@ -223,7 +223,6 @@ USE MODD_EOL_SHARED_IO
USE MODD_FIELD_n USE MODD_FIELD_n
use modd_field, only: NMNHDIM_UNUSED, tfieldmetadata, tfieldlist, NMNHDIM_NI, NMNHDIM_NJ, NMNHDIM_NOTLISTED, & use modd_field, only: NMNHDIM_UNUSED, tfieldmetadata, tfieldlist, NMNHDIM_NI, NMNHDIM_NJ, NMNHDIM_NOTLISTED, &
TYPECHAR, TYPEDATE, TYPEINT, TYPELOG, TYPEREAL TYPECHAR, TYPEDATE, TYPEINT, TYPELOG, TYPEREAL
use mode_field, only: Find_field_id_from_mnhname
USE MODD_FIRE_n USE MODD_FIRE_n
#ifdef MNH_FOREFIRE #ifdef MNH_FOREFIRE
USE MODD_FOREFIRE USE MODD_FOREFIRE
...@@ -1132,30 +1131,30 @@ IF (NSV >= 1 ) THEN ...@@ -1132,30 +1131,30 @@ IF (NSV >= 1 ) THEN
IF ((LORILAM).AND.(CPROGRAM == 'MESONH')) THEN IF ((LORILAM).AND.(CPROGRAM == 'MESONH')) THEN
DO JSV = 1 , NSV_AER DO JSV = 1 , NSV_AER
TZFIELD = TFIELDMETADATA( & TZFIELD = TFIELDMETADATA( &
CMNHNAME = 'FLX_'//TRIM(UPCASE(CAERONAMES(JSV))), & CMNHNAME = 'FLX_'//TRIM(UPCASE(CAERONAMES(JSV))), &
CSTDNAME = '', & CSTDNAME = '', &
CLONGNAME = 'FLX_'//TRIM(UPCASE(CAERONAMES(JSV))), & CLONGNAME = 'FLX_'//TRIM(UPCASE(CAERONAMES(JSV))), &
CUNITS = 'kg m-2 s-1', & CUNITS = 'kg m-2 s-1', &
CDIR = 'XY', & CDIR = 'XY', &
CCOMMENT = 'Aerosols mass flux', & CCOMMENT = 'Aerosol mass flux', &
NGRID = 1, & NGRID = 1, &
NTYPE = TYPEREAL, & NTYPE = TYPEREAL, &
NDIMS = 2, & NDIMS = 2, &
LTIMEDEP = .TRUE. ) LTIMEDEP = .TRUE. )
CALL IO_Field_write(TPFILE,TZFIELD,XFLX_AER(:,:,JSV)) CALL IO_Field_write(TPFILE,TZFIELD,XFLX_AER(:,:,JSV))
TZFIELD = TFIELDMETADATA( & TZFIELD = TFIELDMETADATA( &
CMNHNAME = 'FLXT_'//TRIM(UPCASE(CAERONAMES(JSV))), & CMNHNAME = 'FLXT_'//TRIM(UPCASE(CAERONAMES(JSV))), &
CSTDNAME = '', & CSTDNAME = '', &
CLONGNAME = 'FLXT_'//TRIM(UPCASE(CAERONAMES(JSV))), & CLONGNAME = 'FLXT_'//TRIM(UPCASE(CAERONAMES(JSV))), &
CUNITS = 'kg m-2', & CUNITS = 'kg m-2', &
CDIR = 'XY', & CDIR = 'XY', &
CCOMMENT = 'Integrated aerosols flux since start/restart', & CCOMMENT = 'Integrated aerosol mass flux since start/restart', &
NGRID = 1, & NGRID = 1, &
NTYPE = TYPEREAL, & NTYPE = TYPEREAL, &
NDIMS = 2, & NDIMS = 2, &
LTIMEDEP = .TRUE. ) LTIMEDEP = .TRUE. )
CALL IO_Field_write(TPFILE,TZFIELD,XFLXT_AER(:,:,JSV)) CALL IO_Field_write(TPFILE,TZFIELD,XFLXT_AER(:,:,JSV))
END DO END DO
END IF END IF
...@@ -1164,59 +1163,59 @@ IF ((LSALT).AND.(CPROGRAM == 'MESONH')) THEN ...@@ -1164,59 +1163,59 @@ IF ((LSALT).AND.(CPROGRAM == 'MESONH')) THEN
DO JSV = 1 , NMODE_SLT DO JSV = 1 , NMODE_SLT
WRITE (YFRC1,'(I1.1)') JSV WRITE (YFRC1,'(I1.1)') JSV
TZFIELD = TFIELDMETADATA( & TZFIELD = TFIELDMETADATA( &
CMNHNAME = 'FLX_SLT'//YFRC1, & CMNHNAME = 'FLX_SLT'//YFRC1, &
CSTDNAME = '', & CSTDNAME = '', &
CLONGNAME = 'FLX_SLT'//YFRC1, & CLONGNAME = 'FLX_SLT'//YFRC1, &
CUNITS = 'part m-2 s-1', & CUNITS = 'm-2 s-1', &
CDIR = 'XY', & CDIR = 'XY', &
CCOMMENT = 'Sea salt mass flux', & CCOMMENT = 'Sea salt number flux', &
NGRID = 1, & NGRID = 1, &
NTYPE = TYPEREAL, & NTYPE = TYPEREAL, &
NDIMS = 2, & NDIMS = 2, &
LTIMEDEP = .TRUE. ) LTIMEDEP = .TRUE. )
CALL IO_Field_write(TPFILE,TZFIELD,XFLX_SLT(:,:,JSV)) CALL IO_Field_write(TPFILE,TZFIELD,XFLX_SLT(:,:,JSV))
TZFIELD = TFIELDMETADATA( & TZFIELD = TFIELDMETADATA( &
CMNHNAME = 'FLXT_SLT'//YFRC1, & CMNHNAME = 'FLXT_SLT'//YFRC1, &
CSTDNAME = '', & CSTDNAME = '', &
CLONGNAME = 'FLXT_SLT'//YFRC1, & CLONGNAME = 'FLXT_SLT'//YFRC1, &
CUNITS = 'part m-2', & CUNITS = 'm-2', &
CDIR = 'XY', & CDIR = 'XY', &
CCOMMENT = 'Integrated sea salt flux since start/restart', & CCOMMENT = 'Integrated sea salt number flux since start/restart', &
NGRID = 1, & NGRID = 1, &
NTYPE = TYPEREAL, & NTYPE = TYPEREAL, &
NDIMS = 2, & NDIMS = 2, &
LTIMEDEP = .TRUE. ) LTIMEDEP = .TRUE. )
CALL IO_Field_write(TPFILE,TZFIELD,XFLXT_SLT(:,:,JSV)) CALL IO_Field_write(TPFILE,TZFIELD,XFLXT_SLT(:,:,JSV))
END DO END DO
END IF END IF
IF ((LUSECHEM).AND.(CPROGRAM == 'MESONH')) THEN IF ((LUSECHEM).AND.(CPROGRAM == 'MESONH')) THEN
TZFIELD = TFIELDMETADATA( & TZFIELD = TFIELDMETADATA( &
CMNHNAME = 'FLX_DMS', & CMNHNAME = 'FLX_DMS', &
CSTDNAME = '', & CSTDNAME = '', &
CLONGNAME = 'FLX_DMS', & CLONGNAME = 'FLX_DMS', &
CUNITS = 'kg m-2 s-1', & CUNITS = 'kg m-2 s-1', &
CDIR = 'XY', & CDIR = 'XY', &
CCOMMENT = 'Sea salt mass flux', & CCOMMENT = 'Sea salt mass flux', &
NGRID = 1, & NGRID = 1, &
NTYPE = TYPEREAL, & NTYPE = TYPEREAL, &
NDIMS = 2, & NDIMS = 2, &
LTIMEDEP = .TRUE. ) LTIMEDEP = .TRUE. )
CALL IO_Field_write(TPFILE,TZFIELD,XFLX_DMS(:,:)) CALL IO_Field_write(TPFILE,TZFIELD,XFLX_DMS(:,:))
TZFIELD = TFIELDMETADATA( & TZFIELD = TFIELDMETADATA( &
CMNHNAME = 'FLXT_DMS', & CMNHNAME = 'FLXT_DMS', &
CSTDNAME = '', & CSTDNAME = '', &
CLONGNAME = 'FLXT_DMS', & CLONGNAME = 'FLXT_DMS', &
CUNITS = 'kg m-2', & CUNITS = 'kg m-2', &
CDIR = 'XY', & CDIR = 'XY', &
CCOMMENT = 'Integrated sea salt flux since start/restart', & CCOMMENT = 'Integrated sea salt mass flux since start/restart', &
NGRID = 1, & NGRID = 1, &
NTYPE = TYPEREAL, & NTYPE = TYPEREAL, &
NDIMS = 2, & NDIMS = 2, &
LTIMEDEP = .TRUE. ) LTIMEDEP = .TRUE. )
CALL IO_Field_write(TPFILE,TZFIELD,XFLXT_DMS(:,:)) CALL IO_Field_write(TPFILE,TZFIELD,XFLXT_DMS(:,:))
END IF END IF
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment

Documentation | Mentions légales | CGU | Accessibilité : non conforme