From e8a037cfe770b9c7ef309a5cbffbbee8acf4482a Mon Sep 17 00:00:00 2001
From: Gaelle Tanguy <gaelle.tanguy@meteo.fr>
Date: Fri, 24 Jul 2015 08:31:15 +0000
Subject: [PATCH] Gaelle 24/07/2015: add UM10/VM10 for LCARTESIAN=T cases
---
src/MNH/write_lfifm1_for_diag_supp.f90 | 17 +++++++++++++++--
1 file changed, 15 insertions(+), 2 deletions(-)
diff --git a/src/MNH/write_lfifm1_for_diag_supp.f90 b/src/MNH/write_lfifm1_for_diag_supp.f90
index 49b87ef25..a4a61b436 100644
--- a/src/MNH/write_lfifm1_for_diag_supp.f90
+++ b/src/MNH/write_lfifm1_for_diag_supp.f90
@@ -79,6 +79,7 @@ END MODULE MODI_WRITE_LFIFM1_FOR_DIAG_SUPP
!! change of YCOMMENT
!! October 2011 (C.Lac) FF10MAX : interpolation of 10m wind
!! between 2 Meso-NH levels if 10m is above the first atmospheric level
+!! 2015 : D.Ricard add UM10/VM10 for LCARTESIAN=T cases
!-------------------------------------------------------------------------------
!
!* 0. DECLARATIONS
@@ -727,9 +728,21 @@ IF (CSURF=='EXTE') THEN
IGRID=0
! in this case (argument IGRID=0), input winds are ZONal and MERidien
! and, output ones are in MesoNH grid
- IF (.NOT. LCARTESIAN) &
- CALL UV_TO_ZONAL_AND_MERID(XCURRENT_ZON10M,XCURRENT_MER10M,IGRID, &
+ IF (.NOT. LCARTESIAN) THEN
+ CALL UV_TO_ZONAL_AND_MERID(XCURRENT_ZON10M,XCURRENT_MER10M,IGRID, &
HFMFILE=HFMFILE,HRECU='UM10',HRECV='VM10',HCOMMENT=YCOMMENT)
+ ELSE
+ YRECFM ='UM10'
+ YCOMMENT ='X_Y_UM10 (m/s)'
+ IGRID =1
+ ILENCH =LEN(YCOMMENT)
+ CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',XCURRENT_ZON10M,IGRID,ILENCH,YCOMMENT,IRESP)
+ YRECFM ='VM10'
+ YCOMMENT ='X_Y_VM10 (m/s)'
+ IGRID =1
+ ILENCH =LEN(YCOMMENT)
+ CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',XCURRENT_MER10M,IGRID,ILENCH,YCOMMENT,IRESP)
+ ENDIF
!
IF (SIZE(XTKET)>0) THEN
ZWORK21(:,:)= 0.
--
GitLab