From 487e84f972a0a0b01bdc9a52a80f0d8dd605bb81 Mon Sep 17 00:00:00 2001 From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr> Date: Mon, 5 Feb 2018 17:14:24 +0100 Subject: [PATCH] Philippe 05/02/2018: IO: write integer cluster fields as 3D integer arrays (and not real) --- src/MNH/write_lfifm1_for_diag.f90 | 13 +++++-------- 1 file changed, 5 insertions(+), 8 deletions(-) diff --git a/src/MNH/write_lfifm1_for_diag.f90 b/src/MNH/write_lfifm1_for_diag.f90 index 12431d0d3..db5bf7f61 100644 --- a/src/MNH/write_lfifm1_for_diag.f90 +++ b/src/MNH/write_lfifm1_for_diag.f90 @@ -303,7 +303,7 @@ REAL,DIMENSION(:,:,:), ALLOCATABLE :: ZRHOT, ZTMP ! work array LOGICAL :: GBOTUP ! clustering propagation LOGICAL,DIMENSION(SIZE(XTHT,1),SIZE(XTHT,2),SIZE(XTHT,3)) :: GCLOUD ! mask INTEGER,DIMENSION(SIZE(XTHT,1),SIZE(XTHT,2),SIZE(XTHT,3)) :: ICLUSTERID, ICLUSTERLV -REAL, DIMENSION(SIZE(XTHT,1),SIZE(XTHT,2),SIZE(XTHT,3)) :: ZCLUSTERID, ZCLUSTERLV, ZCLDSIZE +REAL, DIMENSION(SIZE(XTHT,1),SIZE(XTHT,2),SIZE(XTHT,3)) :: ZCLDSIZE !ECRITURE DANS UN FICHIER ASCII DE RESULTATS !INITIALISATION DU NOM DE FICHIER CREE EN PARALLELE AVEC CELUI LFI @@ -2570,9 +2570,6 @@ IF (LCLSTR) THEN CALL CLUSTERING(GBOTUP,GCLOUD,XWT,ICLUSTERID,ICLUSTERLV,ZCLDSIZE) PRINT *,'GOT OUT OF CLUSTERING' ! - ZCLUSTERID=ICLUSTERID - ZCLUSTERLV=ICLUSTERLV - ! TZFIELD%CMNHNAME = 'CLUSTERID' TZFIELD%CSTDNAME = '' TZFIELD%CLONGNAME = 'CLUSTERID' @@ -2580,9 +2577,9 @@ IF (LCLSTR) THEN TZFIELD%CDIR = 'XY' TZFIELD%CCOMMENT = 'X_Y_Z_CLUSTER (ID NUMBER)' TZFIELD%NGRID = 1 - TZFIELD%NTYPE = TYPEREAL + TZFIELD%NTYPE = TYPEINT TZFIELD%NDIMS = 3 - CALL IO_WRITE_FIELD(TPFILE,TZFIELD,ZCLUSTERID) + CALL IO_WRITE_FIELD(TPFILE,TZFIELD,ICLUSTERID) ! TZFIELD%CMNHNAME = 'CLUSTERLV' TZFIELD%CSTDNAME = '' @@ -2591,9 +2588,9 @@ IF (LCLSTR) THEN TZFIELD%CDIR = 'XY' TZFIELD%CCOMMENT = 'X_Y_Z_CLUSTER (BASE OR TOP LEVEL)' TZFIELD%NGRID = 1 - TZFIELD%NTYPE = TYPEREAL + TZFIELD%NTYPE = TYPEINT TZFIELD%NDIMS = 3 - CALL IO_WRITE_FIELD(TPFILE,TZFIELD,ZCLUSTERLV) + CALL IO_WRITE_FIELD(TPFILE,TZFIELD,ICLUSTERLV) ! TZFIELD%CMNHNAME = 'CLDSIZE' TZFIELD%CSTDNAME = '' -- GitLab