Changeset 13284 for NEMO/releases/r4.0/r4.0-HEAD/src/ICE/icedyn_rdgrft.F90
- Timestamp:
- 2020-07-09T17:12:23+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/releases/r4.0/r4.0-HEAD/src/ICE/icedyn_rdgrft.F90
r11732 r13284 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 … … 565 565 oirft2(ji) = oa_i_2d(ji,jl1) * afrft * hi_hrft 566 566 567 IF ( ln_pnd_ H12) THEN567 IF ( ln_pnd_LEV ) THEN 568 568 aprdg1 = a_ip_2d(ji,jl1) * afrdg 569 569 aprdg2(ji) = a_ip_2d(ji,jl1) * afrdg * hi_hrdg(ji,jl1) … … 572 572 aprft2(ji) = a_ip_2d(ji,jl1) * afrft * hi_hrft 573 573 vprft (ji) = v_ip_2d(ji,jl1) * afrft 574 IF ( ln_pnd_lids ) THEN 575 vlrdg (ji) = v_il_2d(ji,jl1) * afrdg 576 vlrft (ji) = v_il_2d(ji,jl1) * afrft 577 ENDIF 574 578 ENDIF 575 579 … … 598 602 sv_i_2d(ji,jl1) = sv_i_2d(ji,jl1) - sirdg1 - sirft(ji) 599 603 oa_i_2d(ji,jl1) = oa_i_2d(ji,jl1) - oirdg1 - oirft1 600 IF ( ln_pnd_ H12) THEN604 IF ( ln_pnd_LEV ) THEN 601 605 a_ip_2d(ji,jl1) = a_ip_2d(ji,jl1) - aprdg1 - aprft1 602 606 v_ip_2d(ji,jl1) = v_ip_2d(ji,jl1) - vprdg(ji) - vprft(ji) 607 IF ( ln_pnd_lids ) THEN 608 v_il_2d(ji,jl1) = v_il_2d(ji,jl1) - vlrdg(ji) - vlrft(ji) 609 ENDIF 603 610 ENDIF 604 611 ENDIF … … 692 699 v_s_2d (ji,jl2) = v_s_2d (ji,jl2) + ( vsrdg (ji) * rn_fsnwrdg * fvol(ji) + & 693 700 & vsrft (ji) * rn_fsnwrft * zswitch(ji) ) 694 IF ( ln_pnd_ H12) THEN701 IF ( ln_pnd_LEV ) THEN 695 702 v_ip_2d (ji,jl2) = v_ip_2d(ji,jl2) + ( vprdg (ji) * rn_fpndrdg * fvol (ji) & 696 703 & + vprft (ji) * rn_fpndrft * zswitch(ji) ) 697 704 a_ip_2d (ji,jl2) = a_ip_2d(ji,jl2) + ( aprdg2(ji) * rn_fpndrdg * farea & 698 705 & + aprft2(ji) * rn_fpndrft * zswitch(ji) ) 706 IF ( ln_pnd_lids ) THEN 707 v_il_2d (ji,jl2) = v_il_2d(ji,jl2) + ( vlrdg(ji) * rn_fpndrdg * fvol (ji) & 708 & + vlrft(ji) * rn_fpndrft * zswitch(ji) ) 709 ENDIF 699 710 ENDIF 700 711 … … 727 738 !---------------- 728 739 ! 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 )740 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 741 ! 731 742 END SUBROUTINE rdgrft_shift … … 839 850 CALL tab_3d_2d( npti, nptidx(1:npti), a_ip_2d(1:npti,1:jpl), a_ip(:,:,:) ) 840 851 CALL tab_3d_2d( npti, nptidx(1:npti), v_ip_2d(1:npti,1:jpl), v_ip(:,:,:) ) 852 CALL tab_3d_2d( npti, nptidx(1:npti), v_il_2d(1:npti,1:jpl), v_il(:,:,:) ) 841 853 DO jl = 1, jpl 842 854 DO jk = 1, nlay_s … … 865 877 CALL tab_2d_3d( npti, nptidx(1:npti), a_ip_2d(1:npti,1:jpl), a_ip(:,:,:) ) 866 878 CALL tab_2d_3d( npti, nptidx(1:npti), v_ip_2d(1:npti,1:jpl), v_ip(:,:,:) ) 879 CALL tab_2d_3d( npti, nptidx(1:npti), v_il_2d(1:npti,1:jpl), v_il(:,:,:) ) 867 880 DO jl = 1, jpl 868 881 DO jk = 1, nlay_s
Note: See TracChangeset
for help on using the changeset viewer.