Changeset 5058 for branches/2015
- Timestamp:
- 2015-02-04T12:25:57+01:00 (10 years ago)
- Location:
- branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_2/ice_2.F90
r5053 r5058 56 56 REAL(wp), PUBLIC :: sangvg, cangvg !: sin and cos of the turning angle for ocean stress 57 57 REAL(wp), PUBLIC :: pstarh !: pstar / 2.0 58 59 ! !!** switch for presence of ice or not 60 REAL(wp), PUBLIC :: rswitch 58 61 59 62 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: ahiu , ahiv !: hor. diffusivity coeff. at ocean U- and V-points (m2/s) -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/ice.F90
r5055 r5058 164 164 REAL(wp), PUBLIC :: r1_rdtice !: = 1. / rdt_ice 165 165 166 ! Clem je sais pas ou les mettre167 INTEGER , PUBLIC :: nn_ hibnd !: thickness category boundaries: tanh function (1), or h^(-alpha) function (2)168 REAL(wp), PUBLIC :: rn_ hibnd !: mean thickness of the domain (used to compute category boundaries, nn_hibnd= 2 only)166 ! !!** ice-thickness distribution namelist (namiceitd) ** 167 INTEGER , PUBLIC :: nn_itdshp !: categories distribution following: tanh function (1), or h^(-alpha) function (2) 168 REAL(wp), PUBLIC :: rn_itmean !: mean thickness of the domain (used to compute the distribution, nn_itdshp = 2 only) 169 169 170 170 ! !!** ice-dynamics namelist (namicedyn) ** -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90
r5055 r5058 423 423 !! ** Purpose : Initializes the ice thickness distribution 424 424 !! ** Method : ... 425 !! Note : hi_max(jpl) is here set up to a value close to 7 m for 426 !! limistate (only) and is changed to 99 m in sbc_lim_init 427 !!------------------------------------------------------------------ 425 !! ** input : Namelist namiceitd 426 !!------------------------------------------------------------------- 427 INTEGER :: ios ! Local integer output status for namelist read 428 NAMELIST/namiceitd/ nn_itdshp, rn_itmean 429 ! 428 430 INTEGER :: jl ! dummy loop index 429 431 REAL(wp) :: zc1, zc2, zc3, zx1 ! local scalars 430 432 REAL(wp) :: zhmax, znum, zden, zalpha ! 431 433 !!------------------------------------------------------------------ 432 433 IF(lwp) WRITE(numout,*) 434 IF(lwp) WRITE(numout,*) 'lim_itd_init : Initialization of ice thickness distribution ' 435 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 436 437 !------------------------------------------------------------------------------! 438 ! 1) Ice thickness distribution parameters initialization 439 !------------------------------------------------------------------------------! 440 IF(lwp) WRITE(numout,*) ' Number of ice categories jpl = ', jpl 441 434 ! 435 REWIND( numnam_ice_ref ) ! Namelist namiceitd in reference namelist : Parameters for ice 436 READ ( numnam_ice_ref, namiceitd, IOSTAT = ios, ERR = 903) 437 903 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namiceitd in reference namelist', lwp ) 438 439 REWIND( numnam_ice_cfg ) ! Namelist namiceitd in configuration namelist : Parameters for ice 440 READ ( numnam_ice_cfg, namiceitd, IOSTAT = ios, ERR = 904 ) 441 904 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namiceitd in configuration namelist', lwp ) 442 IF(lwm) WRITE ( numoni, namiceitd ) 443 ! 444 ! 445 IF(lwp) THEN ! control print 446 WRITE(numout,*) 447 WRITE(numout,*) 'ice_itd : ice cat distribution' 448 WRITE(numout,*) ' ~~~~~~' 449 WRITE(numout,*) ' shape of ice categories distribution nn_itdshp = ', nn_itdshp 450 WRITE(numout,*) ' mean ice thickness in the domain (only active if nn_itdshp=2) rn_itmean = ', rn_itmean 451 ENDIF 452 453 !---------------------------------- 442 454 !- Thickness categories boundaries 443 455 !---------------------------------- 444 ! Clem - je sais pas encore dans quelle namelist les mettre, ca depend des chgts liés à bdy445 nn_hibnd = 2 ! thickness category boundaries: tanh function (1) h^(-alpha) (2)446 rn_hibnd = 2.5 ! mean thickness of the domain (used to compute category boundaries, nn_hibnd = 2 only)456 IF(lwp) WRITE(numout,*) 457 IF(lwp) WRITE(numout,*) 'lim_itd_init : Initialization of ice cat distribution ' 458 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 447 459 448 460 hi_max(:) = 0._wp 449 461 450 SELECT CASE ( nn_ hibnd)462 SELECT CASE ( nn_itdshp ) 451 463 !---------------------- 452 464 CASE (1) ! tanh function (CICE) … … 466 478 zalpha = 0.05 ! exponent of the transform function 467 479 468 zhmax = 3.*rn_ hibnd480 zhmax = 3.*rn_itmean 469 481 470 482 DO jl = 1, jpl … … 479 491 hi_mean(jl) = ( hi_max(jl) + hi_max(jl-1) ) * 0.5_wp 480 492 END DO 493 481 494 ! Set hi_max(jpl) to a big value to ensure that all ice is thinner than hi_max(jpl) 482 495 hi_max(jpl) = 99._wp
Note: See TracChangeset
for help on using the changeset viewer.