Changeset 3978
- Timestamp:
- 2013-07-16T15:56:48+02:00 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limrhg.F90
r3938 r3978 447 447 CALL lbc_lnk( zs12(:,:), 'F', 1. ) 448 448 449 !#if defined key_bdy 450 ! ! clem: change zs1, zs2, zs12 at the boundary for each iteration 451 ! CALL bdy_ice_lim_dyn( 2, zs1, zs2, zs12 ) 452 ! CALL lbc_lnk( zs1 (:,:), 'T', 1. ) 453 ! CALL lbc_lnk( zs2 (:,:), 'T', 1. ) 454 ! CALL lbc_lnk( zs12(:,:), 'F', 1. ) 455 !#endif 456 449 457 ! Ice internal stresses (Appendix C of Hunke and Dukowicz, 2002) 450 458 DO jj = k_j1+1, k_jpj-1 … … 578 586 !#if defined key_bdy 579 587 ! ! clem: change u_ice and v_ice at the boundary for each iteration 580 ! CALL bdy_ice_lim_dyn( )588 ! CALL bdy_ice_lim_dyn( 1 ) 581 589 !#endif 582 590 … … 607 615 !CDIR NOVERRCHK 608 616 DO ji = fs_2, fs_jpim1 609 zindb = MAX( 0.0, SIGN( 1.0, at_i(ji,jj) - 1.0e-6 ) ) 610 zdummy = zindb * vt_i(ji,jj) / MAX(at_i(ji,jj) , 1.0e-06 ) 617 !zindb = MAX( 0.0, SIGN( 1.0, at_i(ji,jj) - 1.0e-6 ) ) 618 !zdummy = zindb * vt_i(ji,jj) / MAX(at_i(ji,jj) , 1.0e-06 ) 619 zdummy = vt_i(ji,jj) 611 620 IF ( zdummy .LE. hminrhg ) THEN 612 621 u_ice(ji,jj) = u_oce(ji,jj) … … 619 628 CALL lbc_lnk( v_ice(:,:), 'V', -1. ) 620 629 630 #if defined key_bdy 621 631 ! clem: change u_ice and v_ice at the boundary 622 #if defined key_bdy 623 CALL bdy_ice_lim_dyn() 632 CALL bdy_ice_lim_dyn( 1 ) 624 633 #endif 625 634 626 635 DO jj = k_j1+1, k_jpj-1 627 636 DO ji = fs_2, fs_jpim1 628 zindb = MAX( 0.0, SIGN( 1.0, at_i(ji,jj) - 1.0e-6 ) ) 629 zdummy = zindb * vt_i(ji,jj) / MAX(at_i(ji,jj) , 1.0e-06 ) 637 !zindb = MAX( 0.0, SIGN( 1.0, at_i(ji,jj) - 1.0e-6 ) ) 638 !zdummy = zindb * vt_i(ji,jj) / MAX(at_i(ji,jj) , 1.0e-06 ) 639 zdummy = vt_i(ji,jj) 630 640 IF ( zdummy .LE. hminrhg ) THEN 631 641 v_ice1(ji,jj) = 0.5*( (v_ice(ji,jj)+v_ice(ji,jj-1))*e1t(ji+1,jj) & … … 650 660 !- zdd(:,:), zdt(:,:): divergence and tension at centre 651 661 !- zds(:,:): shear on northeast corner of grid cells 652 zindb = MAX( 0.0, SIGN( 1.0, at_i(ji,jj) - 1.0e-6 ) )653 zdummy = zindb * vt_i(ji,jj) / MAX(at_i(ji,jj) , 1.0e-06 )654 662 !zindb = MAX( 0.0, SIGN( 1.0, at_i(ji,jj) - 1.0e-6 ) ) 663 !zdummy = zindb * vt_i(ji,jj) / MAX(at_i(ji,jj) , 1.0e-06 ) 664 zdummy = vt_i(ji,jj) 655 665 IF ( zdummy .LE. hminrhg ) THEN 656 666
Note: See TracChangeset
for help on using the changeset viewer.