C.Lac 15/11/2016 Clarify top boundary condition for W
authorGaelle TANGUY <gaelle.tanguy@meteo.fr>
Tue, 15 Nov 2016 09:38:09 +0000 (10:38 +0100)
committerJuan Escobar <escj@aero.obs-mip.fr>
Tue, 22 Nov 2016 14:45:01 +0000 (15:45 +0100)
(cherry picked from commit e387e7c86015172c5f9cf55b772ee95804ed3005)

src/MNH/adv_boundaries.f90

index 63f9001..2dd15d8 100644 (file)
@@ -32,9 +32,10 @@ END MODULE MODI_ADV_BOUNDARIES
 !!
 !!    AUTHOR
 !!    ------
-!!
+!!   V.Masson
 !! Correction :        
-!!   J.Escobar : 15/09/2015 : WENO5 & JPHEXT <> 1 
+!!   J.Escobar : 15/09/2015 : WENO5 & JPHEXT <> 1
+!!   C.Lac     : 10/16 : top BC for W 
 !-------------------------------------------------------------------------------
 !
 !*       0.    DECLARATIONS
@@ -119,6 +120,11 @@ IF( PRESENT(PFIELDI) )  THEN
   END IF
 END IF
 !
+!*       4. TOP BC for W
+!
+ IF (PRESENT(HFIELD)) THEN
+   IF (HFIELD=='W') PFIELD(:,:,IKE+1) = 0.
+ END IF
 !-------------------------------------------------------------------------------
 !
 END SUBROUTINE ADV_BOUNDARIES