From a1a06e22ec61e147db1dc678d45d61f94964816f Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Mon, 12 Mar 2018 13:51:11 +0100
Subject: [PATCH] Philippe 12/03/2018: IO: add CDIR='ZZ' for several vertical
fields
---
src/LIB/SURCOUCHE/src/mode_field.f90 | 6 +++---
src/LIB/SURCOUCHE/src/mode_netcdf.f90 | 4 +++-
2 files changed, 6 insertions(+), 4 deletions(-)
diff --git a/src/LIB/SURCOUCHE/src/mode_field.f90 b/src/LIB/SURCOUCHE/src/mode_field.f90
index fdcda06a5..3950453a2 100644
--- a/src/LIB/SURCOUCHE/src/mode_field.f90
+++ b/src/LIB/SURCOUCHE/src/mode_field.f90
@@ -544,7 +544,7 @@ TFIELDLIST(IDX)%CMNHNAME = 'ZHAT'
TFIELDLIST(IDX)%CSTDNAME = ''
TFIELDLIST(IDX)%CLONGNAME = 'ZHAT'
TFIELDLIST(IDX)%CUNITS = 'm'
-TFIELDLIST(IDX)%CDIR = '--'
+TFIELDLIST(IDX)%CDIR = 'ZZ'
TFIELDLIST(IDX)%CCOMMENT = 'Height level without orography'
TFIELDLIST(IDX)%NGRID = 4
TFIELDLIST(IDX)%NTYPE = TYPEREAL
@@ -1950,7 +1950,7 @@ TFIELDLIST(IDX)%CMNHNAME = 'RHOREFZ'
TFIELDLIST(IDX)%CSTDNAME = ''
TFIELDLIST(IDX)%CLONGNAME = 'RHOREFZ'
TFIELDLIST(IDX)%CUNITS = 'kg m-3'
-TFIELDLIST(IDX)%CDIR = '--'
+TFIELDLIST(IDX)%CDIR = 'ZZ'
TFIELDLIST(IDX)%CCOMMENT = 'rhodz for reference state without orography'
TFIELDLIST(IDX)%NGRID = 4
TFIELDLIST(IDX)%NTYPE = TYPEREAL
@@ -1963,7 +1963,7 @@ TFIELDLIST(IDX)%CMNHNAME = 'THVREFZ'
TFIELDLIST(IDX)%CSTDNAME = ''
TFIELDLIST(IDX)%CLONGNAME = 'THVREFZ'
TFIELDLIST(IDX)%CUNITS = 'K'
-TFIELDLIST(IDX)%CDIR = '--'
+TFIELDLIST(IDX)%CDIR = 'ZZ'
TFIELDLIST(IDX)%CCOMMENT = 'thetavz for reference state without orography'
TFIELDLIST(IDX)%NGRID = 4
TFIELDLIST(IDX)%NTYPE = TYPEREAL
diff --git a/src/LIB/SURCOUCHE/src/mode_netcdf.f90 b/src/LIB/SURCOUCHE/src/mode_netcdf.f90
index fa072c285..d17d039f4 100644
--- a/src/LIB/SURCOUCHE/src/mode_netcdf.f90
+++ b/src/LIB/SURCOUCHE/src/mode_netcdf.f90
@@ -926,8 +926,10 @@ DO JI=1,SIZE(KSHAPE)
IF (JI == 1) THEN
IF ( (YDIR == 'XX' .OR. YDIR == 'XY') .AND. KSHAPE(1)==NCOORDID(1,IGRID)%LEN) THEN
KVDIMS(1) = NCOORDID(1,IGRID)%ID
- ELSE IF ( YDIR == 'YY' .AND. KSHAPE(1)==NCOORDID(1,IGRID)%LEN) THEN
+ ELSE IF ( YDIR == 'YY' .AND. KSHAPE(1)==NCOORDID(2,IGRID)%LEN) THEN
KVDIMS(1) = NCOORDID(2,IGRID)%ID
+ ELSE IF ( YDIR == 'ZZ' .AND. KSHAPE(1)==NCOORDID(3,IGRID)%LEN) THEN
+ KVDIMS(1) = NCOORDID(3,IGRID)%ID
ELSE
PTDIM => GETDIMCDF(TPFILE, KSHAPE(1)); KVDIMS(1) = PTDIM%ID
END IF
--
GitLab