From f88eab07cec4540a9e3de3fe7d0cb150d4df0ea1 Mon Sep 17 00:00:00 2001
From: Guillaume S <gsamson@mercator-ocean.fr>
Date: Tue, 15 Mar 2022 10:58:44 +0100
Subject: [PATCH] use zero wind value if wind speed components are set to 'NOT
USED'
---
src/OCE/SBC/sbcblk.F90 | 2 ++
1 file changed, 2 insertions(+)
diff --git a/src/OCE/SBC/sbcblk.F90 b/src/OCE/SBC/sbcblk.F90
index a68659642..3322532d6 100644
--- a/src/OCE/SBC/sbcblk.F90
+++ b/src/OCE/SBC/sbcblk.F90
@@ -369,6 +369,8 @@ CONTAINS
sf(jfpr)%fnow(:,:,1:ipka) = 101325._wp ! use standard pressure in Pa
ELSEIF( jfpr == jp_prec .OR. jfpr == jp_snow .OR. jfpr == jp_uoatm .OR. jfpr == jp_voatm ) THEN
sf(jfpr)%fnow(:,:,1:ipka) = 0._wp ! no precip or no snow or no surface currents
+ ELSEIF( jfpr == jp_wndi .OR. jfpr == jp_wndj ) THEN
+ sf(jfpr)%fnow(:,:,1:ipka) = 0._wp
ELSEIF( jfpr == jp_hpgi .OR. jfpr == jp_hpgj ) THEN
IF( .NOT. ln_abl ) THEN
DEALLOCATE( sf(jfpr)%fnow ) ! deallocate as not used in this case
--
GitLab