Index: /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/ERA5_arctic/oce+ice/namelist_ecmwf-cstc_cfg
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/ERA5_arctic/oce+ice/namelist_ecmwf-cstc_cfg	(revision 13675)
+++ /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/ERA5_arctic/oce+ice/namelist_ecmwf-cstc_cfg	(revision 13676)
@@ -2,5 +2,5 @@
 !! NEMO/OCE  Configuration namelist : overwrite default values defined in SHARED/namelist_ref
 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-!!             STATION_ASF configuration with ICE, ERA5 forcing
+!!             STATION_ASF configuration with ICE, ERA5_arctic forcing
 !!======================================================================
 !!              ***  Domain & Run management namelists  ***           !!
@@ -26,5 +26,5 @@
 &namrun        !   parameters of the run
 !-----------------------------------------------------------------------
-   cn_exp      =  'STATION_ASF-ECMWF-CSTC'  !  experience name
+   cn_exp      =  'STATION_ASF-ECMWF-CSTC_ERA5_arctic'  !  experience name
    nn_it000    =    1   !  first time step
    nn_itend    =    8760   !  last  time step => 3 years (including 1 leap!) at dt=3600s
@@ -75,5 +75,5 @@
 !!            ***  Surface Boundary Condition namelists  ***          !!
 !!                                                                    !!
-!!                           STATION_ASF case                         !!
+!!                           STATION_ASF test-case                         !!
 !!                                                                    !!
 !!   namsbc          surface boundary condition manager                 (default: NO selection)
Index: /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/ERA5_arctic/oce+ice/namelist_ecmwf-lg15_cfg
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/ERA5_arctic/oce+ice/namelist_ecmwf-lg15_cfg	(revision 13675)
+++ /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/ERA5_arctic/oce+ice/namelist_ecmwf-lg15_cfg	(revision 13676)
@@ -2,5 +2,5 @@
 !! NEMO/OCE  Configuration namelist : overwrite default values defined in SHARED/namelist_ref
 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-!!             STATION_ASF configuration with ICE, ERA5 forcing
+!!             STATION_ASF configuration with ICE, ERA5_arctic forcing
 !!======================================================================
 !!              ***  Domain & Run management namelists  ***           !!
@@ -26,5 +26,5 @@
 &namrun        !   parameters of the run
 !-----------------------------------------------------------------------
-   cn_exp      =  'STATION_ASF-ECMWF-LG15'  !  experience name
+   cn_exp      =  'STATION_ASF-ECMWF-LG15_ERA5_arctic'  !  experience name
    nn_it000    =    1   !  first time step
    nn_itend    =    8760   !  last  time step => 3 years (including 1 leap!) at dt=3600s
Index: /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/ERA5_arctic/oce+ice/namelist_ecmwf-lu12_cfg
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/ERA5_arctic/oce+ice/namelist_ecmwf-lu12_cfg	(revision 13675)
+++ /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/ERA5_arctic/oce+ice/namelist_ecmwf-lu12_cfg	(revision 13676)
@@ -2,5 +2,5 @@
 !! NEMO/OCE  Configuration namelist : overwrite default values defined in SHARED/namelist_ref
 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-!!             STATION_ASF configuration with ICE, ERA5 forcing
+!!             STATION_ASF configuration with ICE, ERA5_arctic forcing
 !!======================================================================
 !!              ***  Domain & Run management namelists  ***           !!
@@ -26,5 +26,5 @@
 &namrun        !   parameters of the run
 !-----------------------------------------------------------------------
-   cn_exp      =  'STATION_ASF-ECMWF-LU12'  !  experience name
+   cn_exp      =  'STATION_ASF-ECMWF-LU12_ERA5_arctic'  !  experience name
    nn_it000    =    1   !  first time step
    nn_itend    =    8760   !  last  time step => 3 years (including 1 leap!) at dt=3600s
Index: /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/ERA5_arctic/oce/namelist_andreas_cfg
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/ERA5_arctic/oce/namelist_andreas_cfg	(revision 13675)
+++ /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/ERA5_arctic/oce/namelist_andreas_cfg	(revision 13676)
@@ -2,5 +2,5 @@
 !! NEMO/OCE  Configuration namelist : overwrite default values defined in SHARED/namelist_ref
 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-!!             STATION_ASF configuration, ERA5 forcing
+!!             STATION_ASF configuration, ERA5_arctic forcing
 !!======================================================================
 !!              ***  Domain & Run management namelists  ***           !!
@@ -26,5 +26,5 @@
 &namrun        !   parameters of the run
 !-----------------------------------------------------------------------
-   cn_exp      =  'STATION_ASF-ANDREAS'  !  experience name
+   cn_exp      =  'STATION_ASF-ANDREAS_ERA5_arctic'  !  experience name
    nn_it000    =    1   !  first time step
    nn_itend    =    8760   !  last  time step => 3 years (including 1 leap!) at dt=3600s
Index: /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/ERA5_arctic/oce/namelist_coare3p6_cfg
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/ERA5_arctic/oce/namelist_coare3p6_cfg	(revision 13675)
+++ /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/ERA5_arctic/oce/namelist_coare3p6_cfg	(revision 13676)
@@ -2,5 +2,5 @@
 !! NEMO/OCE  Configuration namelist : overwrite default values defined in SHARED/namelist_ref
 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-!!             STATION_ASF configuration, ERA5 forcing
+!!             STATION_ASF configuration, ERA5_arctic forcing
 !!======================================================================
 !!              ***  Domain & Run management namelists  ***           !!
@@ -26,5 +26,5 @@
 &namrun        !   parameters of the run
 !-----------------------------------------------------------------------
-   cn_exp      =  'STATION_ASF-COARE3p6'  !  experience name
+   cn_exp      =  'STATION_ASF-COARE3p6_ERA5_arctic'  !  experience name
    nn_it000    =    1   !  first time step
    nn_itend    =    8760   !  last  time step => 3 years (including 1 leap!) at dt=3600s
@@ -115,6 +115,6 @@
       rn_zu        = 10.     !  Wind vector reference height (m)
       nn_iter_algo = 20      !  Number of iterations in bulk param. algo ("stable ABL + weak wind" requires more)
-      ln_skin_cs = .true.    !  use the cool-skin parameterization
-      ln_skin_wl = .true.    !  use the warm-layer parameterization
+      ln_skin_cs   = .true.  !  use the cool-skin parameterization
+      ln_skin_wl   = .true.  !  use the warm-layer parameterization
    !
    ln_humi_sph = .false. !  humidity "sn_humi" is specific humidity  [kg/kg]
Index: /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/ERA5_arctic/oce/namelist_ecmwf_cfg
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/ERA5_arctic/oce/namelist_ecmwf_cfg	(revision 13675)
+++ /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/ERA5_arctic/oce/namelist_ecmwf_cfg	(revision 13676)
@@ -2,5 +2,5 @@
 !! NEMO/OCE  Configuration namelist : overwrite default values defined in SHARED/namelist_ref
 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-!!             STATION_ASF configuration, ERA5 forcing
+!!             STATION_ASF configuration, ERA5_arctic forcing
 !!======================================================================
 !!              ***  Domain & Run management namelists  ***           !!
@@ -26,5 +26,5 @@
 &namrun        !   parameters of the run
 !-----------------------------------------------------------------------
-   cn_exp      =  'STATION_ASF-ECMWF'  !  experience name
+   cn_exp      =  'STATION_ASF-ECMWF_ERA5_arctic'  !  experience name
    nn_it000    =    1   !  first time step
    nn_itend    =    8760   !  last  time step => 3 years (including 1 leap!) at dt=3600s
Index: /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/ERA5_arctic/oce/namelist_ncar_cfg
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/ERA5_arctic/oce/namelist_ncar_cfg	(revision 13675)
+++ /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/ERA5_arctic/oce/namelist_ncar_cfg	(revision 13676)
@@ -2,5 +2,5 @@
 !! NEMO/OCE  Configuration namelist : overwrite default values defined in SHARED/namelist_ref
 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-!!             STATION_ASF configuration, ERA5 forcing
+!!             STATION_ASF configuration, ERA5_arctic forcing
 !!======================================================================
 !!              ***  Domain & Run management namelists  ***           !!
@@ -26,5 +26,5 @@
 &namrun        !   parameters of the run
 !-----------------------------------------------------------------------
-   cn_exp      =  'STATION_ASF-NCAR'  !  experience name
+   cn_exp      =  'STATION_ASF-NCAR_ERA5_arctic'  !  experience name
    nn_it000    =    1   !  first time step
    nn_itend    =    8760   !  last  time step => 3 years (including 1 leap!) at dt=3600s
