New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
usrdef_sbc.F90 in NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/tests/ISOMIP/MY_SRC – NEMO

source: NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/tests/ISOMIP/MY_SRC/usrdef_sbc.F90 @ 10976

Last change on this file since 10976 was 10976, checked in by acc, 4 years ago

2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps : Convert usrdef routines in tests MY_SRC directories to enable SETTE compilation and testing of TEST-CASES (non-AGRIF varieties at least; enables ISOMIP, OVERFLOW and LOCK_EXCHANGE)

  • Property svn:keywords set to Id
File size: 3.9 KB
Line 
1MODULE usrdef_sbc
2   !!======================================================================
3   !!                     ***  MODULE  usrdef_sbc  ***
4   !!
5   !!                  ===  ISOMIP configuration  ===
6   !!
7   !! User defined :   surface forcing of a user configuration
8   !!======================================================================
9   !! History :  4.0   ! 2016-03  (S. Flavoni, G. Madec)   user defined interface
10   !!                  ! 2017-02  (P. Mathiot, S. Flavoni) adapt code to ISOMIP case
11   !!----------------------------------------------------------------------
12
13   !!----------------------------------------------------------------------
14   !!   usr_def_sbc    : user defined surface bounday conditions in ISOMIP case
15   !!----------------------------------------------------------------------
16   USE oce             ! ocean dynamics and tracers
17   USE dom_oce         ! ocean space and time domain
18   USE sbc_oce         ! Surface boundary condition: ocean fields
19   USE sbc_ice         ! Surface boundary condition: ice fields
20   USE phycst          ! physical constants
21   !
22   USE in_out_manager  ! I/O manager
23   USE lib_mpp         ! distribued memory computing library
24   USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
25   USE lib_fortran     ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined)
26
27   IMPLICIT NONE
28   PRIVATE
29
30   PUBLIC   usrdef_sbc_oce      ! routine called in sbcmod module
31   PUBLIC   usrdef_sbc_ice_tau  ! routine called by icestp.F90 for ice dynamics
32   PUBLIC   usrdef_sbc_ice_flx  ! routine called by icestp.F90 for ice thermo
33
34   !! * Substitutions
35#  include "vectopt_loop_substitute.h90"
36   !!----------------------------------------------------------------------
37   !! NEMO/OCE 4.0 , NEMO Consortium (2018)
38   !! $Id$
39   !! Software governed by the CeCILL license (see ./LICENSE)
40   !!----------------------------------------------------------------------
41CONTAINS
42
43   SUBROUTINE usrdef_sbc_oce( kt, Kbb )
44      !!---------------------------------------------------------------------
45      !!                    ***  ROUTINE usr_def_sbc  ***
46      !!             
47      !! ** Purpose :   provide at each time-step the surface boundary
48      !!              condition, i.e. the momentum, heat and freshwater fluxes.
49      !!
50      !! ** Method  :   all 0 fields, for ISOMIP case
51      !!                CAUTION : never mask the surface stress field !
52      !!
53      !! ** Action  : - set to ZERO all the ocean surface boundary condition, i.e.   
54      !!                   utau, vtau, taum, wndm, qns, qsr, emp, sfx
55      !!
56      !!----------------------------------------------------------------------
57      INTEGER, INTENT(in) ::   kt   ! ocean time step
58      INTEGER, INTENT(in) ::   Kbb  ! ocean time index
59      !!---------------------------------------------------------------------
60      !
61      IF( kt == nit000 ) THEN
62         !
63         IF(lwp) WRITE(numout,*)' usr_sbc : ISOMIP case: NO surface forcing'
64         IF(lwp) WRITE(numout,*)' ~~~~~~~~~~~   utau = vtau = taum = wndm = qns = qsr = emp = sfx = 0'
65         !
66         utau(:,:) = 0._wp
67         vtau(:,:) = 0._wp
68         taum(:,:) = 0._wp
69         wndm(:,:) = 0._wp
70         !
71         emp (:,:) = 0._wp
72         sfx (:,:) = 0._wp
73         qns (:,:) = 0._wp
74         qsr (:,:) = 0._wp
75         !         
76      ENDIF
77      !
78   END SUBROUTINE usrdef_sbc_oce
79
80   SUBROUTINE usrdef_sbc_ice_tau( kt )
81      INTEGER, INTENT(in) ::   kt   ! ocean time step
82   END SUBROUTINE usrdef_sbc_ice_tau
83
84
85   SUBROUTINE usrdef_sbc_ice_flx( kt, phs, phi )
86      INTEGER, INTENT(in) ::   kt   ! ocean time step
87      REAL(wp), DIMENSION(:,:,:), INTENT(in)  ::   phs    ! snow thickness
88      REAL(wp), DIMENSION(:,:,:), INTENT(in)  ::   phi    ! ice thickness
89   END SUBROUTINE usrdef_sbc_ice_flx
90
91   !!======================================================================
92END MODULE usrdef_sbc
Note: See TracBrowser for help on using the repository browser.