New URL for NEMO forge!   http://forge.nemo-ocean.eu

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.
#670 (LIM-3: small bug in latent heat of solid precipitation melting) – NEMO

Opened 14 years ago

Closed 14 years ago

Last modified 2 years ago

#670 closed Bug (fixed)

LIM-3: small bug in latent heat of solid precipitation melting

Reported by: gm Owned by: nemo
Priority: low Milestone:
Component: LIM3 Version: v3.2
Severity: Keywords: LIM* forcing ice seas snow v3.2
Cc:

Description

In LIM-3 a betas exponent has been introduced when computing the solid precipitation over the sea ice. Part of the snow that falls over sea-ice is put in the lead areas by the wind. Nevertheless, the latent heat flux due to snow melting is computed with sprecip *(1-at_i) whereas the actual solid precipitation over lead area is sprecip*(1-at_ibetas) (see limsbc.F90). If betas /=1, some latent heat is missing.

In limthd.F90 the following lines:

            zfontn         = sprecip(ji,jj) * lfus              ! energy of melting
            zfnsol         = qns(ji,jj)                         ! total non solar flux
            qldif(ji,jj)   = tms(ji,jj) * ( qsr(ji,jj)                               &
               &                               + zfnsol + fdtcn(ji,jj) - zfontn      &
               &                               + ( 1.0 - zindb ) * fsbbq(ji,jj)  )   &
               &                               * ( 1.0 - at_i(ji,jj) ) * rdt_ice    

should be replace by

            qldif(ji,jj) =  tms(ji,jj) * rdt_ice  &
               &   * (   ( 1.0 - at_i(ji,jj) ) * (   qsr(ji,jj)                          &   ! solar heat
               &                                   + qns(ji,jj)                          &   ! non solar heat
               &                                   + fdtcn(ji,jj)                        &   ! turbulent ice-ocean heat
               &                                   + ( 1.0 - zindb ) * fsbbq(ji,jj)  )   &   ! residual heat from previous step
               &       - sprecip(ji,jj) * lfus * pfrld(ji,jj)**betas                   ) &   ! latent heat of sprecip melting

and remove the declaration of zfontn and zfnsol, not used anymore

Gurvan

Commit History (2)

ChangesetAuthorTimeChangeLog
2387smasson2010-11-15T16:18:00+01:00

ticket #670 lim3: small bug in latent heat of solid precipitation melting

2383gm2010-11-13T16:01:26+01:00

v3.3beta: LIM-3 bug correction ticket #670

Change History (5)

comment:1 Changed 14 years ago by gm

done in v3.3 beta (see changeset 2383)

TO BE DONE in v3.2

comment:2 Changed 14 years ago by smasson

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

done in trunk of v3.2.1, see changeset #2387

comment:3 Changed 8 years ago by nicolasmartin

  • Keywords seas added; sea removed

comment:4 Changed 6 years ago by nemo

  • Keywords LIM* added

comment:5 Changed 2 years ago by nemo

  • Keywords v3.2 added
Note: See TracTickets for help on using tickets.