Index: MO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/namelist_andreas_cfg
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/namelist_andreas_cfg	(revision 13675)
+++ 	(revision )
@@ -1,274 +1,0 @@
-!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-!! NEMO/OCE  Configuration namelist : overwrite default values defined in SHARED/namelist_ref
-!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-!!
-!!             STATION_ASF configuration, IDEALIZED forcing
-!!                   ==> meant to be used with forcing from file:
-!!                       "input_output_VALIDATION_IDEALIZED.nc"
-!!
-!!======================================================================
-!!              ***  Domain & Run management namelists  ***           !!
-!!                                                                    !!
-!!   namrun       parameters of the run
-!!   namdom       space and time domain
-!!   namcfg       parameters of the configuration                       (default: user defined GYRE)
-!!   namwad       Wetting and drying                                    (default: OFF)
-!!   namtsd       data: temperature & salinity                          (default: OFF)
-!!   namcrs       coarsened grid (for outputs and/or TOP)               (ln_crs =T)
-!!   namc1d       1D configuration options                              ("key_c1d")
-!!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d")
-!!   namc1d_uvd   1D data (currents)                                    ("key_c1d")
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&namusr_def    !   STATION_ASF user defined namelist
-!-----------------------------------------------------------------------
-   rn_dept1    =  1.0   ! depth (m) at which the prescribed SST is taken  == depth of first T point! (IMPORTANT for warm-layer impact, the deeper, the more pronounced!)
-   !                    ! [PAPA station = 1 m]
-/
-!-----------------------------------------------------------------------
-&namrun        !   parameters of the run
-!-----------------------------------------------------------------------
-   cn_exp      =  'STATION_ASF-ANDREAS-IDEALIZED'  !  experience name
-   nn_it000    =    1   !  first time step
-   nn_itend    =       104  !  last  time step => 3 years (including 1 leap!) at dt=3600s
-   nn_date0    = 20200101  !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
-   nn_leapy    =       1   !  Leap year calendar (1) or not (0)
-   !
-   nn_stock    =    104
-   nn_write    =    104
-/
-!-----------------------------------------------------------------------
-&namdom        !   time and space domain / STATION_ASF
-!-----------------------------------------------------------------------
-   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time
-   rn_Dt       =  3600.    !  time step for the dynamics and tracer => frequency of atmospheric forcing
-   ln_meshmask = .false.
-/
-!-----------------------------------------------------------------------
-&namcfg        !   parameters of the configuration                      (default: use namusr_def in namelist_cfg)
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namc1d        !   1D configuration options / STATION_ASF               ("key_c1d" default: PAPA station)
-!-----------------------------------------------------------------------
-   rn_lat1d    =      50.1 !  Column latitude    [PAPA station]
-   rn_lon1d    =    215.1  !  Column longitude   [PAPA station]
-/
-!-----------------------------------------------------------------------
-&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d" default: OFF)
-!-----------------------------------------------------------------------
-   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F)
-/
-!-----------------------------------------------------------------------
-&namc1d_uvd    !   data: U & V currents                                 ("key_c1d" default: OFF)
-!-----------------------------------------------------------------------
-   !                       !  =T read U-V fields for:
-   ln_uvd_init   = .false.       !  ocean initialisation
-   ln_uvd_dyndmp = .false.       !  U-V restoring
-
-   cn_dir      = './'      !  root directory for the U-V data location
-   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
-   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
-   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
-   sn_ucur     = 'ucurrent'              ,         -1.       ,'u_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Ume'   , ''
-   sn_vcur     = 'vcurrent'              ,         -1.       ,'v_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Vme'   , ''
-/
-
-!!======================================================================
-!!            ***  Surface Boundary Condition namelists  ***          !!
-!!                                                                    !!
-!!                           STATION_ASF test-case                         !!
-!!                                                                    !!
-!!   namsbc          surface boundary condition manager                 (default: NO selection)
-!!   namsbc_flx      flux               formulation                     (ln_flx     =T)
-!!   namsbc_blk      Bulk formulae formulation                          (ln_blk     =T)
-!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3" )
-!!   namsbc_sas      Stand-Alone Surface module                         (SAS_SRC  only)
-!!   namsbc_iif      Ice-IF: use observed ice cover                     (nn_ice = 1   )
-!!   namtra_qsr      penetrative solar radiation                        (ln_traqsr  =T)
-!!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T)
-!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T)
-!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T)
-!!   namsbc_wave     external fields from wave model                    (ln_wave    =T)
-!!   namberg         iceberg floats                                     (ln_icebergs=T)
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&namsbc        !   Surface Boundary Condition manager                   (default: NO selection)
-!-----------------------------------------------------------------------
-   nn_fsbc     = 1         !!LOLO!  frequency of SBC module call
-   ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk )
-   !
-   nn_ice      = 0  !LOLO  !  =0 no ice boundary condition
-      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif )
-      !                    !  =2 or 3 for SI3 and CICE, respectively
-   !    
-   ln_dm2dc    = .false.   !#lolo  daily mean to diurnal cycle on short wave
-/
-!-----------------------------------------------------------------------
-&namsbc_blk    !   namsbc_blk  generic Bulk formula / STATION_ASF       (ln_blk =T)
-!-----------------------------------------------------------------------
-   !                    !  bulk algorithm :
-   ln_NCAR      = .false.    ! "NCAR"      algorithm   (Large and Yeager 2008)
-   ln_COARE_3p0 = .false.    ! "COARE 3.0" algorithm   (Fairall et al. 2003)
-   ln_COARE_3p6 = .false.    ! "COARE 3.6" algorithm   (Edson et al. 2013)
-   ln_ECMWF     = .false.    ! "ECMWF"     algorithm   (IFS cycle 45r1)
-   ln_ANDREAS   = .true.     ! "ANDREAS"   algorithm   (Andreas et al. 2015)
-      rn_zqt       =  2.     !  Air temperature & humidity reference height (m)
-      rn_zu        = 10.     !  Wind vector reference height (m)
-      nn_iter_algo = 20      !  Number of iterations in bulk param. algo ("stable ABL + weak wind" requires more)
-      ln_skin_cs   = .false. !  use the cool-skin parameterization
-      ln_skin_wl   = .false. !  use the warm-layer parameterization
-   !
-   ln_humi_sph = .false. !  humidity "sn_humi" is specific humidity  [kg/kg]
-   ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K]
-   ln_humi_rlh = .true.  !  humidity "sn_humi" is relative humidity     [%]
-   ln_tpot     = .true.  !!GS: compute potential temperature or not
-   !
-   cn_dir      = './'      !  root directory for the bulk data location
-   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________!
-   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !       weights filename               ! rotation ! land/sea mask !
-   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                                      ! pairing  !    filename   !
-   sn_wndi     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'u10'    ,   .false.   , .true.  , 'yearly'  , ''   ,''   , ''
-   sn_wndj     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'v10'    ,   .false.   , .true.  , 'yearly'  , ''   ,''   , ''
-   sn_qsr      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rad_sw' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_qlw      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rad_lw' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_tair     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  't_air'  ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_humi     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rh_air' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_prec     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'precip' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_snow     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'snow'   ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_slp      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'slp'    ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-/
-!-----------------------------------------------------------------------
-&namsbc_abl    !   Atmospheric Boundary Layer formulation           (ln_abl = T)
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3")
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namsbc_sas    !   Stand-Alone Surface boundary condition / STATION_ASF
-!-----------------------------------------------------------------------
-   l_sasread   = .true.    !  =T Read in file ;  =F set all to 0. (see sbcssm)
-   !
-   cn_dir      = './'      !  root directory for the ocean data location
-   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
-   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
-   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
-   sn_usp  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssu'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   sn_vsp  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssv'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   sn_tem  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'sst'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   sn_sal  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'sss',        .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   sn_ssh  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssh',        .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   !!
-/
-!!======================================================================
-!!               ***  Lateral boundary condition  ***                 !!
-!!                                                                    !!
-!!   namlbc        lateral momentum boundary condition                  (default: NO selection)
-!!   namagrif      agrif nested grid   (read by child model only)       ("key_agrif")
-!!   nam_tide      Tidal forcing                                        (default: OFF)
-!!   nambdy        Unstructured open boundaries                         (default: OFF)
-!!   nambdy_dta    Unstructured open boundaries - external data         (see  nambdy)
-!!   nambdy_tide   tidal forcing at open boundaries                     (default: OFF)
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&namlbc        !   lateral momentum boundary condition                  (default: NO selection)
-!-----------------------------------------------------------------------
-   rn_shlat    =    0.     !  free slip
-/
-!!======================================================================
-!!                ***  Top/Bottom boundary condition  ***             !!
-!!                                                                    !!
-!!   namdrg        top/bottom drag coefficient                          (default: NO selection)
-!!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T)
-!!   namdrg_bot    bottom friction                                      (ln_OFF=F)
-!!   nambbc        bottom temperature boundary condition                (default: OFF)
-!!   nambbl        bottom boundary layer scheme                         (default: OFF)
-!!======================================================================
-!
-!!======================================================================
-!!                        Tracer (T-S) namelists                      !!
-!!                                                                    !!
-!!   nameos        equation of state                                    (default: NO selection)
-!!   namtra_adv    advection scheme                                     (default: NO selection)
-!!   namtra_ldf    lateral diffusion scheme                             (default: NO selection)
-!!   namtra_mle    mixed layer eddy param. (Fox-Kemper param.)          (default: OFF)
-!!   namtra_eiv    eddy induced velocity param.                         (default: OFF)
-!!   namtra_dmp    T & S newtonian damping                              (default: OFF)
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&nameos        !   ocean Equation Of Seawater                           (default: NO selection)
-!-----------------------------------------------------------------------
-   ln_eos80    = .true.         !  = Use EOS80
-/
-!!======================================================================
-!!                      ***  Dynamics namelists  ***                  !!
-!!                                                                    !!
-!!   nam_vvl       vertical coordinate options                          (default: z-star)
-!!   namdyn_adv    formulation of the momentum advection                (default: NO selection)
-!!   namdyn_vor    advection scheme                                     (default: NO selection)
-!!   namdyn_hpg    hydrostatic pressure gradient                        (default: NO selection)
-!!   namdyn_spg    surface pressure gradient                            (default: NO selection)
-!!   namdyn_ldf    lateral diffusion scheme                             (default: NO selection)
-!!   namdta_dyn    offline TOP: dynamics read in files                  (OFF_SRC only)
-!!======================================================================
-!
-!!======================================================================
-!!                     vertical physics namelists                     !!
-!!                                                                    !!
-!!    namzdf        vertical physics manager                            (default: NO selection)
-!!    namzdf_ric    richardson number vertical mixing                   (ln_zdfric=T)
-!!    namzdf_tke    TKE vertical mixing                                 (ln_zdftke=T)
-!!    namzdf_gls    GLS vertical mixing                                 (ln_zdfgls=T)
-!!    namzdf_osm    OSM vertical diffusion                              (ln_zdfosm=T)
-!!    namzdf_iwm    tidal mixing parameterization                       (ln_zdfiwm=T)
-!!======================================================================
-!
-!!======================================================================
-!!                  ***  Diagnostics namelists  ***                   !!
-!!                                                                    !!
-!!   namtrd       dynamics and/or tracer trends                         (default: OFF)
-!!   namptr       Poleward Transport Diagnostics                        (default: OFF)
-!!   namhsb       Heat and salt budgets                                 (default: OFF)
-!!   namdiu       Cool skin and warm layer models                       (default: OFF)
-!!   namdiu       Cool skin and warm layer models                       (default: OFF)
-!!   namflo       float parameters                                      (default: OFF)
-!!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF)
-!!   nam_diadct   transports through some sections                      (default: OFF)
-!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF)
-!!   nam_dia25h   25h Mean Output                                       (default: OFF)
-!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4")
-!!======================================================================
-!
-!!======================================================================
-!!               ***  Observation & Assimilation  ***                 !!
-!!                                                                    !!
-!!   namobs       observation and model comparison                      (default: OFF)
-!!   nam_asminc   assimilation increments                               ('key_asminc')
-!!======================================================================
-!
-!!======================================================================
-!!                  ***  Miscellaneous namelists  ***                 !!
-!!                                                                    !!
-!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi")
-!!   namctl            Control prints                                   (default: OFF)
-!!   namsto            Stochastic parametrization of EOS                (default: OFF)
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi")
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namctl        !   Control prints                                       (default: OFF)
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namsto        ! Stochastic parametrization of EOS                      (default: OFF)
-!-----------------------------------------------------------------------
-/
Index: MO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/namelist_coare3p0_cfg
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/namelist_coare3p0_cfg	(revision 13675)
+++ 	(revision )
@@ -1,274 +1,0 @@
-!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-!! NEMO/OCE  Configuration namelist : overwrite default values defined in SHARED/namelist_ref
-!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-!!
-!!             STATION_ASF configuration, IDEALIZED forcing
-!!                   ==> meant to be used with forcing from file:
-!!                       "input_output_VALIDATION_IDEALIZED.nc"
-!!
-!!======================================================================
-!!              ***  Domain & Run management namelists  ***           !!
-!!                                                                    !!
-!!   namrun       parameters of the run
-!!   namdom       space and time domain
-!!   namcfg       parameters of the configuration                       (default: user defined GYRE)
-!!   namwad       Wetting and drying                                    (default: OFF)
-!!   namtsd       data: temperature & salinity                          (default: OFF)
-!!   namcrs       coarsened grid (for outputs and/or TOP)               (ln_crs =T)
-!!   namc1d       1D configuration options                              ("key_c1d")
-!!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d")
-!!   namc1d_uvd   1D data (currents)                                    ("key_c1d")
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&namusr_def    !   STATION_ASF user defined namelist
-!-----------------------------------------------------------------------
-   rn_dept1    =  1.0   ! depth (m) at which the prescribed SST is taken  == depth of first T point! (IMPORTANT for warm-layer impact, the deeper, the more pronounced!)
-   !                    ! [PAPA station = 1 m]
-/
-!-----------------------------------------------------------------------
-&namrun        !   parameters of the run
-!-----------------------------------------------------------------------
-   cn_exp      =  'STATION_ASF-COARE3p0-IDEALIZED'  !  experience name
-   nn_it000    =    1   !  first time step
-   nn_itend    =       104  !  last  time step => 3 years (including 1 leap!) at dt=3600s
-   nn_date0    = 20200101  !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
-   nn_leapy    =       1   !  Leap year calendar (1) or not (0)
-   !
-   nn_stock    =    104
-   nn_write    =    104
-/
-!-----------------------------------------------------------------------
-&namdom        !   time and space domain / STATION_ASF
-!-----------------------------------------------------------------------
-   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time
-   rn_Dt       =  3600.    !  time step for the dynamics and tracer => frequency of atmospheric forcing
-   ln_meshmask = .false.
-/
-!-----------------------------------------------------------------------
-&namcfg        !   parameters of the configuration                      (default: use namusr_def in namelist_cfg)
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namc1d        !   1D configuration options / STATION_ASF               ("key_c1d" default: PAPA station)
-!-----------------------------------------------------------------------
-   rn_lat1d    =      50.1 !  Column latitude    [PAPA station]
-   rn_lon1d    =    215.1  !  Column longitude   [PAPA station]
-/
-!-----------------------------------------------------------------------
-&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d" default: OFF)
-!-----------------------------------------------------------------------
-   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F)
-/
-!-----------------------------------------------------------------------
-&namc1d_uvd    !   data: U & V currents                                 ("key_c1d" default: OFF)
-!-----------------------------------------------------------------------
-   !                       !  =T read U-V fields for:
-   ln_uvd_init   = .false.       !  ocean initialisation
-   ln_uvd_dyndmp = .false.       !  U-V restoring
-
-   cn_dir      = './'      !  root directory for the U-V data location
-   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
-   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
-   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
-   sn_ucur     = 'ucurrent'              ,         -1.       ,'u_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Ume'   , ''
-   sn_vcur     = 'vcurrent'              ,         -1.       ,'v_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Vme'   , ''
-/
-
-!!======================================================================
-!!            ***  Surface Boundary Condition namelists  ***          !!
-!!                                                                    !!
-!!                           STATION_ASF test-case                    !!
-!!                                                                    !!
-!!   namsbc          surface boundary condition manager                 (default: NO selection)
-!!   namsbc_flx      flux               formulation                     (ln_flx     =T)
-!!   namsbc_blk      Bulk formulae formulation                          (ln_blk     =T)
-!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3" )
-!!   namsbc_sas      Stand-Alone Surface module                         (SAS_SRC  only)
-!!   namsbc_iif      Ice-IF: use observed ice cover                     (nn_ice = 1   )
-!!   namtra_qsr      penetrative solar radiation                        (ln_traqsr  =T)
-!!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T)
-!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T)
-!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T)
-!!   namsbc_wave     external fields from wave model                    (ln_wave    =T)
-!!   namberg         iceberg floats                                     (ln_icebergs=T)
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&namsbc        !   Surface Boundary Condition manager                   (default: NO selection)
-!-----------------------------------------------------------------------
-   nn_fsbc     = 1         !!LOLO!  frequency of SBC module call
-   ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk )
-   !
-   nn_ice      = 0  !LOLO  !  =0 no ice boundary condition
-      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif )
-      !                    !  =2 or 3 for SI3 and CICE, respectively
-   !    
-   ln_dm2dc    = .false.   !#lolo  daily mean to diurnal cycle on short wave
-/
-!-----------------------------------------------------------------------
-&namsbc_blk    !   namsbc_blk  generic Bulk formula / STATION_ASF       (ln_blk =T)
-!-----------------------------------------------------------------------
-   !                    !  bulk algorithm :
-   ln_NCAR      = .false.    ! "NCAR"      algorithm   (Large and Yeager 2008)
-   ln_COARE_3p0 = .true.     ! "COARE 3.0" algorithm   (Fairall et al. 2003)
-   ln_COARE_3p6 = .false.    ! "COARE 3.6" algorithm   (Edson et al. 2013)
-   ln_ECMWF     = .false.    ! "ECMWF"     algorithm   (IFS cycle 45r1)
-   ln_ANDREAS   = .false.    ! "ANDREAS"   algorithm   (Andreas et al. 2015)
-      rn_zqt       =  2.     !  Air temperature & humidity reference height (m)
-      rn_zu        = 10.     !  Wind vector reference height (m)
-      nn_iter_algo = 20      !  Number of iterations in bulk param. algo ("stable ABL + weak wind" requires more)
-      ln_skin_cs   = .true.  !  use the cool-skin parameterization
-      ln_skin_wl   = .true.  !  use the warm-layer parameterization
-   !
-   ln_humi_sph = .false. !  humidity "sn_humi" is specific humidity  [kg/kg]
-   ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K]
-   ln_humi_rlh = .true.  !  humidity "sn_humi" is relative humidity     [%]
-   ln_tpot     = .true.  !!GS: compute potential temperature or not
-   !
-   cn_dir      = './'      !  root directory for the bulk data location
-   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________!
-   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !       weights filename               ! rotation ! land/sea mask !
-   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                                      ! pairing  !    filename   !
-   sn_wndi     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'u10'    ,   .false.   , .true.  , 'yearly'  , ''   ,''   , ''
-   sn_wndj     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'v10'    ,   .false.   , .true.  , 'yearly'  , ''   ,''   , ''
-   sn_qsr      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rad_sw' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_qlw      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rad_lw' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_tair     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  't_air'  ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_humi     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rh_air' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_prec     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'precip' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_snow     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'snow'   ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_slp      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'slp'    ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-/
-!-----------------------------------------------------------------------
-&namsbc_abl    !   Atmospheric Boundary Layer formulation           (ln_abl = T)
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3")
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namsbc_sas    !   Stand-Alone Surface boundary condition / STATION_ASF
-!-----------------------------------------------------------------------
-   l_sasread   = .true.    !  =T Read in file ;  =F set all to 0. (see sbcssm)
-   !
-   cn_dir      = './'      !  root directory for the ocean data location
-   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
-   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
-   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
-   sn_usp  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssu'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   sn_vsp  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssv'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   sn_tem  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'sst'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   sn_sal  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'sss',        .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   sn_ssh  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssh',        .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   !!
-/
-!!======================================================================
-!!               ***  Lateral boundary condition  ***                 !!
-!!                                                                    !!
-!!   namlbc        lateral momentum boundary condition                  (default: NO selection)
-!!   namagrif      agrif nested grid   (read by child model only)       ("key_agrif")
-!!   nam_tide      Tidal forcing                                        (default: OFF)
-!!   nambdy        Unstructured open boundaries                         (default: OFF)
-!!   nambdy_dta    Unstructured open boundaries - external data         (see  nambdy)
-!!   nambdy_tide   tidal forcing at open boundaries                     (default: OFF)
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&namlbc        !   lateral momentum boundary condition                  (default: NO selection)
-!-----------------------------------------------------------------------
-   rn_shlat    =    0.     !  free slip
-/
-!!======================================================================
-!!                ***  Top/Bottom boundary condition  ***             !!
-!!                                                                    !!
-!!   namdrg        top/bottom drag coefficient                          (default: NO selection)
-!!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T)
-!!   namdrg_bot    bottom friction                                      (ln_OFF=F)
-!!   nambbc        bottom temperature boundary condition                (default: OFF)
-!!   nambbl        bottom boundary layer scheme                         (default: OFF)
-!!======================================================================
-!
-!!======================================================================
-!!                        Tracer (T-S) namelists                      !!
-!!                                                                    !!
-!!   nameos        equation of state                                    (default: NO selection)
-!!   namtra_adv    advection scheme                                     (default: NO selection)
-!!   namtra_ldf    lateral diffusion scheme                             (default: NO selection)
-!!   namtra_mle    mixed layer eddy param. (Fox-Kemper param.)          (default: OFF)
-!!   namtra_eiv    eddy induced velocity param.                         (default: OFF)
-!!   namtra_dmp    T & S newtonian damping                              (default: OFF)
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&nameos        !   ocean Equation Of Seawater                           (default: NO selection)
-!-----------------------------------------------------------------------
-   ln_eos80    = .true.         !  = Use EOS80
-/
-!!======================================================================
-!!                      ***  Dynamics namelists  ***                  !!
-!!                                                                    !!
-!!   nam_vvl       vertical coordinate options                          (default: z-star)
-!!   namdyn_adv    formulation of the momentum advection                (default: NO selection)
-!!   namdyn_vor    advection scheme                                     (default: NO selection)
-!!   namdyn_hpg    hydrostatic pressure gradient                        (default: NO selection)
-!!   namdyn_spg    surface pressure gradient                            (default: NO selection)
-!!   namdyn_ldf    lateral diffusion scheme                             (default: NO selection)
-!!   namdta_dyn    offline TOP: dynamics read in files                  (OFF_SRC only)
-!!======================================================================
-!
-!!======================================================================
-!!                     vertical physics namelists                     !!
-!!                                                                    !!
-!!    namzdf        vertical physics manager                            (default: NO selection)
-!!    namzdf_ric    richardson number vertical mixing                   (ln_zdfric=T)
-!!    namzdf_tke    TKE vertical mixing                                 (ln_zdftke=T)
-!!    namzdf_gls    GLS vertical mixing                                 (ln_zdfgls=T)
-!!    namzdf_osm    OSM vertical diffusion                              (ln_zdfosm=T)
-!!    namzdf_iwm    tidal mixing parameterization                       (ln_zdfiwm=T)
-!!======================================================================
-!
-!!======================================================================
-!!                  ***  Diagnostics namelists  ***                   !!
-!!                                                                    !!
-!!   namtrd       dynamics and/or tracer trends                         (default: OFF)
-!!   namptr       Poleward Transport Diagnostics                        (default: OFF)
-!!   namhsb       Heat and salt budgets                                 (default: OFF)
-!!   namdiu       Cool skin and warm layer models                       (default: OFF)
-!!   namdiu       Cool skin and warm layer models                       (default: OFF)
-!!   namflo       float parameters                                      (default: OFF)
-!!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF)
-!!   nam_diadct   transports through some sections                      (default: OFF)
-!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF)
-!!   nam_dia25h   25h Mean Output                                       (default: OFF)
-!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4")
-!!======================================================================
-!
-!!======================================================================
-!!               ***  Observation & Assimilation  ***                 !!
-!!                                                                    !!
-!!   namobs       observation and model comparison                      (default: OFF)
-!!   nam_asminc   assimilation increments                               ('key_asminc')
-!!======================================================================
-!
-!!======================================================================
-!!                  ***  Miscellaneous namelists  ***                 !!
-!!                                                                    !!
-!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi")
-!!   namctl            Control prints                                   (default: OFF)
-!!   namsto            Stochastic parametrization of EOS                (default: OFF)
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi")
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namctl        !   Control prints                                       (default: OFF)
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namsto        ! Stochastic parametrization of EOS                      (default: OFF)
-!-----------------------------------------------------------------------
-/
Index: MO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/namelist_coare3p6_cfg
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/namelist_coare3p6_cfg	(revision 13675)
+++ 	(revision )
@@ -1,274 +1,0 @@
-!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-!! NEMO/OCE  Configuration namelist : overwrite default values defined in SHARED/namelist_ref
-!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-!!
-!!             STATION_ASF configuration, IDEALIZED forcing
-!!                   ==> meant to be used with forcing from file:
-!!                       "input_output_VALIDATION_IDEALIZED.nc"
-!!
-!!======================================================================
-!!              ***  Domain & Run management namelists  ***           !!
-!!                                                                    !!
-!!   namrun       parameters of the run
-!!   namdom       space and time domain
-!!   namcfg       parameters of the configuration                       (default: user defined GYRE)
-!!   namwad       Wetting and drying                                    (default: OFF)
-!!   namtsd       data: temperature & salinity                          (default: OFF)
-!!   namcrs       coarsened grid (for outputs and/or TOP)               (ln_crs =T)
-!!   namc1d       1D configuration options                              ("key_c1d")
-!!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d")
-!!   namc1d_uvd   1D data (currents)                                    ("key_c1d")
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&namusr_def    !   STATION_ASF user defined namelist
-!-----------------------------------------------------------------------
-   rn_dept1    =  1.0   ! depth (m) at which the prescribed SST is taken  == depth of first T point! (IMPORTANT for warm-layer impact, the deeper, the more pronounced!)
-   !                    ! [PAPA station = 1 m]
-/
-!-----------------------------------------------------------------------
-&namrun        !   parameters of the run
-!-----------------------------------------------------------------------
-   cn_exp      =  'STATION_ASF-COARE3p6-IDEALIZED'  !  experience name
-   nn_it000    =    1   !  first time step
-   nn_itend    =       104  !  last  time step => 3 years (including 1 leap!) at dt=3600s
-   nn_date0    = 20200101  !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
-   nn_leapy    =       1   !  Leap year calendar (1) or not (0)
-   !
-   nn_stock    =    104
-   nn_write    =    104
-/
-!-----------------------------------------------------------------------
-&namdom        !   time and space domain / STATION_ASF
-!-----------------------------------------------------------------------
-   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time
-   rn_Dt       =  3600.    !  time step for the dynamics and tracer => frequency of atmospheric forcing
-   ln_meshmask = .false.
-/
-!-----------------------------------------------------------------------
-&namcfg        !   parameters of the configuration                      (default: use namusr_def in namelist_cfg)
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namc1d        !   1D configuration options / STATION_ASF               ("key_c1d" default: PAPA station)
-!-----------------------------------------------------------------------
-   rn_lat1d    =      50.1 !  Column latitude    [PAPA station]
-   rn_lon1d    =    215.1  !  Column longitude   [PAPA station]
-/
-!-----------------------------------------------------------------------
-&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d" default: OFF)
-!-----------------------------------------------------------------------
-   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F)
-/
-!-----------------------------------------------------------------------
-&namc1d_uvd    !   data: U & V currents                                 ("key_c1d" default: OFF)
-!-----------------------------------------------------------------------
-   !                       !  =T read U-V fields for:
-   ln_uvd_init   = .false.       !  ocean initialisation
-   ln_uvd_dyndmp = .false.       !  U-V restoring
-
-   cn_dir      = './'      !  root directory for the U-V data location
-   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
-   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
-   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
-   sn_ucur     = 'ucurrent'              ,         -1.       ,'u_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Ume'   , ''
-   sn_vcur     = 'vcurrent'              ,         -1.       ,'v_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Vme'   , ''
-/
-
-!!======================================================================
-!!            ***  Surface Boundary Condition namelists  ***          !!
-!!                                                                    !!
-!!                           STATION_ASF test-case                    !!
-!!                                                                    !!
-!!   namsbc          surface boundary condition manager                 (default: NO selection)
-!!   namsbc_flx      flux               formulation                     (ln_flx     =T)
-!!   namsbc_blk      Bulk formulae formulation                          (ln_blk     =T)
-!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3" )
-!!   namsbc_sas      Stand-Alone Surface module                         (SAS_SRC  only)
-!!   namsbc_iif      Ice-IF: use observed ice cover                     (nn_ice = 1   )
-!!   namtra_qsr      penetrative solar radiation                        (ln_traqsr  =T)
-!!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T)
-!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T)
-!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T)
-!!   namsbc_wave     external fields from wave model                    (ln_wave    =T)
-!!   namberg         iceberg floats                                     (ln_icebergs=T)
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&namsbc        !   Surface Boundary Condition manager                   (default: NO selection)
-!-----------------------------------------------------------------------
-   nn_fsbc     = 1         !!LOLO!  frequency of SBC module call
-   ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk )
-   !
-   nn_ice      = 0  !LOLO  !  =0 no ice boundary condition
-      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif )
-      !                    !  =2 or 3 for SI3 and CICE, respectively
-   !    
-   ln_dm2dc    = .false.   !#lolo  daily mean to diurnal cycle on short wave
-/
-!-----------------------------------------------------------------------
-&namsbc_blk    !   namsbc_blk  generic Bulk formula / STATION_ASF       (ln_blk =T)
-!-----------------------------------------------------------------------
-   !                    !  bulk algorithm :
-   ln_NCAR      = .false.    ! "NCAR"      algorithm   (Large and Yeager 2008)
-   ln_COARE_3p0 = .false.    ! "COARE 3.0" algorithm   (Fairall et al. 2003)
-   ln_COARE_3p6 = .true.     ! "COARE 3.6" algorithm   (Edson et al. 2013)
-   ln_ECMWF     = .false.    ! "ECMWF"     algorithm   (IFS cycle 45r1)
-   ln_ANDREAS   = .false.    ! "ANDREAS"   algorithm   (Andreas et al. 2015)
-      rn_zqt       =  2.     !  Air temperature & humidity reference height (m)
-      rn_zu        = 10.     !  Wind vector reference height (m)
-      nn_iter_algo = 20      !  Number of iterations in bulk param. algo ("stable ABL + weak wind" requires more)
-      ln_skin_cs   = .true.  !  use the cool-skin parameterization
-      ln_skin_wl   = .true.  !  use the warm-layer parameterization
-   !
-   ln_humi_sph = .false. !  humidity "sn_humi" is specific humidity  [kg/kg]
-   ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K]
-   ln_humi_rlh = .true.  !  humidity "sn_humi" is relative humidity     [%]
-   ln_tpot     = .true.  !!GS: compute potential temperature or not
-   !
-   cn_dir      = './'      !  root directory for the bulk data location
-   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________!
-   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !       weights filename               ! rotation ! land/sea mask !
-   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                                      ! pairing  !    filename   !
-   sn_wndi     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'u10'    ,   .false.   , .true.  , 'yearly'  , ''   ,''   , ''
-   sn_wndj     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'v10'    ,   .false.   , .true.  , 'yearly'  , ''   ,''   , ''
-   sn_qsr      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rad_sw' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_qlw      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rad_lw' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_tair     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  't_air'  ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_humi     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rh_air' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_prec     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'precip' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_snow     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'snow'   ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_slp      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'slp'    ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-/
-!-----------------------------------------------------------------------
-&namsbc_abl    !   Atmospheric Boundary Layer formulation           (ln_abl = T)
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3")
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namsbc_sas    !   Stand-Alone Surface boundary condition / STATION_ASF
-!-----------------------------------------------------------------------
-   l_sasread   = .true.    !  =T Read in file ;  =F set all to 0. (see sbcssm)
-   !
-   cn_dir      = './'      !  root directory for the ocean data location
-   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
-   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
-   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
-   sn_usp  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssu'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   sn_vsp  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssv'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   sn_tem  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'sst'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   sn_sal  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'sss',        .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   sn_ssh  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssh',        .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   !!
-/
-!!======================================================================
-!!               ***  Lateral boundary condition  ***                 !!
-!!                                                                    !!
-!!   namlbc        lateral momentum boundary condition                  (default: NO selection)
-!!   namagrif      agrif nested grid   (read by child model only)       ("key_agrif")
-!!   nam_tide      Tidal forcing                                        (default: OFF)
-!!   nambdy        Unstructured open boundaries                         (default: OFF)
-!!   nambdy_dta    Unstructured open boundaries - external data         (see  nambdy)
-!!   nambdy_tide   tidal forcing at open boundaries                     (default: OFF)
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&namlbc        !   lateral momentum boundary condition                  (default: NO selection)
-!-----------------------------------------------------------------------
-   rn_shlat    =    0.     !  free slip
-/
-!!======================================================================
-!!                ***  Top/Bottom boundary condition  ***             !!
-!!                                                                    !!
-!!   namdrg        top/bottom drag coefficient                          (default: NO selection)
-!!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T)
-!!   namdrg_bot    bottom friction                                      (ln_OFF=F)
-!!   nambbc        bottom temperature boundary condition                (default: OFF)
-!!   nambbl        bottom boundary layer scheme                         (default: OFF)
-!!======================================================================
-!
-!!======================================================================
-!!                        Tracer (T-S) namelists                      !!
-!!                                                                    !!
-!!   nameos        equation of state                                    (default: NO selection)
-!!   namtra_adv    advection scheme                                     (default: NO selection)
-!!   namtra_ldf    lateral diffusion scheme                             (default: NO selection)
-!!   namtra_mle    mixed layer eddy param. (Fox-Kemper param.)          (default: OFF)
-!!   namtra_eiv    eddy induced velocity param.                         (default: OFF)
-!!   namtra_dmp    T & S newtonian damping                              (default: OFF)
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&nameos        !   ocean Equation Of Seawater                           (default: NO selection)
-!-----------------------------------------------------------------------
-   ln_eos80    = .true.         !  = Use EOS80
-/
-!!======================================================================
-!!                      ***  Dynamics namelists  ***                  !!
-!!                                                                    !!
-!!   nam_vvl       vertical coordinate options                          (default: z-star)
-!!   namdyn_adv    formulation of the momentum advection                (default: NO selection)
-!!   namdyn_vor    advection scheme                                     (default: NO selection)
-!!   namdyn_hpg    hydrostatic pressure gradient                        (default: NO selection)
-!!   namdyn_spg    surface pressure gradient                            (default: NO selection)
-!!   namdyn_ldf    lateral diffusion scheme                             (default: NO selection)
-!!   namdta_dyn    offline TOP: dynamics read in files                  (OFF_SRC only)
-!!======================================================================
-!
-!!======================================================================
-!!                     vertical physics namelists                     !!
-!!                                                                    !!
-!!    namzdf        vertical physics manager                            (default: NO selection)
-!!    namzdf_ric    richardson number vertical mixing                   (ln_zdfric=T)
-!!    namzdf_tke    TKE vertical mixing                                 (ln_zdftke=T)
-!!    namzdf_gls    GLS vertical mixing                                 (ln_zdfgls=T)
-!!    namzdf_osm    OSM vertical diffusion                              (ln_zdfosm=T)
-!!    namzdf_iwm    tidal mixing parameterization                       (ln_zdfiwm=T)
-!!======================================================================
-!
-!!======================================================================
-!!                  ***  Diagnostics namelists  ***                   !!
-!!                                                                    !!
-!!   namtrd       dynamics and/or tracer trends                         (default: OFF)
-!!   namptr       Poleward Transport Diagnostics                        (default: OFF)
-!!   namhsb       Heat and salt budgets                                 (default: OFF)
-!!   namdiu       Cool skin and warm layer models                       (default: OFF)
-!!   namdiu       Cool skin and warm layer models                       (default: OFF)
-!!   namflo       float parameters                                      (default: OFF)
-!!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF)
-!!   nam_diadct   transports through some sections                      (default: OFF)
-!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF)
-!!   nam_dia25h   25h Mean Output                                       (default: OFF)
-!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4")
-!!======================================================================
-!
-!!======================================================================
-!!               ***  Observation & Assimilation  ***                 !!
-!!                                                                    !!
-!!   namobs       observation and model comparison                      (default: OFF)
-!!   nam_asminc   assimilation increments                               ('key_asminc')
-!!======================================================================
-!
-!!======================================================================
-!!                  ***  Miscellaneous namelists  ***                 !!
-!!                                                                    !!
-!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi")
-!!   namctl            Control prints                                   (default: OFF)
-!!   namsto            Stochastic parametrization of EOS                (default: OFF)
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi")
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namctl        !   Control prints                                       (default: OFF)
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namsto        ! Stochastic parametrization of EOS                      (default: OFF)
-!-----------------------------------------------------------------------
-/
Index: MO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/namelist_ecmwf_cfg
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/namelist_ecmwf_cfg	(revision 13675)
+++ 	(revision )
@@ -1,274 +1,0 @@
-!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-!! NEMO/OCE  Configuration namelist : overwrite default values defined in SHARED/namelist_ref
-!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-!!
-!!             STATION_ASF configuration, IDEALIZED forcing
-!!                   ==> meant to be used with forcing from file:
-!!                       "input_output_VALIDATION_IDEALIZED.nc"
-!!
-!!======================================================================
-!!              ***  Domain & Run management namelists  ***           !!
-!!                                                                    !!
-!!   namrun       parameters of the run
-!!   namdom       space and time domain
-!!   namcfg       parameters of the configuration                       (default: user defined GYRE)
-!!   namwad       Wetting and drying                                    (default: OFF)
-!!   namtsd       data: temperature & salinity                          (default: OFF)
-!!   namcrs       coarsened grid (for outputs and/or TOP)               (ln_crs =T)
-!!   namc1d       1D configuration options                              ("key_c1d")
-!!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d")
-!!   namc1d_uvd   1D data (currents)                                    ("key_c1d")
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&namusr_def    !   STATION_ASF user defined namelist
-!-----------------------------------------------------------------------
-   rn_dept1    =  1.0   ! depth (m) at which the prescribed SST is taken  == depth of first T point! (IMPORTANT for warm-layer impact, the deeper, the more pronounced!)
-   !                    ! [PAPA station = 1 m]
-/
-!-----------------------------------------------------------------------
-&namrun        !   parameters of the run
-!-----------------------------------------------------------------------
-   cn_exp      =  'STATION_ASF-ECMWF-IDEALIZED'  !  experience name
-   nn_it000    =    1   !  first time step
-   nn_itend    =       104  !  last  time step => 3 years (including 1 leap!) at dt=3600s
-   nn_date0    = 20200101  !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
-   nn_leapy    =       1   !  Leap year calendar (1) or not (0)
-   !
-   nn_stock    =    104
-   nn_write    =    104
-/
-!-----------------------------------------------------------------------
-&namdom        !   time and space domain / STATION_ASF
-!-----------------------------------------------------------------------
-   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time
-   rn_Dt       =  3600.    !  time step for the dynamics and tracer => frequency of atmospheric forcing
-   ln_meshmask = .false.
-/
-!-----------------------------------------------------------------------
-&namcfg        !   parameters of the configuration                      (default: use namusr_def in namelist_cfg)
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namc1d        !   1D configuration options / STATION_ASF               ("key_c1d" default: PAPA station)
-!-----------------------------------------------------------------------
-   rn_lat1d    =      50.1 !  Column latitude    [PAPA station]
-   rn_lon1d    =    215.1  !  Column longitude   [PAPA station]
-/
-!-----------------------------------------------------------------------
-&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d" default: OFF)
-!-----------------------------------------------------------------------
-   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F)
-/
-!-----------------------------------------------------------------------
-&namc1d_uvd    !   data: U & V currents                                 ("key_c1d" default: OFF)
-!-----------------------------------------------------------------------
-   !                       !  =T read U-V fields for:
-   ln_uvd_init   = .false.       !  ocean initialisation
-   ln_uvd_dyndmp = .false.       !  U-V restoring
-
-   cn_dir      = './'      !  root directory for the U-V data location
-   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
-   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
-   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
-   sn_ucur     = 'ucurrent'              ,         -1.       ,'u_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Ume'   , ''
-   sn_vcur     = 'vcurrent'              ,         -1.       ,'v_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Vme'   , ''
-/
-
-!!======================================================================
-!!            ***  Surface Boundary Condition namelists  ***          !!
-!!                                                                    !!
-!!                           STATION_ASF test-case                    !!
-!!                                                                    !!
-!!   namsbc          surface boundary condition manager                 (default: NO selection)
-!!   namsbc_flx      flux               formulation                     (ln_flx     =T)
-!!   namsbc_blk      Bulk formulae formulation                          (ln_blk     =T)
-!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3" )
-!!   namsbc_sas      Stand-Alone Surface module                         (SAS_SRC  only)
-!!   namsbc_iif      Ice-IF: use observed ice cover                     (nn_ice = 1   )
-!!   namtra_qsr      penetrative solar radiation                        (ln_traqsr  =T)
-!!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T)
-!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T)
-!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T)
-!!   namsbc_wave     external fields from wave model                    (ln_wave    =T)
-!!   namberg         iceberg floats                                     (ln_icebergs=T)
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&namsbc        !   Surface Boundary Condition manager                   (default: NO selection)
-!-----------------------------------------------------------------------
-   nn_fsbc     = 1         !!LOLO!  frequency of SBC module call
-   ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk )
-   !
-   nn_ice      = 0  !LOLO  !  =0 no ice boundary condition
-      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif )
-      !                    !  =2 or 3 for SI3 and CICE, respectively
-   !    
-   ln_dm2dc    = .false.   !#lolo  daily mean to diurnal cycle on short wave
-/
-!-----------------------------------------------------------------------
-&namsbc_blk    !   namsbc_blk  generic Bulk formula / STATION_ASF       (ln_blk =T)
-!-----------------------------------------------------------------------
-   !                    !  bulk algorithm :
-   ln_NCAR      = .false.    ! "NCAR"      algorithm   (Large and Yeager 2008)
-   ln_COARE_3p0 = .false.    ! "COARE 3.0" algorithm   (Fairall et al. 2003)
-   ln_COARE_3p6 = .false.    ! "COARE 3.6" algorithm   (Edson et al. 2013)
-   ln_ECMWF     = .true.     ! "ECMWF"     algorithm   (IFS cycle 45r1)
-   ln_ANDREAS   = .false.    ! "ANDREAS"   algorithm   (Andreas et al. 2015)
-      rn_zqt       =  2.     !  Air temperature & humidity reference height (m)
-      rn_zu        = 10.     !  Wind vector reference height (m)
-      nn_iter_algo = 20      !  Number of iterations in bulk param. algo ("stable ABL + weak wind" requires more)
-      ln_skin_cs   = .true.  !  use the cool-skin parameterization
-      ln_skin_wl   = .true.  !  use the warm-layer parameterization
-   !
-   ln_humi_sph = .false. !  humidity "sn_humi" is specific humidity  [kg/kg]
-   ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K]
-   ln_humi_rlh = .true.  !  humidity "sn_humi" is relative humidity     [%]
-   ln_tpot     = .true.  !!GS: compute potential temperature or not
-   !
-   cn_dir      = './'      !  root directory for the bulk data location
-   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________!
-   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !       weights filename               ! rotation ! land/sea mask !
-   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                                      ! pairing  !    filename   !
-   sn_wndi     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'u10'    ,   .false.   , .true.  , 'yearly'  , ''   ,''   , ''
-   sn_wndj     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'v10'    ,   .false.   , .true.  , 'yearly'  , ''   ,''   , ''
-   sn_qsr      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rad_sw' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_qlw      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rad_lw' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_tair     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  't_air'  ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_humi     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rh_air' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_prec     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'precip' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_snow     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'snow'   ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_slp      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'slp'    ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-/
-!-----------------------------------------------------------------------
-&namsbc_abl    !   Atmospheric Boundary Layer formulation           (ln_abl = T)
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3")
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namsbc_sas    !   Stand-Alone Surface boundary condition / STATION_ASF
-!-----------------------------------------------------------------------
-   l_sasread   = .true.    !  =T Read in file ;  =F set all to 0. (see sbcssm)
-   !
-   cn_dir      = './'      !  root directory for the ocean data location
-   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
-   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
-   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
-   sn_usp  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssu'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   sn_vsp  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssv'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   sn_tem  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'sst'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   sn_sal  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'sss',        .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   sn_ssh  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssh',        .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   !!
-/
-!!======================================================================
-!!               ***  Lateral boundary condition  ***                 !!
-!!                                                                    !!
-!!   namlbc        lateral momentum boundary condition                  (default: NO selection)
-!!   namagrif      agrif nested grid   (read by child model only)       ("key_agrif")
-!!   nam_tide      Tidal forcing                                        (default: OFF)
-!!   nambdy        Unstructured open boundaries                         (default: OFF)
-!!   nambdy_dta    Unstructured open boundaries - external data         (see  nambdy)
-!!   nambdy_tide   tidal forcing at open boundaries                     (default: OFF)
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&namlbc        !   lateral momentum boundary condition                  (default: NO selection)
-!-----------------------------------------------------------------------
-   rn_shlat    =    0.     !  free slip
-/
-!!======================================================================
-!!                ***  Top/Bottom boundary condition  ***             !!
-!!                                                                    !!
-!!   namdrg        top/bottom drag coefficient                          (default: NO selection)
-!!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T)
-!!   namdrg_bot    bottom friction                                      (ln_OFF=F)
-!!   nambbc        bottom temperature boundary condition                (default: OFF)
-!!   nambbl        bottom boundary layer scheme                         (default: OFF)
-!!======================================================================
-!
-!!======================================================================
-!!                        Tracer (T-S) namelists                      !!
-!!                                                                    !!
-!!   nameos        equation of state                                    (default: NO selection)
-!!   namtra_adv    advection scheme                                     (default: NO selection)
-!!   namtra_ldf    lateral diffusion scheme                             (default: NO selection)
-!!   namtra_mle    mixed layer eddy param. (Fox-Kemper param.)          (default: OFF)
-!!   namtra_eiv    eddy induced velocity param.                         (default: OFF)
-!!   namtra_dmp    T & S newtonian damping                              (default: OFF)
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&nameos        !   ocean Equation Of Seawater                           (default: NO selection)
-!-----------------------------------------------------------------------
-   ln_eos80    = .true.         !  = Use EOS80
-/
-!!======================================================================
-!!                      ***  Dynamics namelists  ***                  !!
-!!                                                                    !!
-!!   nam_vvl       vertical coordinate options                          (default: z-star)
-!!   namdyn_adv    formulation of the momentum advection                (default: NO selection)
-!!   namdyn_vor    advection scheme                                     (default: NO selection)
-!!   namdyn_hpg    hydrostatic pressure gradient                        (default: NO selection)
-!!   namdyn_spg    surface pressure gradient                            (default: NO selection)
-!!   namdyn_ldf    lateral diffusion scheme                             (default: NO selection)
-!!   namdta_dyn    offline TOP: dynamics read in files                  (OFF_SRC only)
-!!======================================================================
-!
-!!======================================================================
-!!                     vertical physics namelists                     !!
-!!                                                                    !!
-!!    namzdf        vertical physics manager                            (default: NO selection)
-!!    namzdf_ric    richardson number vertical mixing                   (ln_zdfric=T)
-!!    namzdf_tke    TKE vertical mixing                                 (ln_zdftke=T)
-!!    namzdf_gls    GLS vertical mixing                                 (ln_zdfgls=T)
-!!    namzdf_osm    OSM vertical diffusion                              (ln_zdfosm=T)
-!!    namzdf_iwm    tidal mixing parameterization                       (ln_zdfiwm=T)
-!!======================================================================
-!
-!!======================================================================
-!!                  ***  Diagnostics namelists  ***                   !!
-!!                                                                    !!
-!!   namtrd       dynamics and/or tracer trends                         (default: OFF)
-!!   namptr       Poleward Transport Diagnostics                        (default: OFF)
-!!   namhsb       Heat and salt budgets                                 (default: OFF)
-!!   namdiu       Cool skin and warm layer models                       (default: OFF)
-!!   namdiu       Cool skin and warm layer models                       (default: OFF)
-!!   namflo       float parameters                                      (default: OFF)
-!!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF)
-!!   nam_diadct   transports through some sections                      (default: OFF)
-!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF)
-!!   nam_dia25h   25h Mean Output                                       (default: OFF)
-!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4")
-!!======================================================================
-!
-!!======================================================================
-!!               ***  Observation & Assimilation  ***                 !!
-!!                                                                    !!
-!!   namobs       observation and model comparison                      (default: OFF)
-!!   nam_asminc   assimilation increments                               ('key_asminc')
-!!======================================================================
-!
-!!======================================================================
-!!                  ***  Miscellaneous namelists  ***                 !!
-!!                                                                    !!
-!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi")
-!!   namctl            Control prints                                   (default: OFF)
-!!   namsto            Stochastic parametrization of EOS                (default: OFF)
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi")
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namctl        !   Control prints                                       (default: OFF)
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namsto        ! Stochastic parametrization of EOS                      (default: OFF)
-!-----------------------------------------------------------------------
-/
Index: MO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/namelist_ncar_cfg
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/namelist_ncar_cfg	(revision 13675)
+++ 	(revision )
@@ -1,274 +1,0 @@
-!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-!! NEMO/OCE  Configuration namelist : overwrite default values defined in SHARED/namelist_ref
-!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-!!
-!!             STATION_ASF configuration, IDEALIZED forcing
-!!                   ==> meant to be used with forcing from file:
-!!                       "input_output_VALIDATION_IDEALIZED.nc"
-!!
-!!======================================================================
-!!              ***  Domain & Run management namelists  ***           !!
-!!                                                                    !!
-!!   namrun       parameters of the run
-!!   namdom       space and time domain
-!!   namcfg       parameters of the configuration                       (default: user defined GYRE)
-!!   namwad       Wetting and drying                                    (default: OFF)
-!!   namtsd       data: temperature & salinity                          (default: OFF)
-!!   namcrs       coarsened grid (for outputs and/or TOP)               (ln_crs =T)
-!!   namc1d       1D configuration options                              ("key_c1d")
-!!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d")
-!!   namc1d_uvd   1D data (currents)                                    ("key_c1d")
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&namusr_def    !   STATION_ASF user defined namelist
-!-----------------------------------------------------------------------
-   rn_dept1    =  1.0   ! depth (m) at which the prescribed SST is taken  == depth of first T point! (IMPORTANT for warm-layer impact, the deeper, the more pronounced!)
-   !                    ! [PAPA station = 1 m]
-/
-!-----------------------------------------------------------------------
-&namrun        !   parameters of the run
-!-----------------------------------------------------------------------
-   cn_exp      =  'STATION_ASF-NCAR-IDEALIZED'  !  experience name
-   nn_it000    =    1   !  first time step
-   nn_itend    =       104  !  last  time step => 3 years (including 1 leap!) at dt=3600s
-   nn_date0    = 20200101  !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
-   nn_leapy    =       1   !  Leap year calendar (1) or not (0)
-   !
-   nn_stock    =    104
-   nn_write    =    104
-/
-!-----------------------------------------------------------------------
-&namdom        !   time and space domain / STATION_ASF
-!-----------------------------------------------------------------------
-   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time
-   rn_Dt       =  3600.    !  time step for the dynamics and tracer => frequency of atmospheric forcing
-   ln_meshmask = .false.
-/
-!-----------------------------------------------------------------------
-&namcfg        !   parameters of the configuration                      (default: use namusr_def in namelist_cfg)
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namc1d        !   1D configuration options / STATION_ASF               ("key_c1d" default: PAPA station)
-!-----------------------------------------------------------------------
-   rn_lat1d    =      50.1 !  Column latitude    [PAPA station]
-   rn_lon1d    =    215.1  !  Column longitude   [PAPA station]
-/
-!-----------------------------------------------------------------------
-&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d" default: OFF)
-!-----------------------------------------------------------------------
-   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F)
-/
-!-----------------------------------------------------------------------
-&namc1d_uvd    !   data: U & V currents                                 ("key_c1d" default: OFF)
-!-----------------------------------------------------------------------
-   !                       !  =T read U-V fields for:
-   ln_uvd_init   = .false.       !  ocean initialisation
-   ln_uvd_dyndmp = .false.       !  U-V restoring
-
-   cn_dir      = './'      !  root directory for the U-V data location
-   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
-   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
-   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
-   sn_ucur     = 'ucurrent'              ,         -1.       ,'u_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Ume'   , ''
-   sn_vcur     = 'vcurrent'              ,         -1.       ,'v_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Vme'   , ''
-/
-
-!!======================================================================
-!!            ***  Surface Boundary Condition namelists  ***          !!
-!!                                                                    !!
-!!                           STATION_ASF test-case                         !!
-!!                                                                    !!
-!!   namsbc          surface boundary condition manager                 (default: NO selection)
-!!   namsbc_flx      flux               formulation                     (ln_flx     =T)
-!!   namsbc_blk      Bulk formulae formulation                          (ln_blk     =T)
-!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3" )
-!!   namsbc_sas      Stand-Alone Surface module                         (SAS_SRC  only)
-!!   namsbc_iif      Ice-IF: use observed ice cover                     (nn_ice = 1   )
-!!   namtra_qsr      penetrative solar radiation                        (ln_traqsr  =T)
-!!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T)
-!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T)
-!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T)
-!!   namsbc_wave     external fields from wave model                    (ln_wave    =T)
-!!   namberg         iceberg floats                                     (ln_icebergs=T)
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&namsbc        !   Surface Boundary Condition manager                   (default: NO selection)
-!-----------------------------------------------------------------------
-   nn_fsbc     = 1         !!LOLO!  frequency of SBC module call
-   ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk )
-   !
-   nn_ice      = 0  !LOLO  !  =0 no ice boundary condition
-      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif )
-      !                    !  =2 or 3 for SI3 and CICE, respectively
-   !    
-   ln_dm2dc    = .false.   !#lolo  daily mean to diurnal cycle on short wave
-/
-!-----------------------------------------------------------------------
-&namsbc_blk    !   namsbc_blk  generic Bulk formula / STATION_ASF       (ln_blk =T)
-!-----------------------------------------------------------------------
-   !                    !  bulk algorithm :
-   ln_NCAR      = .true.     ! "NCAR"      algorithm   (Large and Yeager 2008)
-   ln_COARE_3p0 = .false.    ! "COARE 3.0" algorithm   (Fairall et al. 2003)
-   ln_COARE_3p6 = .false.    ! "COARE 3.6" algorithm   (Edson et al. 2013)
-   ln_ECMWF     = .false.    ! "ECMWF"     algorithm   (IFS cycle 45r1)
-   ln_ANDREAS   = .false.    ! "ANDREAS"   algorithm   (Andreas et al. 2015)
-      rn_zqt       =  2.     !  Air temperature & humidity reference height (m)
-      rn_zu        = 10.     !  Wind vector reference height (m)
-      nn_iter_algo = 20      !  Number of iterations in bulk param. algo ("stable ABL + weak wind" requires more)
-      ln_skin_cs   = .false. !  use the cool-skin parameterization
-      ln_skin_wl   = .false. !  use the warm-layer parameterization
-   !
-   ln_humi_sph = .false. !  humidity "sn_humi" is specific humidity  [kg/kg]
-   ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K]
-   ln_humi_rlh = .true.  !  humidity "sn_humi" is relative humidity     [%]
-   ln_tpot     = .true.  !!GS: compute potential temperature or not
-   !
-   cn_dir      = './'      !  root directory for the bulk data location
-   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________!
-   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !       weights filename               ! rotation ! land/sea mask !
-   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                                      ! pairing  !    filename   !
-   sn_wndi     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'u10'    ,   .false.   , .true.  , 'yearly'  , ''   ,''   , ''
-   sn_wndj     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'v10'    ,   .false.   , .true.  , 'yearly'  , ''   ,''   , ''
-   sn_qsr      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rad_sw' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_qlw      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rad_lw' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_tair     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  't_air'  ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_humi     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rh_air' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_prec     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'precip' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_snow     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'snow'   ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-   sn_slp      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'slp'    ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
-/
-!-----------------------------------------------------------------------
-&namsbc_abl    !   Atmospheric Boundary Layer formulation           (ln_abl = T)
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3")
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namsbc_sas    !   Stand-Alone Surface boundary condition / STATION_ASF
-!-----------------------------------------------------------------------
-   l_sasread   = .true.    !  =T Read in file ;  =F set all to 0. (see sbcssm)
-   !
-   cn_dir      = './'      !  root directory for the ocean data location
-   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
-   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
-   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
-   sn_usp  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssu'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   sn_vsp  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssv'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   sn_tem  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'sst'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   sn_sal  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'sss',        .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   sn_ssh  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssh',        .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
-   !!
-/
-!!======================================================================
-!!               ***  Lateral boundary condition  ***                 !!
-!!                                                                    !!
-!!   namlbc        lateral momentum boundary condition                  (default: NO selection)
-!!   namagrif      agrif nested grid   (read by child model only)       ("key_agrif")
-!!   nam_tide      Tidal forcing                                        (default: OFF)
-!!   nambdy        Unstructured open boundaries                         (default: OFF)
-!!   nambdy_dta    Unstructured open boundaries - external data         (see  nambdy)
-!!   nambdy_tide   tidal forcing at open boundaries                     (default: OFF)
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&namlbc        !   lateral momentum boundary condition                  (default: NO selection)
-!-----------------------------------------------------------------------
-   rn_shlat    =    0.     !  free slip
-/
-!!======================================================================
-!!                ***  Top/Bottom boundary condition  ***             !!
-!!                                                                    !!
-!!   namdrg        top/bottom drag coefficient                          (default: NO selection)
-!!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T)
-!!   namdrg_bot    bottom friction                                      (ln_OFF=F)
-!!   nambbc        bottom temperature boundary condition                (default: OFF)
-!!   nambbl        bottom boundary layer scheme                         (default: OFF)
-!!======================================================================
-!
-!!======================================================================
-!!                        Tracer (T-S) namelists                      !!
-!!                                                                    !!
-!!   nameos        equation of state                                    (default: NO selection)
-!!   namtra_adv    advection scheme                                     (default: NO selection)
-!!   namtra_ldf    lateral diffusion scheme                             (default: NO selection)
-!!   namtra_mle    mixed layer eddy param. (Fox-Kemper param.)          (default: OFF)
-!!   namtra_eiv    eddy induced velocity param.                         (default: OFF)
-!!   namtra_dmp    T & S newtonian damping                              (default: OFF)
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&nameos        !   ocean Equation Of Seawater                           (default: NO selection)
-!-----------------------------------------------------------------------
-   ln_eos80    = .true.         !  = Use EOS80
-/
-!!======================================================================
-!!                      ***  Dynamics namelists  ***                  !!
-!!                                                                    !!
-!!   nam_vvl       vertical coordinate options                          (default: z-star)
-!!   namdyn_adv    formulation of the momentum advection                (default: NO selection)
-!!   namdyn_vor    advection scheme                                     (default: NO selection)
-!!   namdyn_hpg    hydrostatic pressure gradient                        (default: NO selection)
-!!   namdyn_spg    surface pressure gradient                            (default: NO selection)
-!!   namdyn_ldf    lateral diffusion scheme                             (default: NO selection)
-!!   namdta_dyn    offline TOP: dynamics read in files                  (OFF_SRC only)
-!!======================================================================
-!
-!!======================================================================
-!!                     vertical physics namelists                     !!
-!!                                                                    !!
-!!    namzdf        vertical physics manager                            (default: NO selection)
-!!    namzdf_ric    richardson number vertical mixing                   (ln_zdfric=T)
-!!    namzdf_tke    TKE vertical mixing                                 (ln_zdftke=T)
-!!    namzdf_gls    GLS vertical mixing                                 (ln_zdfgls=T)
-!!    namzdf_osm    OSM vertical diffusion                              (ln_zdfosm=T)
-!!    namzdf_iwm    tidal mixing parameterization                       (ln_zdfiwm=T)
-!!======================================================================
-!
-!!======================================================================
-!!                  ***  Diagnostics namelists  ***                   !!
-!!                                                                    !!
-!!   namtrd       dynamics and/or tracer trends                         (default: OFF)
-!!   namptr       Poleward Transport Diagnostics                        (default: OFF)
-!!   namhsb       Heat and salt budgets                                 (default: OFF)
-!!   namdiu       Cool skin and warm layer models                       (default: OFF)
-!!   namdiu       Cool skin and warm layer models                       (default: OFF)
-!!   namflo       float parameters                                      (default: OFF)
-!!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF)
-!!   nam_diadct   transports through some sections                      (default: OFF)
-!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF)
-!!   nam_dia25h   25h Mean Output                                       (default: OFF)
-!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4")
-!!======================================================================
-!
-!!======================================================================
-!!               ***  Observation & Assimilation  ***                 !!
-!!                                                                    !!
-!!   namobs       observation and model comparison                      (default: OFF)
-!!   nam_asminc   assimilation increments                               ('key_asminc')
-!!======================================================================
-!
-!!======================================================================
-!!                  ***  Miscellaneous namelists  ***                 !!
-!!                                                                    !!
-!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi")
-!!   namctl            Control prints                                   (default: OFF)
-!!   namsto            Stochastic parametrization of EOS                (default: OFF)
-!!======================================================================
-!
-!-----------------------------------------------------------------------
-&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi")
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namctl        !   Control prints                                       (default: OFF)
-!-----------------------------------------------------------------------
-/
-!-----------------------------------------------------------------------
-&namsto        ! Stochastic parametrization of EOS                      (default: OFF)
-!-----------------------------------------------------------------------
-/
Index: /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/oce/namelist_andreas_cfg
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/oce/namelist_andreas_cfg	(revision 13676)
+++ /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/oce/namelist_andreas_cfg	(revision 13676)
@@ -0,0 +1,274 @@
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!! NEMO/OCE  Configuration namelist : overwrite default values defined in SHARED/namelist_ref
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!!
+!!             STATION_ASF configuration, IDEALIZED forcing
+!!                   ==> meant to be used with forcing from file:
+!!                       "input_output_VALIDATION_IDEALIZED.nc"
+!!
+!!======================================================================
+!!              ***  Domain & Run management namelists  ***           !!
+!!                                                                    !!
+!!   namrun       parameters of the run
+!!   namdom       space and time domain
+!!   namcfg       parameters of the configuration                       (default: user defined GYRE)
+!!   namwad       Wetting and drying                                    (default: OFF)
+!!   namtsd       data: temperature & salinity                          (default: OFF)
+!!   namcrs       coarsened grid (for outputs and/or TOP)               (ln_crs =T)
+!!   namc1d       1D configuration options                              ("key_c1d")
+!!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d")
+!!   namc1d_uvd   1D data (currents)                                    ("key_c1d")
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namusr_def    !   STATION_ASF user defined namelist
+!-----------------------------------------------------------------------
+   rn_dept1    =  1.0   ! depth (m) at which the prescribed SST is taken  == depth of first T point! (IMPORTANT for warm-layer impact, the deeper, the more pronounced!)
+   !                    ! [PAPA station = 1 m]
+/
+!-----------------------------------------------------------------------
+&namrun        !   parameters of the run
+!-----------------------------------------------------------------------
+   cn_exp      =  'STATION_ASF-ANDREAS_IDEALIZED'  !  experience name
+   nn_it000    =    1   !  first time step
+   nn_itend    =       104  !  last  time step => 3 years (including 1 leap!) at dt=3600s
+   nn_date0    = 20200101  !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
+   nn_leapy    =       1   !  Leap year calendar (1) or not (0)
+   !
+   nn_stock    =    104
+   nn_write    =    104
+/
+!-----------------------------------------------------------------------
+&namdom        !   time and space domain / STATION_ASF
+!-----------------------------------------------------------------------
+   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time
+   rn_Dt       =  3600.    !  time step for the dynamics and tracer => frequency of atmospheric forcing
+   ln_meshmask = .false.
+/
+!-----------------------------------------------------------------------
+&namcfg        !   parameters of the configuration                      (default: use namusr_def in namelist_cfg)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namc1d        !   1D configuration options / STATION_ASF               ("key_c1d" default: PAPA station)
+!-----------------------------------------------------------------------
+   rn_lat1d    =      50.1 !  Column latitude    [PAPA station]
+   rn_lon1d    =    215.1  !  Column longitude   [PAPA station]
+/
+!-----------------------------------------------------------------------
+&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d" default: OFF)
+!-----------------------------------------------------------------------
+   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F)
+/
+!-----------------------------------------------------------------------
+&namc1d_uvd    !   data: U & V currents                                 ("key_c1d" default: OFF)
+!-----------------------------------------------------------------------
+   !                       !  =T read U-V fields for:
+   ln_uvd_init   = .false.       !  ocean initialisation
+   ln_uvd_dyndmp = .false.       !  U-V restoring
+
+   cn_dir      = './'      !  root directory for the U-V data location
+   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
+   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
+   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
+   sn_ucur     = 'ucurrent'              ,         -1.       ,'u_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Ume'   , ''
+   sn_vcur     = 'vcurrent'              ,         -1.       ,'v_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Vme'   , ''
+/
+
+!!======================================================================
+!!            ***  Surface Boundary Condition namelists  ***          !!
+!!                                                                    !!
+!!                           STATION_ASF test-case                         !!
+!!                                                                    !!
+!!   namsbc          surface boundary condition manager                 (default: NO selection)
+!!   namsbc_flx      flux               formulation                     (ln_flx     =T)
+!!   namsbc_blk      Bulk formulae formulation                          (ln_blk     =T)
+!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3" )
+!!   namsbc_sas      Stand-Alone Surface module                         (SAS_SRC  only)
+!!   namsbc_iif      Ice-IF: use observed ice cover                     (nn_ice = 1   )
+!!   namtra_qsr      penetrative solar radiation                        (ln_traqsr  =T)
+!!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T)
+!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T)
+!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T)
+!!   namsbc_wave     external fields from wave model                    (ln_wave    =T)
+!!   namberg         iceberg floats                                     (ln_icebergs=T)
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namsbc        !   Surface Boundary Condition manager                   (default: NO selection)
+!-----------------------------------------------------------------------
+   nn_fsbc     = 1         !!LOLO!  frequency of SBC module call
+   ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk )
+   !
+   nn_ice      = 0  !LOLO  !  =0 no ice boundary condition
+      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif )
+      !                    !  =2 or 3 for SI3 and CICE, respectively
+   !    
+   ln_dm2dc    = .false.   !#lolo  daily mean to diurnal cycle on short wave
+/
+!-----------------------------------------------------------------------
+&namsbc_blk    !   namsbc_blk  generic Bulk formula / STATION_ASF       (ln_blk =T)
+!-----------------------------------------------------------------------
+   !                    !  bulk algorithm :
+   ln_NCAR      = .false.    ! "NCAR"      algorithm   (Large and Yeager 2008)
+   ln_COARE_3p0 = .false.    ! "COARE 3.0" algorithm   (Fairall et al. 2003)
+   ln_COARE_3p6 = .false.    ! "COARE 3.6" algorithm   (Edson et al. 2013)
+   ln_ECMWF     = .false.    ! "ECMWF"     algorithm   (IFS cycle 45r1)
+   ln_ANDREAS   = .true.     ! "ANDREAS"   algorithm   (Andreas et al. 2015)
+      rn_zqt       =  2.     !  Air temperature & humidity reference height (m)
+      rn_zu        = 10.     !  Wind vector reference height (m)
+      nn_iter_algo = 20      !  Number of iterations in bulk param. algo ("stable ABL + weak wind" requires more)
+      ln_skin_cs   = .false. !  use the cool-skin parameterization
+      ln_skin_wl   = .false. !  use the warm-layer parameterization
+   !
+   ln_humi_sph = .false. !  humidity "sn_humi" is specific humidity  [kg/kg]
+   ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K]
+   ln_humi_rlh = .true.  !  humidity "sn_humi" is relative humidity     [%]
+   ln_tpot     = .true.  !!GS: compute potential temperature or not
+   !
+   cn_dir      = './'      !  root directory for the bulk data location
+   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________!
+   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !       weights filename               ! rotation ! land/sea mask !
+   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                                      ! pairing  !    filename   !
+   sn_wndi     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'u10'    ,   .false.   , .true.  , 'yearly'  , ''   ,''   , ''
+   sn_wndj     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'v10'    ,   .false.   , .true.  , 'yearly'  , ''   ,''   , ''
+   sn_qsr      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rad_sw' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_qlw      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rad_lw' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_tair     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  't_air'  ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_humi     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rh_air' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_prec     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'precip' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_snow     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'snow'   ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_slp      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'slp'    ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+/
+!-----------------------------------------------------------------------
+&namsbc_abl    !   Atmospheric Boundary Layer formulation           (ln_abl = T)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3")
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namsbc_sas    !   Stand-Alone Surface boundary condition / STATION_ASF
+!-----------------------------------------------------------------------
+   l_sasread   = .true.    !  =T Read in file ;  =F set all to 0. (see sbcssm)
+   !
+   cn_dir      = './'      !  root directory for the ocean data location
+   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
+   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
+   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
+   sn_usp  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssu'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   sn_vsp  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssv'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   sn_tem  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'sst'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   sn_sal  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'sss',        .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   sn_ssh  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssh',        .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   !!
+/
+!!======================================================================
+!!               ***  Lateral boundary condition  ***                 !!
+!!                                                                    !!
+!!   namlbc        lateral momentum boundary condition                  (default: NO selection)
+!!   namagrif      agrif nested grid   (read by child model only)       ("key_agrif")
+!!   nam_tide      Tidal forcing                                        (default: OFF)
+!!   nambdy        Unstructured open boundaries                         (default: OFF)
+!!   nambdy_dta    Unstructured open boundaries - external data         (see  nambdy)
+!!   nambdy_tide   tidal forcing at open boundaries                     (default: OFF)
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namlbc        !   lateral momentum boundary condition                  (default: NO selection)
+!-----------------------------------------------------------------------
+   rn_shlat    =    0.     !  free slip
+/
+!!======================================================================
+!!                ***  Top/Bottom boundary condition  ***             !!
+!!                                                                    !!
+!!   namdrg        top/bottom drag coefficient                          (default: NO selection)
+!!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T)
+!!   namdrg_bot    bottom friction                                      (ln_OFF=F)
+!!   nambbc        bottom temperature boundary condition                (default: OFF)
+!!   nambbl        bottom boundary layer scheme                         (default: OFF)
+!!======================================================================
+!
+!!======================================================================
+!!                        Tracer (T-S) namelists                      !!
+!!                                                                    !!
+!!   nameos        equation of state                                    (default: NO selection)
+!!   namtra_adv    advection scheme                                     (default: NO selection)
+!!   namtra_ldf    lateral diffusion scheme                             (default: NO selection)
+!!   namtra_mle    mixed layer eddy param. (Fox-Kemper param.)          (default: OFF)
+!!   namtra_eiv    eddy induced velocity param.                         (default: OFF)
+!!   namtra_dmp    T & S newtonian damping                              (default: OFF)
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&nameos        !   ocean Equation Of Seawater                           (default: NO selection)
+!-----------------------------------------------------------------------
+   ln_eos80    = .true.         !  = Use EOS80
+/
+!!======================================================================
+!!                      ***  Dynamics namelists  ***                  !!
+!!                                                                    !!
+!!   nam_vvl       vertical coordinate options                          (default: z-star)
+!!   namdyn_adv    formulation of the momentum advection                (default: NO selection)
+!!   namdyn_vor    advection scheme                                     (default: NO selection)
+!!   namdyn_hpg    hydrostatic pressure gradient                        (default: NO selection)
+!!   namdyn_spg    surface pressure gradient                            (default: NO selection)
+!!   namdyn_ldf    lateral diffusion scheme                             (default: NO selection)
+!!   namdta_dyn    offline TOP: dynamics read in files                  (OFF_SRC only)
+!!======================================================================
+!
+!!======================================================================
+!!                     vertical physics namelists                     !!
+!!                                                                    !!
+!!    namzdf        vertical physics manager                            (default: NO selection)
+!!    namzdf_ric    richardson number vertical mixing                   (ln_zdfric=T)
+!!    namzdf_tke    TKE vertical mixing                                 (ln_zdftke=T)
+!!    namzdf_gls    GLS vertical mixing                                 (ln_zdfgls=T)
+!!    namzdf_osm    OSM vertical diffusion                              (ln_zdfosm=T)
+!!    namzdf_iwm    tidal mixing parameterization                       (ln_zdfiwm=T)
+!!======================================================================
+!
+!!======================================================================
+!!                  ***  Diagnostics namelists  ***                   !!
+!!                                                                    !!
+!!   namtrd       dynamics and/or tracer trends                         (default: OFF)
+!!   namptr       Poleward Transport Diagnostics                        (default: OFF)
+!!   namhsb       Heat and salt budgets                                 (default: OFF)
+!!   namdiu       Cool skin and warm layer models                       (default: OFF)
+!!   namdiu       Cool skin and warm layer models                       (default: OFF)
+!!   namflo       float parameters                                      (default: OFF)
+!!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF)
+!!   nam_diadct   transports through some sections                      (default: OFF)
+!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF)
+!!   nam_dia25h   25h Mean Output                                       (default: OFF)
+!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4")
+!!======================================================================
+!
+!!======================================================================
+!!               ***  Observation & Assimilation  ***                 !!
+!!                                                                    !!
+!!   namobs       observation and model comparison                      (default: OFF)
+!!   nam_asminc   assimilation increments                               ('key_asminc')
+!!======================================================================
+!
+!!======================================================================
+!!                  ***  Miscellaneous namelists  ***                 !!
+!!                                                                    !!
+!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi")
+!!   namctl            Control prints                                   (default: OFF)
+!!   namsto            Stochastic parametrization of EOS                (default: OFF)
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi")
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namctl        !   Control prints                                       (default: OFF)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namsto        ! Stochastic parametrization of EOS                      (default: OFF)
+!-----------------------------------------------------------------------
+/
Index: /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/oce/namelist_coare3p0_cfg
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/oce/namelist_coare3p0_cfg	(revision 13676)
+++ /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/oce/namelist_coare3p0_cfg	(revision 13676)
@@ -0,0 +1,274 @@
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!! NEMO/OCE  Configuration namelist : overwrite default values defined in SHARED/namelist_ref
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!!
+!!             STATION_ASF configuration, IDEALIZED forcing
+!!                   ==> meant to be used with forcing from file:
+!!                       "input_output_VALIDATION_IDEALIZED.nc"
+!!
+!!======================================================================
+!!              ***  Domain & Run management namelists  ***           !!
+!!                                                                    !!
+!!   namrun       parameters of the run
+!!   namdom       space and time domain
+!!   namcfg       parameters of the configuration                       (default: user defined GYRE)
+!!   namwad       Wetting and drying                                    (default: OFF)
+!!   namtsd       data: temperature & salinity                          (default: OFF)
+!!   namcrs       coarsened grid (for outputs and/or TOP)               (ln_crs =T)
+!!   namc1d       1D configuration options                              ("key_c1d")
+!!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d")
+!!   namc1d_uvd   1D data (currents)                                    ("key_c1d")
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namusr_def    !   STATION_ASF user defined namelist
+!-----------------------------------------------------------------------
+   rn_dept1    =  1.0   ! depth (m) at which the prescribed SST is taken  == depth of first T point! (IMPORTANT for warm-layer impact, the deeper, the more pronounced!)
+   !                    ! [PAPA station = 1 m]
+/
+!-----------------------------------------------------------------------
+&namrun        !   parameters of the run
+!-----------------------------------------------------------------------
+   cn_exp      =  'STATION_ASF-COARE3p0_IDEALIZED'  !  experience name
+   nn_it000    =    1   !  first time step
+   nn_itend    =       104  !  last  time step => 3 years (including 1 leap!) at dt=3600s
+   nn_date0    = 20200101  !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
+   nn_leapy    =       1   !  Leap year calendar (1) or not (0)
+   !
+   nn_stock    =    104
+   nn_write    =    104
+/
+!-----------------------------------------------------------------------
+&namdom        !   time and space domain / STATION_ASF
+!-----------------------------------------------------------------------
+   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time
+   rn_Dt       =  3600.    !  time step for the dynamics and tracer => frequency of atmospheric forcing
+   ln_meshmask = .false.
+/
+!-----------------------------------------------------------------------
+&namcfg        !   parameters of the configuration                      (default: use namusr_def in namelist_cfg)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namc1d        !   1D configuration options / STATION_ASF               ("key_c1d" default: PAPA station)
+!-----------------------------------------------------------------------
+   rn_lat1d    =      50.1 !  Column latitude    [PAPA station]
+   rn_lon1d    =    215.1  !  Column longitude   [PAPA station]
+/
+!-----------------------------------------------------------------------
+&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d" default: OFF)
+!-----------------------------------------------------------------------
+   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F)
+/
+!-----------------------------------------------------------------------
+&namc1d_uvd    !   data: U & V currents                                 ("key_c1d" default: OFF)
+!-----------------------------------------------------------------------
+   !                       !  =T read U-V fields for:
+   ln_uvd_init   = .false.       !  ocean initialisation
+   ln_uvd_dyndmp = .false.       !  U-V restoring
+
+   cn_dir      = './'      !  root directory for the U-V data location
+   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
+   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
+   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
+   sn_ucur     = 'ucurrent'              ,         -1.       ,'u_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Ume'   , ''
+   sn_vcur     = 'vcurrent'              ,         -1.       ,'v_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Vme'   , ''
+/
+
+!!======================================================================
+!!            ***  Surface Boundary Condition namelists  ***          !!
+!!                                                                    !!
+!!                           STATION_ASF test-case                    !!
+!!                                                                    !!
+!!   namsbc          surface boundary condition manager                 (default: NO selection)
+!!   namsbc_flx      flux               formulation                     (ln_flx     =T)
+!!   namsbc_blk      Bulk formulae formulation                          (ln_blk     =T)
+!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3" )
+!!   namsbc_sas      Stand-Alone Surface module                         (SAS_SRC  only)
+!!   namsbc_iif      Ice-IF: use observed ice cover                     (nn_ice = 1   )
+!!   namtra_qsr      penetrative solar radiation                        (ln_traqsr  =T)
+!!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T)
+!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T)
+!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T)
+!!   namsbc_wave     external fields from wave model                    (ln_wave    =T)
+!!   namberg         iceberg floats                                     (ln_icebergs=T)
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namsbc        !   Surface Boundary Condition manager                   (default: NO selection)
+!-----------------------------------------------------------------------
+   nn_fsbc     = 1         !!LOLO!  frequency of SBC module call
+   ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk )
+   !
+   nn_ice      = 0  !LOLO  !  =0 no ice boundary condition
+      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif )
+      !                    !  =2 or 3 for SI3 and CICE, respectively
+   !    
+   ln_dm2dc    = .false.   !#lolo  daily mean to diurnal cycle on short wave
+/
+!-----------------------------------------------------------------------
+&namsbc_blk    !   namsbc_blk  generic Bulk formula / STATION_ASF       (ln_blk =T)
+!-----------------------------------------------------------------------
+   !                    !  bulk algorithm :
+   ln_NCAR      = .false.    ! "NCAR"      algorithm   (Large and Yeager 2008)
+   ln_COARE_3p0 = .true.     ! "COARE 3.0" algorithm   (Fairall et al. 2003)
+   ln_COARE_3p6 = .false.    ! "COARE 3.6" algorithm   (Edson et al. 2013)
+   ln_ECMWF     = .false.    ! "ECMWF"     algorithm   (IFS cycle 45r1)
+   ln_ANDREAS   = .false.    ! "ANDREAS"   algorithm   (Andreas et al. 2015)
+      rn_zqt       =  2.     !  Air temperature & humidity reference height (m)
+      rn_zu        = 10.     !  Wind vector reference height (m)
+      nn_iter_algo = 20      !  Number of iterations in bulk param. algo ("stable ABL + weak wind" requires more)
+      ln_skin_cs   = .true.  !  use the cool-skin parameterization
+      ln_skin_wl   = .true.  !  use the warm-layer parameterization
+   !
+   ln_humi_sph = .false. !  humidity "sn_humi" is specific humidity  [kg/kg]
+   ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K]
+   ln_humi_rlh = .true.  !  humidity "sn_humi" is relative humidity     [%]
+   ln_tpot     = .true.  !!GS: compute potential temperature or not
+   !
+   cn_dir      = './'      !  root directory for the bulk data location
+   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________!
+   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !       weights filename               ! rotation ! land/sea mask !
+   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                                      ! pairing  !    filename   !
+   sn_wndi     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'u10'    ,   .false.   , .true.  , 'yearly'  , ''   ,''   , ''
+   sn_wndj     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'v10'    ,   .false.   , .true.  , 'yearly'  , ''   ,''   , ''
+   sn_qsr      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rad_sw' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_qlw      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rad_lw' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_tair     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  't_air'  ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_humi     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rh_air' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_prec     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'precip' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_snow     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'snow'   ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_slp      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'slp'    ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+/
+!-----------------------------------------------------------------------
+&namsbc_abl    !   Atmospheric Boundary Layer formulation           (ln_abl = T)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3")
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namsbc_sas    !   Stand-Alone Surface boundary condition / STATION_ASF
+!-----------------------------------------------------------------------
+   l_sasread   = .true.    !  =T Read in file ;  =F set all to 0. (see sbcssm)
+   !
+   cn_dir      = './'      !  root directory for the ocean data location
+   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
+   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
+   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
+   sn_usp  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssu'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   sn_vsp  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssv'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   sn_tem  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'sst'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   sn_sal  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'sss',        .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   sn_ssh  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssh',        .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   !!
+/
+!!======================================================================
+!!               ***  Lateral boundary condition  ***                 !!
+!!                                                                    !!
+!!   namlbc        lateral momentum boundary condition                  (default: NO selection)
+!!   namagrif      agrif nested grid   (read by child model only)       ("key_agrif")
+!!   nam_tide      Tidal forcing                                        (default: OFF)
+!!   nambdy        Unstructured open boundaries                         (default: OFF)
+!!   nambdy_dta    Unstructured open boundaries - external data         (see  nambdy)
+!!   nambdy_tide   tidal forcing at open boundaries                     (default: OFF)
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namlbc        !   lateral momentum boundary condition                  (default: NO selection)
+!-----------------------------------------------------------------------
+   rn_shlat    =    0.     !  free slip
+/
+!!======================================================================
+!!                ***  Top/Bottom boundary condition  ***             !!
+!!                                                                    !!
+!!   namdrg        top/bottom drag coefficient                          (default: NO selection)
+!!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T)
+!!   namdrg_bot    bottom friction                                      (ln_OFF=F)
+!!   nambbc        bottom temperature boundary condition                (default: OFF)
+!!   nambbl        bottom boundary layer scheme                         (default: OFF)
+!!======================================================================
+!
+!!======================================================================
+!!                        Tracer (T-S) namelists                      !!
+!!                                                                    !!
+!!   nameos        equation of state                                    (default: NO selection)
+!!   namtra_adv    advection scheme                                     (default: NO selection)
+!!   namtra_ldf    lateral diffusion scheme                             (default: NO selection)
+!!   namtra_mle    mixed layer eddy param. (Fox-Kemper param.)          (default: OFF)
+!!   namtra_eiv    eddy induced velocity param.                         (default: OFF)
+!!   namtra_dmp    T & S newtonian damping                              (default: OFF)
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&nameos        !   ocean Equation Of Seawater                           (default: NO selection)
+!-----------------------------------------------------------------------
+   ln_eos80    = .true.         !  = Use EOS80
+/
+!!======================================================================
+!!                      ***  Dynamics namelists  ***                  !!
+!!                                                                    !!
+!!   nam_vvl       vertical coordinate options                          (default: z-star)
+!!   namdyn_adv    formulation of the momentum advection                (default: NO selection)
+!!   namdyn_vor    advection scheme                                     (default: NO selection)
+!!   namdyn_hpg    hydrostatic pressure gradient                        (default: NO selection)
+!!   namdyn_spg    surface pressure gradient                            (default: NO selection)
+!!   namdyn_ldf    lateral diffusion scheme                             (default: NO selection)
+!!   namdta_dyn    offline TOP: dynamics read in files                  (OFF_SRC only)
+!!======================================================================
+!
+!!======================================================================
+!!                     vertical physics namelists                     !!
+!!                                                                    !!
+!!    namzdf        vertical physics manager                            (default: NO selection)
+!!    namzdf_ric    richardson number vertical mixing                   (ln_zdfric=T)
+!!    namzdf_tke    TKE vertical mixing                                 (ln_zdftke=T)
+!!    namzdf_gls    GLS vertical mixing                                 (ln_zdfgls=T)
+!!    namzdf_osm    OSM vertical diffusion                              (ln_zdfosm=T)
+!!    namzdf_iwm    tidal mixing parameterization                       (ln_zdfiwm=T)
+!!======================================================================
+!
+!!======================================================================
+!!                  ***  Diagnostics namelists  ***                   !!
+!!                                                                    !!
+!!   namtrd       dynamics and/or tracer trends                         (default: OFF)
+!!   namptr       Poleward Transport Diagnostics                        (default: OFF)
+!!   namhsb       Heat and salt budgets                                 (default: OFF)
+!!   namdiu       Cool skin and warm layer models                       (default: OFF)
+!!   namdiu       Cool skin and warm layer models                       (default: OFF)
+!!   namflo       float parameters                                      (default: OFF)
+!!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF)
+!!   nam_diadct   transports through some sections                      (default: OFF)
+!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF)
+!!   nam_dia25h   25h Mean Output                                       (default: OFF)
+!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4")
+!!======================================================================
+!
+!!======================================================================
+!!               ***  Observation & Assimilation  ***                 !!
+!!                                                                    !!
+!!   namobs       observation and model comparison                      (default: OFF)
+!!   nam_asminc   assimilation increments                               ('key_asminc')
+!!======================================================================
+!
+!!======================================================================
+!!                  ***  Miscellaneous namelists  ***                 !!
+!!                                                                    !!
+!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi")
+!!   namctl            Control prints                                   (default: OFF)
+!!   namsto            Stochastic parametrization of EOS                (default: OFF)
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi")
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namctl        !   Control prints                                       (default: OFF)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namsto        ! Stochastic parametrization of EOS                      (default: OFF)
+!-----------------------------------------------------------------------
+/
Index: /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/oce/namelist_coare3p6_cfg
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/oce/namelist_coare3p6_cfg	(revision 13676)
+++ /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/oce/namelist_coare3p6_cfg	(revision 13676)
@@ -0,0 +1,274 @@
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!! NEMO/OCE  Configuration namelist : overwrite default values defined in SHARED/namelist_ref
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!!
+!!             STATION_ASF configuration, IDEALIZED forcing
+!!                   ==> meant to be used with forcing from file:
+!!                       "input_output_VALIDATION_IDEALIZED.nc"
+!!
+!!======================================================================
+!!              ***  Domain & Run management namelists  ***           !!
+!!                                                                    !!
+!!   namrun       parameters of the run
+!!   namdom       space and time domain
+!!   namcfg       parameters of the configuration                       (default: user defined GYRE)
+!!   namwad       Wetting and drying                                    (default: OFF)
+!!   namtsd       data: temperature & salinity                          (default: OFF)
+!!   namcrs       coarsened grid (for outputs and/or TOP)               (ln_crs =T)
+!!   namc1d       1D configuration options                              ("key_c1d")
+!!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d")
+!!   namc1d_uvd   1D data (currents)                                    ("key_c1d")
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namusr_def    !   STATION_ASF user defined namelist
+!-----------------------------------------------------------------------
+   rn_dept1    =  1.0   ! depth (m) at which the prescribed SST is taken  == depth of first T point! (IMPORTANT for warm-layer impact, the deeper, the more pronounced!)
+   !                    ! [PAPA station = 1 m]
+/
+!-----------------------------------------------------------------------
+&namrun        !   parameters of the run
+!-----------------------------------------------------------------------
+   cn_exp      =  'STATION_ASF-COARE3p6_IDEALIZED'  !  experience name
+   nn_it000    =    1   !  first time step
+   nn_itend    =       104  !  last  time step => 3 years (including 1 leap!) at dt=3600s
+   nn_date0    = 20200101  !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
+   nn_leapy    =       1   !  Leap year calendar (1) or not (0)
+   !
+   nn_stock    =    104
+   nn_write    =    104
+/
+!-----------------------------------------------------------------------
+&namdom        !   time and space domain / STATION_ASF
+!-----------------------------------------------------------------------
+   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time
+   rn_Dt       =  3600.    !  time step for the dynamics and tracer => frequency of atmospheric forcing
+   ln_meshmask = .false.
+/
+!-----------------------------------------------------------------------
+&namcfg        !   parameters of the configuration                      (default: use namusr_def in namelist_cfg)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namc1d        !   1D configuration options / STATION_ASF               ("key_c1d" default: PAPA station)
+!-----------------------------------------------------------------------
+   rn_lat1d    =      50.1 !  Column latitude    [PAPA station]
+   rn_lon1d    =    215.1  !  Column longitude   [PAPA station]
+/
+!-----------------------------------------------------------------------
+&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d" default: OFF)
+!-----------------------------------------------------------------------
+   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F)
+/
+!-----------------------------------------------------------------------
+&namc1d_uvd    !   data: U & V currents                                 ("key_c1d" default: OFF)
+!-----------------------------------------------------------------------
+   !                       !  =T read U-V fields for:
+   ln_uvd_init   = .false.       !  ocean initialisation
+   ln_uvd_dyndmp = .false.       !  U-V restoring
+
+   cn_dir      = './'      !  root directory for the U-V data location
+   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
+   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
+   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
+   sn_ucur     = 'ucurrent'              ,         -1.       ,'u_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Ume'   , ''
+   sn_vcur     = 'vcurrent'              ,         -1.       ,'v_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Vme'   , ''
+/
+
+!!======================================================================
+!!            ***  Surface Boundary Condition namelists  ***          !!
+!!                                                                    !!
+!!                           STATION_ASF test-case                    !!
+!!                                                                    !!
+!!   namsbc          surface boundary condition manager                 (default: NO selection)
+!!   namsbc_flx      flux               formulation                     (ln_flx     =T)
+!!   namsbc_blk      Bulk formulae formulation                          (ln_blk     =T)
+!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3" )
+!!   namsbc_sas      Stand-Alone Surface module                         (SAS_SRC  only)
+!!   namsbc_iif      Ice-IF: use observed ice cover                     (nn_ice = 1   )
+!!   namtra_qsr      penetrative solar radiation                        (ln_traqsr  =T)
+!!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T)
+!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T)
+!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T)
+!!   namsbc_wave     external fields from wave model                    (ln_wave    =T)
+!!   namberg         iceberg floats                                     (ln_icebergs=T)
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namsbc        !   Surface Boundary Condition manager                   (default: NO selection)
+!-----------------------------------------------------------------------
+   nn_fsbc     = 1         !!LOLO!  frequency of SBC module call
+   ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk )
+   !
+   nn_ice      = 0  !LOLO  !  =0 no ice boundary condition
+      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif )
+      !                    !  =2 or 3 for SI3 and CICE, respectively
+   !    
+   ln_dm2dc    = .false.   !#lolo  daily mean to diurnal cycle on short wave
+/
+!-----------------------------------------------------------------------
+&namsbc_blk    !   namsbc_blk  generic Bulk formula / STATION_ASF       (ln_blk =T)
+!-----------------------------------------------------------------------
+   !                    !  bulk algorithm :
+   ln_NCAR      = .false.    ! "NCAR"      algorithm   (Large and Yeager 2008)
+   ln_COARE_3p0 = .false.    ! "COARE 3.0" algorithm   (Fairall et al. 2003)
+   ln_COARE_3p6 = .true.     ! "COARE 3.6" algorithm   (Edson et al. 2013)
+   ln_ECMWF     = .false.    ! "ECMWF"     algorithm   (IFS cycle 45r1)
+   ln_ANDREAS   = .false.    ! "ANDREAS"   algorithm   (Andreas et al. 2015)
+      rn_zqt       =  2.     !  Air temperature & humidity reference height (m)
+      rn_zu        = 10.     !  Wind vector reference height (m)
+      nn_iter_algo = 20      !  Number of iterations in bulk param. algo ("stable ABL + weak wind" requires more)
+      ln_skin_cs   = .true.  !  use the cool-skin parameterization
+      ln_skin_wl   = .true.  !  use the warm-layer parameterization
+   !
+   ln_humi_sph = .false. !  humidity "sn_humi" is specific humidity  [kg/kg]
+   ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K]
+   ln_humi_rlh = .true.  !  humidity "sn_humi" is relative humidity     [%]
+   ln_tpot     = .true.  !!GS: compute potential temperature or not
+   !
+   cn_dir      = './'      !  root directory for the bulk data location
+   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________!
+   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !       weights filename               ! rotation ! land/sea mask !
+   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                                      ! pairing  !    filename   !
+   sn_wndi     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'u10'    ,   .false.   , .true.  , 'yearly'  , ''   ,''   , ''
+   sn_wndj     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'v10'    ,   .false.   , .true.  , 'yearly'  , ''   ,''   , ''
+   sn_qsr      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rad_sw' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_qlw      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rad_lw' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_tair     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  't_air'  ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_humi     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rh_air' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_prec     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'precip' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_snow     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'snow'   ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_slp      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'slp'    ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+/
+!-----------------------------------------------------------------------
+&namsbc_abl    !   Atmospheric Boundary Layer formulation           (ln_abl = T)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3")
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namsbc_sas    !   Stand-Alone Surface boundary condition / STATION_ASF
+!-----------------------------------------------------------------------
+   l_sasread   = .true.    !  =T Read in file ;  =F set all to 0. (see sbcssm)
+   !
+   cn_dir      = './'      !  root directory for the ocean data location
+   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
+   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
+   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
+   sn_usp  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssu'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   sn_vsp  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssv'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   sn_tem  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'sst'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   sn_sal  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'sss',        .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   sn_ssh  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssh',        .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   !!
+/
+!!======================================================================
+!!               ***  Lateral boundary condition  ***                 !!
+!!                                                                    !!
+!!   namlbc        lateral momentum boundary condition                  (default: NO selection)
+!!   namagrif      agrif nested grid   (read by child model only)       ("key_agrif")
+!!   nam_tide      Tidal forcing                                        (default: OFF)
+!!   nambdy        Unstructured open boundaries                         (default: OFF)
+!!   nambdy_dta    Unstructured open boundaries - external data         (see  nambdy)
+!!   nambdy_tide   tidal forcing at open boundaries                     (default: OFF)
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namlbc        !   lateral momentum boundary condition                  (default: NO selection)
+!-----------------------------------------------------------------------
+   rn_shlat    =    0.     !  free slip
+/
+!!======================================================================
+!!                ***  Top/Bottom boundary condition  ***             !!
+!!                                                                    !!
+!!   namdrg        top/bottom drag coefficient                          (default: NO selection)
+!!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T)
+!!   namdrg_bot    bottom friction                                      (ln_OFF=F)
+!!   nambbc        bottom temperature boundary condition                (default: OFF)
+!!   nambbl        bottom boundary layer scheme                         (default: OFF)
+!!======================================================================
+!
+!!======================================================================
+!!                        Tracer (T-S) namelists                      !!
+!!                                                                    !!
+!!   nameos        equation of state                                    (default: NO selection)
+!!   namtra_adv    advection scheme                                     (default: NO selection)
+!!   namtra_ldf    lateral diffusion scheme                             (default: NO selection)
+!!   namtra_mle    mixed layer eddy param. (Fox-Kemper param.)          (default: OFF)
+!!   namtra_eiv    eddy induced velocity param.                         (default: OFF)
+!!   namtra_dmp    T & S newtonian damping                              (default: OFF)
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&nameos        !   ocean Equation Of Seawater                           (default: NO selection)
+!-----------------------------------------------------------------------
+   ln_eos80    = .true.         !  = Use EOS80
+/
+!!======================================================================
+!!                      ***  Dynamics namelists  ***                  !!
+!!                                                                    !!
+!!   nam_vvl       vertical coordinate options                          (default: z-star)
+!!   namdyn_adv    formulation of the momentum advection                (default: NO selection)
+!!   namdyn_vor    advection scheme                                     (default: NO selection)
+!!   namdyn_hpg    hydrostatic pressure gradient                        (default: NO selection)
+!!   namdyn_spg    surface pressure gradient                            (default: NO selection)
+!!   namdyn_ldf    lateral diffusion scheme                             (default: NO selection)
+!!   namdta_dyn    offline TOP: dynamics read in files                  (OFF_SRC only)
+!!======================================================================
+!
+!!======================================================================
+!!                     vertical physics namelists                     !!
+!!                                                                    !!
+!!    namzdf        vertical physics manager                            (default: NO selection)
+!!    namzdf_ric    richardson number vertical mixing                   (ln_zdfric=T)
+!!    namzdf_tke    TKE vertical mixing                                 (ln_zdftke=T)
+!!    namzdf_gls    GLS vertical mixing                                 (ln_zdfgls=T)
+!!    namzdf_osm    OSM vertical diffusion                              (ln_zdfosm=T)
+!!    namzdf_iwm    tidal mixing parameterization                       (ln_zdfiwm=T)
+!!======================================================================
+!
+!!======================================================================
+!!                  ***  Diagnostics namelists  ***                   !!
+!!                                                                    !!
+!!   namtrd       dynamics and/or tracer trends                         (default: OFF)
+!!   namptr       Poleward Transport Diagnostics                        (default: OFF)
+!!   namhsb       Heat and salt budgets                                 (default: OFF)
+!!   namdiu       Cool skin and warm layer models                       (default: OFF)
+!!   namdiu       Cool skin and warm layer models                       (default: OFF)
+!!   namflo       float parameters                                      (default: OFF)
+!!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF)
+!!   nam_diadct   transports through some sections                      (default: OFF)
+!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF)
+!!   nam_dia25h   25h Mean Output                                       (default: OFF)
+!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4")
+!!======================================================================
+!
+!!======================================================================
+!!               ***  Observation & Assimilation  ***                 !!
+!!                                                                    !!
+!!   namobs       observation and model comparison                      (default: OFF)
+!!   nam_asminc   assimilation increments                               ('key_asminc')
+!!======================================================================
+!
+!!======================================================================
+!!                  ***  Miscellaneous namelists  ***                 !!
+!!                                                                    !!
+!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi")
+!!   namctl            Control prints                                   (default: OFF)
+!!   namsto            Stochastic parametrization of EOS                (default: OFF)
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi")
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namctl        !   Control prints                                       (default: OFF)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namsto        ! Stochastic parametrization of EOS                      (default: OFF)
+!-----------------------------------------------------------------------
+/
Index: /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/oce/namelist_ecmwf_cfg
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/oce/namelist_ecmwf_cfg	(revision 13676)
+++ /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/oce/namelist_ecmwf_cfg	(revision 13676)
@@ -0,0 +1,274 @@
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!! NEMO/OCE  Configuration namelist : overwrite default values defined in SHARED/namelist_ref
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!!
+!!             STATION_ASF configuration, IDEALIZED forcing
+!!                   ==> meant to be used with forcing from file:
+!!                       "input_output_VALIDATION_IDEALIZED.nc"
+!!
+!!======================================================================
+!!              ***  Domain & Run management namelists  ***           !!
+!!                                                                    !!
+!!   namrun       parameters of the run
+!!   namdom       space and time domain
+!!   namcfg       parameters of the configuration                       (default: user defined GYRE)
+!!   namwad       Wetting and drying                                    (default: OFF)
+!!   namtsd       data: temperature & salinity                          (default: OFF)
+!!   namcrs       coarsened grid (for outputs and/or TOP)               (ln_crs =T)
+!!   namc1d       1D configuration options                              ("key_c1d")
+!!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d")
+!!   namc1d_uvd   1D data (currents)                                    ("key_c1d")
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namusr_def    !   STATION_ASF user defined namelist
+!-----------------------------------------------------------------------
+   rn_dept1    =  1.0   ! depth (m) at which the prescribed SST is taken  == depth of first T point! (IMPORTANT for warm-layer impact, the deeper, the more pronounced!)
+   !                    ! [PAPA station = 1 m]
+/
+!-----------------------------------------------------------------------
+&namrun        !   parameters of the run
+!-----------------------------------------------------------------------
+   cn_exp      =  'STATION_ASF-ECMWF_IDEALIZED'  !  experience name
+   nn_it000    =    1   !  first time step
+   nn_itend    =       104  !  last  time step => 3 years (including 1 leap!) at dt=3600s
+   nn_date0    = 20200101  !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
+   nn_leapy    =       1   !  Leap year calendar (1) or not (0)
+   !
+   nn_stock    =    104
+   nn_write    =    104
+/
+!-----------------------------------------------------------------------
+&namdom        !   time and space domain / STATION_ASF
+!-----------------------------------------------------------------------
+   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time
+   rn_Dt       =  3600.    !  time step for the dynamics and tracer => frequency of atmospheric forcing
+   ln_meshmask = .false.
+/
+!-----------------------------------------------------------------------
+&namcfg        !   parameters of the configuration                      (default: use namusr_def in namelist_cfg)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namc1d        !   1D configuration options / STATION_ASF               ("key_c1d" default: PAPA station)
+!-----------------------------------------------------------------------
+   rn_lat1d    =      50.1 !  Column latitude    [PAPA station]
+   rn_lon1d    =    215.1  !  Column longitude   [PAPA station]
+/
+!-----------------------------------------------------------------------
+&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d" default: OFF)
+!-----------------------------------------------------------------------
+   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F)
+/
+!-----------------------------------------------------------------------
+&namc1d_uvd    !   data: U & V currents                                 ("key_c1d" default: OFF)
+!-----------------------------------------------------------------------
+   !                       !  =T read U-V fields for:
+   ln_uvd_init   = .false.       !  ocean initialisation
+   ln_uvd_dyndmp = .false.       !  U-V restoring
+
+   cn_dir      = './'      !  root directory for the U-V data location
+   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
+   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
+   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
+   sn_ucur     = 'ucurrent'              ,         -1.       ,'u_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Ume'   , ''
+   sn_vcur     = 'vcurrent'              ,         -1.       ,'v_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Vme'   , ''
+/
+
+!!======================================================================
+!!            ***  Surface Boundary Condition namelists  ***          !!
+!!                                                                    !!
+!!                           STATION_ASF test-case                    !!
+!!                                                                    !!
+!!   namsbc          surface boundary condition manager                 (default: NO selection)
+!!   namsbc_flx      flux               formulation                     (ln_flx     =T)
+!!   namsbc_blk      Bulk formulae formulation                          (ln_blk     =T)
+!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3" )
+!!   namsbc_sas      Stand-Alone Surface module                         (SAS_SRC  only)
+!!   namsbc_iif      Ice-IF: use observed ice cover                     (nn_ice = 1   )
+!!   namtra_qsr      penetrative solar radiation                        (ln_traqsr  =T)
+!!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T)
+!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T)
+!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T)
+!!   namsbc_wave     external fields from wave model                    (ln_wave    =T)
+!!   namberg         iceberg floats                                     (ln_icebergs=T)
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namsbc        !   Surface Boundary Condition manager                   (default: NO selection)
+!-----------------------------------------------------------------------
+   nn_fsbc     = 1         !!LOLO!  frequency of SBC module call
+   ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk )
+   !
+   nn_ice      = 0  !LOLO  !  =0 no ice boundary condition
+      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif )
+      !                    !  =2 or 3 for SI3 and CICE, respectively
+   !    
+   ln_dm2dc    = .false.   !#lolo  daily mean to diurnal cycle on short wave
+/
+!-----------------------------------------------------------------------
+&namsbc_blk    !   namsbc_blk  generic Bulk formula / STATION_ASF       (ln_blk =T)
+!-----------------------------------------------------------------------
+   !                    !  bulk algorithm :
+   ln_NCAR      = .false.    ! "NCAR"      algorithm   (Large and Yeager 2008)
+   ln_COARE_3p0 = .false.    ! "COARE 3.0" algorithm   (Fairall et al. 2003)
+   ln_COARE_3p6 = .false.    ! "COARE 3.6" algorithm   (Edson et al. 2013)
+   ln_ECMWF     = .true.     ! "ECMWF"     algorithm   (IFS cycle 45r1)
+   ln_ANDREAS   = .false.    ! "ANDREAS"   algorithm   (Andreas et al. 2015)
+      rn_zqt       =  2.     !  Air temperature & humidity reference height (m)
+      rn_zu        = 10.     !  Wind vector reference height (m)
+      nn_iter_algo = 20      !  Number of iterations in bulk param. algo ("stable ABL + weak wind" requires more)
+      ln_skin_cs   = .true.  !  use the cool-skin parameterization
+      ln_skin_wl   = .true.  !  use the warm-layer parameterization
+   !
+   ln_humi_sph = .false. !  humidity "sn_humi" is specific humidity  [kg/kg]
+   ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K]
+   ln_humi_rlh = .true.  !  humidity "sn_humi" is relative humidity     [%]
+   ln_tpot     = .true.  !!GS: compute potential temperature or not
+   !
+   cn_dir      = './'      !  root directory for the bulk data location
+   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________!
+   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !       weights filename               ! rotation ! land/sea mask !
+   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                                      ! pairing  !    filename   !
+   sn_wndi     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'u10'    ,   .false.   , .true.  , 'yearly'  , ''   ,''   , ''
+   sn_wndj     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'v10'    ,   .false.   , .true.  , 'yearly'  , ''   ,''   , ''
+   sn_qsr      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rad_sw' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_qlw      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rad_lw' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_tair     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  't_air'  ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_humi     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rh_air' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_prec     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'precip' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_snow     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'snow'   ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_slp      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'slp'    ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+/
+!-----------------------------------------------------------------------
+&namsbc_abl    !   Atmospheric Boundary Layer formulation           (ln_abl = T)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3")
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namsbc_sas    !   Stand-Alone Surface boundary condition / STATION_ASF
+!-----------------------------------------------------------------------
+   l_sasread   = .true.    !  =T Read in file ;  =F set all to 0. (see sbcssm)
+   !
+   cn_dir      = './'      !  root directory for the ocean data location
+   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
+   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
+   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
+   sn_usp  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssu'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   sn_vsp  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssv'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   sn_tem  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'sst'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   sn_sal  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'sss',        .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   sn_ssh  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssh',        .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   !!
+/
+!!======================================================================
+!!               ***  Lateral boundary condition  ***                 !!
+!!                                                                    !!
+!!   namlbc        lateral momentum boundary condition                  (default: NO selection)
+!!   namagrif      agrif nested grid   (read by child model only)       ("key_agrif")
+!!   nam_tide      Tidal forcing                                        (default: OFF)
+!!   nambdy        Unstructured open boundaries                         (default: OFF)
+!!   nambdy_dta    Unstructured open boundaries - external data         (see  nambdy)
+!!   nambdy_tide   tidal forcing at open boundaries                     (default: OFF)
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namlbc        !   lateral momentum boundary condition                  (default: NO selection)
+!-----------------------------------------------------------------------
+   rn_shlat    =    0.     !  free slip
+/
+!!======================================================================
+!!                ***  Top/Bottom boundary condition  ***             !!
+!!                                                                    !!
+!!   namdrg        top/bottom drag coefficient                          (default: NO selection)
+!!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T)
+!!   namdrg_bot    bottom friction                                      (ln_OFF=F)
+!!   nambbc        bottom temperature boundary condition                (default: OFF)
+!!   nambbl        bottom boundary layer scheme                         (default: OFF)
+!!======================================================================
+!
+!!======================================================================
+!!                        Tracer (T-S) namelists                      !!
+!!                                                                    !!
+!!   nameos        equation of state                                    (default: NO selection)
+!!   namtra_adv    advection scheme                                     (default: NO selection)
+!!   namtra_ldf    lateral diffusion scheme                             (default: NO selection)
+!!   namtra_mle    mixed layer eddy param. (Fox-Kemper param.)          (default: OFF)
+!!   namtra_eiv    eddy induced velocity param.                         (default: OFF)
+!!   namtra_dmp    T & S newtonian damping                              (default: OFF)
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&nameos        !   ocean Equation Of Seawater                           (default: NO selection)
+!-----------------------------------------------------------------------
+   ln_eos80    = .true.         !  = Use EOS80
+/
+!!======================================================================
+!!                      ***  Dynamics namelists  ***                  !!
+!!                                                                    !!
+!!   nam_vvl       vertical coordinate options                          (default: z-star)
+!!   namdyn_adv    formulation of the momentum advection                (default: NO selection)
+!!   namdyn_vor    advection scheme                                     (default: NO selection)
+!!   namdyn_hpg    hydrostatic pressure gradient                        (default: NO selection)
+!!   namdyn_spg    surface pressure gradient                            (default: NO selection)
+!!   namdyn_ldf    lateral diffusion scheme                             (default: NO selection)
+!!   namdta_dyn    offline TOP: dynamics read in files                  (OFF_SRC only)
+!!======================================================================
+!
+!!======================================================================
+!!                     vertical physics namelists                     !!
+!!                                                                    !!
+!!    namzdf        vertical physics manager                            (default: NO selection)
+!!    namzdf_ric    richardson number vertical mixing                   (ln_zdfric=T)
+!!    namzdf_tke    TKE vertical mixing                                 (ln_zdftke=T)
+!!    namzdf_gls    GLS vertical mixing                                 (ln_zdfgls=T)
+!!    namzdf_osm    OSM vertical diffusion                              (ln_zdfosm=T)
+!!    namzdf_iwm    tidal mixing parameterization                       (ln_zdfiwm=T)
+!!======================================================================
+!
+!!======================================================================
+!!                  ***  Diagnostics namelists  ***                   !!
+!!                                                                    !!
+!!   namtrd       dynamics and/or tracer trends                         (default: OFF)
+!!   namptr       Poleward Transport Diagnostics                        (default: OFF)
+!!   namhsb       Heat and salt budgets                                 (default: OFF)
+!!   namdiu       Cool skin and warm layer models                       (default: OFF)
+!!   namdiu       Cool skin and warm layer models                       (default: OFF)
+!!   namflo       float parameters                                      (default: OFF)
+!!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF)
+!!   nam_diadct   transports through some sections                      (default: OFF)
+!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF)
+!!   nam_dia25h   25h Mean Output                                       (default: OFF)
+!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4")
+!!======================================================================
+!
+!!======================================================================
+!!               ***  Observation & Assimilation  ***                 !!
+!!                                                                    !!
+!!   namobs       observation and model comparison                      (default: OFF)
+!!   nam_asminc   assimilation increments                               ('key_asminc')
+!!======================================================================
+!
+!!======================================================================
+!!                  ***  Miscellaneous namelists  ***                 !!
+!!                                                                    !!
+!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi")
+!!   namctl            Control prints                                   (default: OFF)
+!!   namsto            Stochastic parametrization of EOS                (default: OFF)
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi")
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namctl        !   Control prints                                       (default: OFF)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namsto        ! Stochastic parametrization of EOS                      (default: OFF)
+!-----------------------------------------------------------------------
+/
Index: /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/oce/namelist_ncar_cfg
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/oce/namelist_ncar_cfg	(revision 13676)
+++ /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/IDEALIZED/oce/namelist_ncar_cfg	(revision 13676)
@@ -0,0 +1,274 @@
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!! NEMO/OCE  Configuration namelist : overwrite default values defined in SHARED/namelist_ref
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!!
+!!             STATION_ASF configuration, IDEALIZED forcing
+!!                   ==> meant to be used with forcing from file:
+!!                       "input_output_VALIDATION_IDEALIZED.nc"
+!!
+!!======================================================================
+!!              ***  Domain & Run management namelists  ***           !!
+!!                                                                    !!
+!!   namrun       parameters of the run
+!!   namdom       space and time domain
+!!   namcfg       parameters of the configuration                       (default: user defined GYRE)
+!!   namwad       Wetting and drying                                    (default: OFF)
+!!   namtsd       data: temperature & salinity                          (default: OFF)
+!!   namcrs       coarsened grid (for outputs and/or TOP)               (ln_crs =T)
+!!   namc1d       1D configuration options                              ("key_c1d")
+!!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d")
+!!   namc1d_uvd   1D data (currents)                                    ("key_c1d")
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namusr_def    !   STATION_ASF user defined namelist
+!-----------------------------------------------------------------------
+   rn_dept1    =  1.0   ! depth (m) at which the prescribed SST is taken  == depth of first T point! (IMPORTANT for warm-layer impact, the deeper, the more pronounced!)
+   !                    ! [PAPA station = 1 m]
+/
+!-----------------------------------------------------------------------
+&namrun        !   parameters of the run
+!-----------------------------------------------------------------------
+   cn_exp      =  'STATION_ASF-NCAR_IDEALIZED'  !  experience name
+   nn_it000    =    1   !  first time step
+   nn_itend    =       104  !  last  time step => 3 years (including 1 leap!) at dt=3600s
+   nn_date0    = 20200101  !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
+   nn_leapy    =       1   !  Leap year calendar (1) or not (0)
+   !
+   nn_stock    =    104
+   nn_write    =    104
+/
+!-----------------------------------------------------------------------
+&namdom        !   time and space domain / STATION_ASF
+!-----------------------------------------------------------------------
+   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time
+   rn_Dt       =  3600.    !  time step for the dynamics and tracer => frequency of atmospheric forcing
+   ln_meshmask = .false.
+/
+!-----------------------------------------------------------------------
+&namcfg        !   parameters of the configuration                      (default: use namusr_def in namelist_cfg)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namc1d        !   1D configuration options / STATION_ASF               ("key_c1d" default: PAPA station)
+!-----------------------------------------------------------------------
+   rn_lat1d    =      50.1 !  Column latitude    [PAPA station]
+   rn_lon1d    =    215.1  !  Column longitude   [PAPA station]
+/
+!-----------------------------------------------------------------------
+&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d" default: OFF)
+!-----------------------------------------------------------------------
+   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F)
+/
+!-----------------------------------------------------------------------
+&namc1d_uvd    !   data: U & V currents                                 ("key_c1d" default: OFF)
+!-----------------------------------------------------------------------
+   !                       !  =T read U-V fields for:
+   ln_uvd_init   = .false.       !  ocean initialisation
+   ln_uvd_dyndmp = .false.       !  U-V restoring
+
+   cn_dir      = './'      !  root directory for the U-V data location
+   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
+   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
+   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
+   sn_ucur     = 'ucurrent'              ,         -1.       ,'u_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Ume'   , ''
+   sn_vcur     = 'vcurrent'              ,         -1.       ,'v_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Vme'   , ''
+/
+
+!!======================================================================
+!!            ***  Surface Boundary Condition namelists  ***          !!
+!!                                                                    !!
+!!                           STATION_ASF test-case                         !!
+!!                                                                    !!
+!!   namsbc          surface boundary condition manager                 (default: NO selection)
+!!   namsbc_flx      flux               formulation                     (ln_flx     =T)
+!!   namsbc_blk      Bulk formulae formulation                          (ln_blk     =T)
+!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3" )
+!!   namsbc_sas      Stand-Alone Surface module                         (SAS_SRC  only)
+!!   namsbc_iif      Ice-IF: use observed ice cover                     (nn_ice = 1   )
+!!   namtra_qsr      penetrative solar radiation                        (ln_traqsr  =T)
+!!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T)
+!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T)
+!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T)
+!!   namsbc_wave     external fields from wave model                    (ln_wave    =T)
+!!   namberg         iceberg floats                                     (ln_icebergs=T)
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namsbc        !   Surface Boundary Condition manager                   (default: NO selection)
+!-----------------------------------------------------------------------
+   nn_fsbc     = 1         !!LOLO!  frequency of SBC module call
+   ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk )
+   !
+   nn_ice      = 0  !LOLO  !  =0 no ice boundary condition
+      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif )
+      !                    !  =2 or 3 for SI3 and CICE, respectively
+   !    
+   ln_dm2dc    = .false.   !#lolo  daily mean to diurnal cycle on short wave
+/
+!-----------------------------------------------------------------------
+&namsbc_blk    !   namsbc_blk  generic Bulk formula / STATION_ASF       (ln_blk =T)
+!-----------------------------------------------------------------------
+   !                    !  bulk algorithm :
+   ln_NCAR      = .true.     ! "NCAR"      algorithm   (Large and Yeager 2008)
+   ln_COARE_3p0 = .false.    ! "COARE 3.0" algorithm   (Fairall et al. 2003)
+   ln_COARE_3p6 = .false.    ! "COARE 3.6" algorithm   (Edson et al. 2013)
+   ln_ECMWF     = .false.    ! "ECMWF"     algorithm   (IFS cycle 45r1)
+   ln_ANDREAS   = .false.    ! "ANDREAS"   algorithm   (Andreas et al. 2015)
+      rn_zqt       =  2.     !  Air temperature & humidity reference height (m)
+      rn_zu        = 10.     !  Wind vector reference height (m)
+      nn_iter_algo = 20      !  Number of iterations in bulk param. algo ("stable ABL + weak wind" requires more)
+      ln_skin_cs   = .false. !  use the cool-skin parameterization
+      ln_skin_wl   = .false. !  use the warm-layer parameterization
+   !
+   ln_humi_sph = .false. !  humidity "sn_humi" is specific humidity  [kg/kg]
+   ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K]
+   ln_humi_rlh = .true.  !  humidity "sn_humi" is relative humidity     [%]
+   ln_tpot     = .true.  !!GS: compute potential temperature or not
+   !
+   cn_dir      = './'      !  root directory for the bulk data location
+   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________!
+   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !       weights filename               ! rotation ! land/sea mask !
+   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                                      ! pairing  !    filename   !
+   sn_wndi     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'u10'    ,   .false.   , .true.  , 'yearly'  , ''   ,''   , ''
+   sn_wndj     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'v10'    ,   .false.   , .true.  , 'yearly'  , ''   ,''   , ''
+   sn_qsr      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rad_sw' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_qlw      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rad_lw' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_tair     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  't_air'  ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_humi     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'rh_air' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_prec     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'precip' ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_snow     = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'snow'   ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+   sn_slp      = 'input_output_VALIDATION_IDEALIZED',  1.  ,  'slp'    ,   .false.   , .true.  , 'yearly'  , ''  , ''   , ''
+/
+!-----------------------------------------------------------------------
+&namsbc_abl    !   Atmospheric Boundary Layer formulation           (ln_abl = T)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3")
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namsbc_sas    !   Stand-Alone Surface boundary condition / STATION_ASF
+!-----------------------------------------------------------------------
+   l_sasread   = .true.    !  =T Read in file ;  =F set all to 0. (see sbcssm)
+   !
+   cn_dir      = './'      !  root directory for the ocean data location
+   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
+   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
+   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
+   sn_usp  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssu'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   sn_vsp  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssv'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   sn_tem  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'sst'     ,   .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   sn_sal  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'sss',        .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   sn_ssh  = 'input_output_VALIDATION_IDEALIZED', 1.         , 'ssh',        .false.   , .true.  , 'yearly'  ,    ''            ,    ''    ,    ''
+   !!
+/
+!!======================================================================
+!!               ***  Lateral boundary condition  ***                 !!
+!!                                                                    !!
+!!   namlbc        lateral momentum boundary condition                  (default: NO selection)
+!!   namagrif      agrif nested grid   (read by child model only)       ("key_agrif")
+!!   nam_tide      Tidal forcing                                        (default: OFF)
+!!   nambdy        Unstructured open boundaries                         (default: OFF)
+!!   nambdy_dta    Unstructured open boundaries - external data         (see  nambdy)
+!!   nambdy_tide   tidal forcing at open boundaries                     (default: OFF)
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namlbc        !   lateral momentum boundary condition                  (default: NO selection)
+!-----------------------------------------------------------------------
+   rn_shlat    =    0.     !  free slip
+/
+!!======================================================================
+!!                ***  Top/Bottom boundary condition  ***             !!
+!!                                                                    !!
+!!   namdrg        top/bottom drag coefficient                          (default: NO selection)
+!!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T)
+!!   namdrg_bot    bottom friction                                      (ln_OFF=F)
+!!   nambbc        bottom temperature boundary condition                (default: OFF)
+!!   nambbl        bottom boundary layer scheme                         (default: OFF)
+!!======================================================================
+!
+!!======================================================================
+!!                        Tracer (T-S) namelists                      !!
+!!                                                                    !!
+!!   nameos        equation of state                                    (default: NO selection)
+!!   namtra_adv    advection scheme                                     (default: NO selection)
+!!   namtra_ldf    lateral diffusion scheme                             (default: NO selection)
+!!   namtra_mle    mixed layer eddy param. (Fox-Kemper param.)          (default: OFF)
+!!   namtra_eiv    eddy induced velocity param.                         (default: OFF)
+!!   namtra_dmp    T & S newtonian damping                              (default: OFF)
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&nameos        !   ocean Equation Of Seawater                           (default: NO selection)
+!-----------------------------------------------------------------------
+   ln_eos80    = .true.         !  = Use EOS80
+/
+!!======================================================================
+!!                      ***  Dynamics namelists  ***                  !!
+!!                                                                    !!
+!!   nam_vvl       vertical coordinate options                          (default: z-star)
+!!   namdyn_adv    formulation of the momentum advection                (default: NO selection)
+!!   namdyn_vor    advection scheme                                     (default: NO selection)
+!!   namdyn_hpg    hydrostatic pressure gradient                        (default: NO selection)
+!!   namdyn_spg    surface pressure gradient                            (default: NO selection)
+!!   namdyn_ldf    lateral diffusion scheme                             (default: NO selection)
+!!   namdta_dyn    offline TOP: dynamics read in files                  (OFF_SRC only)
+!!======================================================================
+!
+!!======================================================================
+!!                     vertical physics namelists                     !!
+!!                                                                    !!
+!!    namzdf        vertical physics manager                            (default: NO selection)
+!!    namzdf_ric    richardson number vertical mixing                   (ln_zdfric=T)
+!!    namzdf_tke    TKE vertical mixing                                 (ln_zdftke=T)
+!!    namzdf_gls    GLS vertical mixing                                 (ln_zdfgls=T)
+!!    namzdf_osm    OSM vertical diffusion                              (ln_zdfosm=T)
+!!    namzdf_iwm    tidal mixing parameterization                       (ln_zdfiwm=T)
+!!======================================================================
+!
+!!======================================================================
+!!                  ***  Diagnostics namelists  ***                   !!
+!!                                                                    !!
+!!   namtrd       dynamics and/or tracer trends                         (default: OFF)
+!!   namptr       Poleward Transport Diagnostics                        (default: OFF)
+!!   namhsb       Heat and salt budgets                                 (default: OFF)
+!!   namdiu       Cool skin and warm layer models                       (default: OFF)
+!!   namdiu       Cool skin and warm layer models                       (default: OFF)
+!!   namflo       float parameters                                      (default: OFF)
+!!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF)
+!!   nam_diadct   transports through some sections                      (default: OFF)
+!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF)
+!!   nam_dia25h   25h Mean Output                                       (default: OFF)
+!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4")
+!!======================================================================
+!
+!!======================================================================
+!!               ***  Observation & Assimilation  ***                 !!
+!!                                                                    !!
+!!   namobs       observation and model comparison                      (default: OFF)
+!!   nam_asminc   assimilation increments                               ('key_asminc')
+!!======================================================================
+!
+!!======================================================================
+!!                  ***  Miscellaneous namelists  ***                 !!
+!!                                                                    !!
+!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi")
+!!   namctl            Control prints                                   (default: OFF)
+!!   namsto            Stochastic parametrization of EOS                (default: OFF)
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi")
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namctl        !   Control prints                                       (default: OFF)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namsto        ! Stochastic parametrization of EOS                      (default: OFF)
+!-----------------------------------------------------------------------
+/
Index: /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/PAPA/oce/namelist_andreas_cfg
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/PAPA/oce/namelist_andreas_cfg	(revision 13675)
+++ /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/PAPA/oce/namelist_andreas_cfg	(revision 13676)
@@ -26,5 +26,5 @@
 &namrun        !   parameters of the run
 !-----------------------------------------------------------------------
