w3ogcmmd.F90

Go to the documentation of this file.

#include "w3macros.h"
!/ ------------------------------------------------------------------- /
MODULE w3ogcmmd
  !/
  !/                  +-----------------------------------+
  !/                  | WAVEWATCH III           NOAA/NCEP |
  !/                  |           A. Thevenin             |
  !/                  |                        FORTRAN 90 |
  !/                  | Last update :         22-Mar-2021 |
  !/                  +-----------------------------------+
  !/
  !/    Jul-2013 : Origination.                       ( version 4.18 )
  !/               For upgrades see subroutines.
  !/    Apr-2016 : Add comments (J. Pianezze)         ( version 5.07 )
  !/ 22-Mar-2021 : Add extra coupling variables       ( version 7.13 )
  !/
  !/    Copyright 2009-2012 National Weather Service (NWS),
  !/       National Oceanic and Atmospheric Administration.  All rights
  !/       reserved.  WAVEWATCH III is a trademark of the NWS.
  !/       No unauthorized use without permission.
  !/
  !  1. Purpose :
  !
  !     Module used for coupling applications between oceanic model and WW3 with OASIS3-MCT
  !
  !  2. Variables and types :
  !
  !  3. Subroutines and functions :
  !
  !      Name                   Type  Scope    Description
  !     ----------------------------------------------------------------
  !      SND_FIELDS_TO_OCEAN    Subr. Public   Send fields to ocean model
  !      RCV_FIELDS_FROM_OCEAN  Subr. Public   Receive fields from ocean model
  !     ----------------------------------------------------------------
  !
  !  4. Subroutines and functions used :
  !
  !      Name                Type  Module    Description
  !     ----------------------------------------------------------------
  !      CPL_OASIS_SEND       Subr.   W3OACPMD   Send fields
  !      CPL_OASIS_RECV       Subr.   W3OACPMD   Receive fields
  !     ----------------------------------------------------------------
  !
  !  5. Remarks
  !  6. Switches :
  !  7. Source code :
  !
  !/ ------------------------------------------------------------------- /
  !
  IMPLICIT NONE
  !
  include "mpif.h"
  !
  PRIVATE
  !
  ! * Accessibility
  PUBLIC snd_fields_to_ocean
  PUBLIC rcv_fields_from_ocean
  !
