Philippe 01/12/2023: stations: add temperature field (field already existing...

Philippe 01/12/2023: stations: add temperature field (field already existing but not computed nor written)

src/MNH/stationn.f90

@@ -77,7 +77,7 @@ END MODULE MODI_STATION_n
 !  P. Wautelet  13/09/2019: budget: simplify and modernize date/time management
 !  R. Schoetter    11/2019: use LCARTESIAN instead of LSTATLAT for multiproc in cartesian
 !  P. Wautelet     04/2022: restructure stations for better performance, reduce memory usage and correct some problems/bugs
-!
+!  P. Wautelet  01/12/2023: compute temperature (XT field)
 ! --------------------------------------------------------------------------
 !
 !*      0. DECLARATIONS
@@ -85,7 +85,7 @@ END MODULE MODI_STATION_n
 !
 USE MODD_ALLSTATION_n,  ONLY: LDIAG_SURFRAD_STAT
 USE MODD_CONF,          ONLY: LCARTESIAN
-USE MODD_CST,           ONLY: XPI
+USE MODD_CST,           ONLY: XCPD, XP00, XPI, XRD
 USE MODD_GRID,          ONLY: XBETA, XLON0, XRPK
 USE MODD_PARAMETERS,    ONLY: JPVEXT
 USE MODD_PARAM_n,       ONLY: CRAD
@@ -158,6 +158,7 @@ STATION: DO JS = 1, NUMBSTAT_LOC
   TSTATIONS(JS)%XW (1,IN) = TSTATIONS(JS)%INTERP_HOR_FROM_MASSPOINT( PW(:,:,JK) )
   TSTATIONS(JS)%XTH(1,IN) = TSTATIONS(JS)%INTERP_HOR_FROM_MASSPOINT( PTH(:,:,JK) )
   TSTATIONS(JS)%XP (1,IN) = TSTATIONS(JS)%INTERP_HOR_FROM_MASSPOINT( PP(:,:,JK) )
+  TSTATIONS(JS)%XT (1,IN) = TSTATIONS(JS)%XTH(1,IN) * ( TSTATIONS(JS)%XP (1,IN) / XP00 ) **( XRD / XCPD )
 
   DO JSV=1,SIZE(PR,4)
     TSTATIONS(JS)%XR(1,IN,JSV) = TSTATIONS(JS)%INTERP_HOR_FROM_MASSPOINT( PR(:,:,JK,JSV) )

src/MNH/write_stationn.f90

@@ -13,6 +13,7 @@
 !  P. Wautelet 03/03/2021: budgets: add tbudiachrometadata type (useful to pass more information to Write_diachro)
 !  P. Wautelet 04/02/2022: use TSVLIST to manage metadata of scalar variables
 !  P. Wautelet    04/2022: restructure stations for better performance, reduce memory usage and correct some problems/bugs
+!  P. Wautelet 01/12/2023: write temperature field
 ! --------------------------------------------------------------------------
 !      ###########################
 MODULE MODE_WRITE_STATION_n
@@ -70,7 +71,7 @@ type(tfieldmetadata_base), dimension(:), allocatable :: tzfields
 !
 !----------------------------------------------------------------------------
 !
-IPROC = 6 + SIZE(TPSTATION%XR,3) + SIZE(TPSTATION%XSV,3)
+IPROC = 7 + SIZE(TPSTATION%XR,3) + SIZE(TPSTATION%XSV,3)
 
 IF ( CTURB == 'TKEL' ) IPROC = IPROC + 1
 IF (LDIAG_SURFRAD_STAT) THEN
@@ -99,8 +100,9 @@ else
   call Add_point( 'MER_WIND', 'Meridional wind', 'm s-1', tpstation%xmer(1,:) )
 end if
 
-call Add_point( 'W',  'Air vertical speed',    'm s-1', tpstation%xw(1,:)  )
+call Add_point( 'W',  'Air vertical speed',    'm s-1', tpstation%xw (1,:) )
 call Add_point( 'Th', 'Potential temperature', 'K',     tpstation%xth(1,:) )
+call Add_point( 'T',  'Temperature',           'K',     tpstation%xt (1,:) )
 
 if ( ldiag_surfrad_stat ) call Add_diag_surfrad_data( tpstation )