-   cn_exp      =  'STATION_ASF-ANDREAS'  !  experience name
+   cn_exp      =  'STATION_ASF-ANDREAS_PAPA'  !  experience name
    nn_it000    =    1   !  first time step
    nn_itend    =    8760   !  last  time step => 3 years (including 1 leap!) at dt=3600s
@@ -117,5 +117,5 @@
       ln_skin_cs   = .false. !  use the cool-skin parameterization
       ln_skin_wl   = .false. !  use the warm-layer parameterization
-      !
+   !
    ln_humi_sph = .false. !  humidity "sn_humi" is specific humidity  [kg/kg]
    ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K]
Index: /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/PAPA/oce/namelist_coare3p6_cfg
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/PAPA/oce/namelist_coare3p6_cfg	(revision 13675)
+++ /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/PAPA/oce/namelist_coare3p6_cfg	(revision 13676)
@@ -26,5 +26,5 @@
 &namrun        !   parameters of the run
 !-----------------------------------------------------------------------
-   cn_exp      =  'STATION_ASF-COARE3p6'  !  experience name
+   cn_exp      =  'STATION_ASF-COARE3p6_PAPA'  !  experience name
    nn_it000    =    1   !  first time step
    nn_itend    =    8760   !  last  time step => 3 years (including 1 leap!) at dt=3600s
