- Timestamp:
- 2020-04-08T18:54:44+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icedyn_rdgrft.F90
r11732 r12720 494 494 REAL(wp) :: airdg1, oirdg1, aprdg1, virdg1, sirdg1 495 495 REAL(wp) :: airft1, oirft1, aprft1 496 REAL(wp), DIMENSION(jpij) :: airdg2, oirdg2, aprdg2, virdg2, sirdg2, vsrdg, vprdg ! area etc of new ridges497 REAL(wp), DIMENSION(jpij) :: airft2, oirft2, aprft2, virft , sirft , vsrft, vprft ! area etc of rafted ice496 REAL(wp), DIMENSION(jpij) :: airdg2, oirdg2, aprdg2, virdg2, sirdg2, vsrdg, vprdg, vlrdg ! area etc of new ridges 497 REAL(wp), DIMENSION(jpij) :: airft2, oirft2, aprft2, virft , sirft , vsrft, vprft, vlrft ! area etc of rafted ice 498 498 ! 499 499 REAL(wp), DIMENSION(jpij) :: ersw ! enth of water trapped into ridges … … 569 569 aprdg2(ji) = a_ip_2d(ji,jl1) * afrdg * hi_hrdg(ji,jl1) 570 570 vprdg (ji) = v_ip_2d(ji,jl1) * afrdg 571 vlrdg (ji) = v_il_2d(ji,jl1) * afrdg 571 572 aprft1 = a_ip_2d(ji,jl1) * afrft 572 573 aprft2(ji) = a_ip_2d(ji,jl1) * afrft * hi_hrft 573 574 vprft (ji) = v_ip_2d(ji,jl1) * afrft 575 vlrft (ji) = v_il_2d(ji,jl1) * afrft 574 576 ENDIF 575 577 … … 601 603 a_ip_2d(ji,jl1) = a_ip_2d(ji,jl1) - aprdg1 - aprft1 602 604 v_ip_2d(ji,jl1) = v_ip_2d(ji,jl1) - vprdg(ji) - vprft(ji) 605 v_il_2d(ji,jl1) = v_il_2d(ji,jl1) - vlrdg(ji) - vlrft(ji) 603 606 ENDIF 604 607 ENDIF … … 697 700 a_ip_2d (ji,jl2) = a_ip_2d(ji,jl2) + ( aprdg2(ji) * rn_fpndrdg * farea & 698 701 & + aprft2(ji) * rn_fpndrft * zswitch(ji) ) 702 v_il_2d (ji,jl2) = v_il_2d(ji,jl2) + ( vlrdg (ji) * rn_fpndrdg * fvol (ji) & 703 & + vlrft (ji) * rn_fpndrft * zswitch(ji) ) 699 704 ENDIF 700 705 … … 727 732 !---------------- 728 733 ! In case ridging/rafting lead to very small negative values (sometimes it happens) 729 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 )734 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 ) 730 735 ! 731 736 END SUBROUTINE rdgrft_shift … … 839 844 CALL tab_3d_2d( npti, nptidx(1:npti), a_ip_2d(1:npti,1:jpl), a_ip(:,:,:) ) 840 845 CALL tab_3d_2d( npti, nptidx(1:npti), v_ip_2d(1:npti,1:jpl), v_ip(:,:,:) ) 846 CALL tab_3d_2d( npti, nptidx(1:npti), v_il_2d(1:npti,1:jpl), v_il(:,:,:) ) 841 847 DO jl = 1, jpl 842 848 DO jk = 1, nlay_s … … 865 871 CALL tab_2d_3d( npti, nptidx(1:npti), a_ip_2d(1:npti,1:jpl), a_ip(:,:,:) ) 866 872 CALL tab_2d_3d( npti, nptidx(1:npti), v_ip_2d(1:npti,1:jpl), v_ip(:,:,:) ) 873 CALL tab_2d_3d( npti, nptidx(1:npti), v_il_2d(1:npti,1:jpl), v_il(:,:,:) ) 867 874 DO jl = 1, jpl 868 875 DO jk = 1, nlay_s
Note: See TracChangeset
for help on using the changeset viewer.