New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
ticket/1613_NOC3_vvl_by_default (diff) – NEMO

Changes between Version 17 and Version 18 of ticket/1613_NOC3_vvl_by_default


Ignore:
Timestamp:
2015-10-31T09:11:33+01:00 (8 years ago)
Author:
gm
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • ticket/1613_NOC3_vvl_by_default

    v17 v18  
    2525 
    2626== 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] 
    2828 
    2929'''   [[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 !!!! [[BR]]  '''  The error is thus probably in dynldf_lap_blp.F90[[BR]]      [[BR]]      other test to do: try no non-linear term ;  vvl off ; filtered free surf ; mld already off (!) ; 
     
    3939 
    4040---- 
    41 == Misc. == 
     41== Misc. issues also solved == 
    4242 * domhpg.F90  : replace everywhere  /e1u (/e2v) by *r1_e1u (*r1_e2v)   ==>> put in comment !!gm to get same results 
    4343 * use nn_print test to bypass some control print. (domcfg.F90 ; ...)