- Timestamp:
- 2018-01-19T18:56:15+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icedyn_rdgrft.F90
r9169 r9271 570 570 REAL(wp) :: airdg1, oirdg1, aprdg1, virdg1, sirdg1 571 571 REAL(wp) :: airft1, oirft1, aprft1 572 REAL(wp), DIMENSION(jpij) :: airdg2, oirdg2, aprdg2, virdg2, sirdg2, vsrdg, vprdg , esrdg! area etc of new ridges573 REAL(wp), DIMENSION(jpij) :: airft2, oirft2, aprft2, virft , sirft , vsrft, vprft , esrft! area etc of rafted ice572 REAL(wp), DIMENSION(jpij) :: airdg2, oirdg2, aprdg2, virdg2, sirdg2, vsrdg, vprdg ! area etc of new ridges 573 REAL(wp), DIMENSION(jpij) :: airft2, oirft2, aprft2, virft , sirft , vsrft, vprft ! area etc of rafted ice 574 574 ! 575 575 REAL(wp), DIMENSION(jpij) :: ersw ! enth of water trapped into ridges … … 577 577 REAL(wp), DIMENSION(jpij) :: z1_ai ! 1 / a 578 578 ! 579 REAL(wp), DIMENSION(jpij,nlay_i) :: eirft ! ice energy of rafting ice 579 REAL(wp), DIMENSION(jpij,nlay_s) :: esrft ! snow energy of rafting ice 580 REAL(wp), DIMENSION(jpij,nlay_i) :: eirft ! ice energy of rafting ice 581 REAL(wp), DIMENSION(jpij,nlay_s) :: esrdg ! enth*volume of new ridges 580 582 REAL(wp), DIMENSION(jpij,nlay_i) :: eirdg ! enth*volume of new ridges 581 583 !!------------------------------------------------------------------- … … 633 635 oirdg1 = oa_i_2d(ji,jl1) * afrdg 634 636 oirdg2(ji) = oa_i_2d(ji,jl1) * afrdg * hi_hrdg(ji,jl1) 635 esrdg(ji) = ze_s_2d(ji,1,jl1) * afrdg636 637 637 638 virft(ji) = v_i_2d (ji,jl1) * afrft … … 640 641 oirft1 = oa_i_2d(ji,jl1) * afrft 641 642 oirft2(ji) = oa_i_2d(ji,jl1) * afrft * hi_hrft 642 esrft(ji) = ze_s_2d(ji,1,jl1) * afrft643 643 644 644 IF ( ln_pnd_H12 ) THEN … … 663 663 wfx_snw_dyn_1d(ji) = wfx_snw_dyn_1d(ji) + ( rhosn * vsrdg(ji) * ( 1._wp - rn_fsnwrdg ) & ! fresh water source for ocean 664 664 & + rhosn * vsrft(ji) * ( 1._wp - rn_fsnwrft ) ) * r1_rdtice 665 666 hfx_dyn_1d(ji) = hfx_dyn_1d(ji) + ( - esrdg(ji) * ( 1._wp - rn_fsnwrdg ) & ! heat sink for ocean (<0, W.m-2)667 & - esrft(ji) * ( 1._wp - rn_fsnwrft ) ) * r1_rdtice668 665 669 666 ! Put the melt pond water into the ocean … … 697 694 v_ip_2d(ji,jl1) = v_ip_2d(ji,jl1) - vprdg(ji) - vprft(ji) 698 695 ENDIF 699 ze_s_2d(ji,1,jl1) = ze_s_2d(ji,1,jl1) - esrdg (ji) - esrft (ji)700 696 ENDIF 701 697 702 698 END DO ! ji 703 699 700 ! special loop for e_s because of layers jk 701 DO jk = 1, nlay_s 702 DO ji = 1, npti 703 IF( apartf(ji,jl1) > 0._wp .AND. closing_gross(ji) > 0._wp ) THEN 704 ! Compute ridging /rafting fractions 705 afrdg = aridge(ji,jl1) * closing_gross(ji) * rdt_ice * z1_ai(ji) 706 afrft = araft (ji,jl1) * closing_gross(ji) * rdt_ice * z1_ai(ji) 707 ! Compute ridging ice and new ridges for es 708 esrdg(ji,jk) = ze_s_2d (ji,jk,jl1) * afrdg 709 esrft(ji,jk) = ze_s_2d (ji,jk,jl1) * afrft 710 ! Put the snow lost by ridging into the ocean 711 hfx_dyn_1d(ji) = hfx_dyn_1d(ji) + ( - esrdg(ji,jk) * ( 1._wp - rn_fsnwrdg ) & ! heat sink for ocean (<0, W.m-2) 712 & - esrft(ji,jk) * ( 1._wp - rn_fsnwrft ) ) * r1_rdtice 713 ! Update jl1 714 ze_s_2d(ji,jk,jl1) = ze_s_2d(ji,jk,jl1) * ( 1._wp - afrdg - afrft ) 715 ENDIF 716 END DO 717 END DO 718 704 719 ! special loop for e_i because of layers jk 705 720 DO jk = 1, nlay_i … … 754 769 sv_i_2d(ji,jl2) = sv_i_2d(ji,jl2) + ( sirdg2(ji) * fvol(ji) + sirft (ji) * zswitch(ji) ) 755 770 v_s_2d (ji,jl2) = v_s_2d (ji,jl2) + ( vsrdg (ji) * rn_fsnwrdg * fvol(ji) + & 756 & 771 & vsrft (ji) * rn_fsnwrft * zswitch(ji) ) 757 772 IF ( ln_pnd_H12 ) THEN 758 773 v_ip_2d (ji,jl2) = v_ip_2d(ji,jl2) + ( vprdg (ji) * rn_fpndrdg * fvol (ji) & … … 761 776 & + aprft2(ji) * rn_fpndrft * zswitch(ji) ) 762 777 ENDIF 763 ze_s_2d(ji,1,jl2) = ze_s_2d(ji,1,jl2) + ( esrdg (ji) * rn_fsnwrdg * fvol(ji) + &764 & esrft (ji) * rn_fsnwrft * zswitch(ji) )765 778 766 779 ENDIF 767 780 768 781 END DO 769 782 ! for snow enthalpy 783 DO jk = 1, nlay_s 784 DO ji = 1, npti 785 IF( apartf(ji,jl1) > 0._wp .AND. closing_gross(ji) > 0._wp ) & 786 & ze_s_2d(ji,jk,jl2) = ze_s_2d(ji,jk,jl2) + ( esrdg(ji,jk) * rn_fsnwrdg * fvol(ji) + & 787 & esrft(ji,jk) * rn_fsnwrft * zswitch(ji) ) 788 END DO 789 END DO 790 ! for ice enthalpy 770 791 DO jk = 1, nlay_i 771 792 DO ji = 1, npti
Note: See TracChangeset
for help on using the changeset viewer.