- Timestamp:
- 2020-02-12T17:57:19+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/iceitd.F90
r11081 r12379 409 409 CALL tab_3d_2d( npti, nptidx(1:npti), a_ip_2d(1:npti,1:jpl), a_ip ) 410 410 CALL tab_3d_2d( npti, nptidx(1:npti), v_ip_2d(1:npti,1:jpl), v_ip ) 411 CALL tab_3d_2d( npti, nptidx(1:npti), lh_ip_2d(1:npti,1:jpl), lh_ip ) 411 412 CALL tab_3d_2d( npti, nptidx(1:npti), t_su_2d(1:npti,1:jpl), t_su ) 412 413 DO jl = 1, jpl … … 481 482 v_ip_2d(ji,jl1) = v_ip_2d(ji,jl1) - ztrans 482 483 v_ip_2d(ji,jl2) = v_ip_2d(ji,jl2) + ztrans 484 ! 485 ztrans = lh_ip_2d(ji,jl1) * zworka(ji) ! Pond lid thickness 486 lh_ip_2d(ji,jl1) = lh_ip_2d(ji,jl1) - ztrans 487 lh_ip_2d(ji,jl2) = lh_ip_2d(ji,jl2) + ztrans 483 488 ENDIF 484 489 ! … … 525 530 ! clem: The transfer between one category to another can lead to very small negative values (-1.e-20) 526 531 ! because of truncation error ( i.e. 1. - 1. /= 0 ) 527 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 )532 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, lh_ip_2d, ze_s_2d, ze_i_2d ) 528 533 529 534 ! at_i must be <= rn_amax … … 553 558 CALL tab_2d_3d( npti, nptidx(1:npti), a_ip_2d(1:npti,1:jpl), a_ip ) 554 559 CALL tab_2d_3d( npti, nptidx(1:npti), v_ip_2d(1:npti,1:jpl), v_ip ) 560 CALL tab_2d_3d( npti, nptidx(1:npti), lh_ip_2d(1:npti,1:jpl), lh_ip ) 555 561 CALL tab_2d_3d( npti, nptidx(1:npti), t_su_2d(1:npti,1:jpl), t_su ) 556 562 DO jl = 1, jpl
Note: See TracChangeset
for help on using the changeset viewer.