Changes between Version 21 and Version 22 of ticket/1613_NOC3_vvl_by_default
- Timestamp:
- 2015-11-03T13:38:12+01:00 (9 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
ticket/1613_NOC3_vvl_by_default
v21 v22 25 25 26 26 == I. Before/Now/After vertical coordinate system by default == 27 ===>>> runs ORCA2_LIM3 trunk and this branch identical to the last digit ===>>> ''' OK''' see revision [https://forge.ipsl.jussieu.fr/nemo/log/branches/2015/dev_r5836_NOC3_vvl_by_default?rev= 5845]27 ===>>> runs ORCA2_LIM3 trunk and this branch identical to the last digit ===>>> ''' OK''' see revision [https://forge.ipsl.jussieu.fr/nemo/log/branches/2015/dev_r5836_NOC3_vvl_by_default?rev= 5845] [[BR]]===>>> small bug correction in in_out_manager.F90 ===>>> ''' OK''' see revision [https://forge.ipsl.jussieu.fr/nemo/log/branches/2015/dev_r5836_NOC3_vvl_by_default?rev= 5858] 28 28 29 29 ''' [[BR]] Nevertheless''' for both the trunk & this branch, ORCA2_LIM3 blow up atter 138 time step....[[BR]] NB: to obtain the same results 1 things : revert domhpg.F90 optimization (use of r1_e...) [[BR]] =========>>>>>>>>>>>> try to find the origine of the model explosion:[[BR]] (1) dynvor.F90 reverted to the v3.6_stable version, i.e. with use of fmask at the proper place, the trunk still explode round the same time-step... So the problem does not come from dynvor.[[BR]] (2) iso neutral diffusion switch off (as well as div): it is not the issue the model still blow up (pdt 139)[[BR]] (3) partial cells switch off (ln_zco=T): it is not the issue the model still blow up (pdt 142)[[BR]] (4) tra_adv : use 2nd order centered scheme: blow up pdt 134[[BR]] (5) just set hat as a constant= 2000 (instead of grid depend coeff nn_aht_ijk = 20) : blow up pdt 134[[BR]] (6) dynldf off (ln_dynldf_lap=F) & use of UBS (ln_dynadv_vec= F ; ln_dynad_ubs=T): ''' it is working !!!! '''(I stop the run after 3531)[[BR]]The error is thus probably in dynldf_lap_blp.F90. Note that if the code does not blow up it is very noisy at mid latitudes with such a setting... (probably logic with UBS on such a low resolution...). … … 35 35 36 36 == [[BR]]II. Non-linear free surface by default == 37 * replace lk_vvl by .NOT.ln_linssh everywhere 37 * replace lk_vvl by .NOT.ln_linssh everywhere 38 38 * remove key_vvl and introduce ln_linssh in nam_zgr namelist 39 39