- Timestamp:
- 2020-07-04T10:22:08+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/iceupdate.F90
r12785 r13249 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 … … 321 322 REAL(wp) :: zat_u, zutau_ice, zu_t, zmodt ! local scalar 322 323 REAL(wp) :: zat_v, zvtau_ice, zv_t, zrhoco ! - - 324 REAL(wp) :: zflagi ! - - 323 325 !!--------------------------------------------------------------------- 324 326 IF( ln_timing ) CALL timing_start('ice_update_tau') … … 353 355 ! 354 356 ! !== every ocean time-step ==! 357 IF ( ln_drgice_imp ) THEN 358 ! Save drag with right sign to update top drag in the ocean implicit friction 359 rCdU_ice(:,:) = -r1_rau0 * tmod_io(:,:) * at_i(:,:) * tmask(:,:,1) 360 zflagi = 0._wp 361 ELSE 362 zflagi = 1._wp 363 ENDIF 355 364 ! 356 365 DO jj = 2, jpjm1 !* update the stress WITHOUT an ice-ocean rotation angle … … 362 371 & / MAX( 1.0_wp , tmask(ji,jj,1) + tmask(ji ,jj+1,1) ) 363 372 ! ! linearized quadratic drag formulation 364 zutau_ice = 0.5_wp * ( tmod_io(ji,jj) + tmod_io(ji+1,jj) ) * ( u_ice(ji,jj) - pu_oce(ji,jj) )365 zvtau_ice = 0.5_wp * ( tmod_io(ji,jj) + tmod_io(ji,jj+1) ) * ( v_ice(ji,jj) - pv_oce(ji,jj) )373 zutau_ice = 0.5_wp * ( tmod_io(ji,jj) + tmod_io(ji+1,jj) ) * ( u_ice(ji,jj) - zflagi * pu_oce(ji,jj) ) 374 zvtau_ice = 0.5_wp * ( tmod_io(ji,jj) + tmod_io(ji,jj+1) ) * ( v_ice(ji,jj) - zflagi * pv_oce(ji,jj) ) 366 375 ! ! stresses at the ocean surface 367 376 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.