From 4484abb18309a5e51440e3f962a43920c5362689 Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Thu, 3 Feb 2022 15:18:00 +0100
Subject: [PATCH] Philippe 03/02/2022: ini_nsv: correct some units
---
src/MNH/ini_nsv.f90 | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/src/MNH/ini_nsv.f90 b/src/MNH/ini_nsv.f90
index 01d0547f9..b117f3edb 100644
--- a/src/MNH/ini_nsv.f90
+++ b/src/MNH/ini_nsv.f90
@@ -351,7 +351,6 @@ ELSE
END IF
!
#ifdef MNH_FOREFIRE
-
! ForeFire tracers
IF (LFOREFIRE .AND. NFFSCALARS .GT. 0) THEN
NSV_FF_A(KMI) = NFFSCALARS
@@ -903,6 +902,7 @@ DO JSV = NSV_LIMA_BEG_A(KMI), NSV_LIMA_END_A(KMI)
TSVLIST_A(JSV, KMI)%CMNHNAME = TRIM( CLIMA_COLD_NAMES(5) )
ELSE IF ( JSV == NSV_LIMA_SPRO_A(KMI) ) THEN
CSVNAMES_A(JSV,KMI) = TRIM( CLIMA_WARM_NAMES(5) )
+ TSVLIST_A(JSV, KMI)%CUNITS = '1'
TSVLIST_A(JSV, KMI)%CMNHNAME = TRIM( CLIMA_WARM_NAMES(5) )
ELSE
CALL Print_msg( NVERB_FATAL, 'GEN', 'INI_NSV', 'invalid index for LIMA' )
@@ -915,11 +915,11 @@ DO JSV = NSV_ELECBEG_A(KMI), NSV_ELECEND_A(KMI)
CSVNAMES_A(JSV,KMI) = TRIM( CELECNAMES(JSV-NSV_ELECBEG_A(KMI)+1) )
IF ( JSV > NSV_ELECBEG .AND. JSV < NSV_ELECEND ) THEN
- YUNITS = 'C m-3'
+ YUNITS = 'C kg-1'
WRITE( YCOMMENT, '( A6, A3, I3.3 )' ) 'X_Y_Z_', 'SVT', JSV
ELSE
- YUNITS = 'm-3'
- WRITE( YCOMMENT, '( A6, A3, I3.3, A8 )' ) 'X_Y_Z_', 'SVT', JSV, ' (nb ions/m3)'
+ YUNITS = 'kg-1'
+ WRITE( YCOMMENT, '( A6, A3, I3.3, A8 )' ) 'X_Y_Z_', 'SVT', JSV, ' (nb ions/kg)'
END IF
TSVLIST_A(JSV, KMI) = TFIELDMETADATA( &
--
GitLab