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.
#1361 (initialization of TKE when missing in the restart) – NEMO

Opened 10 years ago

Closed 9 years ago

#1361 closed Bug (fixed)

initialization of TKE when missing in the restart

Reported by: avidard Owned by: gm
Priority: low Milestone:
Component: OCE Version: v3.4
Severity: Keywords:
Cc:

Description

from 3.4 onward (probably before)

In the unlikely event you want to use tke while it is not in the restart, the code will do a spinup like

ELSE                                     ! No TKE array found: initialisation

IF(lwp) WRITE(numout,*) ' ===>>>> : previous run without tke scheme, en computed by iterative loop'
en (:,:,:) = rn_emin * tmask(:,:,:)
CALL tke_avn                               ! recompute avt, avm, avmu, avmv and dissl (approximation)
DO jit = nit000 + 1, nit000 + 10   ;   CALL zdf_tke( jit )   ;   END DO

ENDIF

however zdf_tke starts by

IF( kt /= nit000 ) THEN   ! restore before value to compute tke

avt (:,:,:) = avt_k (:,:,:)
avm (:,:,:) = avm_k (:,:,:)
avmu(:,:,:) = avmu_k(:,:,:)
avmv(:,:,:) = avmv_k(:,:,:)

ENDIF

by starting the loop of the spinup at nit000+1 the av*_k are not initialized

Commit History (2)

ChangesetAuthorTimeChangeLog
5113cbricaud2015-03-02T17:04:43+01:00

bugfix for ticket #1361

5112cbricaud2015-03-02T17:03:24+01:00

bugfix for ticket #1361

Change History (2)

comment:1 Changed 10 years ago by gm

  • Owner changed from NEMO team to gm

comment:2 Changed 9 years ago by cbricaud

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

avt_k, avm_k ,avmu_k, avmv_k need to be initialised by avt, avm,avmu, avmv before call to tke_avn.

bugfix done for zdftke and zdfgls in NEMO_3.4_1:
https://forge.ipsl.jussieu.fr/nemo/changeset/5113/branches/2012/dev_v3_4_STABLE_2012

and zdftke for NEMO_3.6_beta ( zdfgls was ok) :
https://forge.ipsl.jussieu.fr/nemo/changeset/5112/trunk

Note: See TracTickets for help on using tickets.