source: NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/tests/ISOMIP+/MY_SRC/usrdef_sbc.F90 @ 13159

Last change on this file since 13159 was 12353, checked in by acc, 9 months ago

Branch 2019/dev_r11943_MERGE_2019. Additions to the do loop macro implementation: converted a few loops previously missed because they used jpi-1 instead of jpim1 etc.; changed internal macro names in do_loop_substitute.h90 to strings that are much more unlikely to appear in any future code elsewhere and removed the key_vectopt_loop option (and all related code) since the do loop macros have suppressed this option. These changes have been fully SETTE-tested and this branch should now be ready to go back to the trunk.

File size: 3.7 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   !!----------------------------------------------------------------------
35   !! NEMO/OCE 4.0 , NEMO Consortium (2018)
36   !! $Id: usrdef_sbc.F90 10074 2018-08-28 16:15:49Z nicolasmartin $
37   !! Software governed by the CeCILL license (see ./LICENSE)
38   !!----------------------------------------------------------------------
39CONTAINS
40
41   SUBROUTINE usrdef_sbc_oce( kt, Kbb )
42      !!---------------------------------------------------------------------
43      !!                    ***  ROUTINE usr_def_sbc  ***
44      !!             
45      !! ** Purpose :   provide at each time-step the surface boundary
46      !!              condition, i.e. the momentum, heat and freshwater fluxes.
47      !!
48      !! ** Method  :   all 0 fields, for ISOMIP case
49      !!                CAUTION : never mask the surface stress field !
50      !!
51      !! ** Action  : - set to ZERO all the ocean surface boundary condition, i.e.   
52      !!                   utau, vtau, taum, wndm, qns, qsr, emp, sfx
53      !!
54      !!----------------------------------------------------------------------
55      INTEGER, INTENT(in) ::   kt   ! ocean time step
56      INTEGER, INTENT(in) ::   Kbb  ! ocean time index
57      !!---------------------------------------------------------------------
58      !
59      !
60      IF (kt == nit000) THEN
61         IF(lwp) WRITE(numout,*)' usr_sbc : ISOMIP case: NO surface forcing'
62         IF(lwp) WRITE(numout,*)' ~~~~~~~~~~~   utau = vtau = taum = wndm = qns = qsr = emp = sfx = 0'
63      END IF
64      !
65      utau(:,:) = 0._wp
66      vtau(:,:) = 0._wp
67      taum(:,:) = 0._wp
68      wndm(:,:) = 0._wp
69      !
70      emp (:,:) = 0._wp
71      sfx (:,:) = 0._wp
72      qns (:,:) = 0._wp
73      qsr (:,:) = 0._wp
74      !         
75   END SUBROUTINE usrdef_sbc_oce
76
77   SUBROUTINE usrdef_sbc_ice_tau( kt )
78      INTEGER, INTENT(in) ::   kt   ! ocean time step
79   END SUBROUTINE usrdef_sbc_ice_tau
80
81   SUBROUTINE usrdef_sbc_ice_flx( kt )
82      INTEGER, INTENT(in) ::   kt   ! ocean time step
83   END SUBROUTINE usrdef_sbc_ice_flx
84
85   !!======================================================================
86END MODULE usrdef_sbc
Note: See TracBrowser for help on using the repository browser.