#2174 closed Enhancement (fixed)

Update iceload calculation in dyn_hpg_init

Reported by: jamesharle Owned by: systeam
Priority: low Milestone:
Component: DYN Version: trunk
Severity: minor Keywords: ISF, HPG
Cc: Review:
MP ready?:
Progress:

Description

Context

In dyn_hpg_init the iceload calculation uses gdept_1d. For sigma style vertical coordinates this needs to be replaced by gdept_0.

Proposal

  • dynhpg.F90

     
    244244                     &                              * (1._wp - tmask(ji,jj,jk)) 
    245245               END DO 
    246246               IF (ikt  >=  2) ziceload(ji,jj) = ziceload(ji,jj) + (2._wp * znad + zrhdtop_isf(ji,jj) + zrhd(ji,jj,ikt-1)) & 
    247                   &                                              * ( risfdep(ji,jj) - gdept_1d(ikt-1) ) 
     247                  &                                              * ( risfdep(ji,jj) - gdept_0(ji,jj,ikt-1) ) 
    248248            END DO 
    249249         END DO 
    250250         riceload(:,:) = ziceload(:,:)  ! need to be saved for diaar5 

Commit History (0)

(No commits)

Change History (1)

comment:1 Changed 15 months ago by mathiot

  • Resolution set to fixed
  • Status changed from new to closed

Fix in r10491

Note: See TracTickets for help on using tickets.