Opened 3 years ago
Closed 3 years ago
#2746 closed Defect (fixed)
correction specific to BIO in 1D configuration: density has to be computed before the first time step of BIO
Reported by: | cperruche | Owned by: | systeam |
---|---|---|---|
Priority: | low | Milestone: | |
Component: | DOM | Version: | v4.0.* |
Severity: | minor | Keywords: | BIO, C1D, |
Cc: |
Description
#!correction specific to BIO in 1D configuration
Context
the density is not computed before the first time step of the BIO, but some BIO variables needs the density, so the code uses rhop=0.00.
This impacts DIC during several months.
Analysis
In stp_c1d, trc_stp needs the density, but trc_stp is called before before eos. So rhop=0 is used at the first time step.
Recommendation
I suggest to add a call to eos in istate_init (istate.F90), at the end of the loop:
IF( ln_rstart ) THEN ! Restart from a file
…
ELSE ! Start from rest
…
CALL eos( tsn, rhd, rhop, gdept_0(:,:,:) ) ! now potential density for zdfmxl
ENDIF
Commit History (1)
Changeset | Author | Time | ChangeLog |
---|---|---|---|
15594 | cetlod | 2021-12-13T12:05:33+01:00 | v40-HEAD : bugfix in 1D configuration, see ticket #2746 |
Change History (2)
comment:1 Changed 3 years ago by cetlod
comment:2 Changed 3 years ago by cetlod
- Resolution set to fixed
- Status changed from new to closed
In 15594: