In adjoint-based application, I have made some correction of get_ngfldr.F, get_2dfldr.F, get_3dfldr.F in version 424.
Then it works well.
The correction is listed as follows:
IF (tdays(ng).lt.Tmax) THEN
Tmono=Tend
ELSE
! Tmono=Tend+(tdays(ng)-MOD(tdays(ng),Clength)) !original code
Tmono=Tend ...
Search found 2 matches
- Thu Feb 25, 2010 11:30 am
- Forum: ROMS Bugs
- Topic: Possible bug in get_2dfldr.F and get_3dfldr.F
- Replies: 3
- Views: 3159
- Thu Jan 14, 2010 12:43 pm
- Forum: ROMS Bugs
- Topic: Possible bugs in tl_uv3dmix2_geo.h
- Replies: 1
- Views: 2468
Possible bugs in tl_uv3dmix2_geo.h
I have been working on a strong constraint data assimilation test with version 424. If # define MIX_GEO_UV and # define UV_VIS2
are selected, the model may blow up suddenly with NaN error. The first NaN error came from the variable dVdz in tl_uv3dmix2_geo.h.
Furtherly, I found that the Loop index ...
are selected, the model may blow up suddenly with NaN error. The first NaN error came from the variable dVdz in tl_uv3dmix2_geo.h.
Furtherly, I found that the Loop index ...