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