@@ -117,5 +117,5 @@
       ln_skin_cs   = .true.  !  use the cool-skin parameterization
       ln_skin_wl   = .true.  !  use the warm-layer parameterization
-      !
+   !
    ln_humi_sph = .false. !  humidity "sn_humi" is specific humidity  [kg/kg]
    ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K]
Index: /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/PAPA/oce/namelist_ecmwf_cfg
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/PAPA/oce/namelist_ecmwf_cfg	(revision 13675)
+++ /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/PAPA/oce/namelist_ecmwf_cfg	(revision 13676)
@@ -26,5 +26,5 @@
 &namrun        !   parameters of the run
 !-----------------------------------------------------------------------
-   cn_exp      =  'STATION_ASF-ECMWF'  !  experience name
+   cn_exp      =  'STATION_ASF-ECMWF_PAPA'  !  experience name
    nn_it000    =    1   !  first time step
    nn_itend    =    8760   !  last  time step => 3 years (including 1 leap!) at dt=3600s
@@ -115,7 +115,7 @@
       rn_zu        =  4.     !  Wind vector reference height (m)                [PAPA station]
       nn_iter_algo = 20      !  Number of iterations in bulk param. algo ("stable ABL + weak wind" requires more)
-      ln_skin_cs   = .true.  !  use the cool-skin parameterization
-      ln_skin_wl   = .true.  !  use the warm-layer parameterization
-      !
+      ln_skin_cs = .true.    !  use the cool-skin parameterization
+      ln_skin_wl = .true.    !  use the warm-layer parameterization
+   !
    ln_humi_sph = .false. !  humidity "sn_humi" is specific humidity  [kg/kg]
    ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K]
