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 12928 for NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/tests/CANAL/MY_SRC/usrdef_sbc.F90 – NEMO

Ignore:
Timestamp:
2020-05-14T21:46:00+02:00 (4 years ago)
Author:
smueller
Message:

Synchronizing with /NEMO/trunk@12925 (ticket #2170)

Location:
NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser

    • Property svn:externals
      •  

        old new  
        66^/vendors/FCM@HEAD            ext/FCM 
        77^/vendors/IOIPSL@HEAD         ext/IOIPSL 
         8 
         9# SETTE 
         10^/utils/CI/sette@HEAD         sette 
  • NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/tests/CANAL/MY_SRC/usrdef_sbc.F90

    r10074 r12928  
    3131   PUBLIC   usrdef_sbc_ice_flx  ! routine called by icestp.F90 for ice thermo 
    3232 
    33    !! * Substitutions 
    34 #  include "vectopt_loop_substitute.h90" 
    3533   !!---------------------------------------------------------------------- 
    3634   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
     
    4038CONTAINS 
    4139 
    42    SUBROUTINE usrdef_sbc_oce( kt ) 
     40   SUBROUTINE usrdef_sbc_oce( kt, Kbb ) 
    4341      !!--------------------------------------------------------------------- 
    4442      !!                    ***  ROUTINE usr_def_sbc  *** 
     
    5452      !! 
    5553      !!---------------------------------------------------------------------- 
    56       INTEGER, INTENT(in) ::   kt   ! ocean time step 
    57        
     54      INTEGER, INTENT(in) ::   kt        ! ocean time step 
     55      INTEGER, INTENT(in) ::   Kbb       ! ocean time index 
    5856      INTEGER  ::   ji, jj               ! dummy loop indices 
    5957      REAL(wp) :: zrhoair = 1.22     ! approximate air density [Kg/m3] 
     
    8886          
    8987         WHERE( ABS(gphit) <= rn_windszy/2. ) 
    90             zwndrel(:,:) = rn_u10 - rn_uofac * un(:,:,1) 
     88            zwndrel(:,:) = rn_u10 - rn_uofac * uu(:,:,1,Kbb) 
    9189         ELSEWHERE 
    92             zwndrel(:,:) =        - rn_uofac * un(:,:,1) 
     90            zwndrel(:,:) =        - rn_uofac * uu(:,:,1,Kbb) 
    9391         END WHERE 
    9492         utau(:,:) = zrhocd * zwndrel(:,:) * zwndrel(:,:) 
    9593 
    96          zwndrel(:,:) = - rn_uofac * vn(:,:,1) 
     94         zwndrel(:,:) = - rn_uofac * vv(:,:,1,Kbb) 
    9795         vtau(:,:) = zrhocd * zwndrel(:,:) * zwndrel(:,:) 
    9896 
     
    105103   END SUBROUTINE usrdef_sbc_ice_tau 
    106104 
    107    SUBROUTINE usrdef_sbc_ice_flx( kt ) 
     105 
     106   SUBROUTINE usrdef_sbc_ice_flx( kt, phs, phi ) 
    108107      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 
    109110   END SUBROUTINE usrdef_sbc_ice_flx 
    110111 
Note: See TracChangeset for help on using the changeset viewer.