attempt to diagnose early blow up with Ash namelist

EXP_Apr19/file_def_nemo.xml

@@ -12,10 +12,12 @@
     
       <file_group id="1ts" output_freq="1ts"  output_level="10" enabled=".TRUE."/> <!-- 1 time step files -->
 
-      <file_group id="1h" output_freq="1h"  output_level="10" enabled=".TRUE."> <!-- 1h files -->
+      <file_group id="1ts" output_freq="1ts"  output_level="10" enabled=".TRUE."> <!-- 1h files -->
   	<file id="file19" name_suffix="_SSH" description="ocean T grid variables" >
     	 <field field_ref="ssh"          name="zos"   />
          <field field_ref="sst"          name="tos"   />
+         <field field_ref="ssu"          name="uos"  grid_ref="grid_U_2D" />
+         <field field_ref="ssv"          name="vos"  grid_ref="grid_V_2D" />
   	</file>
       </file_group>
 

EXP_Apr19/namelist_cfg

@@ -6,7 +6,7 @@
 !-----------------------------------------------------------------------
    cn_exp      =  "BoBEAS"  !  experience name
    nn_it000    = 1 ! 57601  !  first time step ! CHNAGE THIS TO END OF LAST RUN + 1
-   nn_itend    = 14400 ! 218880 crashed in Apr16 ! 138240 !  last  time step (std 1 day = 144) (for dt = 6 min, so 240*dt = 1 day) ! CHANGE THIS TO THIS + 56 DAYS
+   nn_itend    = 100! 218880 crashed in Apr16 ! 138240 !  last  time step (std 1 day = 144) (for dt = 6 min, so 240*dt = 1 day) ! CHANGE THIS TO THIS + 56 DAYS
    nn_date0    =  20190401  !20110717 !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
    nn_time0    =       0   !  initial time of day in hhmm
    nn_leapy    =       1   !  Leap year calendar (1) or not (0)
@@ -23,8 +23,8 @@
 !   nn_euler    =       1   !  = 0 : start with forward time step if ln_rstart=.true.
 !   cn_ocerst_in   = "ammsurge_restart_oce"   !  suffix of ocean restart name (input)
 !   cn_ocerst_out  = "restart_oce_out"   !  suffix of ocean restart name (input)
-   nn_istate   =       0   !  output the initial state (1) or not (0)
-   nn_stock    =  14400 ! 40320 =28days ! frequency of creation of a restart file (modulo referenced to 1) ! 28 DAYS
+   nn_istate   =       1   !  output the initial state (1) or not (0)
+   nn_stock    =  10 ! 40320 =28days ! frequency of creation of a restart file (modulo referenced to 1) ! 28 DAYS
    nn_write    =  10 ! 86400    !  frequency of write in the output file   (modulo referenced to nit000) 
 /
 !-----------------------------------------------------------------------
@@ -52,7 +52,7 @@
 &namdom        !   space and time domain (bathymetry, mesh, timestep)
 !-----------------------------------------------------------------------
 !   ln_2d        = .false.  !  (=T) run in 2D barotropic mode (no tracer processes or vertical diffusion)
-   rn_rdt      =   60.    !  time step for the dynamics (and tracer if nn_acc=0)
+   rn_rdt      =   30.    !  time step for the dynamics (and tracer if nn_acc=0)
 /
 
 !-----------------------------------------------------------------------
@@ -68,7 +68,7 @@
    !
    cn_dir      = './ICS/'      !  root directory for the location of the runoff files
    ln_tsd_init = .true.    !  Initialisation of ocean T & S with T & S input data (T) or not (F)
-   ln_tsd_interp = .false.    !  Interpolation of T & S in the verticalinput data (T) or not (F)
+   ln_tsd_interp = .true.    !  Interpolation of T & S in the verticalinput data (T) or not (F)
    ln_tsd_tradmp = .false.  !  damping of ocean T & S toward T & S input data (T) or not (F)
 /
 !-----------------------------------------------------------------------
@@ -76,8 +76,8 @@
 !-----------------------------------------------------------------------
    nn_fsbc     = 1         !  frequency of surface boundary condition computation
                            !     (also = the frequency of sea-ice model call)
-   ln_usr =  .false.
-   ln_blk =  .true.
+   ln_usr =  .true.
+   ln_blk =  .false.
    ln_apr_dyn  = .false.    !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr )
    nn_ice      = 0         !  =0 no ice boundary condition   ,
    ln_rnf      = .true.   !  Runoffs                                   (T => fill namsbc_rnf) 
@@ -174,7 +174,7 @@
    nb_bdy         = 1                    !  number of open boundary sets
    cn_coords_file = 'coordinates.bdy.nc' !  bdy coordinates files
    cn_dyn2d       = 'flather'            ! was flather
-   nn_dyn2d_dta   =  0                   !  = 0, bdy data are equal to the initial state
+   nn_dyn2d_dta   =  2                   !  = 0, bdy data are equal to the initial state
                                          !  = 1, bdy data are read in 'bdydata   .nc' files
                                          !  = 2, use tidal harmonic forcing data from files
                                          !  = 3, use external data AND tidal harmonic forcing
@@ -271,13 +271,30 @@
    ln_hpg_sco  = .true.    !  s-coordinate (Standard Jacobian scheme)
 /
 !-----------------------------------------------------------------------
-&namdyn_spg    !   surface pressure gradient   (CPP key only)
+&namdyn_spg    !   surface pressure gradient                            (default: NO)
 !-----------------------------------------------------------------------
-   ln_dynspg_ts = .true.    ! split-explicit free surface
-   ln_bt_auto =    .true.           !  Set nn_baro automatically to be just below
-                                       !  a user defined maximum courant number (rn_bt_cmax)
+  ln_dynspg_exp  = .false.   ! explicit free surface
+  ln_dynspg_ts   = .true.   ! split-explicit free surface
+      ln_bt_fw      = .true.     ! Forward integration of barotropic Eqs.
+      ln_bt_av      = .true.     ! Time filtering of barotropic variables
+        nn_bt_flt     = 1          ! Time filter choice  = 0 None
+         !                          !                     = 1 Boxcar over   nn_baro sub-steps
+        !                          !                     = 2 Boxcar over 2*nn_baro  "    "
+      ln_bt_auto    = .true.     ! Number of sub-step defined from:
+       rn_bt_cmax   =  0.5        ! =T : the Maximum Courant Number allowed
+        nn_baro      = 30          ! =F : the number of sub-step in rn_rdt seconds
+      ln_ulimit     =.true.
+       cn_ulimit     = 0.9
+       cnn_ulimit    = 0.54
 /
 !-----------------------------------------------------------------------
+!&namdyn_spg    !   surface pressure gradient   (CPP key only)
+!-----------------------------------------------------------------------
+!   ln_dynspg_ts = .true.    ! split-explicit free surface
+!   ln_bt_auto =    .true.           !  Set nn_baro automatically to be just below
+!                                       !  a user defined maximum courant number (rn_bt_cmax)
+!/
+!-----------------------------------------------------------------------
 &namdyn_ldf    !   lateral diffusion on momentum
 !-----------------------------------------------------------------------
    !                       !  Type of the operator :