Index: /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/PAPA/oce/namelist_ncar_cfg
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/PAPA/oce/namelist_ncar_cfg	(revision 13675)
+++ /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/PAPA/oce/namelist_ncar_cfg	(revision 13676)
@@ -26,5 +26,5 @@
 &namrun        !   parameters of the run
 !-----------------------------------------------------------------------
-   cn_exp      =  'STATION_ASF-NCAR'  !  experience name
+   cn_exp      =  'STATION_ASF-NCAR_PAPA'  !  experience name
    nn_it000    =    1   !  first time step
    nn_itend    =    8760   !  last  time step => 3 years (including 1 leap!) at dt=3600s
Index: /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/analyze_output.py
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/analyze_output.py	(revision 13676)
+++ /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/analyze_output.py	(revision 13676)
@@ -0,0 +1,267 @@
+#!/usr/bin/env python3
+# -*- Mode: Python; coding: utf-8; indent-tabs-mode: nil; tab-width: 4 -*-
+#
+#######################################################################################
+# This script analyzes the output of STATION_ASF test-case with IDEALIZED forcing
+# in order to test the validity of computed fluxes and bulk transfer coefficient.
+# Beside an explicit standard output message, a result file is spawned:
+#    * SBCBLK.success => the test passed !
+#    * SBCBLK.fail    => the test failed !
+#
+#   Brodeau, 2020
+#
+########################################################################################
+
+import sys
+from os import path as path
+import math
+import numpy as nmp
+from netCDF4 import Dataset,num2date
+import matplotlib as mpl
+mpl.use('Agg')
+import matplotlib.pyplot as plt
+import matplotlib.dates as mdates
+
+l_t_shift = False  ; # because time interp. is set to FALSE into "&namsbc_blk" of NEMO...
+#                    #  ==> so time array is shifted by 30 minutes but fluxes are the same (persitence of input fields ???)
+
+l_alg = [ 'ECMWF' , 'NCAR' , 'COARE3p0', 'COARE3p6', 'ANDREAS' ]
+nb_alg     =    len(l_alg)
+
+
+# Variables to check:
+l_var_rf = [  'Qsen'   ,   'Qlat'   ,  'Qlw'   ,   'Tau',   'Cd'  ,   'Ce'   ] ; # In forcing file "IDEALIZED/input_output_VALIDATION_IDEALIZED.nc"
+l_var_ot = [ 'qsb_oce' , 'qla_oce'  , 'qlw_oce',  'taum', 'Cd_oce', 'Ce_oce' ] ; # names in the NEMO/STATION_ASF output file (check file_def_oce.xml)
+nb_var  = len(l_var_rf)
+
+dir_figs='.'
+size_fig=(13,8)
+fig_ext='png'
+
+
+rDPI=100.
+
+
+class fclrs:
+    OKGR = '\033[92m'
+    FAIL = '\033[91m'
+    ENDC = '\033[0m'
+
+
+
+# Getting arguments:
+narg = len(sys.argv)
+if not narg in [3,4]:
+    print('Usage: '+sys.argv[0]+' <forcing_+_validation_file> <NEMO-STATION_ASF_output_directory> (<m> for more/debug)'); sys.exit(0)
+cf_rf = sys.argv[1]
+cdir_out = sys.argv[2]
+l_more = False
+if narg==4:
+    l_more = ( sys.argv[3] in ['m','M'] )
+
+
+    
+# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+# Populating and checking existence of files to be read
+# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+def chck4f(cf):
+    cmesg = 'ERROR: File '+cf+' does not exist !!!'
+    if not path.exists(cf): print(cmesg) ; sys.exit(0)
+
+print('\n')
+
+# Input forcing/valid file:
+chck4f(cf_rf)
+id_rf = Dataset(cf_rf)
+vt = id_rf.variables['time'][:]
+cunit_t = id_rf.variables['time'].units
+id_rf.close()
+Nt_rf = len(vt)
+vtime_rf = nmp.zeros(Nt_rf); vtime_rf[:]  = vt[:]
+del vt
+print(' *** in forcing/valid file, "time" is in "'+cunit_t+'", Nt = '+str(Nt_rf)+'\n')
+
+
+# STATION_ASF output files (1 file per algorithm):
+cf_nemo = []
+for ja in range(nb_alg):
+    cfi = cdir_out+'/STATION_ASF-'+l_alg[ja]+'_IDEALIZED_1h_20200101_20200105_gridT.nc'
+    chck4f(cfi)
+    cf_nemo.append(cfi)
+print('\n *** NEMO/STATION_ASF output files we are goin to check:')
+for ja in range(nb_alg): print(cf_nemo[ja])
+print('\n')
+#-----------------------------------------------------------------
+
+# Getting time array from the first file:
+id_nm = Dataset(cf_nemo[0])
+vt = id_nm.variables['time_counter'][:]
+cunit_t = id_nm.variables['time_counter'].units ; print(' "time_counter" is in "'+cunit_t+'"')
+id_nm.close()
+Nt = len(vt)
+vtime_nm = nmp.zeros(Nt); vtime_nm[:]  = vt[:]
+del vt
+
+if Nt != Nt_rf-1: print('ERROR: the two files do not agree in terms of record lengrth: '+Nt_rf-1+' vs '+Nt)
+
+print('\n *** Excellent! We are going to look at surface fluxes under '+str(Nt)+' different scenarios of air-sea stability/wind-speed conditions...\n')
+
+# 30 minute shift, just like NEMO
+vtime = nmp.zeros(Nt)
+if l_t_shift:
+    vtime[:] = 0.5*(vtime_rf[:-1] + vtime_rf[1:])
+else:
+     vtime[:] =     vtime_rf[:-1]
+# Debug:
+#for jt in range(3):
+#    print(' Ref.  , Nemo "', vtime_rf[jt], vtime_nm[jt], vtime[jt])
+#sys.exit(0)
+
+
+
+
+##
+IREPORT = nmp.zeros((nb_alg,nb_var), dtype=int)
+
+
+
+
+# Loop on the fields to control...
+###################################
+
+jv=0
+for cv in l_var_rf:
+    cv_rf_m = cv+'_mean'
+    cv_rf_t = cv+'_tol'
+    cv_nemo = l_var_ot[jv]
+    print('\n\n ==== Checking flux '+cv_rf_m+' ('+cv_nemo+') !')
+
+    F_rf_m = nmp.zeros( Nt )
+    F_rf_t = nmp.zeros( Nt )
+    F_nemo = nmp.zeros((Nt,nb_alg))
+
+    wnd_rf = nmp.zeros( Nt ) ; #DEBUG
+    tzt_rf = nmp.zeros( Nt ) ; #DEBUG
+    qzt_rf = nmp.zeros( Nt ) ; #DEBUG
+    sst_rf = nmp.zeros( Nt ) ; #DEBUG
+    
+    id_rf = Dataset(cf_rf)
+    trf_m   = id_rf.variables[cv_rf_m][:] ; # Nt+1
+    trf_t   = id_rf.variables[cv_rf_t][:] ; # Nt+1
+    trf_wnd = id_rf.variables['wndspd'][:,1,1] ; # DEBUG
+    trf_tzt = id_rf.variables['t_air'][:,1,1] ; # DEBUG
+    trf_qzt = id_rf.variables['rh_air'][:,1,1] ; # DEBUG
+    trf_sst = id_rf.variables['sst'][:,1,1] ; # DEBUG
+    id_rf.close()
+    
+    if l_t_shift:
+        # 30 minute shift, just like NEMO
+        F_rf_m[:] = 0.5 * (trf_m[:-1] + trf_m[1:])
+        F_rf_t[:] = 0.5 * (trf_t[:-1] + trf_t[1:])
+    else:
+        F_rf_m[:] =        trf_m[:-1]
+        F_rf_t[:] =        trf_t[:-1]
+
+        wnd_rf[:] =        trf_wnd[:-1]        
+        tzt_rf[:] =        trf_tzt[:-1]
+        qzt_rf[:] =        trf_qzt[:-1]
+        sst_rf[:] =        trf_sst[:-1]        
+    
+    
+    for ja in range(nb_alg):
+        
+        calgo = l_alg[ja]
+
+        print(' *** '+calgo+' => '+cf_nemo[ja])
+
+        id_nemo = Dataset(cf_nemo[ja])
+        F_nemo[:,ja] = id_nemo.variables[cv_nemo][:,1,1] ; # it's 3x3 spatial domain, taking middle point !
+        id_nemo.close()
+
+
+        if l_more:
+            #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+            # How does all this look on on a figure?
+            cfig = l_var_rf[jv]+'_'+calgo+'.'+fig_ext
+            print(' *** will plot '+cfig)
+            fig = plt.figure(num=1, figsize=size_fig, facecolor='w', edgecolor='k')
+            ax1 = plt.axes([0.08, 0.25, 0.9, 0.7])
+            #
+            plt.plot(vtime, F_nemo[:,ja], label='NEMO['+calgo+']', zorder=1)
+            plt.plot(vtime, F_rf_m[:], color='k', label='MEAN REF!', zorder=10)
+            # +- rtol enveloppe:
+            plt.fill_between(vtime, F_rf_m[:]-F_rf_t[:], F_rf_m[:]+F_rf_t[:], alpha=0.2)
+            #
+            ax1.grid(color='k', linestyle='-', linewidth=0.3)
+            plt.legend(loc='best', ncol=1, shadow=True, fancybox=True)
+            plt.savefig(cfig, dpi=int(rDPI), transparent=False)
+            plt.close(1)
+            print('')
+            #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+
+        # Does the field look okay with respect to reference +- tolerance ?
+
+        l_overshoot_a = nmp.any( F_nemo[:,ja] > F_rf_m[:]+F_rf_t[:] )
+        l_overshoot_b = nmp.any( F_nemo[:,ja] < F_rf_m[:]-F_rf_t[:] )
+
+        if l_overshoot_a: print(fclrs.FAIL+'\n ***** BAD overshoot + for '+calgo+' for variable '+cv+' !!!\n'+fclrs.ENDC )
+            
+        if l_overshoot_b: print(fclrs.FAIL+'\n ***** BAD overshoot - for '+calgo+' for variable '+cv+' !!!'+fclrs.ENDC )
+
+        if l_overshoot_a or l_overshoot_b:
+            print(fclrs.FAIL+'\n ***** TEST NOT PASSED FOR '+calgo+' for variable '+cv+' !!!\n'+fclrs.ENDC )
+        else:
+            # We're all good !
+            IREPORT[ja,jv] = 1
+            if l_more: print(fclrs.OKGR+'\n ***** TEST PASSED FOR '+calgo+' for variable '+cv+' :D !!!\n'+fclrs.ENDC )
+                    
+    jv=jv+1
+
+
+
+l_ok =  nmp.sum(IREPORT[:,:]) == nb_var*nb_alg
+
+
+
+if l_ok:
+    ctxt      = 'PASSED'
+    ccol      = fclrs.OKGR
+    cf_report = 'SBCBLK.success'
+else:
+    ctxt      = 'FAILED'
+    ccol      = fclrs.FAIL
+    cf_report = 'SBCBLK.fail'
+    
+
+print(ccol+'\n\n ############ FINAL REPORT ############\n'+fclrs.ENDC)
+
+f = open(cf_report, 'w')
+
+f.write("### Sanity-check report for SBCBLK generated via 'STATION_ASF/EXP00/sbcblk_sanity_check.sh'\n\n")
+
+for ja in range(nb_alg):
+    calgo = l_alg[ja]
+    
+    if nmp.sum(IREPORT[ja,:]) == nb_var:
+        # Success for this algo
+        cbla = ' ***** Algorithm "'+calgo+'" PASSED sanity check !!!\n\n'
+        print(fclrs.OKGR+cbla+fclrs.ENDC ) ; f.write(cbla)
+    else:
+        # Algo FAILS!
+        cbla = ' ***** Algorithm "'+calgo+'" FAILED sanity check !!!\n'
+        print(fclrs.FAIL+cbla+fclrs.ENDC ) ; f.write(cbla)
+        (idx_fail,) = nmp.where(IREPORT[ja,:]==0)
+        for jv in idx_fail:
+            cbla = '   ==> on variable '+l_var_rf[jv]+' !\n\n'
+            print(fclrs.FAIL+cbla+fclrs.ENDC ) ; f.write(cbla)
+
+# Conclusion:
+
+cbla = '    ####################################\n'  +\
+       '    ###   TEST '+ctxt+' FOR SBCBLK !   ###\n'+\
+       '    ####################################\n\n'
+
+print(ccol+cbla+fclrs.ENDC ) ; f.write(cbla)
+
+f.close()
Index: /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/file_def_nemo-oce.xml
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/file_def_nemo-oce.xml	(revision 13675)
+++ /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/file_def_nemo-oce.xml	(revision 13676)
@@ -13,5 +13,4 @@
       <field field_ref="ssv_m"        name="ssv"      />
       <!--    -->
