@@ -507,8 +507,7 @@ end subroutine atmo_boundary_const
507
507
!
508
508
! changes: Andrew Roberts, NPS (RASM/CESM coupling and documentation)
509
509
510
- subroutine neutral_drag_coeffs (apnd , hpnd , &
511
- ipnd , &
510
+ subroutine neutral_drag_coeffs (apondn , &
512
511
alvl , vlvl , &
513
512
aice , vice , &
514
513
vsno , aicen , &
@@ -526,9 +525,7 @@ subroutine neutral_drag_coeffs (apnd, hpnd, &
526
525
use icepack_tracers, only: tr_pond
527
526
528
527
real (kind= dbl_kind), dimension (:), intent (in ) :: &
529
- apnd ,& ! melt pond fraction of sea ice
530
- hpnd ,& ! mean melt pond depth over sea ice
531
- ipnd ,& ! mean ice pond depth over sea ice in cat n
528
+ apondn ,& ! melt pond fraction of sea ice category
532
529
alvl ,& ! level ice area fraction (of grid cell ?)
533
530
vlvl ! level ice mean thickness
534
531
@@ -666,7 +663,7 @@ subroutine neutral_drag_coeffs (apnd, hpnd, &
666
663
if (tr_pond) then
667
664
do n = 1 ,ncat
668
665
! area of pond per unit area of grid cell
669
- apond = apond+ apnd (n)* aicen(n)
666
+ apond = apond+ apondn (n)* aicen(n)
670
667
enddo
671
668
endif
672
669
0 commit comments