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.
#318 (PISCES model : change the computation of xnegtr) – NEMO

Opened 15 years ago

Closed 15 years ago

#318 closed Enhancement (fixed)

PISCES model : change the computation of xnegtr

Reported by: cetlod Owned by: nemo
Priority: low Milestone:
Component: TOP Version: v3.0
Severity: Keywords:
Cc:

Description (last modified by cetlod)

xnegtr - array used to indicate negative values of tracer concentration - is computes such as his values must be between [0,1].


Older version:


xnegtr(:,:,:) = 1.e0
DO jn = jp_pcs0, jp_pcs1

DO jk = 1, jpk

DO jj = 1, jpj

DO ji = 1, jpi

IF( ( trn(ji,jj,jk,jn) + tra(ji,jj,jk,jn) ) < 0.0)

xnegtr(ji,jj,jk) = 0.e0

END DO

END DO

END DO

END DO


New version


xnegtr(:,:,:) = 1.e0
DO jn = jp_pcs0, jp_pcs1

DO jk = 1, jpk

DO jj = 1, jpj

DO ji = 1, jpi

IF( ( trn(ji,jj,jk,jn) + tra(ji,jj,jk,jn) ) < 0.e0 ) THEN

ztra = ABS( ( trn(ji,jj,jk,jn) - rtrn ) &

/ ( tra(ji,jj,jk,jn) + rtrn ) )

xnegtr(ji,jj,jk) = MIN( xnegtr(ji,jj,jk), ztra )

ENDIF

END DO

END DO

END DO

END DO


DO jn = jp_pcs0, jp_pcs1

trn(:,:,:,jn) = trn(:,:,:,jn) + xnegtr(:,:,:) * tra(:,:,:,jn)

END DO

Changes made by Olivier Aumont

Commit History (0)

(No commits)

Change History (1)

comment:1 Changed 15 years ago by cetlod

  • Description modified (diff)
  • Resolution set to fixed
  • Status changed from new to closed
Note: See TracTickets for help on using tickets.