-      <field field_ref="ssq"          name="ssq"      />
       <field field_ref="rho_air"      name="rho_air"  />
       <field field_ref="dt_skin"      name="dt_skin"  />
@@ -27,5 +26,5 @@
       <field field_ref="evap_oce"     name="evap_oce" />
       <field field_ref="precip"       name="precip" />
-      <!-- <field field_ref="empmr"        name="empmr" /> -->
+      <!-- <field field_ref="empmr"   name="empmr" /> -->
       <!-- -->
       <field field_ref="taum"         name="taum"   />
@@ -36,5 +35,10 @@
       <field field_ref="Ch_oce"       name="Ch_oce" />
       <!-- -->
-      <field field_ref="ice_cover"    name="ice_cover" />
+      <field field_ref="theta_zt"     name="theta_zt" />
+      <field field_ref="q_zt"         name="q_zt" />
+      <field field_ref="theta_zu"     name="theta_zu" />
+      <field field_ref="q_zu"         name="q_zu" />
+      <field field_ref="ssq"          name="ssq" />
+      <field field_ref="wspd_blk"     name="wspd_blk" />      
     </file>
 
Index: /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/launch_sasf.sh
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/launch_sasf.sh	(revision 13675)
+++ /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/launch_sasf.sh	(revision 13676)
@@ -12,12 +12,14 @@
 
 CONFIG="STATION_ASF" ; # general name of the configuration
