From 1d101fb9a5f5c811714362b3bbb8741362ba4f1b Mon Sep 17 00:00:00 2001
From: Juan Escobar <juan.escobar@aero.obs-mip.fr>
Date: Wed, 19 Mar 2014 16:02:44 +0000
Subject: [PATCH] Juan 01/2014 : correction for // initialisation geostrophic
ZUF,ZVF,ZWF
---
src/MNH/forcing.f90 | 13 +++++++++++++
1 file changed, 13 insertions(+)
diff --git a/src/MNH/forcing.f90 b/src/MNH/forcing.f90
index 4af45421e..3aac75017 100644
--- a/src/MNH/forcing.f90
+++ b/src/MNH/forcing.f90
@@ -144,6 +144,7 @@ END MODULE MODI_FORCING
!! 01/2004 V. Masson surface externalization, removes SST
!! forcing
!! 06/2012 V. Masson Adds tendency of geostrophic wind itself to wind tendency
+!! 01/2014 J. escobar correction for // initialisation geostrophic ZUF,ZVF,ZWF
!-------------------------------------------------------------------------------
!
!* 0. DECLARATIONS
@@ -166,6 +167,8 @@ USE MODI_TEMPORAL_LT
USE MODI_TEMPORAL_DIST
USE MODI_BUDGET
!
+USE MODI_GET_HALO
+!
IMPLICIT NONE
!
!* 0.1 Declarations of dummy arguments :
@@ -459,6 +462,8 @@ ELSE
!
ZZF(:,:,:) = MXM( ZZA(:,:,:) )
ZZF(1,:,:) = 2.0*ZZA(1,:,:) - ZZF(2,:,:)
+!
+ ZUF = 0.0
!
DO JL=1,IKU-1
ZDZHAT_INV = 1.0 / ( PZHAT(JL+1)-PZHAT(JL) )
@@ -476,9 +481,12 @@ ELSE
END DO
END DO
END DO
+ CALL GET_HALO(ZUF)
!
ZZF(:,:,:) = MYM( ZZA(:,:,:) )
ZZF(:,1,:) = 2.0*ZZA(:,1,:)-ZZF(:,2,:)
+!
+ ZVF = 0.0
!
DO JL=1,IKU-1
ZDZHAT_INV = 1.0 / ( PZHAT(JL+1)-PZHAT(JL) )
@@ -496,6 +504,10 @@ ELSE
END DO
END DO
END DO
+ CALL GET_HALO(ZVF)
+!
+ ZWF = 0.0
+!
DO JL=1,IKU-1
ZDZHAT_INV = 1.0 / ( PZHAT(JL+1)-PZHAT(JL) )
DO JK=1,IKU
@@ -512,6 +524,7 @@ ELSE
END DO
END DO
END DO
+ CALL GET_HALO(ZWF)
!
ZZF(:,:,:) = MZF(1,IKU,1, PZZ(:,:,:) )
ZZF(:,:,IKU) = 2.0*PZZ(:,:,IKU)-ZZF(:,:,IKU-1)
--
GitLab