- Timestamp:
- 2020-09-16T18:37:46+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0.1_FKOSM_m11715/src/ICE/iceupdate.F90
r13454 r13478 26 26 USE icectl ! sea-ice: control prints 27 27 USE bdy_oce , ONLY : ln_bdy 28 USE zdfdrg , ONLY: ln_drgice_imp29 28 ! 30 29 USE in_out_manager ! I/O manager … … 324 323 REAL(wp) :: zat_u, zutau_ice, zu_t, zmodt ! local scalar 325 324 REAL(wp) :: zat_v, zvtau_ice, zv_t, zrhoco ! - - 326 REAL(wp) :: zflagi ! - -327 325 !!--------------------------------------------------------------------- 328 326 IF( ln_timing ) CALL timing_start('ice_update_tau') … … 357 355 ! 358 356 ! !== every ocean time-step ==! 359 IF ( ln_drgice_imp ) THEN360 ! Save drag with right sign to update top drag in the ocean implicit361 ! friction362 rCdU_ice(:,:) = -r1_rau0 * tmod_io(:,:) * at_i(:,:) * tmask(:,:,1)363 zflagi = 0._wp364 ELSE365 zflagi = 1._wp366 ENDIF367 357 ! 368 358 DO jj = 2, jpjm1 !* update the stress WITHOUT an ice-ocean rotation angle … … 374 364 & / MAX( 1.0_wp , tmask(ji,jj,1) + tmask(ji ,jj+1,1) ) 375 365 ! ! linearized quadratic drag formulation 376 zutau_ice = 0.5_wp * ( tmod_io(ji,jj) + tmod_io(ji+1,jj) ) * ( u_ice(ji,jj) - zflagi *pu_oce(ji,jj) )377 zvtau_ice = 0.5_wp * ( tmod_io(ji,jj) + tmod_io(ji,jj+1) ) * ( v_ice(ji,jj) - zflagi *pv_oce(ji,jj) )366 zutau_ice = 0.5_wp * ( tmod_io(ji,jj) + tmod_io(ji+1,jj) ) * ( u_ice(ji,jj) - pu_oce(ji,jj) ) 367 zvtau_ice = 0.5_wp * ( tmod_io(ji,jj) + tmod_io(ji,jj+1) ) * ( v_ice(ji,jj) - pv_oce(ji,jj) ) 378 368 ! ! stresses at the ocean surface 379 369 utau(ji,jj) = ( 1._wp - zat_u ) * utau_oce(ji,jj) + zat_u * zutau_ice
Note: See TracChangeset
for help on using the changeset viewer.