From aeb7ed585caba7b5ff042e9e25adeab040992aa2 Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Fri, 18 Aug 2017 15:01:53 +0200
Subject: [PATCH] Philippe 18/08/2017: IO: correction in fieldlist (errors in
 metadata + add more fields when DIAG)

---
 src/LIB/SURCOUCHE/src/mode_field.f90 | 19 +++++++------------
 1 file changed, 7 insertions(+), 12 deletions(-)

diff --git a/src/LIB/SURCOUCHE/src/mode_field.f90 b/src/LIB/SURCOUCHE/src/mode_field.f90
index deb862876..6f96aad63 100644
--- a/src/LIB/SURCOUCHE/src/mode_field.f90
+++ b/src/LIB/SURCOUCHE/src/mode_field.f90
@@ -439,7 +439,7 @@ TFIELDLIST(IDX)%CLONGNAME  = 'MesoNH: XHAT'
 TFIELDLIST(IDX)%CUNITS     = 'm'
 TFIELDLIST(IDX)%CDIR       = 'XX'
 TFIELDLIST(IDX)%CCOMMENT   = 'Position x in the conformal or cartesian plane'
-TFIELDLIST(IDX)%NGRID      = 2
+TFIELDLIST(IDX)%NGRID      = 4
 TFIELDLIST(IDX)%NTYPE      = TYPEREAL
 TFIELDLIST(IDX)%NDIMS      = 1
 ALLOCATE(TFIELDLIST(IDX)%TFIELD_X1D(IMODEL))
@@ -453,7 +453,7 @@ TFIELDLIST(IDX)%CLONGNAME  = 'MesoNH: YHAT'
 TFIELDLIST(IDX)%CUNITS     = 'm'
 TFIELDLIST(IDX)%CDIR       = 'YY'
 TFIELDLIST(IDX)%CCOMMENT   = 'Position y in the conformal or cartesian plane'
-TFIELDLIST(IDX)%NGRID      = 3
+TFIELDLIST(IDX)%NGRID      = 4
 TFIELDLIST(IDX)%NTYPE      = TYPEREAL
 TFIELDLIST(IDX)%NDIMS      = 1
 ALLOCATE(TFIELDLIST(IDX)%TFIELD_X1D(IMODEL))
@@ -481,7 +481,7 @@ TFIELDLIST(IDX)%CLONGNAME  = 'MesoNH: DXHAT'
 TFIELDLIST(IDX)%CUNITS     = 'm'
 TFIELDLIST(IDX)%CDIR       = 'XX'
 TFIELDLIST(IDX)%CCOMMENT   = 'Horizontal stretching in x'
-TFIELDLIST(IDX)%NGRID      = 2
+TFIELDLIST(IDX)%NGRID      = 4
 TFIELDLIST(IDX)%NTYPE      = TYPEREAL
 TFIELDLIST(IDX)%NDIMS      = 1
 ALLOCATE(TFIELDLIST(IDX)%TFIELD_X1D(IMODEL))
@@ -495,7 +495,7 @@ TFIELDLIST(IDX)%CLONGNAME  = 'MesoNH: DYHAT'
 TFIELDLIST(IDX)%CUNITS     = 'm'
 TFIELDLIST(IDX)%CDIR       = 'YY'
 TFIELDLIST(IDX)%CCOMMENT   = 'Horizontal stretching in y'
-TFIELDLIST(IDX)%NGRID      = 3
+TFIELDLIST(IDX)%NGRID      = 4
 TFIELDLIST(IDX)%NTYPE      = TYPEREAL
 TFIELDLIST(IDX)%NDIMS      = 1
 ALLOCATE(TFIELDLIST(IDX)%TFIELD_X1D(IMODEL))
@@ -1786,7 +1786,7 @@ TFIELDLIST(IDX)%NDIMS      = 0
 IDX = IDX+1
 !
 !
-IF (CPROGRAM == 'MESONH') THEN
+IF (TRIM(CPROGRAM) == 'MESONH' .OR. TRIM(CPROGRAM) == 'DIAG') THEN
 !
 IF(IDX>MAXFIELDS) CALL ERR_INI_FIELD_LIST()
 TFIELDLIST(IDX)%CMNHNAME   = 'US_PRES'
@@ -1961,11 +1961,6 @@ TFIELDLIST(IDX)%NDIMS      = 3
 ALLOCATE(TFIELDLIST(IDX)%TFIELD_X3D(IMODEL))
 IDX = IDX+1
 !
-END IF ! CPROGRAM=MESONH
-!
-!
-IF (CPROGRAM == 'MESONH' .OR. TRIM(CPROGRAM) == 'DIAG') THEN
-!
 IF(IDX>MAXFIELDS) CALL ERR_INI_FIELD_LIST()
 TFIELDLIST(IDX)%CMNHNAME   = 'CIT'
 TFIELDLIST(IDX)%CSTDNAME   = ''
@@ -2009,7 +2004,7 @@ ALLOCATE(TFIELDLIST(IDX)%TFIELD_X3D(IMODEL))
 IDX = IDX+1
 !
 !
-IF (CPROGRAM=='MESONH') THEN
+IF (TRIM(CPROGRAM)=='MESONH' .OR. TRIM(CPROGRAM)=='DIAG') THEN
 !
 IF(IDX>MAXFIELDS) CALL ERR_INI_FIELD_LIST()
 TFIELDLIST(IDX)%CMNHNAME   = 'DTHRAD'
@@ -2167,7 +2162,7 @@ TFIELDLIST(IDX)%NDIMS      = 2
 ALLOCATE(TFIELDLIST(IDX)%TFIELD_X2D(IMODEL))
 IDX = IDX+1
 !
-END IF !CPROGRAM=MESONH
+END IF !CPROGRAM=MESONH .OR. DIAG
 !
 !
 IF (TRIM(CPROGRAM) /= 'PGD' .AND. TRIM(CPROGRAM) /= 'NESPGD' .AND. TRIM(CPROGRAM) /= 'SPAWN') THEN
-- 
GitLab