-CONFIG_BLD="STATION_ASF" ; # name of config as build in NEMO... (directory inside "tests" actually contains the compiled test-case?)
+CONFIG_BLD="STATION_ASF2" ; # name of config as build in NEMO... (directory inside "tests" actually contains the compiled test-case?)
 
 # Root directory NEMOGCM reference depository where to fetch compiled STATION_ASF nemo.exe + default namelists:
 NEMO_REF_DIR="/home/laurent/DEV/NEMO/NEMOGCM_dev_r13648_ASINTER-04_laurent_bulk_ice"
-NEMO_WRK_DIR="/home/laurent/NEMO/NEMOvdev_r13648_ASINTER-04_laurent_bulk_ice"
 
 # NEMOGCM root directory where to fetch compiled STATION_ASF nemo.exe:
-SASF_WRK_DIR="${NEMO_WRK_DIR}/tests/${CONFIG_BLD}"
+SASF_WRK_DIR="${NEMO_REF_DIR}/tests/${CONFIG_BLD}"
+
+# DATA_IN_DIR => Directory containing sea-surface + atmospheric forcings:
+DATA_IN_DIR="${NEMO_REF_DIR}/tests/${CONFIG}/input_data"
 
 # Directory where to run the simulation:
@@ -30,7 +32,6 @@
 i_sea_ice=0 ; # default does not compute fluxes over sea-ice
 #
-NFORC="ERA5_arctic" ; i_sea_ice=1 ; FORC="ERA5_arctic_surface_81N_36p75E_1h" ; # string sufficient to copy relevant files as in "*${FORC}*.nc"
-
-#NFORC="PAPA" ; FORC="Station_PAPA_50N-145W" ; # string sufficient to copy relevant files as in "*${FORC}*.nc"
+#NFORC="ERA5_arctic" ; i_sea_ice=1 ; FORC="ERA5_arctic_surface_81N_36p75E_1h" ; # string sufficient to copy relevant files as in "*${FORC}*.nc"
+NFORC="PAPA" ; FORC="Station_PAPA_50N-145W" ; # string sufficient to copy relevant files as in "*${FORC}*.nc"
 
 
@@ -40,4 +41,13 @@
 
 #================================================================================
+
+if [ "`which ncks`" = "" ]; then
+    echo
+    echo "WARNING: you do not seem to have NCO installed here... You should!"
+    echo "      => anyway! will do without, output fields will remaine 3x3 in space :( "
+    echo
+    sleep 4
+fi
+
 HERE=`pwd`
 
@@ -68,10 +78,10 @@
     LIST_OA_ALGOS="ECMWF" ;          # list of air-sea algorithms to test
     LIST_IA_ALGOS="LG15 LU12 CSTC" ; # list of air-ice algorithms to test
-    DIR_NL=${NFORC}/ice ; # where to get the namelists from...
+    DIR_NL=${NFORC}/ice ; # where to fetch the namelists from...
 else
     # Only OPEN-OCEAN/AIR flux computation
-    LIST_OA_ALGOS="NCAR ECMWF" ; # "NCAR ECMWF COARE3p6"; do   ; # list of air-sea algorithms to test
+    LIST_OA_ALGOS="NCAR ECMWF COARE3p6 ANDREAS"; # list of air-sea algorithms to test
     LIST_IA_ALGOS=""
-    DIR_NL=${NFORC}/oce ; # where to get the namelists from...
+    DIR_NL=${NFORC}/oce ; # where to fetch the namelists from...
 fi
 if [ ! -d ${DIR_NL} ]; then echo " Mhhh, seems like forcing ${NFORC} is not meant to be used with sea-ice/air flux computation !"; exit; fi
@@ -86,17 +96,4 @@
 
 
-# DATA_IN_DIR => Directory containing sea-surface + atmospheric forcings
-#             (get it there https://drive.google.com/file/d/1MxNvjhRHmMrL54y6RX7WIaM9-LGl--ZP/):
-if [ `hostname` = "merlat"        ]; then
-    DATA_IN_DIR="/MEDIA/data/STATION_ASF/input_data"
-elif [ `hostname` = "luitel"        ]; then
-    DATA_IN_DIR="/data/gcm_setup/STATION_ASF/input_data"
-elif [ `hostname` = "ige-meom-cal1" ]; then
-    DATA_IN_DIR="/mnt/meom/workdir/brodeau/STATION_ASF/input_data"
-elif [ `hostname` = "salvelinus" ]; then
-    DATA_IN_DIR="/opt/data/STATION_ASF/input_data"
-else
-    echo "Boo!"; exit
-fi
 if [ ! -d ${DATA_IN_DIR} ]; then echo "PROBLEM!!! => did not find directory 'input_data' with input forcing..."; exit; fi
 
@@ -164,17 +161,17 @@
 if [ "${LIST_IA_ALGOS}" = "" ]; then LIST_IA_ALGOS="none"; fi
 
-echo " {LIST_OA_ALGOS} = ${LIST_OA_ALGOS} !"
-echo " {LIST_IA_ALGOS} = ${LIST_IA_ALGOS} !"
+echo
 
 for OA_A in ${LIST_OA_ALGOS}; do
     for IA_A in ${LIST_IA_ALGOS}; do
-        
+
+        echo        
         CASE="${OA_A}"
         if [ ${i_sea_ice} -eq 1 ] && [ "${IA_A}" != "none" ]; then CASE="${OA_A}-${IA_A}"; fi
 
         echo ; echo
-        echo "===================================="
-        echo " Going for ${CASE} experiment"
-        echo "===================================="
+        echo "======================================================================="
+        echo " Going for experiment: ${CASE} bulk param. with ${NFORC} forcing "
+        echo "======================================================================="
         echo
 
@@ -204,6 +201,6 @@
         # Moving output files:
         mkdir -p output
-        mv -f ${CONFIG}-${CASE}_*_grid*.nc output/
-        if [ ${i_sea_ice} -eq 1 ]; then mv -f ${CONFIG}-${CASE}_*_icemod.nc output/; fi
+        mv -f ${CONFIG}-${CASE}_${NFORC}_*_grid*.nc output/
+        if [ ${i_sea_ice} -eq 1 ]; then mv -f ${CONFIG}-${CASE}_${NFORC}_*_icemod.nc output/; fi
 
         # Saving logs:
@@ -213,13 +210,10 @@
         # Making 3x3 to 1 !
         cd output/
-        mkdir -p 3x3
-        list=`\ls ${CONFIG}*${CASE}*.nc | grep -v '_restart_'`
+        list=`\ls ${CONFIG}*${CASE}_${NFORC}_*.nc | grep -v '_restart_'| grep -v '_1p.nc'`
         for ff in ${list}; do
             echo
-            fn=`echo ${ff} | sed -e "s|.nc|_${NFORC}_1p.nc|g"`
-            fb=`echo ${ff} | sed -e "s|.nc|_${NFORC}.nc|g"`
+            fn=`echo ${ff} | sed -e "s|.nc|_1p.nc|g"`
             CMD="ncks -O -d x,1 -d y,1 ${ff} -o ${fn}"
             echo " *** ${CMD}"; ${CMD}
-            mv -f ${ff} 3x3/${fb}
             echo
         done
