New URL for NEMO forge!

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
#1637 (Mismatching units in new contributions to qns in sbcblk_core.F90) – NEMO

Opened 8 years ago

Closed 8 years ago

#1637 closed Bug (invalid)

Mismatching units in new contributions to qns in sbcblk_core.F90

Reported by: royf3 Owned by: nemo
Priority: low Milestone:
Component: OCE Version: v3.6
Severity: Keywords:


In sbcblk_core.F90

I think there are errors in these new contributions to qns (units are not giving W/m2). Also, should temperatures in Kelvin be used when multiplying by cpic?
The other question I have is for CICE interface users, maybe some of these contributions are already dealt with in sbcice_cice.F90.

qns(:,:) = ...

& + ( sf(jp_prec)%fnow(:,:,1) - sf(jp_snow)%fnow(:,:,1) ) * rn_pfac & ! add liquid precip heat content at Tair
& * ( sf(jp_tair)%fnow(:,:,1) - rt0 ) * rcp &
& + sf(jp_snow)%fnow(:,:,1) * rn_pfac & ! add solid precip heat content at min(Tair,Tsnow)
& * ( MIN( sf(jp_tair)%fnow(:,:,1), rt0_snow ) - rt0 ) * cpic * tmask(:,:,1)


Commit History (0)

(No commits)

Change History (2)

comment:1 Changed 8 years ago by royf3

In phycst.F90

units of rcp should be J/Kelvin/kg, use of temperature in Celsius ok as the advective units of T is Celsius. So finally I think the code is ok if rcp is indeed in J/Kelvin/kg


comment:2 Changed 8 years ago by gm

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

Historically ocean and ice models always use Celsius , not Kelvin. All associated thermodynamical quantities and the equation of state are set assuming Celsius.  

===>>> no bug there.

sbcblk_core compute the fluxes between the atoms and the ocean, and the atoms and the ice. 

This does not depend on what is done in sbcice_cice or sbcice_lim(_2) (depending of sea-ice model used) which compute the fluxes between ice and ocean. The heat content of mass exchanged between ice and ocean are taken into account in sbcice_cice or sbcice_lim and this is properly done.

Note: See TracTickets for help on using tickets.