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.
Changeset 9403 for branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/CANAL/MY_SRC/usrdef_sbc.F90 – NEMO

Ignore:
Timestamp:
2018-03-15T15:57:42+01:00 (6 years ago)
Author:
smasson
Message:

dev_merge_2017: update test_cases CANAL

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/CANAL/MY_SRC/usrdef_sbc.F90

    r9302 r9403  
    33   !!                       ***  MODULE  usrdef_sbc  *** 
    44   !!  
    5    !!                      ===  VORTEX configuration  === 
     5   !!                      ===  CANAL configuration  === 
    66   !! 
    77   !! User defined :   surface forcing of a user configuration 
     
    1717   USE sbc_oce         ! Surface boundary condition: ocean fields 
    1818   USE phycst          ! physical constants 
    19    USE usrdef_nam, ONLY : rn_u10, rn_uofac, rn_ujetszy  
     19   USE usrdef_nam, ONLY : rn_u10, rn_uofac, rn_windszy  
    2020   ! 
    2121   USE in_out_manager  ! I/O manager 
     
    4747      !!              condition, i.e. the momentum, heat and freshwater fluxes. 
    4848      !! 
    49       !! ** Method  :   all 0 fields, for VORTEX case 
     49      !! ** Method  :   all 0 fields, for CANAL case 
    5050      !!                CAUTION : never mask the surface stress field ! 
    5151      !! 
     
    7171         ! 
    7272         utau(:,:) = 0._wp 
    73          IF( rn_u10 /= 0. .AND. rn_ujetszy > 0. ) THEN 
    74             WHERE( ABS(gphit) <= rn_ujetszy/2. ) utau(:,:) = zrhocd * rn_u10 * rn_u10 
     73         IF( rn_u10 /= 0. .AND. rn_windszy > 0. ) THEN 
     74            WHERE( ABS(gphit) <= rn_windszy/2. ) utau(:,:) = zrhocd * rn_u10 * rn_u10 
    7575         ENDIF 
    7676         vtau(:,:) = 0._wp 
     
    8787      IF( rn_uofac /= 0. ) THEN 
    8888          
    89          WHERE( ABS(gphit) <= rn_ujetszy/2. ) 
     89         WHERE( ABS(gphit) <= rn_windszy/2. ) 
    9090            zwndrel(:,:) = rn_u10 - rn_uofac * un(:,:,1) 
    9191         ELSEWHERE 
Note: See TracChangeset for help on using the changeset viewer.