- Timestamp:
- 2020-09-08T15:50:53+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
r11715 r13454 26 26 USE icectl ! sea-ice: control prints 27 27 USE bdy_oce , ONLY : ln_bdy 28 USE zdfdrg , ONLY: ln_drgice_imp 28 29 ! 29 30 USE in_out_manager ! I/O manager … … 323 324 REAL(wp) :: zat_u, zutau_ice, zu_t, zmodt ! local scalar 324 325 REAL(wp) :: zat_v, zvtau_ice, zv_t, zrhoco ! - - 326 REAL(wp) :: zflagi ! - - 325 327 !!--------------------------------------------------------------------- 326 328 IF( ln_timing ) CALL timing_start('ice_update_tau') … … 355 357 ! 356 358 ! !== every ocean time-step ==! 359 IF ( ln_drgice_imp ) THEN 360 ! Save drag with right sign to update top drag in the ocean implicit 361 ! friction 362 rCdU_ice(:,:) = -r1_rau0 * tmod_io(:,:) * at_i(:,:) * tmask(:,:,1) 363 zflagi = 0._wp 364 ELSE 365 zflagi = 1._wp 366 ENDIF 357 367 ! 358 368 DO jj = 2, jpjm1 !* update the stress WITHOUT an ice-ocean rotation angle … … 364 374 & / MAX( 1.0_wp , tmask(ji,jj,1) + tmask(ji ,jj+1,1) ) 365 375 ! ! linearized quadratic drag formulation 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) )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) ) 368 378 ! ! stresses at the ocean surface 369 379 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.