From 8abd4fd1579c5018db789fb0a60212cbdaa0491c Mon Sep 17 00:00:00 2001
From: Juan Escobar <juan.escobar@aero.obs-mip.fr>
Date: Thu, 27 Mar 2014 16:17:01 +0000
Subject: [PATCH] Juan 27/03/2014 : write LAT/LON only in not CARTESIAN case

---
 src/MNH/write_lfifm1_for_diag.f90 | 25 ++++++++++++++-----------
 1 file changed, 14 insertions(+), 11 deletions(-)

diff --git a/src/MNH/write_lfifm1_for_diag.f90 b/src/MNH/write_lfifm1_for_diag.f90
index dfe2e2435..4d52fa7e3 100644
--- a/src/MNH/write_lfifm1_for_diag.f90
+++ b/src/MNH/write_lfifm1_for_diag.f90
@@ -135,6 +135,7 @@ END MODULE MODI_WRITE_LFIFM1_FOR_DIAG
 !!       S.Bielli 12/2012 : add latitude and longitude
 !!       F. Duffourg 02/2013 : add new fields
 !!      J.Escobar 21/03/2013: for HALOK get correctly local array dim/bound
+!!       J. escobar 27/03/2014 : write LAT/LON only in not CARTESIAN case
 !-------------------------------------------------------------------------------
 !
 !*       0.    DECLARATIONS
@@ -701,17 +702,19 @@ END IF
 !
 !*   Latitude and Longitude arrays
 !
-YRECFM='LAT'
-YCOMMENT='X_Y_latitude (degrees)'
-IGRID=1
-ILENCH=LEN(YCOMMENT)
-CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',XLAT,IGRID,ILENCH,YCOMMENT,IRESP)
-!
-YRECFM='LON'
-YCOMMENT='X_Y_longitude (degrees)'
-IGRID=1
-ILENCH=LEN(YCOMMENT)
-CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',XLON,IGRID,ILENCH,YCOMMENT,IRESP)
+IF (.NOT.LCARTESIAN) THEN
+   YRECFM='LAT'
+   YCOMMENT='X_Y_latitude (degrees)'
+   IGRID=1
+   ILENCH=LEN(YCOMMENT)
+   CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',XLAT,IGRID,ILENCH,YCOMMENT,IRESP)
+   !
+   YRECFM='LON'
+   YCOMMENT='X_Y_longitude (degrees)'
+   IGRID=1
+   ILENCH=LEN(YCOMMENT)
+   CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',XLON,IGRID,ILENCH,YCOMMENT,IRESP)
+END IF
 !
 !
 !-------------------------------------------------------------------------------
-- 
GitLab