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.
#1606 (Bug in analytical initial temperature profile) – NEMO

Opened 9 years ago

Closed 7 years ago

#1606 closed Bug (wontfix)

Bug in analytical initial temperature profile

Reported by: hadcv Owned by: flavoni
Priority: low Milestone:
Component: OCE Version: v3.6
Severity: Keywords: DOM
Cc: flavoni

Description (last modified by nicolasmartin)

The analytical temperature profile in istate_t_s seems to be incorrect:

      DO jk = 1, jpk
         tsn(:,:,jk,jp_tem) = (  ( ( 7.5 - 0. * ABS( gphit(:,:) )/30. ) * ( 1.-TANH((fsdept(:,:,jk)-80.)/30.) )   &
            &                + 10. * ( 5000. - fsdept(:,:,jk) ) /5000.)  ) * tmask(:,:,jk)
         tsb(:,:,jk,jp_tem) = tsn(:,:,jk,jp_tem)
      END DO

The latitude in the tsn equation is multiplied by 0; the profile will have a ~25 degree SST everywhere.

Commit History (0)

(No commits)

Change History (5)

comment:1 Changed 8 years ago by nicolasmartin

  • Keywords istate removed

comment:2 Changed 7 years ago by clevy

  • Owner changed from nemo to flavoni

comment:3 Changed 7 years ago by clevy

  • Cc flavoni added
  • Status changed from new to assigned

comment:4 Changed 7 years ago by nicolasmartin

  • Description modified (diff)

comment:5 Changed 7 years ago by flavoni

  • Resolution set to wontfix
  • Status changed from assigned to closed

This is not really a bug; in this particular case it was not used, because it is multiplied by 0, but it is here just like an example, to show how it can be used an initial temperature profile that varies with latitude.

Note: See TracTickets for help on using tickets.