source: NEMO/branches/2020/dev_r12512_HPC-04_mcastril_Mixed_Precision_implementation/tests/CANAL/MY_SRC/usrdef_sbc.F90 @ 13135

Last change on this file since 13135 was 13135, checked in by orioltp, 3 months ago

dev_r12512_HPC-04_mcastril_Mixed_Precision_implementation: merge with trunk@13134, see #2364

  • Property svn:keywords set to Id
File size: 4.6 KB
Line 
1MODULE usrdef_sbc
2   !!======================================================================
3   !!                       ***  MODULE  usrdef_sbc  ***
4   !!
5   !!                      ===  CANAL configuration  ===
6   !!
7   !! User defined :   surface forcing of a user configuration
8   !!======================================================================
9   !! History :  4.0   ! 2017-11  (J.Chanut)  user defined interface
10   !!----------------------------------------------------------------------
11
12   !!----------------------------------------------------------------------
13   !!   usr_def_sbc    : user defined surface bounday conditions in OVERFLOW case
14   !!----------------------------------------------------------------------
15   USE oce             ! ocean dynamics and tracers
16   USE dom_oce         ! ocean space and time domain
17   USE sbc_oce         ! Surface boundary condition: ocean fields
18   USE phycst          ! physical constants
19   USE usrdef_nam, ONLY : rn_u10, rn_uofac, rn_windszy 
20   !
21   USE in_out_manager  ! I/O manager
22   USE lib_mpp         ! distribued memory computing library
23   USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
24   USE lib_fortran     ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined)
25
26   IMPLICIT NONE
27   PRIVATE
28
29   PUBLIC   usrdef_sbc_oce      ! routine called in sbcmod module
30   PUBLIC   usrdef_sbc_ice_tau  ! routine called by icestp.F90 for ice dynamics
31   PUBLIC   usrdef_sbc_ice_flx  ! routine called by icestp.F90 for ice thermo
32
33   !!----------------------------------------------------------------------
34   !! NEMO/OCE 4.0 , NEMO Consortium (2018)
35   !! $Id$
36   !! Software governed by the CeCILL license (see ./LICENSE)
37   !!----------------------------------------------------------------------
38CONTAINS
39
40   SUBROUTINE usrdef_sbc_oce( kt, Kbb )
41      !!---------------------------------------------------------------------
42      !!                    ***  ROUTINE usr_def_sbc  ***
43      !!             
44      !! ** Purpose :   provide at each time-step the surface boundary
45      !!              condition, i.e. the momentum, heat and freshwater fluxes.
46      !!
47      !! ** Method  :   all 0 fields, for CANAL case
48      !!                CAUTION : never mask the surface stress field !
49      !!
50      !! ** Action  : - set to ZERO all the ocean surface boundary condition, i.e.   
51      !!                   utau, vtau, taum, wndm, qns, qsr, emp, sfx
52      !!
53      !!----------------------------------------------------------------------
54      INTEGER, INTENT(in) ::   kt        ! ocean time step
55      INTEGER, INTENT(in) ::   Kbb       ! ocean time index
56      INTEGER  ::   ji, jj               ! dummy loop indices
57      REAL(wp) :: zrhoair = 1.22     ! approximate air density [Kg/m3]
58      REAL(wp) :: zcd = 1.13e-3      ! approximate drag coefficient
59      REAL(wp) :: zrhocd             ! Rho * Cd
60      REAL(wp), DIMENSION(jpi,jpj) :: zwndrel   ! relative wind
61      !!---------------------------------------------------------------------
62      !
63      zrhocd = zrhoair * zcd
64     
65      IF( kt == nit000 ) THEN
66         !
67         IF(lwp) WRITE(numout,*)' usr_sbc : EW_CANAL case: surface forcing'
68         IF(lwp) WRITE(numout,*)' ~~~~~~~~~~~   vtau = taum = wndm = qns = qsr = emp = sfx = 0'
69         !
70         utau(:,:) = 0._wp
71         IF( rn_u10 /= 0. .AND. rn_windszy > 0. ) THEN
72            WHERE( ABS(gphit) <= rn_windszy/2. ) utau(:,:) = zrhocd * rn_u10 * rn_u10
73         ENDIF
74         vtau(:,:) = 0._wp
75         taum(:,:) = 0._wp
76         wndm(:,:) = 0._wp
77         !
78         emp (:,:) = 0._wp
79         sfx (:,:) = 0._wp
80         qns (:,:) = 0._wp
81         qsr (:,:) = 0._wp
82         !         
83      ENDIF
84
85      IF( rn_uofac /= 0. ) THEN
86         
87         WHERE( ABS(gphit) <= rn_windszy/2. )
88            zwndrel(:,:) = rn_u10 - rn_uofac * uu(:,:,1,Kbb)
89         ELSEWHERE
90            zwndrel(:,:) =        - rn_uofac * uu(:,:,1,Kbb)
91         END WHERE
92         utau(:,:) = zrhocd * zwndrel(:,:) * zwndrel(:,:)
93
94         zwndrel(:,:) = - rn_uofac * vv(:,:,1,Kbb)
95         vtau(:,:) = zrhocd * zwndrel(:,:) * zwndrel(:,:)
96
97      ENDIF
98      !
99   END SUBROUTINE usrdef_sbc_oce
100
101   SUBROUTINE usrdef_sbc_ice_tau( kt )
102      INTEGER, INTENT(in) ::   kt   ! ocean time step
103   END SUBROUTINE usrdef_sbc_ice_tau
104
105
106   SUBROUTINE usrdef_sbc_ice_flx( kt, phs, phi )
107      INTEGER, INTENT(in) ::   kt   ! ocean time step
108      REAL(wp), DIMENSION(:,:,:), INTENT(in)  ::   phs    ! snow thickness
109      REAL(wp), DIMENSION(:,:,:), INTENT(in)  ::   phi    ! ice thickness
110   END SUBROUTINE usrdef_sbc_ice_flx
111
112   !!======================================================================
113END MODULE usrdef_sbc
Note: See TracBrowser for help on using the repository browser.