Skip to content

Home

Jump to bottom
jdha edited this page Jan 11, 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.

    • ln_full_vel was switched off as it was causing the SZT model to blow up on the boundary. The purpose of this switch is to decompose the BT and BC velocities from the full 3D BDY velocity file. As a Neumann condition was being used for the BC velocities I've just switch ln_full_vel = .false. so only the 2D barotropic BDY files are read in.

  • 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 and also clarify the method for the spin-up of the ice

Clone this wiki locally