Skip to content

Home

Jump to bottom
Adam Blaker edited this page Jan 7, 2021 · 26 revisions

ORCHESTRA

Testing the SZT configuration

In my stripped down version it all runs fine. Copying and modifying Dave's ZPS configuration, things start to break down (always at ts=25). A quick difference of my namelist_cfg with Dave's shows a few candidates. A summary of a few sensitivity tests wrt the namelist_cfg:

  • nn_tra_dta = 0 (fails ts=25; original value 1)
  • ln_isf = false (fails ts=25; original value true)
  • rn_rdt = 240 (fails ts=31; original value 300)
  • nn_dyn3d_dta = 0 (fails ts=25; original value 1)
  • ln_vol = false (fails ts=25; original value true)
  • ln_bdy = false (fails ts=25; original value true)

I've also switched off runoff and geothermal heating switch, but it makes no difference either.

Switching out the initial conditions [which incidentally had different names, but appear to be the same data], with no success, I turned my attention to the cpp keys. By eliminating keys one by one I've found that key_trabbl appears to be the one that is causing the SZT configuration to crash.

Solution: you can still compile with key_trabbl, just set the namelist parameters accordingly (i.e. for SZT set bbl diffusion and advection to 0. This way the same exe can be used for ZPS and SZT configurations). i.e.

!-----------------------------------------------------------------------
&nambbl        !   bottom boundary layer scheme
!-----------------------------------------------------------------------
   nn_bbl_ldf  =  0      !  diffusive bbl (=1)   or not (=0)
   nn_bbl_adv  =  0      !  advective bbl (=1/2) or not (=0)

Why should trabbl be an issue in the SZT configuration?

Queries SZT vs ZPS namelist choices

Reading through the namelists the following choices warrant further investigation/understanding:

  • ln_dm2dc = .false. in DM's JRA forced simulations. Was this the case in the CORE2 ones?

  • ln_full_vel = .true. in DM's simulations. Presumably this is the full velocities being supplied at the northern boundaries? It is currently set to .false. in this repository, but .true. would seem a better choice unless there are sensible reasons otherwise.

  • nn_bbl_ldf and nn_bbl_adv < see above, N006 used ...ldf = 1 and ...adv = 0

  • ln_traldf_iso switched to ln_traldf_hor |
    ln_dynldf_hor switched to ln_dynldf_lev | Are these required changes for the SZT?

  • rn_aht_0 = 1500. ? This seems an extremely high value, given that eORCA1 simulations typically use 1000. The N006 global ORCA12 simulation used rn_aht_0 = 125.

  • Check with DM what the benefits of the MUSCL scheme are - should this be switched for SZT?

Also check with DM about how to handle leap years in the simulations so that e.g. the 5-day means remain in sync.

Clone this wiki locally