- Timestamp:
- 2017-06-06T09:39:43+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90
r8125 r8142 509 509 ! MV MP 2016 510 510 REAL(wp), POINTER, DIMENSION(:) :: vprdg ! pond volume of ridging ice 511 REAL(wp), POINTER, DIMENSION(:) :: aprdg1 ! pond area of ridging ice 512 REAL(wp), POINTER, DIMENSION(:) :: aprdg2 ! pond area of ridging ice 511 513 ! END MV MP 2016 512 514 REAL(wp), POINTER, DIMENSION(:) :: dhr , dhr2 ! hrmax - hrmin & hrmax^2 - hrmin^2 … … 525 527 ! MV MP 2016 526 528 REAL(wp), POINTER, DIMENSION(:) :: vprft ! pond volume of rafting ice 529 REAL(wp), POINTER, DIMENSION(:) :: aprft1 ! pond area of rafted ice 530 REAL(wp), POINTER, DIMENSION(:) :: aprft2 ! pond area of new rafted ice 527 531 ! END MV MP 2016 528 532 REAL(wp), POINTER, DIMENSION(:) :: esrft , smrft ! snow energy & salinity of rafting ice … … 539 543 ! MV MP 2016 540 544 !CALL wrk_alloc( jpij, afrac, ardg1, ardg2, vsrdg, esrdg, dhr, dhr2 ) 541 CALL wrk_alloc( jpij, afrac, ardg1, ardg2, vsrdg, esrdg, vprdg, dhr, dhr2 )545 CALL wrk_alloc( jpij, afrac, ardg1, ardg2, vsrdg, esrdg, vprdg, aprdg1, aprdg2, dhr, dhr2 ) 542 546 ! END MV MP 2016 543 547 CALL wrk_alloc( jpij, vrdg1, vrdg2, vsw , srdg1, srdg2, smsw, oirdg1, oirdg2 ) 544 548 ! MV MP 2016 545 549 !CALL wrk_alloc( jpij, afrft, arft1, arft2, virft, vsrft, esrft, smrft, oirft1, oirft2 ) 546 CALL wrk_alloc( jpij, afrft, arft1, arft2, virft, vsrft, esrft, vprft, smrft, oirft1, oirft2 ) 550 CALL wrk_alloc( jpij, afrft, arft1, arft2, virft, vsrft, esrft, aprft1, aprft2) 551 CALL wrk_alloc ( jpij, vprft, smrft, oirft1, oirft2 ) 547 552 ! END MV MP 2016 548 553 CALL wrk_alloc( jpij,nlay_i, eirft, erdg1, erdg2, ersw ) … … 598 603 599 604 !-------------------------------------------------------------------------- 600 ! 3.4) Sub tract area, volume, and energy from ridging605 ! 3.4) Substract area, volume, and energy from ridging 601 606 ! / rafting category n1. 602 607 !-------------------------------------------------------------------------- … … 609 614 !MV MP 2016 610 615 IF ( nn_pnd_scheme > 0 ) THEN 611 vprdg(ij) = v_ip(ji,jj, jl1) * afrac(ij) 616 aprdg1(ij) = a_ip(ji,jj, jl1) * afrac(ij) 617 aprdg2(ij) = a_ip(ji,jj, jl1) * afrac(ij) * krdg(ji,jj,jl1) 618 vprdg(ij) = v_ip(ji,jj, jl1) * afrac(ij) 612 619 ENDIF 613 620 ! END MV MP 2016 … … 621 628 !MV MP 2016 622 629 IF ( nn_pnd_scheme > 0 ) THEN 623 vprft(ij) = v_ip(ji,jj,jl1) * afrft(ij) 630 aprft1(ij) = a_ip (ji,jj, jl1) * afrft(ij) 631 aprft2(ij) = a_ip (ji,jj, jl1) * afrft(ij) * kraft 632 vprft(ij) = v_ip(ji,jj,jl1) * afrft(ij) 624 633 ENDIF 625 634 ! END MV MP 2016 … … 690 699 IF ( nn_pnd_scheme > 0 ) THEN 691 700 v_ip (ji,jj,jl1) = v_ip (ji,jj,jl1) - vprdg (ij) - vprft (ij) 701 a_ip (ji,jj,jl1) = a_ip (ji,jj,jl1) - aprdg1(ij) - aprft1(ij) 692 702 ENDIF 693 703 ! END MV MP 2016 … … 758 768 ! MV MP 2016 759 769 IF ( nn_pnd_scheme > 0 ) THEN 760 v_ip (ji,jj,jl2) = v_ip (ji,jj,jl2) + ( vprdg (ij) * rn_fpondrdg * fvol(ij) + & 761 & vprft (ij) * rn_fpondrft * zswitch(ij) ) 770 v_ip (ji,jj,jl2) = v_ip (ji,jj,jl2) + ( vprdg (ij) * rn_fpondrdg * fvol(ij) + & 771 & vprft (ij) * rn_fpondrft * zswitch(ij) ) 772 a_ip (ji,jj,jl2) = a_ip(ji,jj,jl2) + ( aprdg2(ij) * rn_fpondrdg * farea + & 773 & aprft2(ij) * rn_fpondrft * zswitch(ji) ) 762 774 ENDIF 763 775 ! END MV MP 2016 … … 782 794 ! MV MP 2016 783 795 !CALL wrk_dealloc( jpij, afrac, ardg1, ardg2, vsrdg, esrdg, dhr, dhr2 ) 784 CALL wrk_dealloc( jpij, afrac, ardg1, ardg2, vsrdg, esrdg, vprdg, dhr, dhr2 )796 CALL wrk_dealloc( jpij, afrac, ardg1, ardg2, vsrdg, esrdg, vprdg, aprdg1, aprdg2, dhr, dhr2 ) 785 797 ! END MV MP 2016 786 798 CALL wrk_dealloc( jpij, vrdg1, vrdg2, vsw , srdg1, srdg2, smsw, oirdg1, oirdg2 ) 787 799 ! MV MP 2016 788 800 !CALL wrk_dealloc( jpij, afrft, arft1, arft2, virft, vsrft, esrft, smrft, oirft1, oirft2 ) 789 CALL wrk_dealloc( jpij, afrft, arft1, arft2, virft, vsrft, esrft, vprft, smrft, oirft1, oirft2 ) 801 CALL wrk_dealloc( jpij, afrft, arft1, arft2, virft, vsrft, esrft, aprft1, aprft2, vprft ) 802 CALL wrk_dealloc( jpij, smrft, oirft1, oirft2 ) 790 803 CALL wrk_dealloc( jpij,nlay_i, eirft, erdg1, erdg2, ersw ) 791 804 !
Note: See TracChangeset
for help on using the changeset viewer.