w3ogcmmd Module Reference

Functions/Subroutines
subroutine, public	snd_fields_to_ocean ()
subroutine, public	rcv_fields_from_ocean (ID_LCOMM, IDFLD, FXN, FYN, FAN)

Function/Subroutine Documentation

rcv_fields_from_ocean()

subroutine, public w3ogcmmd::rcv_fields_from_ocean	(	integer, intent(in)	ID_LCOMM,
		character(len=3), intent(in)	IDFLD,
		real, dimension(:,:), intent(inout)	FXN,
		real, dimension(:,:), intent(inout)	FYN,
		real, dimension(:,:), intent(inout)	FAN
	)

Definition at line 407 of file w3ogcmmd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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
    !
    !/ ------------------------------------------------------------------- /

References w3oacpmd::cpl_oasis_rcv(), w3odatmd::iaproc, w3oacpmd::id_oasis_time, w3oacpmd::il_nb_rcv, w3gdatmd::mapsf, w3odatmd::naproc, w3gdatmd::nsea, w3gdatmd::nseal, w3gdatmd::nx, w3gdatmd::ny, w3oacpmd::rcv_fld, and w3servmd::w3s2xy().

Referenced by w3fldsmd::w3fldg().

snd_fields_to_ocean()

subroutine, public w3ogcmmd::snd_fields_to_ocean

Definition at line 63 of file w3ogcmmd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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
    !/ ------------------------------------------------------------------- /

References w3adatmd::bhd, w3adatmd::charn, w3oacpmd::cpl_oasis_snd(), constants::dera, w3adatmd::hs, w3odatmd::iaproc, w3oacpmd::id_oasis_time, w3oacpmd::il_nb_snd, w3gdatmd::mapsf, w3gdatmd::mapsta, w3odatmd::naproc, w3gdatmd::nseal, w3adatmd::phibbl, w3adatmd::phioc, constants::pi, w3oacpmd::snd_fld, w3adatmd::t01, w3adatmd::t0m1, w3adatmd::taubbl, w3adatmd::tauocx, w3adatmd::tauocy, w3adatmd::tauox, w3adatmd::tauoy, w3adatmd::tauwix, w3adatmd::tauwiy, w3adatmd::thm, w3adatmd::tusx, w3adatmd::tusy, w3adatmd::uba, w3adatmd::ubd, w3adatmd::ussx, w3adatmd::ussy, w3adatmd::wlm, and w3adatmd::wnmean.

Referenced by w3shel(), and w3wavemd::w3wave().