Changeset 14961 for NEMO/trunk/src

Timestamp:

2021-06-09T11:38:51+02:00 (3 years ago)

Author:

amoulin

Message:

read wave number field if ln_zdfswm activated -ticket #2685

File:

• NEMO/trunk/src/OCE/SBC/sbcwave.F90 (modified) (3 diffs)

NEMO/trunk/src/OCE/SBC/sbcwave.F90

@@ -24 +24 @@
    USE bdy_oce        ! open boundary condition variables
    USE domvvl         ! domain: variable volume layers
+   USE zdf_oce,  ONLY : ln_zdfswm ! Qiao wave enhanced mixing 
    !
    USE iom            ! I/O manager library
@@ -310 +311 @@
             IF( jp_vsd > 0 )   vt0sd(:,:) = sf_sd(jp_vsd)%fnow(:,:,1) * tmask(:,:,1)  ! 2D meridional Stokes Drift at T point
          ENDIF
+
+         ! Read also wave number if needed, so that it is available in
+         ! coupling routines
+         IF( ln_zdfswm .AND. .NOT. cpl_wnum ) THEN     !==wavenumber==!
+            CALL fld_read( kt, nn_fsbc, sf_wn )             ! read wave parameters from external forcing
+            wnum(:,:) = sf_wn(1)%fnow(:,:,1) * tmask(:,:,1)
+         ENDIF
+ 
          !
          IF( jpfld == 4 .OR. ln_wave_test )   &
@@ -506 +515 @@
             ENDIF
             !
-            ! 3. Wave number (only needed for Qiao parametrisation, ln_zdfqiao=T)
-            IF( .NOT. cpl_wnum ) THEN
+            ! 3. Wave number (only needed for Qiao parametrisation, ln_zdfswm=T)
+            IF( ln_zdfswm .AND. .NOT. cpl_wnum ) THEN
                ALLOCATE( sf_wn(1), STAT=ierror )           !* allocate and fill sf_wave with sn_wnum
                IF( ierror > 0 )   CALL ctl_stop( 'STOP', 'sbc_wave_init: unable to allocate sf_wn structure' )