Index: MO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/launch_sasf_old.sh
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/launch_sasf_old.sh	(revision 13675)
+++ 	(revision )
@@ -1,98 +1,0 @@
-#!/bin/bash
-
-################################################################
-#
-# Script to launch a set of STATION_ASF simulations
-#
-# L. Brodeau, 2020
-#
-################################################################
-
-# What directory inside "tests" actually contains the compiled "nemo.exe" for STATION_ASF ?
-TC_DIR="STATION_ASF2"
-
-expdir=`basename ${PWD}`; # we expect "EXPREF" or "EXP00" normally...
-
-# NEMOGCM root directory:
-NEMO_ROOT_DIR=`pwd | sed -e "s|/tests/STATION_ASF/${expdir}||g"`
-
-# NEMOGCM root directory where to fetch compiled STATION_ASF nemo.exe:
-SASF_WRK_DIR="${NEMO_ROOT_DIR}/tests/${TC_DIR}"
-
-# Directory where to run the simulation:
-PROD_DIR="${HOME}/tmp/STATION_ASF"
-
-
-####### End of normal user configurable section #######
-
-#================================================================================
-
-SASF_REF_DIR="${NEMO_ROOT_DIR}/tests/STATION_ASF"
-if [ ! -d ${SASF_REF_DIR} ]; then echo " Mhhh, no EXPREF directory ${SASF_REF_DIR} !"; exit; fi
-
-# NEMO executable to use is:
-NEMO_EXE="${SASF_WRK_DIR}/BLD/bin/nemo.exe"
-if [ ! -f ${NEMO_EXE} ]; then echo " Mhhh, no compiled 'nemo.exe' found into `dirname ${NEMO_EXE}` !"; exit; fi
-
-DATA_IN_DIR="${SASF_REF_DIR}/input_data" ; # Directory containing sea-surface + atmospheric input data
-if [ ! -d ${DATA_IN_DIR} ]; then echo "PROBLEM!!! => did not find directory 'input_data' with input forcing..."; exit; fi
-
-SASF_EXPREF=${SASF_REF_DIR}/${expdir}  ; # STATION_ASF EXPREF directory from which to use namelists and XIOS xml files...
-if [ ! -d ${SASF_EXPREF} ]; then echo " Mhhh, no ${expdir} directory ${SASF_EXPREF} !"; exit; fi
-
-
-echo "###########################################################"
-echo "#        S T A T I O N   A i r  -  S e a   F l u x        #"
-echo "###########################################################"
-echo
-echo "  * NEMO reference root directory is: ${NEMO_ROOT_DIR}"
-echo "  * STATION_ASF work directory is: ${SASF_WRK_DIR}"
-echo "       ==> NEMO EXE to use: ${NEMO_EXE}"
-echo
-echo "  * Input forcing data into: ${DATA_IN_DIR}"
-echo "  * Production will be done into: ${PROD_DIR}"
-echo "  * Directory in which namelists and xml files are fetched:"
-echo "       ==> ${SASF_EXPREF}"
-echo
-
-mkdir -p ${PROD_DIR}
-
-rsync -avP ${NEMO_EXE}          ${PROD_DIR}/
-
-for ff in "context_nemo.xml" "domain_def_nemo.xml" "field_def_nemo-oce.xml" "file_def_nemo-oce.xml" "grid_def_nemo.xml" "iodef.xml" "namelist_ref"; do
-    if [ ! -f ${SASF_EXPREF}/${ff} ]; then echo " Mhhh, ${ff} not found into ${SASF_EXPREF} !"; exit; fi
-    rsync -avPL ${SASF_EXPREF}/${ff} ${PROD_DIR}/
-done
-
-# Copy forcing to work directory:
-rsync -avP ${DATA_IN_DIR}/Station_PAPA_50N-145W*.nc ${PROD_DIR}/
-
-for CASE in "ECMWF" "COARE3p6" "NCAR" "ECMWF-noskin" "COARE3p6-noskin"; do
-
-    echo ; echo
-    echo "============================="
-    echo " Going for ${CASE} experiment"
-    echo "============================="
-    echo
-
-    scase=`echo "${CASE}" | tr '[:upper:]' '[:lower:]'`
-
-    rm -f ${PROD_DIR}/namelist_cfg
-    rsync -avPL ${SASF_EXPREF}/namelist_${scase}_cfg ${PROD_DIR}/namelist_cfg
-
-    cd ${PROD_DIR}/
-    echo
-    echo "Launching NEMO !"
-    ./nemo.exe 1>out_nemo.out 2>err_nemo.err
-    echo "Done!"
-    echo
-
-    # Moving output files:
-    mkdir -p output
-    mv -f STATION_ASF-${CASE}_*_grid*.nc output/
-    
-    # Saving logs:
-    mkdir -p ${CASE}_log
-    mv -f *.out *.err ocean.output output.namelist.dyn ${CASE}_log/
-
-done
Index: /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/plot_station_asf_OCE.py
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/plot_station_asf_OCE.py	(revision 13676)
+++ /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/plot_station_asf_OCE.py	(revision 13676)
@@ -0,0 +1,233 @@
+#!/usr/bin/env python3
+# -*- Mode: Python; coding: utf-8; indent-tabs-mode: nil; tab-width: 4 -*-
+
+# Post-diagnostic of STATION_ASF /  L. Brodeau, 2019
+
+import sys
+from os import path as path
+import math
+import numpy as nmp
+from netCDF4 import Dataset,num2date
+import matplotlib as mpl
+mpl.use('Agg')
+import matplotlib.pyplot as plt
+import matplotlib.dates as mdates
+
+
+cforcing = 'PAPA' ; # name of forcing ('PAPA', 'ERA5_arctic', etc...)
+
+cy1     = '2018' ; # First year
+cy2     = '2018' ; # Last year
+
+jt0 = 0
+
+dir_figs='.'
+size_fig=(13,8)
+size_fig0=(12,10)
+fig_ext='svg'
+
+clr_red = '#AD0000'
+clr_sat = '#ffed00'
+clr_mod = '#008ab8'
+
+rDPI=100.
+
+L_ALGOS = [ 'COARE3p6' , 'ECMWF'   , 'NCAR' , 'ANDREAS' ]
+l_color = [  '#ffed00' , '#008ab8' , '0.4'  , '0.8'     ] ; # colors to differentiate algos on the plot
+l_width = [     3      ,    2      ,  1     ,  3        ] ; # line-width to differentiate algos on the plot
+l_style = [    '-'     ,   '-'     , '--'   , '-'       ] ; # line-style
+
+
+# Variables to compare for A GIVEN algorithm
+###############################################
+#L_VNEM0 = [ 
+
+
+
+# Variables to compare between algorithms
+############################################
+L_VNEM  = [   'Cd_oce'  ,   'Ce_oce'  ,   'qla_oce' , 'qsb_oce'     ,     'qt_oce' ,   'qlw_oce' ,  'taum'     ,    'dt_skin'         ]
+L_VARO  = [     'Cd'    ,     'Ce'    ,   'Qlat'    ,    'Qsen'     ,     'Qnet'   ,   'Qlw'     ,  'Tau'      ,    'dT_skin'         ] ; # name of variable on figure
+L_VARL  = [ r'$C_{D}$'  , r'$C_{E}$'  , r'$Q_{lat}$', r'$Q_{sens}$' , r'$Q_{net}$' , r'$Q_{lw}$' , r'$|\tau|$' , r'$\Delta T_{skin}$' ] ; # name of variable in latex mode
+L_VUNT  = [     ''      ,     ''      , r'$W/m^2$'  , r'$W/m^2$'    , r'$W/m^2$'   , r'$W/m^2$'  , r'$N/m^2$'  ,      'K'             ]
+L_VMAX  = [    0.0075   ,    0.005    ,     75.     ,     75.       ,    800.      ,     25.     ,    1.2      ,       0.7            ]
+L_VMIN  = [    0.0005   ,    0.0005   ,   -250.     ,   -125.       ,   -400.      ,   -150.     ,    0.       ,      -0.7            ]
+L_ANOM  = [   False     ,   False     ,   True      ,    True       ,    True      ,    True     ,   True      ,      False           ]
+
+
+nb_algos = len(L_ALGOS) ; print(nb_algos)
+
+# Getting arguments:
+narg = len(sys.argv)
+if narg != 2:
+    print('Usage: '+sys.argv[0]+' <DIR_OUT_SASF>'); sys.exit(0)
+cdir_data = sys.argv[1]
+
+
+
+# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+# Populating and checking existence of files to be read
+# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+def chck4f(cf):
+    cmesg = 'ERROR: File '+cf+' does not exist !!!'
+    if not path.exists(cf): print(cmesg) ; sys.exit(0)
+
+###cf_in = nmp.empty((), dtype="S10")
+cf_in = []
+for ja in range(nb_algos):
+    cfi = cdir_data+'/output/'+'STATION_ASF-'+L_ALGOS[ja]+'_'+cforcing+'_1h_'+cy1+'0101_'+cy2+'1231_gridT.nc'
+    chck4f(cfi)
+    cf_in.append(cfi)
+print('Files we are goin to use:')
+for ja in range(nb_algos): print(cf_in[ja])
+#-----------------------------------------------------------------
+
+
+# Getting time array from the first file:
+id_in = Dataset(cf_in[0])
+vt = id_in.variables['time_counter'][jt0:]
+cunit_t = id_in.variables['time_counter'].units ; print(' "time_counter" is in "'+cunit_t+'"')
+id_in.close()
+nbr = len(vt)
+
+vtime = num2date(vt, units=cunit_t) ; # something understandable!                                                                 
+vtime = vtime.astype(dtype='datetime64[D]')
+
+ii=nbr/300
+ib=max(ii-ii%10,1)
+xticks_d=int(30*ib)
+
+
+rat = 100./float(rDPI)
+params = { 'font.family':'Open Sans',
+           'font.size':       int(15.*rat),
+           'legend.fontsize': int(15.*rat),
+           'xtick.labelsize': int(15.*rat),
+           'ytick.labelsize': int(15.*rat),
+           'axes.labelsize':  int(16.*rat)
+}
+mpl.rcParams.update(params)
+font_inf = { 'fontname':'Open Sans', 'fontweight':'normal', 'fontsize':18.*rat }
+font_x   = { 'fontname':'Open Sans', 'fontweight':'normal', 'fontsize':15.*rat }
+
+
+
+
+
+# First for each algorithm we compare some input vs out put variables:                                                                       
+
+# t_skin
+vtemp_in = [ 'sst'            ,               'theta_zu'                 ,               't_skin'                   ] ; #,               'theta_zt'                
+vtemp_lb = [ 'SST (bulk SST)' , r'$\theta_{zu}$ (pot. air temp. at 10m)' ,           'Skin temperature'             ] ; #, r'$\theta_{zt}$ (pot. air temp. at 2m)' 
+vtemp_cl = [  'k'             ,                clr_mod                   ,                clr_red                   ] ; #,                clr_mod                  
+vtemp_lw = [   3              ,                   1.3                    ,                   0.7                    ] ; #,                 2                       
+vtemp_ls = [     '-'          ,                   '-'                    ,                   '-'                    ] ; #,                 '-'                     
+
+ntemp = len(vtemp_in)
+
+xxx = nmp.zeros((nbr,ntemp))
+
+for ja in range(nb_algos):
+    #
+    # Temperatures...
+    id_in = Dataset(cf_in[ja])
+    for jv in range(ntemp):
+        xxx[:,jv] = id_in.variables[vtemp_in[jv]][jt0:,1,1] # only the center point of the 3x3 spatial domain!
+    id_in.close()
+
+    fig = plt.figure(num = 1, figsize=size_fig0, facecolor='w', edgecolor='k')
+    ax1 = plt.axes([0.07, 0.2, 0.9, 0.75])
+    ax1.set_xticks(vtime[::xticks_d])
+    ax1.xaxis.set_major_formatter(mdates.DateFormatter('%Y-%m-%d %H:%M:%S'))
+    plt.xticks(rotation='60', **font_x)
+
+    for jv in range(ntemp):
+        plt.plot(vtime, xxx[:,jv], '-', color=vtemp_cl[jv], linestyle=vtemp_ls[jv], linewidth=vtemp_lw[jv], label=vtemp_lb[jv], zorder=10)
+
+    ax1.set_ylim(0., 17.) ; ax1.set_xlim(vtime[0],vtime[nbr-1])
+    plt.ylabel(r'Temperature [$^{\circ}$C]')
+
+    ax1.grid(color='k', linestyle='-', linewidth=0.3)
+    plt.legend(loc='best', ncol=1, shadow=True, fancybox=True)
+    ax1.annotate('Algo: '+L_ALGOS[ja]+', station: '+cforcing, xy=(0.4, 1.), xycoords='axes fraction',  bbox={'facecolor':'w', 'alpha':1., 'pad':10}, zorder=50, **font_inf)
+    plt.savefig('01_temperatures_'+L_ALGOS[ja]+'.'+fig_ext, dpi=int(rDPI), transparent=False)
+    plt.close(1)
+
+
+del xxx
+
+
+
+# Now we compare output variables from bulk algorithms between them:
+
+nb_var = len(L_VNEM)
+
+xF  = nmp.zeros((nbr,nb_algos))
+xFa = nmp.zeros((nbr,nb_algos))
+
+
+for jv in range(nb_var):
+    print('\n *** Treating variable: '+L_VARO[jv]+' !')
+
+    for ja in range(nb_algos):
+        #
+        id_in = Dataset(cf_in[ja])
+        xF[:,ja] = id_in.variables[L_VNEM[jv]][jt0:,1,1] # only the center point of the 3x3 spatial domain!
+        if ja == 0: cvar_lnm = id_in.variables[L_VNEM[jv]].long_name
+        id_in.close()
+
+    fig = plt.figure(num = jv, figsize=size_fig, facecolor='w', edgecolor='k')
+    ax1 = plt.axes([0.08, 0.25, 0.9, 0.7])
+    ax1.set_xticks(vtime[::xticks_d])
+    ax1.xaxis.set_major_formatter(mdates.DateFormatter('%Y-%m-%d %H:%M:%S'))
+    plt.xticks(rotation='60', **font_x)
+
+    for ja in range(nb_algos):
+        plt.plot(vtime, xF[:,ja], '-', color=l_color[ja], linestyle=l_style[ja], linewidth=l_width[ja], label=L_ALGOS[ja], zorder=10+ja)
+
+    ax1.set_ylim(L_VMIN[jv], L_VMAX[jv]) ; ax1.set_xlim(vtime[0],vtime[nbr-1])
+    plt.ylabel(L_VARL[jv]+' ['+L_VUNT[jv]+']')
+
+    ax1.grid(color='k', linestyle='-', linewidth=0.3)
+    plt.legend(loc='best', ncol=1, shadow=True, fancybox=True)
+    ax1.annotate(cvar_lnm+', station: '+cforcing, xy=(0.3, 1.), xycoords='axes fraction',  bbox={'facecolor':'w', 'alpha':1., 'pad':10}, zorder=50, **font_inf)
+    plt.savefig(L_VARO[jv]+'.'+fig_ext, dpi=int(rDPI), transparent=False)
+    plt.close(jv)
+
+    if L_ANOM[jv]:
+
+        for ja in range(nb_algos): xFa[:,ja] = xF[:,ja] - nmp.mean(xF,axis=1)
+
+        if nmp.sum(xFa[:,:]) == 0.0:
+            print('     Well! Seems that for variable '+L_VARO[jv]+', choice of algo has no impact a all!')
+            print('          ==> skipping anomaly plot...')
+
+        else:
+
+            # Want a symetric y-range that makes sense for the anomaly we're looking at:
+            rmax = nmp.max(xFa) ; rmin = nmp.min(xFa)
+            rmax = max( abs(rmax) , abs(rmin) )
+            romagn = math.floor(math.log(rmax, 10)) ; # order of magnitude of the anomaly  we're dealing with
+            rmlt = 10.**(int(romagn)) / 2.
+            yrng = math.copysign( math.ceil(abs(rmax)/rmlt)*rmlt , rmax)
+
+            fig = plt.figure(num = 10+jv, figsize=size_fig, facecolor='w', edgecolor='k')
+            ax1 = plt.axes([0.08, 0.25, 0.9, 0.7])
+            ax1.set_xticks(vtime[::xticks_d])
+            ax1.xaxis.set_major_formatter(mdates.DateFormatter('%Y-%m-%d %H:%M:%S'))
+            plt.xticks(rotation='60', **font_x)
+
+            for ja in range(nb_algos):
+                plt.plot(vtime, xFa[:,ja], '-', color=l_color[ja], linewidth=l_width[ja], label=L_ALGOS[ja], zorder=10+ja)
+
+            ax1.set_ylim(-yrng,yrng) ; ax1.set_xlim(vtime[0],vtime[nbr-1])
+            plt.ylabel(L_VARL[jv]+' ['+L_VUNT[jv]+']')
+            ax1.grid(color='k', linestyle='-', linewidth=0.3)
+            plt.legend(bbox_to_anchor=(0.45, 0.2), ncol=1, shadow=True, fancybox=True)
+            ax1.annotate('Anomaly of '+cvar_lnm, xy=(0.3, 0.97), xycoords='axes fraction',  bbox={'facecolor':'w', 'alpha':1., 'pad':10}, zorder=50, **font_inf)
+            plt.savefig(L_VARO[jv]+'_anomaly.'+fig_ext, dpi=int(rDPI), transparent=False)
+            plt.close(10+jv)
+
+
+
+
Index: MO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/plot_station_asf_old.py
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/plot_station_asf_old.py	(revision 13675)
+++ 	(revision )
@@ -1,242 +1,0 @@
-#!/usr/bin/env python3
-# -*- Mode: Python; coding: utf-8; indent-tabs-mode: nil; tab-width: 4 -*-
-
-# Post-diagnostic of STATION_ASF /  L. Brodeau, 2020
-
-import sys
-from os import path as path
-import math
-import numpy as nmp
-from netCDF4 import Dataset
-import matplotlib as mpl
-mpl.use('Agg')
-import matplotlib.pyplot as plt
-import matplotlib.dates as mdates
-
-cy1     = '2018' ; # First year
-cy2     = '2018' ; # Last year
-
-dir_figs='.'
-size_fig=(13,7)
-fig_ext='png'
-
-clr_red = '#AD0000'
-clr_sat = '#ffed00'
-clr_mod = '#008ab8'
-
-rDPI=200.
-
-L_ALGOS = [ 'COARE3p6' , 'ECMWF'   , 'NCAR' ]
-l_xtrns = [ '-noskin'  , '-noskin' ,  ''    ] ; # string to add to algo name (L_ALGOS) to get version without skin params turned on
-l_color = [  '#ffed00' , '#008ab8' , '0.4'  ] ; # colors to differentiate algos on the plot
-l_width = [     3      ,    2      ,  1     ] ; # line-width to differentiate algos on the plot
-l_style = [    '-'     ,   '-'     , '--'   ] ; # line-style
-
-L_VNEM  = [   'qla'     ,     'qsb'     ,     'qt'     ,   'qlw'     ,  'taum'     ,    'dt_skin'         ]
-L_VARO  = [   'Qlat'    ,    'Qsen'     ,     'Qnet'   ,   'Qlw'     ,  'Tau'      ,    'dT_skin'         ] ; # name of variable on figure
-L_VARL  = [ r'$Q_{lat}$', r'$Q_{sens}$' , r'$Q_{net}$' , r'$Q_{lw}$' , r'$|\tau|$' , r'$\Delta T_{skin}$' ] ; # name of variable in latex mode
-L_VUNT  = [ r'$W/m^2$'  , r'$W/m^2$'    , r'$W/m^2$'   , r'$W/m^2$'  , r'$N/m^2$'  ,      'K'             ]
-L_VMAX  = [     75.     ,     75.       ,    800.      ,     25.     ,    1.2      ,       0.7            ]
-L_VMIN  = [   -250.     ,   -125.       ,   -400.      ,   -150.     ,    0.       ,      -0.7            ]
-L_ANOM  = [   True      ,    True       ,    True      ,    True     ,   True      ,      False           ]
-
-#L_VNEM  = [   'qlw'      ]
-#L_VARO  = [   'Qlw'      ] ; # name of variable on figure
-#L_VARL  = [ r'$Q_{lw}$'  ] ; # name of variable in latex mode
-#L_VUNT  = [ r'$W/m^2$'   ]
-#L_VMAX  = [     25.      ]
-#L_VMIN  = [   -150.      ]
-#L_ANOM  = [    True      ]
-
-
-
-nb_algos = len(L_ALGOS) ; print(nb_algos)
-
-# Getting arguments:
-narg = len(sys.argv)
-if narg != 2:
-    print('Usage: '+sys.argv[0]+' <DIR_OUT_SASF>'); sys.exit(0)
-cdir_data = sys.argv[1]
-
-
-
-# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-# Populating and checking existence of files to be read
-# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-def chck4f(cf):
-    cmesg = 'ERROR: File '+cf+' does not exist !!!'
-    if not path.exists(cf): print(cmesg); sys.exit(0)
-
-###cf_in = nmp.empty((), dtype="S10")
-cf_in = [] ; cf_in_ns = []
-for ja in range(nb_algos):
-    cfi = cdir_data+'/output/'+'STATION_ASF-'+L_ALGOS[ja]+'_1h_'+cy1+'0101_'+cy2+'1231_gridT.nc'
-    chck4f(cfi)
-    cf_in.append(cfi)
-# Same but without skin params:
-for ja in range(nb_algos):
-    cfi = cdir_data+'/output/'+'STATION_ASF-'+L_ALGOS[ja]+l_xtrns[ja]+'_1h_'+cy1+'0101_'+cy2+'1231_gridT.nc'
-    chck4f(cfi)
-    cf_in_ns.append(cfi)
-print('Files we are goin to use:')
-for ja in range(nb_algos): print(cf_in[ja])
-print('   --- same without cool-skin/warm-layer:')
-for ja in range(nb_algos): print(cf_in_ns[ja])
-#-----------------------------------------------------------------
-
-
-# Getting time array from the first file:
-id_in = Dataset(cf_in[0])
-vt = id_in.variables['time_counter'][:]
-cunit_t = id_in.variables['time_counter'].units ; print(' "time_counter" is in "'+cunit_t+'"')
-id_in.close()
-nbr = len(vt)
-
-
-
-
-vtime = nmp.zeros(nbr)
-
-vt = vt + 1036800. + 30.*60. # BUG!??? don't get why false in epoch to date conversion, and yet ncview gets it right!
-for jt in range(nbr): vtime[jt] = mdates.epoch2num(vt[jt])
-
-ii=nbr/300
-ib=max(ii-ii%10,1)
-xticks_d=int(30*ib)
-
-font_inf = { 'fontname':'Open Sans', 'fontweight':'normal', 'fontsize':14 }
-
-nb_var = len(L_VNEM)
-
-xF  = nmp.zeros((nbr,nb_algos))
-xFa = nmp.zeros((nbr,nb_algos))
-
-
-for ctest in ['skin','noskin']:
-
-    for jv in range(nb_var):
-        print('\n *** Treating variable: '+L_VARO[jv]+' ('+ctest+') !')
-
-        for ja in range(nb_algos):
-            #
-            if ctest == 'skin':   id_in = Dataset(cf_in[ja])
-            if ctest == 'noskin': id_in = Dataset(cf_in_ns[ja])
-            xF[:,ja] = id_in.variables[L_VNEM[jv]][:,1,1] # only the center point of the 3x3 spatial domain!
-            if ja == 0: cvar_lnm = id_in.variables[L_VNEM[jv]].long_name
-            id_in.close()
-
-        fig = plt.figure(num = jv, figsize=size_fig, facecolor='w', edgecolor='k')
-
-        ax1 = plt.axes([0.07, 0.22, 0.9, 0.75])
-
-        ax1.set_xticks(vtime[::xticks_d])
-        ax1.xaxis.set_major_formatter(mdates.DateFormatter('%Y-%m-%d %H:%M:%S'))
-        plt.xticks(rotation='60')
-
-        for ja in range(nb_algos):
-            plt.plot(vtime, xF[:,ja], '-', color=l_color[ja], linestyle=l_style[ja], linewidth=l_width[ja], label=L_ALGOS[ja], zorder=10+ja)
-
-        ax1.set_ylim(L_VMIN[jv], L_VMAX[jv]) ; ax1.set_xlim(vtime[0],vtime[nbr-1])
-        plt.ylabel(L_VARL[jv]+' ['+L_VUNT[jv]+']')
-
-        ax1.grid(color='k', linestyle='-', linewidth=0.3)
-        plt.legend(bbox_to_anchor=(0.45, 0.2), ncol=1, shadow=True, fancybox=True)
-        ax1.annotate(cvar_lnm+' ('+ctest+')', xy=(0.3, 0.97), xycoords='axes fraction',  bbox={'facecolor':'w', 'alpha':1., 'pad':10}, zorder=50, **font_inf)
-        plt.savefig(L_VARO[jv]+'_'+ctest+'.'+fig_ext, dpi=int(rDPI), transparent=False)
-        plt.close(jv)
-
-
-
-        if L_ANOM[jv]:
-
-            for ja in range(nb_algos): xFa[:,ja] = xF[:,ja] - nmp.mean(xF,axis=1)
-
-            if nmp.sum(xFa[:,:]) == 0.0:
-                print('     Well! Seems that for variable '+L_VARO[jv]+', choice of algo has no impact a all!')
-                print('          ==> skipping anomaly plot...')
-
-            else:
-
-                # Want a symetric y-range that makes sense for the anomaly we're looking at:
-                rmax = nmp.max(xFa) ; rmin = nmp.min(xFa)
-                rmax = max( abs(rmax) , abs(rmin) )
-                romagn = math.floor(math.log(rmax, 10)) ; # order of magnitude of the anomaly  we're dealing with
-                rmlt = 10.**(int(romagn)) / 2.
-                yrng = math.copysign( math.ceil(abs(rmax)/rmlt)*rmlt , rmax)
-
-                fig = plt.figure(num = 10+jv, figsize=size_fig, facecolor='w', edgecolor='k')
-                ax1 = plt.axes([0.07, 0.22, 0.9, 0.75])
-
-                ax1.set_xticks(vtime[::xticks_d])
-                ax1.xaxis.set_major_formatter(mdates.DateFormatter('%Y-%m-%d %H:%M:%S'))
-                plt.xticks(rotation='60')
-
-                for ja in range(nb_algos):
-                    plt.plot(vtime, xFa[:,ja], '-', color=l_color[ja], linewidth=l_width[ja], label=L_ALGOS[ja], zorder=10+ja)
-
-                ax1.set_ylim(-yrng,yrng) ; ax1.set_xlim(vtime[0],vtime[nbr-1])
-                plt.ylabel(L_VARL[jv]+' ['+L_VUNT[jv]+']')
-                ax1.grid(color='k', linestyle='-', linewidth=0.3)
-                plt.legend(bbox_to_anchor=(0.45, 0.2), ncol=1, shadow=True, fancybox=True)
-                ax1.annotate('Anomaly of '+cvar_lnm+' ('+ctest+')', xy=(0.3, 0.97), xycoords='axes fraction',  bbox={'facecolor':'w', 'alpha':1., 'pad':10}, zorder=50, **font_inf)
-                plt.savefig(L_VARO[jv]+'_'+ctest+'_anomaly.'+fig_ext, dpi=int(rDPI), transparent=False)
-                plt.close(10+jv)
-
-
-
-
-# Difference skin vs noskin:
-xFns = nmp.zeros((nbr,nb_algos))
-
-for jv in range(nb_var-1):
-    print('\n *** Treating variable: '+L_VARO[jv]+' ('+ctest+') !')
-
-    for ja in range(nb_algos-1):
-        id_in = Dataset(cf_in[ja])
-        xF[:,ja]   = id_in.variables[L_VNEM[jv]][:,1,1] # only the center point of the 3x3 spatial domain!
-        if ja == 0: cvar_lnm = id_in.variables[L_VNEM[jv]].long_name
-        id_in.close()
-        #
-        id_in = Dataset(cf_in_ns[ja])
-        xFns[:,ja] = id_in.variables[L_VNEM[jv]][:,1,1] # only the center point of the 3x3 spatial domain!
-        if ja == 0: cvar_lnm = id_in.variables[L_VNEM[jv]].long_name
-        id_in.close()
-
-        xFa[:,ja] = xF[:,ja] - xFns[:,ja]  ; # difference!
-
-
-    # Want a symetric y-range that makes sense for the anomaly we're looking at:
-    rmax = nmp.max(xFa) ; rmin = nmp.min(xFa)
-    rmax = max( abs(rmax) , abs(rmin) )
-    romagn = math.floor(math.log(rmax, 10)) ; # order of magnitude of the anomaly  we're dealing with
-    rmlt = 10.**(int(romagn)) / 2.
-    yrng = math.copysign( math.ceil(abs(rmax)/rmlt)*rmlt , rmax)
-
-    
-    for ja in range(nb_algos-1):
-
-        calgo = L_ALGOS[ja]
-
-        if nmp.sum(xFa[:,ja]) == 0.0:
-            print('     Well! Seems that for variable '+L_VARO[jv]+', and algo '+calgo+', skin param has no impact')
-            print('          ==> skipping difference plot...')
-
-        else:
-
-            fig = plt.figure(num = jv, figsize=size_fig, facecolor='w', edgecolor='k')
-            ax1 = plt.axes([0.07, 0.22, 0.9, 0.75])
-
-            ax1.set_xticks(vtime[::xticks_d])
-            ax1.xaxis.set_major_formatter(mdates.DateFormatter('%Y-%m-%d %H:%M:%S'))
-            plt.xticks(rotation='60')
-
-            plt.plot(vtime, xFa[:,ja], '-', color=l_color[ja], linestyle=l_style[ja], linewidth=l_width[ja], label=None, zorder=10+ja)
-
-            ax1.set_ylim(-yrng,yrng) ; ax1.set_xlim(vtime[0],vtime[nbr-1])
-            plt.ylabel(L_VARL[jv]+' ['+L_VUNT[jv]+']')
-
-            ax1.grid(color='k', linestyle='-', linewidth=0.3)
-            #plt.legend(bbox_to_anchor=(0.45, 0.2), ncol=1, shadow=True, fancybox=True)
-            ax1.annotate(cvar_lnm+' ('+ctest+')', xy=(0.3, 0.97), xycoords='axes fraction',  bbox={'facecolor':'w', 'alpha':1., 'pad':10}, zorder=50, **font_inf)
-            plt.savefig('diff_skin-noskin_'+L_VARO[jv]+'_'+calgo+'_'+ctest+'.'+fig_ext, dpi=int(rDPI), transparent=False)
-            plt.close(jv)
Index: /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/sbcblk_sanity_check.sh
===================================================================
--- /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/sbcblk_sanity_check.sh	(revision 13676)
+++ /NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/sbcblk_sanity_check.sh	(revision 13676)
@@ -0,0 +1,195 @@
+#!/bin/bash
+
+################################################################
+#
+# Script to launch a set of STATION_ASF simulations
+#
+#  Ocean-only version....
+#
+# L. Brodeau, 2020
+#
+#
+# Note: the current setup only performs an ocean-only test (no sea-ice),
+#       even though the current script already recongnises if the nemo.exe
+#       has been compiled with "key_si3". For now it will just test the
+#       fluxes over open water.
+#
+#
+#
+#
+#
+################################################################
+
+CONFIG="STATION_ASF" ; # general name of the configuration
+CONFIG_BLD="STATION_ASF2" ; # name of config as build in NEMO... (directory inside "tests" actually contains the compiled test-case?)
+
+# Root directory NEMOGCM reference depository where to fetch compiled STATION_ASF nemo.exe + default namelists:
+NEMO_REF_DIR="/home/laurent/DEV/NEMO/NEMOGCM_dev_r13648_ASINTER-04_laurent_bulk_ice"
+
+# NEMOGCM root directory where to fetch compiled STATION_ASF nemo.exe:
+SASF_WRK_DIR="${NEMO_REF_DIR}/tests/${CONFIG_BLD}"
+
+# DATA_IN_DIR => Directory containing sea-surface + atmospheric forcings:
+DATA_IN_DIR="${NEMO_REF_DIR}/tests/${CONFIG}/input_data"
+
+# Directory where to run the simulation:
+PROD_DIR="${HOME}/tmp/${CONFIG}"
+
+
+
+# Atmo + SSX forcing to use and wheter to compute fluxes over sea-ice as well (i_sea_ice=1)
+#############################################################################################
+i_sea_ice=0 ; # default does not compute fluxes over sea-ice
+NFORC="IDEALIZED" ; i_sea_ice=0 ; FORC="input_output_VALIDATION_IDEALIZED.nc"
+
+
+####### End of normal user configurable section #######
+
+#================================================================================
+
+rm -f SBCBLK.success SBCBLK.fail
+
+HERE=`pwd`
+
+echo
+
+# Is the nemo.exe of STATION_ASF compiled with sea-ice support (SI3) ?
+i_si3=0
+FCPP="${SASF_WRK_DIR}/cpp_${CONFIG_BLD}.fcm"
+if [ ! -f ${FCPP} ]; then echo " Mhhh, we did not find 'cpp_STATION_ASF.fcm' into `dirname ${FCPP}` !"; exit; fi
+ca=`cat ${FCPP} | grep 'key_si3'`
+
+if [ "${ca}" = "" ]; then
+    echo " *** NEMO was NOT compiled with sea-ice support (SI3) !" ; echo
+    if [ ${i_sea_ice} -eq 1 ]; then
+        echo " ===> so you cannot request ice-air flux computation !"
+        echo "     ===> please set i_sea_ice=0 or compile STATION_ASF with CPP key 'key_si3' !"
+        echo ; exit
+    fi
+else
+    echo " *** NEMO was apparently compiled with sea-ice support (SI3) !" ; echo
+    i_si3=1
+fi
+
+
+DIR_NL=${NFORC}/oce ; # directory where to find the namelists...
+
+# NEMO executable to use is:
+NEMO_EXE="${SASF_WRK_DIR}/BLD/bin/nemo.exe"
+if [ ! -f ${NEMO_EXE} ]; then echo " Mhhh, no compiled 'nemo.exe' found into `dirname ${NEMO_EXE}` !"; exit; fi
+
+SASF_EXPREF=`pwd`  ; # STATION_ASF EXPREF directory from which to use namelists and XIOS xml files...
+
+
+if [ ! -d ${DATA_IN_DIR} ]; then echo "PROBLEM!!! => did not find directory 'input_data' with input forcing..."; exit; fi
+
+cdt_cmpl="`\ls -l ${NEMO_EXE} | cut -d' ' -f 6,7,8`"
+
+echo "###########################################################"
+echo "#        S T A T I O N   A i r  -  S e a   F l u x        #"
+echo "###########################################################"
+echo
+echo "  * NEMO reference root directory is: ${NEMO_REF_DIR}"
+echo "  * STATION_ASF work directory is: ${SASF_WRK_DIR}"
+echo "       ==> NEMO EXE to use: ${NEMO_EXE}"
+echo "           ==> compiled: ${cdt_cmpl} !"
+echo
+echo "  * Input forcing data into: ${DATA_IN_DIR}"
+echo "  * Production will be done into: ${PROD_DIR}"
+echo "  * Directory in which namelists and xml files are fetched:"
+echo "       ==> ${SASF_EXPREF}"
+echo
+sleep 2
+
+mkdir -p ${PROD_DIR}
+
+rsync -avP ${NEMO_EXE}          ${PROD_DIR}/
+ln -sf     ${NEMO_EXE}          ${PROD_DIR}/nemo.exe.lnk
+
+
+# XIOS xml file
+################
+
+list_xml="iodef.xml file_def_nemo-oce.xml field_def_nemo-oce.xml"
+fcntxt="context_nemo_OCE.xml"
+if [ ${i_sea_ice} -eq 1 ]; then
+    list_xml+=" file_def_nemo-ice.xml field_def_nemo-ice.xml"
+    fcntxt="context_nemo_OCE+ICE.xml"
+fi
+
+# The "context_nemo.xml" file:
+if [ ! -f ${SASF_EXPREF}/${fcntxt} ]; then echo " Mhhh, ${fcntxt} not found into ${SASF_EXPREF} !"; exit; fi
+rsync -avPL ${SASF_EXPREF}/${fcntxt} ${PROD_DIR}/context_nemo.xml
+
+# All remaining "*.xml" files:
+for ff in ${list_xml} ; do
+    if [ ! -f ${SASF_EXPREF}/${ff} ]; then echo " Mhhh, ${ff} not found into ${SASF_EXPREF} !"; exit; fi
+    rsync -avPL ${SASF_EXPREF}/${ff} ${PROD_DIR}/
+done
+
+
+# Getting reference/defaults files from reference NEMO distro:
+rdir="${NEMO_REF_DIR}/cfgs/SHARED"
+for ff in "domain_def_nemo.xml" "grid_def_nemo.xml"; do
+    if [ ! -f ${rdir}/${ff} ]; then echo " Mhhh, ${ff} not found into ${rdir} !"; exit; fi
+    ln -sf ${rdir}/${ff} ${PROD_DIR}/.
+done
+
+
+# Copy forcing to work directory:
+echo ; echo "Forcing files to use:"
+\ls -l ${DATA_IN_DIR}/${FORC}
+echo
+rsync -avP ${DATA_IN_DIR}/${FORC} ${PROD_DIR}/
+echo; echo
+
+for CASE in "NCAR" "ECMWF" "COARE3p0" "COARE3p6" "ANDREAS"; do
+
+        echo ; echo
+        echo "======================================================================="
+        echo " Going for experiment: ${CASE} bulk param. with ${NFORC} forcing "
+        echo "======================================================================="
+        echo
+
+        scase=`echo "${CASE}" | tr '[:upper:]' '[:lower:]'`
+
+        rm -f ${PROD_DIR}/namelist_cfg
+        fnml="${SASF_EXPREF}/${DIR_NL}/namelist_${scase}_cfg"
+        if [ ! -f ${fnml} ]; then echo " Mhhh, test namelist ${fnml} not found !"; exit; fi
+
+        echo "   ===> namelist to use is: ${DIR_NL}/namelist_${scase}_cfg !"; echo
+        
+        # The namelists:
+        rsync -avPL ${SASF_EXPREF}/${DIR_NL}/namelist_${scase}_cfg ${PROD_DIR}/namelist_cfg
+        rsync -avPL ${SASF_EXPREF}/namelist_ref                   ${PROD_DIR}/namelist_ref
+    
+        cd ${PROD_DIR}/
+        echo
+        echo "Launching NEMO !"
+        ./nemo.exe 1>out_nemo.out 2>err_nemo.err
+        echo "Done!"
+        echo
+
+        # Moving output files:
+        mkdir -p output
+        mv -f ${CONFIG}-${CASE}_${NFORC}_*_grid*.nc output/
+
+        # Saving logs:
+        mkdir -p ${CASE}_${NFORC}_log
+        mv -f *.out *.err ocean.output output.namelist.dyn ${CASE}_${NFORC}_log/
+
+done
+
+
+
+
+# Now we can compare the results with sanity range!
+
+echo ; echo ; echo
+echo " *** Now time for sanity-check of heat flux and wind stress components ***"
+echo
+
+cd ${HERE}
+python3 ./analyze_output.py ${DATA_IN_DIR}/${FORC} ${PROD_DIR}/output
+
+
