#1549 closed Bug (fixed)
qsr_oce not initialised in lim_thd
| Reported by: | ufla | Owned by: | nemo |
|---|---|---|---|
| Priority: | high | Milestone: | |
| Component: | LIM3 | Version: | trunk |
| Severity: | Keywords: | LIM* | |
| Cc: |
Description
Sorry if this is the same all over again, but here is another one from limthd.F90:
! --- Energy received in the lead, zqld is defined everywhere (J.m-2) --- !
zqld = tmask(ji,jj,1) * rdt_ice * &
& ( pfrld(ji,jj) * qsr_oce(ji,jj) * frq_m(ji,jj) + pfrld(ji,jj) * qns_oce(ji,jj) + qemp_oce(ji,jj) )
qsr_oce is not initialised, at least not in the first execution of the routine. Leads to floating point exception and crash.
Coupled configuration and confirmed by initialising qsr_oce to HUGE() in sbc_ice.F90.
Commit History (0)
(No commits)
Change History (2)
comment:1 Changed 9 years ago by clem
- Priority changed from major to critical
- Resolution set to fixed
- Status changed from new to closed
comment:2 Changed 6 years ago by nemo
- Keywords LIM* added
Note: See
TracTickets for help on using
tickets.
This is definitely not the same all over again. You point out a major bug in coupled mode with LIM3. From what I see, qsr_oce is only defined in forced mode and therefore will always be 0 while coupling with the atmosphere! Whoops. I just corrected it at revision 5486.
For coupled simulations, revision 5486 is absolutely needed. It seems to me that all previous simulations can unfortunately be trashed.
Thanks again Uwe. You are being very helpful.