- Timestamp:
- 2020-09-15T09:27:47+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/temporary_r4_trunk/src/ICE/iceitd.F90
r11732 r13466 47 47 LOGICAL :: ln_cat_usr ! ice categories are defined by rn_catbnd 48 48 REAL(wp), DIMENSION(0:100) :: rn_catbnd ! ice categories bounds 49 REAL(wp) :: rn_himax ! maximum ice thickness allowed 49 50 ! 50 51 !!---------------------------------------------------------------------- … … 304 305 IF ( a_i_1d(ji) > epsi10 .AND. h_i_1d(ji) < rn_himin ) THEN 305 306 a_i_1d(ji) = a_i_1d(ji) * h_i_1d(ji) / rn_himin 306 IF( ln_pnd_ H12) a_ip_1d(ji) = a_ip_1d(ji) * h_i_1d(ji) / rn_himin307 IF( ln_pnd_LEV ) a_ip_1d(ji) = a_ip_1d(ji) * h_i_1d(ji) / rn_himin 307 308 h_i_1d(ji) = rn_himin 308 309 ENDIF … … 410 411 CALL tab_3d_2d( npti, nptidx(1:npti), a_ip_2d(1:npti,1:jpl), a_ip ) 411 412 CALL tab_3d_2d( npti, nptidx(1:npti), v_ip_2d(1:npti,1:jpl), v_ip ) 413 CALL tab_3d_2d( npti, nptidx(1:npti), v_il_2d(1:npti,1:jpl), v_il ) 412 414 CALL tab_3d_2d( npti, nptidx(1:npti), t_su_2d(1:npti,1:jpl), t_su ) 413 415 DO jl = 1, jpl … … 474 476 zaTsfn(ji,jl2) = zaTsfn(ji,jl2) + ztrans 475 477 ! 476 IF ( ln_pnd_ H12) THEN478 IF ( ln_pnd_LEV ) THEN 477 479 ztrans = a_ip_2d(ji,jl1) * zworka(ji) ! Pond fraction 478 480 a_ip_2d(ji,jl1) = a_ip_2d(ji,jl1) - ztrans … … 482 484 v_ip_2d(ji,jl1) = v_ip_2d(ji,jl1) - ztrans 483 485 v_ip_2d(ji,jl2) = v_ip_2d(ji,jl2) + ztrans 486 ! 487 IF ( ln_pnd_lids ) THEN ! Pond lid volume 488 ztrans = v_il_2d(ji,jl1) * zworka(ji) 489 v_il_2d(ji,jl1) = v_il_2d(ji,jl1) - ztrans 490 v_il_2d(ji,jl2) = v_il_2d(ji,jl2) + ztrans 491 ENDIF 484 492 ENDIF 485 493 ! … … 526 534 ! clem: The transfer between one category to another can lead to very small negative values (-1.e-20) 527 535 ! because of truncation error ( i.e. 1. - 1. /= 0 ) 528 CALL ice_var_roundoff( a_i_2d, v_i_2d, v_s_2d, sv_i_2d, oa_i_2d, a_ip_2d, v_ip_2d, ze_s_2d, ze_i_2d )536 CALL ice_var_roundoff( a_i_2d, v_i_2d, v_s_2d, sv_i_2d, oa_i_2d, a_ip_2d, v_ip_2d, v_il_2d, ze_s_2d, ze_i_2d ) 529 537 530 538 ! at_i must be <= rn_amax … … 554 562 CALL tab_2d_3d( npti, nptidx(1:npti), a_ip_2d(1:npti,1:jpl), a_ip ) 555 563 CALL tab_2d_3d( npti, nptidx(1:npti), v_ip_2d(1:npti,1:jpl), v_ip ) 564 CALL tab_2d_3d( npti, nptidx(1:npti), v_il_2d(1:npti,1:jpl), v_il ) 556 565 CALL tab_2d_3d( npti, nptidx(1:npti), t_su_2d(1:npti,1:jpl), t_su ) 557 566 DO jl = 1, jpl … … 683 692 REAL(wp) :: zhmax, znum, zden, zalpha ! - - 684 693 ! 685 NAMELIST/namitd/ ln_cat_hfn, rn_himean, ln_cat_usr, rn_catbnd, rn_himin 694 NAMELIST/namitd/ ln_cat_hfn, rn_himean, ln_cat_usr, rn_catbnd, rn_himin, rn_himax 686 695 !!------------------------------------------------------------------ 687 696 ! … … 702 711 WRITE(numout,*) ' mean ice thickness in the domain rn_himean = ', rn_himean 703 712 WRITE(numout,*) ' Ice categories are defined by rn_catbnd ln_cat_usr = ', ln_cat_usr 704 WRITE(numout,*) ' minimum ice thickness rn_himin = ', rn_himin 713 WRITE(numout,*) ' minimum ice thickness allowed rn_himin = ', rn_himin 714 WRITE(numout,*) ' maximum ice thickness allowed rn_himax = ', rn_himax 705 715 ENDIF 706 716 ! … … 739 749 END DO 740 750 ! 741 hi_max(jpl) = 99._wp! set to a big value to ensure that all ice is thinner than hi_max(jpl)751 hi_max(jpl) = rn_himax ! set to a big value to ensure that all ice is thinner than hi_max(jpl) 742 752 ! 743 753 IF(lwp) WRITE(numout,*)
Note: See TracChangeset
for help on using the changeset viewer.