Changeset 3934
- Timestamp:
- 2013-06-21T11:04:14+02:00 (11 years ago)
- Location:
- branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/CONFIG/ORCA2_LIM3
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/namelist_ice_lim3
r3931 r3934 62 62 telast =9600.0 ! timescale for elastic waves, SB, 720.0 63 63 alphaevp = 1.0 ! coefficient for the solution of internal ice stresses 64 hminrhg = 0.05 ! ice thickness (m) below which ice velocity equal ocean velocity 64 65 / 65 66 !----------------------------------------------------------------------- -
branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/CONFIG/ORCA2_LIM3/MY_SRC/ice.F90
r3931 r3934 187 187 REAL(wp), PUBLIC :: alphaevp = 1._wp !: coeficient of the internal stresses !SB 188 188 REAL(wp), PUBLIC :: unit_fac = 1.e+09_wp !: conversion factor for ice / snow enthalpy 189 REAL(wp), PUBLIC :: hminrhg = 0.05_wp !: clem : ice thickness (in m) below which ice velocity is set to ocean velocity 189 190 190 191 ! !!** ice-salinity namelist (namicesal) ** -
branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/CONFIG/ORCA2_LIM3/MY_SRC/limdyn.F90
r3931 r3934 276 276 & dm, nbiter, nbitdr, om, resl, cw, angvg, pstar, & 277 277 & c_rhg, etamn, creepl, ecc, ahi0, & 278 & nevp, telast, alphaevp 278 & nevp, telast, alphaevp, hminrhg 279 279 !!------------------------------------------------------------------- 280 280 … … 304 304 WRITE(numout,*) ' timescale for elastic waves telast = ', telast 305 305 WRITE(numout,*) ' coefficient for the solution of int. stresses alphaevp = ', alphaevp 306 WRITE(numout,*) ' min ice thickness for rheology calculations hminrhg = ', hminrhg 306 307 ENDIF 307 308 ! -
branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/CONFIG/ORCA2_LIM3/MY_SRC/limrhg.F90
r3931 r3934 36 36 #endif 37 37 USE lib_fortran ! to use key_nosignedzero 38 39 #if defined key_bdy 40 USE bdyice_lim 41 #endif 38 42 39 43 IMPLICIT NONE … … 571 575 572 576 ENDIF 577 578 !#if defined key_bdy 579 ! ! clem: change u_ice and v_ice at the boundary for each iteration 580 ! CALL bdy_ice_lim_dyn() 581 !#endif 573 582 574 583 IF(ln_ctl) THEN … … 582 591 ENDIF 583 592 584 ! 593 ! ! ==================== ! 585 594 END DO ! end loop over jter ! 586 595 ! ! ==================== ! 587 588 596 ! 589 597 !------------------------------------------------------------------------------! 590 598 ! 4) Prevent ice velocities when the ice is thin 591 599 !------------------------------------------------------------------------------! 592 !clem : change zdummy from 5.e-2 to 1.e-2593 ! 594 ! If the ice thickness is below 1cmthen ice velocity should equal the600 !clem : add hminrhg in the namelist 601 ! 602 ! If the ice thickness is below hminrhg (5cm) then ice velocity should equal the 595 603 ! ocean velocity, 596 604 ! This prevents high velocity when ice is thin … … 601 609 zindb = MAX( 0.0, SIGN( 1.0, at_i(ji,jj) - 1.0e-6 ) ) 602 610 zdummy = zindb * vt_i(ji,jj) / MAX(at_i(ji,jj) , 1.0e-06 ) 603 IF ( zdummy .LE. 5.0e-2) THEN611 IF ( zdummy .LE. hminrhg ) THEN 604 612 u_ice(ji,jj) = u_oce(ji,jj) 605 613 v_ice(ji,jj) = v_oce(ji,jj) … … 610 618 CALL lbc_lnk( u_ice(:,:), 'U', -1. ) 611 619 CALL lbc_lnk( v_ice(:,:), 'V', -1. ) 620 621 ! clem: change u_ice and v_ice at the boundary 622 #if defined key_bdy 623 CALL bdy_ice_lim_dyn() 624 #endif 612 625 613 626 DO jj = k_j1+1, k_jpj-1 … … 615 628 zindb = MAX( 0.0, SIGN( 1.0, at_i(ji,jj) - 1.0e-6 ) ) 616 629 zdummy = zindb * vt_i(ji,jj) / MAX(at_i(ji,jj) , 1.0e-06 ) 617 IF ( zdummy .LE. 5.0e-2) THEN630 IF ( zdummy .LE. hminrhg ) THEN 618 631 v_ice1(ji,jj) = 0.5*( (v_ice(ji,jj)+v_ice(ji,jj-1))*e1t(ji+1,jj) & 619 632 & +(v_ice(ji+1,jj)+v_ice(ji+1,jj-1))*e1t(ji,jj)) & … … 640 653 zdummy = zindb * vt_i(ji,jj) / MAX(at_i(ji,jj) , 1.0e-06 ) 641 654 642 IF ( zdummy .LE. 5.0e-2) THEN655 IF ( zdummy .LE. hminrhg ) THEN 643 656 644 657 zdd(ji,jj) = ( e2u(ji,jj)*u_ice(ji,jj) &
Note: See TracChangeset
for help on using the changeset viewer.