Changeset 12377 for NEMO/trunk/tests/CANAL/MY_SRC/usrdef_sbc.F90
- Timestamp:
- 2020-02-12T15:39:06+01:00 (4 years ago)
- Location:
- NEMO/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk
- Property svn:externals
-
old new 3 3 ^/utils/build/mk@HEAD mk 4 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev @HEAD ext/AGRIF5 ^/vendors/AGRIF/dev_r11615_ENHANCE-04_namelists_as_internalfiles_agrif@HEAD ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL
-
- Property svn:externals
-
NEMO/trunk/tests/CANAL/MY_SRC/usrdef_sbc.F90
r10074 r12377 31 31 PUBLIC usrdef_sbc_ice_flx ! routine called by icestp.F90 for ice thermo 32 32 33 !! * Substitutions34 # include "vectopt_loop_substitute.h90"35 33 !!---------------------------------------------------------------------- 36 34 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 40 38 CONTAINS 41 39 42 SUBROUTINE usrdef_sbc_oce( kt )40 SUBROUTINE usrdef_sbc_oce( kt, Kmm, Kbb ) 43 41 !!--------------------------------------------------------------------- 44 42 !! *** ROUTINE usr_def_sbc *** … … 54 52 !! 55 53 !!---------------------------------------------------------------------- 56 INTEGER, INTENT(in) :: kt ! ocean time step57 54 INTEGER, INTENT(in) :: kt ! ocean time step 55 INTEGER, INTENT(in) :: Kbb, Kmm ! ocean time index 58 56 INTEGER :: ji, jj ! dummy loop indices 59 57 REAL(wp) :: zrhoair = 1.22 ! approximate air density [Kg/m3] … … 88 86 89 87 WHERE( ABS(gphit) <= rn_windszy/2. ) 90 zwndrel(:,:) = rn_u10 - rn_uofac * u n(:,:,1)88 zwndrel(:,:) = rn_u10 - rn_uofac * uu(:,:,1,Kmm) 91 89 ELSEWHERE 92 zwndrel(:,:) = - rn_uofac * u n(:,:,1)90 zwndrel(:,:) = - rn_uofac * uu(:,:,1,Kmm) 93 91 END WHERE 94 92 utau(:,:) = zrhocd * zwndrel(:,:) * zwndrel(:,:) 95 93 96 zwndrel(:,:) = - rn_uofac * v n(:,:,1)94 zwndrel(:,:) = - rn_uofac * vv(:,:,1,Kmm) 97 95 vtau(:,:) = zrhocd * zwndrel(:,:) * zwndrel(:,:) 98 96 … … 105 103 END SUBROUTINE usrdef_sbc_ice_tau 106 104 107 SUBROUTINE usrdef_sbc_ice_flx( kt ) 105 106 SUBROUTINE usrdef_sbc_ice_flx( kt, phs, phi ) 108 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 109 110 END SUBROUTINE usrdef_sbc_ice_flx 110 111
Note: See TracChangeset
for help on using the changeset viewer.