- Timestamp:
- 2021-04-07T19:16:18+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2021/dev_r14273_HPC-02_Daley_Tiling/src/OCE/SBC/sbcrnf.F90
r14072 r14680 206 206 IF( ln_rnf_depth .OR. ln_rnf_depth_ini ) THEN !== runoff distributed over several levels ==! 207 207 IF( ln_linssh ) THEN !* constant volume case : just apply the runoff input flow 208 DO_2D( 1, 1, 1, 1 ) 208 ! [comm_cleanup] ! DO_2D( 1, 1, 1, 1 ) 209 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 209 210 DO jk = 1, nk_rnf(ji,jj) 210 211 phdivn(ji,jj,jk) = phdivn(ji,jj,jk) - ( rnf(ji,jj) + rnf_b(ji,jj) ) * zfact * r1_rho0 / h_rnf(ji,jj) … … 212 213 END_2D 213 214 ELSE !* variable volume case 214 DO_2D( 1, 1, 1, 1 ) ! update the depth over which runoffs are distributed 215 ! [comm_cleanup] ! DO_2D( 1, 1, 1, 1 ) ! update the depth over which runoffs are distributed 216 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) ! update the depth over which runoffs are distributed 215 217 h_rnf(ji,jj) = 0._wp 216 218 DO jk = 1, nk_rnf(ji,jj) ! recalculates h_rnf to be the depth in metres … … 358 360 ! 359 361 nk_rnf(:,:) = 0 ! set the number of level over which river runoffs are applied 360 DO_2D( 1, 1, 1, 1 ) 362 ! [comm_cleanup] ! DO_2D( 1, 1, 1, 1 ) 363 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 361 364 IF( h_rnf(ji,jj) > 0._wp ) THEN 362 365 jk = 2 … … 371 374 ENDIF 372 375 END_2D 373 DO_2D( 1, 1, 1, 1 ) ! set the associated depth 376 ! [comm_cleanup] ! DO_2D( 1, 1, 1, 1 ) ! set the associated depth 377 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) ! set the associated depth 374 378 h_rnf(ji,jj) = 0._wp 375 379 DO jk = 1, nk_rnf(ji,jj) … … 401 405 WHERE( zrnfcl(:,:,1) > 0._wp ) h_rnf(:,:) = zacoef * zrnfcl(:,:,1) ! compute depth for all runoffs 402 406 ! 403 DO_2D( 1, 1, 1, 1 ) ! take in account min depth of ocean rn_hmin 407 ! [comm_cleanup] ! DO_2D( 1, 1, 1, 1 ) ! take in account min depth of ocean rn_hmin 408 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) ! take in account min depth of ocean rn_hmin 404 409 IF( zrnfcl(ji,jj,1) > 0._wp ) THEN 405 410 jk = mbkt(ji,jj) … … 409 414 ! 410 415 nk_rnf(:,:) = 0 ! number of levels on which runoffs are distributed 411 DO_2D( 1, 1, 1, 1 ) 416 ! [comm_cleanup] ! DO_2D( 1, 1, 1, 1 ) 417 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 412 418 IF( zrnfcl(ji,jj,1) > 0._wp ) THEN 413 419 jk = 2 … … 420 426 END_2D 421 427 ! 422 DO_2D( 1, 1, 1, 1 ) ! set the associated depth 428 ! [comm_cleanup] ! DO_2D( 1, 1, 1, 1 ) ! set the associated depth 429 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) ! set the associated depth 423 430 h_rnf(ji,jj) = 0._wp 424 431 DO jk = 1, nk_rnf(ji,jj)
Note: See TracChangeset
for help on using the changeset viewer.