Changes between Version 23 and Version 24 of ticket/1613_NOC3_vvl_by_default

Timestamp:

2015-11-05T12:53:58+01:00 (9 years ago)

Author:

gm

Comment:

--

ticket/1613_NOC3_vvl_by_default

@@ -25 +25 @@
 
 == I. Before/Now/After vertical coordinate system by default ==
-  ===>>> 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][[BR]]===>>> Not working in non vvl case (blow up in 10 time step !)  ===>>> 
+  ===>>> 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] [[BR]] ===>>> Not working in non vvl case (blow up in 10 time step !)  ===>>>
 
 '''   [[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...).
@@ -40 +40 @@
 ----
 == Misc. issues also solved ==
+ * dynldf_lap_blp.F90 & sbcwaves.F90: correct the 2 bugs found in trunk (see revision [https://forge.ipsl.jussieu.fr/nemo/log/trunk/NEMOGCM?rev= ​5860])
+ * dynldf_lap_blp.F90 : add the correction on zdiv in AGRIF case
  * domhpg.F90  : replace everywhere  /e1u (/e2v) by *r1_e1u (*r1_e2v)   ==>> put in comment !!gm to get same results
- * dynspg_flt.F90 : replace everywhere /e1u (/e2v) by *r1_e1u (*r1_e2v)
+ * dynspg_flt.F90 : replace everywhere /e1u (/e2v) by *r1_e1u (*r1_e2v)
  * use nn_print test to bypass some control print. (domcfg.F90 ; ...)
  * remove the use of prihre, prihin, prizre  ==>>> remove lib_print.f90