Changeset 14766 for NEMO/releases/r4.0/r4.0-HEAD/src/OCE/BDY/bdyice.F90
- Timestamp:
- 2021-04-30T11:08:19+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/releases/r4.0/r4.0-HEAD/src/OCE/BDY/bdyice.F90
r13589 r14766 153 153 h_i (ji,jj, jl) = ( h_i (ji,jj, jl) * zwgt1 + dta%h_i(i_bdy,jl) * zwgt ) * tmask(ji,jj,1) ! Ice depth 154 154 h_s (ji,jj, jl) = ( h_s (ji,jj, jl) * zwgt1 + dta%h_s(i_bdy,jl) * zwgt ) * tmask(ji,jj,1) ! Snow depth 155 t_i (ji,jj,:,jl) = ( t_i (ji,jj,:,jl) * zwgt1 + dta%t_i(i_bdy,jl) * zwgt )* tmask(ji,jj,1) ! Ice temperature156 t_s (ji,jj,:,jl) = ( t_s (ji,jj,:,jl) * zwgt1 + dta%t_s(i_bdy,jl) * zwgt )* tmask(ji,jj,1) ! Snow temperature157 t_su(ji,jj, jl) = ( t_su(ji,jj, jl) * zwgt1 + dta%tsu(i_bdy,jl) * zwgt )* tmask(ji,jj,1) ! Surf temperature158 s_i (ji,jj, jl) = ( s_i (ji,jj, jl) * zwgt1 + dta%s_i(i_bdy,jl) * zwgt )* tmask(ji,jj,1) ! Ice salinity155 t_i (ji,jj,:,jl) = dta%t_i(i_bdy,jl) * tmask(ji,jj,1) ! Ice temperature 156 t_s (ji,jj,:,jl) = dta%t_s(i_bdy,jl) * tmask(ji,jj,1) ! Snow temperature 157 t_su(ji,jj, jl) = dta%tsu(i_bdy,jl) * tmask(ji,jj,1) ! Surf temperature 158 s_i (ji,jj, jl) = dta%s_i(i_bdy,jl) * tmask(ji,jj,1) ! Ice salinity 159 159 a_ip(ji,jj, jl) = ( a_ip(ji,jj, jl) * zwgt1 + dta%aip(i_bdy,jl) * zwgt ) * tmask(ji,jj,1) ! Ice pond concentration 160 160 h_ip(ji,jj, jl) = ( h_ip(ji,jj, jl) * zwgt1 + dta%hip(i_bdy,jl) * zwgt ) * tmask(ji,jj,1) ! Ice pond depth … … 363 363 ! i-1 i i | ! i i i+1 | ! i i i+1 | 364 364 ! > ice > | ! o > ice | ! o > o | 365 ! => set at u_ice(i-1) ! => set to O! => unchanged365 ! => set at u_ice(i-1) ! => set to u_oce ! => unchanged 366 366 IF( zflag == -1. .AND. ji > 1 .AND. ji < jpi ) THEN 367 367 IF ( vt_i(ji ,jj) > 0. ) THEN ; u_ice(ji,jj) = u_ice(ji-1,jj) 368 ELSEIF( vt_i(ji+1,jj) > 0. ) THEN ; u_ice(ji,jj) = 0._wp368 ELSEIF( vt_i(ji+1,jj) > 0. ) THEN ; u_ice(ji,jj) = u_oce(ji,jj) 369 369 END IF 370 370 END IF 371 371 ! | i i+1 i+1 ! | i i i+1 ! | i i i+1 372 372 ! | > ice > ! | ice > o ! | o > o 373 ! => set at u_ice(i+1) ! => set to O! => unchanged373 ! => set at u_ice(i+1) ! => set to u_oce ! => unchanged 374 374 IF( zflag == 1. .AND. ji+1 < jpi+1 ) THEN 375 375 IF ( vt_i(ji+1,jj) > 0. ) THEN ; u_ice(ji,jj) = u_ice(ji+1,jj) 376 ELSEIF( vt_i(ji ,jj) > 0. ) THEN ; u_ice(ji,jj) = 0._wp376 ELSEIF( vt_i(ji ,jj) > 0. ) THEN ; u_ice(ji,jj) = u_oce(ji,jj) 377 377 END IF 378 378 END IF … … 395 395 ! ice (jj ) ! o (jj ) ! o (jj ) 396 396 ! ^ (jj-1) ! ! 397 ! => set to u_ice(jj-1) ! => set to 0! => unchanged397 ! => set to u_ice(jj-1) ! => set to v_oce ! => unchanged 398 398 IF( zflag == -1. .AND. jj > 1 .AND. jj < jpj ) THEN 399 399 IF ( vt_i(ji,jj ) > 0. ) THEN ; v_ice(ji,jj) = v_ice(ji,jj-1) 400 ELSEIF( vt_i(ji,jj+1) > 0. ) THEN ; v_ice(ji,jj) = 0._wp400 ELSEIF( vt_i(ji,jj+1) > 0. ) THEN ; v_ice(ji,jj) = v_oce(ji,jj) 401 401 END IF 402 402 END IF … … 405 405 ! ^ (jj ) ! ^ (jj ) ! ^ (jj ) 406 406 ! ________________ ! ____ice___(jj )_ ! _____o____(jj ) 407 ! => set to u_ice(jj+1) ! => set to 0! => unchanged407 ! => set to u_ice(jj+1) ! => set to v_oce ! => unchanged 408 408 IF( zflag == 1. .AND. jj < jpj ) THEN 409 409 IF ( vt_i(ji,jj+1) > 0. ) THEN ; v_ice(ji,jj) = v_ice(ji,jj+1) 410 ELSEIF( vt_i(ji,jj ) > 0. ) THEN ; v_ice(ji,jj) = 0._wp410 ELSEIF( vt_i(ji,jj ) > 0. ) THEN ; v_ice(ji,jj) = v_oce(ji,jj) 411 411 END IF 412 412 END IF
Note: See TracChangeset
for help on using the changeset viewer.