CONTAINS
  !/ ------------------------------------------------------------------- /
  SUBROUTINE snd_fields_to_ocean()
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           A. Thevenin             |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         22-Mar-2021 |
    !/                  +-----------------------------------+
    !/
    !/    Jul-2013 : Origination.                    ( version 4.18 )
    !/    Apr-2016 : Add comments (J. Pianezze)      ( version 5.07 )
    !/ 22-Mar-2021 : Add extra coupling variables    ( version 7.13 )
    !/
    !  1. Purpose :
    !
    !     Send coupling fields to oceanic model
    !
    !  2. Method :
    !  3. Parameters :
    !  4. Subroutines used :
    !
    !     Name             Type    Module     Description
    !     -------------------------------------------------------------------
    !     CPL_OASIS_SND    Subr.   W3OACPMD   Send field to atmos/ocean model
    !     -------------------------------------------------------------------
    !
    !  5. Called by :
    !
    !      Name            Type    Module     Description
    !     ------------------------------------------------------------------
    !     W3WAVE           Subr.   W3WAVEMD   Wave model
    !     ------------------------------------------------------------------
    !
    !  6. Error messages :
    !  7. Remarks :
    !
    !     According to the present implementation, fields are sent at each coupling time step to OASIS
    !     Consequently, OASIS cannot estimate any time average
    !     For such an application, one must estimate the fields at each time step
    !     (or a time step smaller than the coupling time step).
    !     In such conditions, OASIS get the information every time step
    !     but only send the information to the other code when the time matches the coupling time.
    !
    !  8. Structure :
    !  9. Switches :
    ! 10. Source code :
    !
    !/ ------------------------------------------------------------------- /
    !
    USE w3oacpmd,  ONLY: id_oasis_time, il_nb_snd, snd_fld, cpl_oasis_snd
    USE w3gdatmd,  ONLY: nseal, mapsta, mapsf
    USE w3adatmd,  ONLY: hs, t0m1, t01, thm, bhd, tauox, tauoy, phioc,&
         uba, ubd, tauwix, tauwiy, tusx, tusy, ussx,  &
         ussy, wlm, phibbl,taubbl, charn, tauocx,     &
         tauocy, wnmean
    USE w3odatmd,  ONLY: naproc, iaproc, undef
    USE constants, ONLY: pi, dera
    !
    !/ ------------------------------------------------------------------- /
    !/ Parameter list
    !/
    !
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !/
    INTEGER                          :: i, isea, ix, iy
    INTEGER, DIMENSION(NSEAL)        :: mask
    REAL(kind=8), dimension(nseal,1) :: rla_oasis_snd
    INTEGER                          :: ib_do
    LOGICAL                          :: ll_action
    REAL(kind=8), dimension(nseal)   :: tmp
    !
    !----------------------------------------------------------------------
    ! * Executable part
    !
    DO i = 1, nseal
      isea = iaproc + (i-1)*naproc
      ix = mapsf(isea,1)
      iy = mapsf(isea,2)
      ! Get the mask : 1 - sea 0 - open boundary cells dried cells
      mask(i) = mod(mapsta(iy,ix),2)
    END DO
    !
    DO ib_do = 1, il_nb_snd
      !
      ! Mask - wet-drying
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3_ODRY') THEN
        tmp(1:nseal) = 0.0
        WHERE(mask(1:nseal) /= undef) tmp(1:nseal)=mask(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Mean wave period (tmn in s) (m0,-1)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3_T0M1') THEN
        tmp(1:nseal) = 0.0
        WHERE(t0m1(1:nseal) /= undef) tmp(1:nseal)=t0m1(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Mean wave period (tmn in s) (m0,1)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3__T01') THEN
        tmp(1:nseal) = 0.0
        WHERE(t01(1:nseal) /= undef) tmp(1:nseal)=t01(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Mean wave number (wnm in m-1)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3__WNM') THEN
        tmp(1:nseal) = 0.0
        WHERE(wnmean(1:nseal) /= undef) tmp(1:nseal)=wnmean(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Charnock coefficient  (-)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3_OCHA') THEN
        tmp(1:nseal) = 0.0
        WHERE(charn(1:nseal) /= undef) tmp(1:nseal)=charn(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Wave height (hs in m)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3__OHS') THEN
        tmp(1:nseal) = 0.0
        WHERE(hs(1:nseal) /= undef) tmp(1:nseal)=hs(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Cosinus of Mean wave direction (cos(theta) in radians)
      ! ---------------------------------------------------------------------
      ! dir : nautical convention (GRIDDED files) - 0 degree from north, 90 from east
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3_CDIR') THEN
        tmp(1:nseal) = 0.0
        WHERE(thm(1:nseal) /= undef) tmp(1:nseal)=cos(thm(1:nseal))
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Sinus of Mean wave direction (sin(theta) in radians)
      ! ---------------------------------------------------------------------
      ! dir : nautical convention (GRIDDED files) - 0 degree from north, 90 from east
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3_SDIR') THEN
        tmp(1:nseal) = 0.0
        WHERE(thm(1:nseal) /= undef) tmp(1:nseal)=sin(thm(1:nseal))
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Mean wave direction theta in radians
      ! ---------------------------------------------------------------------
      ! dir : nautical convention (GRIDDED files) - 0 degree from north, 90 from east
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3__DIR') THEN
        tmp(1:nseal) = 0.0
        WHERE(thm /= undef) tmp(1:nseal)=thm(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Wave-induced Bernoulli head pressure (bhd in N.m-1) (J term, Smith JPO 2006)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3__BHD') THEN
        tmp(1:nseal) = 0.0
        WHERE(bhd(1:nseal) /= undef) tmp(1:nseal)=bhd(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Wave-ocean momentum flux (tauox in m2.s-2)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3_TWOX') THEN
        tmp(1:nseal) = 0.0
        WHERE(tauox(1:nseal) /= undef) tmp(1:nseal)=tauox(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Wave-ocean momentum flux (tauoy in m2.s-2)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3_TWOY') THEN
        tmp(1:nseal) = 0.0
        WHERE(tauoy(1:nseal) /= undef) tmp(1:nseal)=tauoy(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Wave-ocean total momentum flux (tauocx in Pa)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3_TOCX') THEN
        tmp(1:nseal) = 0.0
        WHERE(tauocx(1:nseal) /= undef) tmp(1:nseal)=tauocx(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Wave-ocean total momentum flux (tauocy in Pa)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3_TOCY') THEN
        tmp(1:nseal) = 0.0
        WHERE(tauocy(1:nseal) /= undef) tmp(1:nseal)=tauocy(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Wave-to-ocean TKE flux (phioc in W.m-2)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3__FOC') THEN
        tmp(1:nseal) = 0.0
        WHERE(phioc(1:nseal) /= undef) tmp(1:nseal)=phioc(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Momentum flux due to bottom friction (taubblx in m2.s-2)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3_TBBX') THEN
        tmp(1:nseal) = 0.0
        WHERE(taubbl(1:nseal,1) /= undef) tmp(1:nseal)=taubbl(1:nseal,1)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Momentum flux due to bottom friction (taubbly in m2.s-2)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3_TBBY') THEN
        tmp(1:nseal) = 0.0
        WHERE(taubbl(1:nseal,2) /= undef) tmp(1:nseal)=taubbl(1:nseal,2)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Energy flux due to bottom friction (phibbl in W.m-2)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3__FBB') THEN
        tmp(1:nseal) = 0.0
        WHERE(phibbl(1:nseal) /= undef) tmp(1:nseal)=phibbl(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! rms amplitude of orbital velocity of the waves (ubr in m.s-1)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3__UBR') THEN
        tmp(1:nseal) = 0.0
        WHERE(uba(1:nseal) /= undef) tmp(1:nseal)=uba(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! x component of the near-bottom rms wave velocity (in m.s-1)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3_UBRX') THEN
        tmp(1:nseal) = 0.0
        WHERE(uba(1:nseal) /= undef) tmp(1:nseal)=uba(1:nseal)*cos(ubd(1:nseal))
        rla_oasis_snd(:,1) = tmp(1:nseal)
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! y component of the near-bottom rms wave velocity (in m.s-1)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3_UBRY') THEN
        tmp(1:nseal) = 0.0
        WHERE(uba(1:nseal) /= undef) tmp(1:nseal)=uba(1:nseal)*sin(ubd(1:nseal))
        rla_oasis_snd(:,1) = tmp(1:nseal)
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Net wave-supported stress, u component (tauwix in m2.s-2)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3_TAWX') THEN
        tmp(1:nseal) = 0.0
        WHERE(tauwix(1:nseal) /= undef) tmp(1:nseal)=tauwix(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Net wave-supported stress, v component (tauwix in m2.s-2)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3_TAWY') THEN
        tmp(1:nseal) = 0.0
        WHERE(tauwiy(1:nseal) /= undef) tmp(1:nseal)=tauwiy(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Volume transport associated to Stokes drift, u component (tusx in m2.s-1)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3_TUSX') THEN
        tmp(1:nseal) = 0.0
        WHERE(tusx(1:nseal) /= undef) tmp(1:nseal)=tusx(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Volume transport associated to Stokes drift, v component (tusy in m2.s-1)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3_TUSY') THEN
        tmp(1:nseal) = 0.0
        WHERE(tusy(1:nseal) /= undef) tmp(1:nseal)=tusy(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Surface Stokes drift, u component (ussx in m.s-1)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3_USSX') THEN
        tmp(1:nseal) = 0.0
        WHERE(ussx(1:nseal) /= undef) tmp(1:nseal)=ussx(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Surface Stokes drift, v component (ussy in m.s-1)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3_USSY') THEN
        tmp(1:nseal) = 0.0
        WHERE(ussy(1:nseal) /= undef) tmp(1:nseal)=ussy(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
      ! Mean wave length (wlm in m)
      ! ---------------------------------------------------------------------
      IF (snd_fld(ib_do)%CL_FIELD_NAME == 'WW3___LM') THEN
        tmp(1:nseal) = 0.0
        WHERE(wlm(1:nseal) /= undef) tmp(1:nseal)=wlm(1:nseal)
        rla_oasis_snd(:,1) = dble(tmp(1:nseal))
        CALL cpl_oasis_snd(ib_do, id_oasis_time, rla_oasis_snd, ll_action)
      ENDIF
      !
    ENDDO
    !/ ------------------------------------------------------------------- /
  END SUBROUTINE snd_fields_to_ocean
  !/ ------------------------------------------------------------------- /
  SUBROUTINE rcv_fields_from_ocean(ID_LCOMM, IDFLD, FXN, FYN, FAN)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           A. Thevenin             |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :          April-2016 |
    !/                  +-----------------------------------+
    !/
    !/    Jul-2013 : Origination.                               ( version 4.18 )
    !/    Apr-2014 : Add IDFLD, FXN, FYX and FAN (M. Accensi)   ( version 5.07 )
    !/    Apr-2016 : Add comments (J. Pianezze)                 ( version 5.07 )
    !/
    !  1. Purpose :
    !
    !     Receive coupling fields from oceanic model
    !
    !  2. Method :
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !     ID_LCOMM          Char.     I     MPI communicator
    !     IDFLD             Int.      I     Name of the exchange fields
    !     FXN               Int.     I/O    First exchange field
    !     FYN               Int.     I/O    Second exchange field
    !     FAN               Int.     I/O    Third exchange field
    !     ----------------------------------------------------------------
    !
    !  4. Subroutines used :
    !
    !     Name             Type    Module     Description
    !     -------------------------------------------------------------------
    !     CPL_OASIS_RCV    Subr.   W3OACPMD   Receive fields from atmos/ocean model
    !     W3S2XY           Subr.   W3SERVMD   Convert from storage (NSEA) to spatial grid (NX, NY)
    !     -------------------------------------------------------------------
    !
    !  5. Called by :
    !
    !     Name            Type    Module     Description
    !     ------------------------------------------------------------------
    !     W3FLDG          Subr.   W3FLDSMD   Manage input fields of depth,
    !                                        current, wind and ice concentration
    !     ------------------------------------------------------------------
    !
    !  6. Error messages :
    !  7. Remarks :
    !
    !     IDFLD   C*3  I/O ID string for field type, valid are: 'LEV', 'CUR' (J=1,2)
    !
    !  8. Structure :
    !  9. Switches :
    ! 10. Source code :
    !
    !/ ------------------------------------------------------------------- /
    !
    USE w3oacpmd, ONLY: id_oasis_time, il_nb_rcv, rcv_fld, cpl_oasis_rcv
    USE w3gdatmd, ONLY: nx, ny, nseal, nsea, mapsf
    USE w3odatmd, ONLY: naproc, iaproc
    USE w3servmd, ONLY: w3s2xy
    !
    !/ ------------------------------------------------------------------- /
    !/ Parameter list
    !/
    INTEGER, INTENT(IN)              :: id_lcomm
    CHARACTER(LEN=3), INTENT(IN)     :: idfld
    REAL, INTENT(INOUT)              :: fxn(:,:), fyn(:,:), fan(:,:)
    !
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !/
    LOGICAL                          :: ll_action
    INTEGER                          :: ib_do, ib_i, ib_j, il_err
    INTEGER, SAVE                    :: id_oasis_time_wetdryonlyonce = -1
    REAL(kind=8), dimension(nseal,1) :: rla_oasis_rcv
    REAL(kind=8), dimension(nseal)   :: tmp, maskt, masku, maskv
    REAL, DIMENSION(1:NSEA)          :: snd_buff,rcv_buff
    !
    !----------------------------------------------------------------------
    ! * Executable part
    !
    maskt(:)=1.
    masku(:)=1.
    maskv(:)=1.
    rla_oasis_rcv(:,:) = 0.0
    !
    ! ---------------------------------------------------------------------
    ! Perform mask variables
    ! ---------------------------------------------------------------------
    !
    ! For the same coupling time, W3FLDG is called for the level and current variables.
    ! As RCV_FIELDS_FROM_OCEAN is called from W3FLDG, the following test prevents to
    ! exchange the wet-dry variables more than once per coupling time.
    !cval well but it cannot work because MASKT,MASKU,MASKV variable are not global variable
    !cval Anyway we will give up the exchange of mask, it is not a good idea at all
 
    IF (id_oasis_time > id_oasis_time_wetdryonlyonce) THEN
      !
      DO ib_do = 1, il_nb_rcv
        !
        ! Land mask - u
        ! ---------------------------------------------------------------------
        IF (rcv_fld(ib_do)%CL_FIELD_NAME == 'WW3_OWDH') THEN
 
          CALL cpl_oasis_rcv(ib_do, id_oasis_time, rla_oasis_rcv, ll_action)
          IF (ll_action) THEN
            maskt(1:nseal)  = rla_oasis_rcv(1:nseal,1)
          ENDIF
        ENDIF
        !
        ! Land mask - h
        ! ---------------------------------------------------------------------
        IF (rcv_fld(ib_do)%CL_FIELD_NAME == 'WW3_OWDU') THEN
          CALL cpl_oasis_rcv(ib_do, id_oasis_time, rla_oasis_rcv, ll_action)
          IF (ll_action) THEN
            masku(1:nseal)  = rla_oasis_rcv(1:nseal,1)
          ENDIF
        ENDIF
        !
        ! Land mask - v
        ! ---------------------------------------------------------------------
        IF (rcv_fld(ib_do)%CL_FIELD_NAME == 'WW3_OWDV') THEN
          CALL cpl_oasis_rcv(ib_do, id_oasis_time, rla_oasis_rcv, ll_action)
          IF (ll_action) THEN
            maskv(1:nseal)  = rla_oasis_rcv(1:nseal,1)
          ENDIF
        ENDIF
        !
      ENDDO
      !
    ENDIF
    !
    ! ---------------------------------------------------------------------
    ! Treatment of the dynamical variables
    ! ---------------------------------------------------------------------
    DO ib_do = 1, il_nb_rcv
      !
      ! Sea surface Height (m)
      ! ---------------------------------------------------------------------
      IF (idfld == 'LEV') THEN
        !
        IF (rcv_fld(ib_do)%CL_FIELD_NAME == 'WW3__SSH') THEN
          CALL cpl_oasis_rcv(ib_do, id_oasis_time, rla_oasis_rcv, ll_action)
          IF (ll_action) THEN
            tmp(1:nseal) = rla_oasis_rcv(1:nseal,1) * maskt(1:nseal)
            snd_buff(1:nsea) = 0.0
            DO ib_i = 1, nseal
              ib_j = iaproc + (ib_i-1)*naproc
              snd_buff(ib_j) = tmp(ib_i)
            ENDDO
            !
            CALL mpi_allreduce(snd_buff(1:nsea), &
                 rcv_buff(1:nsea), &
                 nsea,     &
                 mpi_real, &
                 mpi_sum,  &
                 id_lcomm, &
                 il_err)
            !
            ! Convert from storage (NSEA) to spatial grid (NX, NY)
            CALL w3s2xy(nsea,nsea,nx,ny,rcv_buff(1:nsea),mapsf,fan)
            !
          ENDIF
        ENDIF
      ENDIF
      !
      ! Ocean sea surface current (m.s-1)
      ! ---------------------------------------------------------------------
      IF (idfld == 'CUR') THEN
        !
        ! u-component
        ! ---------------------------------------------------------------------
        IF (rcv_fld(ib_do)%CL_FIELD_NAME == 'WW3_OSSU') THEN
          CALL cpl_oasis_rcv(ib_do, id_oasis_time, rla_oasis_rcv, ll_action)
          IF (ll_action) THEN
            tmp(1:nseal) = rla_oasis_rcv(1:nseal,1) * masku(1:nseal)
            snd_buff(1:nsea) = 0.0
            DO ib_i = 1, nseal
              ib_j = iaproc + (ib_i-1)*naproc
              snd_buff(ib_j) = tmp(ib_i)
            ENDDO
            !
            CALL mpi_allreduce(snd_buff(1:nsea),       &
                 rcv_buff(1:nsea),       &
                 nsea,     &
                 mpi_real, &
                 mpi_sum,  &
                 id_lcomm, &
                 il_err)
            !
            ! Convert from storage (NSEA) to spatial grid (NX, NY)
            CALL w3s2xy(nsea,nsea,nx,ny,rcv_buff(1:nsea),mapsf,fxn)
            !
          ENDIF
        ENDIF
        !
        ! v-component
        ! ---------------------------------------------------------------------
        IF (rcv_fld(ib_do)%CL_FIELD_NAME == 'WW3_OSSV') THEN
          CALL cpl_oasis_rcv(ib_do, id_oasis_time, rla_oasis_rcv, ll_action)
          IF (ll_action) THEN
            tmp(1:nseal) = rla_oasis_rcv(1:nseal,1) * maskv(1:nseal)
            snd_buff(1:nsea) = 0.0
            DO ib_i = 1, nseal
              ib_j = iaproc + (ib_i-1)*naproc
              snd_buff(ib_j) = tmp(ib_i)
            ENDDO
            !
            CALL mpi_allreduce(snd_buff(1:nsea),       &
                 rcv_buff(1:nsea),       &
                 nsea,     &
                 mpi_real, &
                 mpi_sum,  &
                 id_lcomm, &
                 il_err)
            !
            ! Convert from storage (NSEA) to spatial grid (NX, NY)
            CALL w3s2xy(nsea,nsea,nx,ny,rcv_buff(1:nsea),mapsf,fyn)
            !
          ENDIF
        ENDIF
      ENDIF
    ENDDO
    !
    id_oasis_time_wetdryonlyonce = id_oasis_time
    !
    !/ ------------------------------------------------------------------- /
  END SUBROUTINE rcv_fields_from_ocean
  !/ ------------------------------------------------------------------- /
  !/
END MODULE w3ogcmmd
!/
!/ ------------------------------------------------------------------- /