From d2ee3de57a7b676778950fd5d0cd5d6ec57a6fba Mon Sep 17 00:00:00 2001
From: Gaelle Tanguy <gaelle.tanguy@meteo.fr>
Date: Wed, 23 Apr 2014 07:53:51 +0000
Subject: [PATCH] S.Bielli 23/04/2014 : supress writing of LAt and LON in
NETCDF case
---
src/MNH/prep_pgd.f90 | 61 ++++++++++++++++++++++----------------------
1 file changed, 31 insertions(+), 30 deletions(-)
diff --git a/src/MNH/prep_pgd.f90 b/src/MNH/prep_pgd.f90
index 475696c04..4aea3bc85 100644
--- a/src/MNH/prep_pgd.f90
+++ b/src/MNH/prep_pgd.f90
@@ -58,6 +58,7 @@
!! J.Escobar 4/04/2008 Improve checking --> add STATUS=OLD in open_ll(PRE_PGD1.nam,...
!!
!! Modification 30/03/2012 Add NAM_NCOUT for netcdf output (S.Bielli)
+!! S.Bielli 23/04/2014 supress writing of LAt and LON in NETCDF case
!----------------------------------------------------------------------------
!
!* 0. DECLARATION
@@ -229,22 +230,22 @@ CALL FMWRIT(COUTFMFILE,'PACK ',CLUOUT0,'--',LPACK,0,1,' ',IRESP)
#ifdef MNH_NCWRIT
NC_WRITE = LNETCDF
CALL WRITE_PGD_SURF_ATM_n('MESONH')
-IF (LNETCDF.AND..NOT.LCARTESIAN) THEN
- LLFIFM = .FALSE.
-!!!!! WRITE LAT and LON
- ALLOCATE(ZWORK(NIMAX_ll+NHALO*2,NJMAX_ll+NHALO*2))
- ALLOCATE(ZWORK_LAT(NIMAX_ll+2,NJMAX_ll+2))
- ALLOCATE(ZWORK_LON(NIMAX_ll+2,NJMAX_ll+2))
- ZWORK=RESHAPE(XLAT, (/ (NIMAX_ll+NHALO*2),(NJMAX_ll+NHALO*2) /) )
- ZWORK_LAT=ZWORK(NHALO:(NIMAX_ll+NHALO+1),NHALO:(NJMAX_ll+NHALO+1))
-!
-CALL FMWRIT(COUTFMFILE,'LAT',CLUOUT0,'XY',ZWORK_LAT,1,21,'X_Y_latitude (degree)',IRESP)
- ZWORK=RESHAPE(XLON, (/ NIMAX_ll+NHALO*2,NJMAX_ll+NHALO*2 /) )
- ZWORK_LON=ZWORK(NHALO:(NIMAX_ll+NHALO+1),NHALO:(NJMAX_ll+NHALO+1))
-CALL FMWRIT(COUTFMFILE,'LON',CLUOUT0,'XY',ZWORK_LON,1,22,'X_Y_longitude (degree)',IRESP)
- DEALLOCATE(ZWORK)
- LLFIFM = .TRUE.
-END IF
+!IF (LNETCDF.AND..NOT.LCARTESIAN) THEN
+! LLFIFM = .FALSE.
+!!!! WRITE LAT and LON
+! ALLOCATE(ZWORK(NIMAX_ll+NHALO*2,NJMAX_ll+NHALO*2))
+! ALLOCATE(ZWORK_LAT(NIMAX_ll+2,NJMAX_ll+2))
+! ALLOCATE(ZWORK_LON(NIMAX_ll+2,NJMAX_ll+2))
+! ZWORK=RESHAPE(XLAT, (/ (NIMAX_ll+NHALO*2),(NJMAX_ll+NHALO*2) /) )
+! ZWORK_LAT=ZWORK(NHALO:(NIMAX_ll+NHALO+1),NHALO:(NJMAX_ll+NHALO+1))
+!!
+!CALL FMWRIT(COUTFMFILE,'LAT',CLUOUT0,'XY',ZWORK_LAT,1,21,'X_Y_latitude (degree)',IRESP)
+! ZWORK=RESHAPE(XLON, (/ NIMAX_ll+NHALO*2,NJMAX_ll+NHALO*2 /) )
+! ZWORK_LON=ZWORK(NHALO:(NIMAX_ll+NHALO+1),NHALO:(NJMAX_ll+NHALO+1))
+!CALL FMWRIT(COUTFMFILE,'LON',CLUOUT0,'XY',ZWORK_LON,1,22,'X_Y_longitude (degree)',IRESP)
+! DEALLOCATE(ZWORK)
+! LLFIFM = .TRUE.
+!END IF
!* 4. Computes and writes smooth orography for SLEVE coordinate
! ---------------------------------------------------------
CALL ZSMT_PGD(COUTFMFILE,NZSFILTER,NSLEVE,XSMOOTH_ZS)
@@ -252,20 +253,20 @@ CALL ZSMT_PGD(COUTFMFILE,NZSFILTER,NSLEVE,XSMOOTH_ZS)
IF ( LNETCDF ) THEN
DEF_NC=.FALSE.
CALL WRITE_PGD_SURF_ATM_n('MESONH')
- IF (LNETCDF.AND..NOT.LCARTESIAN) THEN
- LLFIFM = .FALSE.
-!!!! WRITE LAT and LON
- CALL FMWRIT(COUTFMFILE,'LAT',CLUOUT0,'XY',ZWORK_LAT,1,21,'X_Y_latitude (degree)',IRESP)
- CALL FMWRIT(COUTFMFILE,'LON',CLUOUT0,'XY',ZWORK_LON,1,22,'X_Y_longitude (degree)',IRESP)
- END IF
- ALLOCATE(ZZS(NIMAX_ll+2,NJMAX_ll+2))
-!!!! writes smooth orography for SLEVE coordinate in netcdf
- YRECFM = 'ZS '
- CALL FMREAD(COUTFMFILE,YRECFM,CLUOUT0,'XY',ZZS,IGRID,ILENCH,YCOMMENT,IRESP)
- CALL FMWRIT(COUTFMFILE,'ZS',CLUOUT0,'XY',ZZS,IGRID,ILENCH,YCOMMENT,IRESP)
- YRECFM = 'ZSMT '
- CALL FMREAD(COUTFMFILE,YRECFM,CLUOUT0,'XY',ZZS,IGRID,ILENCH,YCOMMENT,IRESP)
- CALL FMWRIT(COUTFMFILE,'ZSMT',CLUOUT0,'XY',ZZS,IGRID,ILENCH,YCOMMENT,IRESP)
+! IF (LNETCDF.AND..NOT.LCARTESIAN) THEN
+! LLFIFM = .FALSE.
+!!!!! WRITE LAT and LON
+! CALL FMWRIT(COUTFMFILE,'LAT',CLUOUT0,'XY',ZWORK_LAT,1,21,'X_Y_latitude (degree)',IRESP)
+! CALL FMWRIT(COUTFMFILE,'LON',CLUOUT0,'XY',ZWORK_LON,1,22,'X_Y_longitude (degree)',IRESP)
+! END IF
+! ALLOCATE(ZZS(NIMAX_ll+2,NJMAX_ll+2))
+!!!!! writes smooth orography for SLEVE coordinate in netcdf
+! YRECFM = 'ZS '
+! CALL FMREAD(COUTFMFILE,YRECFM,CLUOUT0,'XY',ZZS,IGRID,ILENCH,YCOMMENT,IRESP)
+! CALL FMWRIT(COUTFMFILE,'ZS',CLUOUT0,'XY',ZZS,IGRID,ILENCH,YCOMMENT,IRESP)
+! YRECFM = 'ZSMT '
+! CALL FMREAD(COUTFMFILE,YRECFM,CLUOUT0,'XY',ZZS,IGRID,ILENCH,YCOMMENT,IRESP)
+! CALL FMWRIT(COUTFMFILE,'ZSMT',CLUOUT0,'XY',ZZS,IGRID,ILENCH,YCOMMENT,IRESP)
DEF_NC=.TRUE.
NC_WRITE = .FALSE.
END IF
--
GitLab