Skip to content
Snippets Groups Projects
Commit 1d101fb9 authored by ESCOBAR MUNOZ Juan's avatar ESCOBAR MUNOZ Juan Committed by WAUTELET Philippe
Browse files

Juan 01/2014 : correction for // initialisation geostrophic ZUF,ZVF,ZWF

parent 9197911a
No related branches found
No related tags found
No related merge requests found
...@@ -144,6 +144,7 @@ END MODULE MODI_FORCING ...@@ -144,6 +144,7 @@ END MODULE MODI_FORCING
!! 01/2004 V. Masson surface externalization, removes SST !! 01/2004 V. Masson surface externalization, removes SST
!! forcing !! forcing
!! 06/2012 V. Masson Adds tendency of geostrophic wind itself to wind tendency !! 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 !* 0. DECLARATIONS
...@@ -166,6 +167,8 @@ USE MODI_TEMPORAL_LT ...@@ -166,6 +167,8 @@ USE MODI_TEMPORAL_LT
USE MODI_TEMPORAL_DIST USE MODI_TEMPORAL_DIST
USE MODI_BUDGET USE MODI_BUDGET
! !
USE MODI_GET_HALO
!
IMPLICIT NONE IMPLICIT NONE
! !
!* 0.1 Declarations of dummy arguments : !* 0.1 Declarations of dummy arguments :
...@@ -459,6 +462,8 @@ ELSE ...@@ -459,6 +462,8 @@ ELSE
! !
ZZF(:,:,:) = MXM( ZZA(:,:,:) ) ZZF(:,:,:) = MXM( ZZA(:,:,:) )
ZZF(1,:,:) = 2.0*ZZA(1,:,:) - ZZF(2,:,:) ZZF(1,:,:) = 2.0*ZZA(1,:,:) - ZZF(2,:,:)
!
ZUF = 0.0
! !
DO JL=1,IKU-1 DO JL=1,IKU-1
ZDZHAT_INV = 1.0 / ( PZHAT(JL+1)-PZHAT(JL) ) ZDZHAT_INV = 1.0 / ( PZHAT(JL+1)-PZHAT(JL) )
...@@ -476,9 +481,12 @@ ELSE ...@@ -476,9 +481,12 @@ ELSE
END DO END DO
END DO END DO
END DO END DO
CALL GET_HALO(ZUF)
! !
ZZF(:,:,:) = MYM( ZZA(:,:,:) ) ZZF(:,:,:) = MYM( ZZA(:,:,:) )
ZZF(:,1,:) = 2.0*ZZA(:,1,:)-ZZF(:,2,:) ZZF(:,1,:) = 2.0*ZZA(:,1,:)-ZZF(:,2,:)
!
ZVF = 0.0
! !
DO JL=1,IKU-1 DO JL=1,IKU-1
ZDZHAT_INV = 1.0 / ( PZHAT(JL+1)-PZHAT(JL) ) ZDZHAT_INV = 1.0 / ( PZHAT(JL+1)-PZHAT(JL) )
...@@ -496,6 +504,10 @@ ELSE ...@@ -496,6 +504,10 @@ ELSE
END DO END DO
END DO END DO
END DO END DO
CALL GET_HALO(ZVF)
!
ZWF = 0.0
!
DO JL=1,IKU-1 DO JL=1,IKU-1
ZDZHAT_INV = 1.0 / ( PZHAT(JL+1)-PZHAT(JL) ) ZDZHAT_INV = 1.0 / ( PZHAT(JL+1)-PZHAT(JL) )
DO JK=1,IKU DO JK=1,IKU
...@@ -512,6 +524,7 @@ ELSE ...@@ -512,6 +524,7 @@ ELSE
END DO END DO
END DO END DO
END DO END DO
CALL GET_HALO(ZWF)
! !
ZZF(:,:,:) = MZF(1,IKU,1, PZZ(:,:,:) ) ZZF(:,:,:) = MZF(1,IKU,1, PZZ(:,:,:) )
ZZF(:,:,IKU) = 2.0*PZZ(:,:,IKU)-ZZF(:,:,IKU-1) ZZF(:,:,IKU) = 2.0*PZZ(:,:,IKU)-ZZF(:,:,IKU-1)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment