Opened 5 years ago
Closed 4 years ago
#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: |
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
244 244 & * (1._wp - tmask(ji,jj,jk)) 245 245 END DO 246 246 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) ) 248 248 END DO 249 249 END DO 250 250 riceload(:,:) = ziceload(:,:) ! need to be saved for diaar5
Commit History (0)
(No commits)
Change History (1)
comment:1 Changed 4 years ago by mathiot
- Resolution set to fixed
- Status changed from new to closed
Note: See
TracTickets for help on using
tickets.
Fix in r10491