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

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 11835 – NEMO

Changeset 11835


Ignore:
Timestamp:
2019-10-29T19:21:00+01:00 (4 years ago)
Author:
mathiot
Message:

ENHANCE-02_ISF_nemo: change required to pass sette test

Location:
NEMO/branches/2019/ENHANCE-02_ISF_nemo
Files:
16 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/cfgs/SHARED/field_def_nemo-oce.xml

    r11521 r11835  
    264264 
    265265          <!-- * variable related to ice shelf forcing * --> 
    266           <field id="isfdebug3d"                                                                                grid_ref="grid_T_3D" /> 
     266          <field id="isfdebug3d2"                                                                                grid_ref="grid_T_3D" /> 
     267          <field id="isfdebug3d1"                                                                                grid_ref="grid_T_3D" /> 
    267268          <field id="isfdebug2d"                                                                                /> 
    268269          <field id="fwfisf_cav"      long_name="Ice shelf melting"                             unit="kg/m2/s"  /> 
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/cfgs/SHARED/namelist_ref

    r11529 r11835  
    11981198&nam_diatmb    !  Top Middle Bottom Output                              (default: OFF) 
    11991199!----------------------------------------------------------------------- 
    1200    ln_diatmb   = .false.   !  Choose Top Middle and Bottom output or not 
     1200   ln_diatmb   = .true.   !  Choose Top Middle and Bottom output or not 
    12011201/ 
    12021202!----------------------------------------------------------------------- 
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/ISOMIP/EXPREF/namelist_cfg

    r10075 r11835  
    4646   rn_rdt      = 1800.     !  time step for the dynamics (and tracer if nn_acc=0) 
    4747/ 
     48!----------------------------------------------------------------------- 
     49&namcfg        !   parameters of the configuration                      (default: use namusr_def in namelist_cfg) 
     50!----------------------------------------------------------------------- 
     51/ 
     52!----------------------------------------------------------------------- 
     53&namtsd        !    Temperature & Salinity Data  (init/dmp)             (default: OFF) 
     54!----------------------------------------------------------------------- 
     55    
     56/ 
     57!----------------------------------------------------------------------- 
     58&namwad        !   Wetting and Drying (WaD)                             (default: OFF) 
     59!----------------------------------------------------------------------- 
     60/ 
     61!----------------------------------------------------------------------- 
     62&namcrs        !   coarsened grid (for outputs and/or TOP)              (ln_crs =T) 
     63!----------------------------------------------------------------------- 
     64/ 
     65!----------------------------------------------------------------------- 
     66&namc1d        !   1D configuration options                             ("key_c1d" default: PAPA station) 
     67!----------------------------------------------------------------------- 
     68/ 
     69!----------------------------------------------------------------------- 
     70&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d" default: OFF) 
     71!----------------------------------------------------------------------- 
     72/ 
     73!----------------------------------------------------------------------- 
     74&namc1d_uvd    !   data: U & V currents                                 ("key_c1d" default: OFF) 
     75!----------------------------------------------------------------------- 
     76 
     77/ 
     78 
    4879!!====================================================================== 
    4980!!            ***  Surface Boundary Condition namelists  ***          !! 
     
    5990!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T) 
    6091!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T) 
    61 !!   namsbc_isf      ice shelf melting/freezing                         (ln_isfcav  =T : read (ln_read_cfg=T) or set or usr_def_zgr ) 
    62 !!   namsbc_iscpl    coupling option between land ice model and ocean   (ln_isfcav  =T) 
    6392!!   namsbc_wave     external fields from wave model                    (ln_wave    =T) 
    6493!!   namberg         iceberg floats                                     (ln_icebergs=T) 
     
    6695! 
    6796!----------------------------------------------------------------------- 
    68 &namsbc        !   Surface Boundary Condition (surface module) 
    69 !----------------------------------------------------------------------- 
    70    nn_fsbc     = 1         !  frequency of surface boundary condition computation 
     97&namsbc        !   Surface Boundary Condition manager                   (default: NO selection) 
     98!----------------------------------------------------------------------- 
     99   nn_fsbc     = 1         !  frequency of SBC module call 
     100      !                    !  (control sea-ice & iceberg model call) 
    71101   ln_usr      = .true.    !  user defined formulation                  (T => check usrdef_sbc) 
    72    ln_isf      = .true.    !  ice shelf melting/freezing                (T => fill namsbc_isf) 
    73 / 
    74 !----------------------------------------------------------------------- 
    75 &namsbc_isf    !  Top boundary layer (ISF)                              (ln_isfcav =T : read (ln_read_cfg=T)  
    76 !-----------------------------------------------------------------------             or set or usr_def_zgr ) 
    77 !              ! file name ! frequency (hours) ! variable ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    78 !              !           !  (if <0  months)  !   name   !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
    79 ! nn_isf == 4 
    80    sn_fwfisf   = 'rnfisf'  ,         -12       ,'sowflisf',     .false.    , .true.  , 'yearly'  ,  ''      ,   ''     , '' 
    81 ! nn_isf == 3 
    82    sn_rnfisf   = 'rnfisf'  ,         -12       ,'sofwfisf',     .false.    , .true.  , 'yearly'  ,  ''      ,   ''     , '' 
    83 ! nn_isf == 2 and 3 
    84    sn_depmax_isf = 'rnfisf' ,        -12       ,'sozisfmax' ,   .false.    , .true.  , 'yearly'  ,  ''      ,   ''     , '' 
    85    sn_depmin_isf = 'rnfisf' ,        -12       ,'sozisfmin' ,   .false.    , .true.  , 'yearly'  ,  ''      ,   ''     , '' 
    86 ! nn_isf == 2 
    87    sn_Leff_isf = 'rnfisf'  ,         -12       ,'Leff'    ,     .false.    , .true.  , 'yearly'  ,  ''      ,   ''     , '' 
    88 ! for all case 
    89    nn_isf      = 1         !  ice shelf melting/freezing 
    90                            !  1 = presence of ISF    2 = bg03 parametrisation  
    91                            !  3 = rnf file for isf   4 = ISF fwf specified 
    92                            !  option 1 and 4 need ln_isfcav = .true. (domzgr) 
    93 ! only for nn_isf = 1 or 2 
    94    rn_gammat0  = 1.0e-4   ! gammat coefficient used in blk formula 
    95    rn_gammas0  = 1.0e-4   ! gammas coefficient used in blk formula 
    96 ! only for nn_isf = 1 or 4 
    97    rn_hisf_tbl =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
    98                           ! 0 => thickness of the tbl = thickness of the first wet cell 
    99 ! only for nn_isf = 1 
    100    nn_isfblk   = 1        ! 1 ISOMIP  like: 2 equations formulation (Hunter et al., 2006) 
    101                           ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) 
    102    nn_gammablk = 0        ! 0 = cst Gammat (= gammat/s) 
    103                           ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
    104                           ! 2 = velocity and stability dependent Gamma    (Holland et al. 1999) 
    105 / 
    106 !----------------------------------------------------------------------- 
    107 &namsbc_iscpl  !   land ice / ocean coupling option                      
    108 !----------------------------------------------------------------------- 
    109 / 
     102/ 
     103!----------------------------------------------------------------------- 
     104&namsbc_flx    !   surface boundary condition : flux formulation        (ln_flx =T) 
     105!----------------------------------------------------------------------- 
     106/ 
     107!----------------------------------------------------------------------- 
     108&namsbc_blk    !   namsbc_blk  generic Bulk formula                     (ln_blk =T) 
     109!----------------------------------------------------------------------- 
     110 
     111/ 
     112!----------------------------------------------------------------------- 
     113&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3") 
     114!----------------------------------------------------------------------- 
     115 
     116/ 
     117!----------------------------------------------------------------------- 
     118&namsbc_sas    !   Stand-Alone Surface module: ocean data               (SAS_SRC  only) 
     119!----------------------------------------------------------------------- 
     120 
     121/ 
     122!----------------------------------------------------------------------- 
     123&namsbc_iif    !   Ice-IF : use observed ice cover                      (nn_ice = 1) 
     124!----------------------------------------------------------------------- 
     125/ 
     126!----------------------------------------------------------------------- 
     127&namtra_qsr    !   penetrative solar radiation                          (ln_traqsr =T) 
     128!----------------------------------------------------------------------- 
     129    
     130/ 
     131!----------------------------------------------------------------------- 
     132&namsbc_ssr    !   surface boundary condition : sea surface restoring   (ln_ssr =T) 
     133!----------------------------------------------------------------------- 
     134 
     135/ 
     136!----------------------------------------------------------------------- 
     137&namsbc_rnf    !   runoffs                                              (ln_rnf =T) 
     138!----------------------------------------------------------------------- 
     139 
     140/ 
     141!----------------------------------------------------------------------- 
     142&namsbc_apr    !   Atmospheric pressure used as ocean forcing           (ln_apr_dyn =T) 
     143!----------------------------------------------------------------------- 
     144 
     145/ 
     146!----------------------------------------------------------------------- 
     147&namisf       !  Top boundary layer (ISF)                               (default: OFF) 
     148!----------------------------------------------------------------------- 
     149   ! 
     150   ! ---------------- ice shelf load ------------------------------- 
     151   ! 
     152   cn_isfload = 'isomip'      ! scheme to compute ice shelf load (ln_isfcav = .true. in domain_cfg.nc) 
     153   ! 
     154   ! ---------------- ice shelf melt formulation ------------------------------- 
     155   ! 
     156   ln_isf = .true.           ! activate ice shelf module 
     157      cn_isfdir = './'           ! directory for all ice shelf input file 
     158      ! 
     159      ! ---------------- cavities opened ------------------------------- 
     160      ! 
     161      ln_isfcav_mlt = .true.    ! ice shelf melting into the cavity (need ln_isfcav = .true. in domain_cfg.nc) 
     162         cn_isfcav_mlt = '2eq'   ! ice shelf melting formulation (spe/2eq/3eq/oasis) 
     163         !                       ! spe = fwfisf is read from a forcing field 
     164         !                       ! 2eq = ISOMIP  like: 2 equations formulation (Hunter et al., 2006) 
     165         !                       ! 3eq = ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) 
     166         !                       ! oasis = fwfisf is given by oasis and pattern by file sn_isfcav_fwf 
     167         !              !  cn_isfcav_mlt = 2eq or 3eq cases: 
     168         cn_gammablk = 'spe'    ! scheme to compute gammat/s (spe,ad15,hj99) 
     169         !                       ! ad15 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
     170         !                       ! hj99 = velocity and stability dependent Gamma    (Holland et al. 1999) 
     171         rn_gammat0  = 1.e-4     ! gammat coefficient used in blk formula 
     172         rn_gammas0  = 1.e-4     ! gammas coefficient used in blk formula 
     173         ! 
     174         rn_htbl     =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
     175         !                       ! 0 => thickness of the tbl = thickness of the first wet cell 
     176         ! 
     177/ 
     178!----------------------------------------------------------------------- 
     179&namsbc_wave   ! External fields from wave model                        (ln_wave=T) 
     180!----------------------------------------------------------------------- 
     181/ 
     182!----------------------------------------------------------------------- 
     183&namberg       !   iceberg parameters                                   (default: OFF) 
     184!----------------------------------------------------------------------- 
     185 
     186/ 
     187 
    110188!!====================================================================== 
    111189!!               ***  Lateral boundary condition  ***                 !! 
     
    124202   rn_shlat    =    0.     !  free slip 
    125203/ 
     204!----------------------------------------------------------------------- 
     205&namagrif      !  AGRIF zoom                                            ("key_agrif") 
     206!----------------------------------------------------------------------- 
     207/ 
     208!----------------------------------------------------------------------- 
     209&nam_tide      !   tide parameters                                      (default: OFF) 
     210!----------------------------------------------------------------------- 
     211/ 
     212!----------------------------------------------------------------------- 
     213&nambdy        !  unstructured open boundaries                          (default: OFF) 
     214!----------------------------------------------------------------------- 
     215/ 
     216!----------------------------------------------------------------------- 
     217&nambdy_dta    !  open boundaries - external data                       (see nam_bdy) 
     218!----------------------------------------------------------------------- 
     219 
     220/ 
     221!----------------------------------------------------------------------- 
     222&nambdy_tide   !  tidal forcing at open boundaries                      (default: OFF) 
     223!----------------------------------------------------------------------- 
     224/ 
     225 
    126226!!====================================================================== 
    127227!!                ***  Top/Bottom boundary condition  ***             !! 
     
    145245/ 
    146246!----------------------------------------------------------------------- 
    147 &namdrg_top    !   TOP friction                                         (ln_isfcav=T) 
     247&namdrg_top    !   TOP friction                                         (ln_OFF =F & ln_isfcav=T) 
    148248!----------------------------------------------------------------------- 
    149249   rn_Cd0     =  2.5e-3    !  drag coefficient [-] 
     
    156256/ 
    157257!----------------------------------------------------------------------- 
    158 &namdrg_bot    !   BOTTOM friction                                   
     258&namdrg_bot    !   BOTTOM friction                                      (ln_OFF =F) 
    159259!----------------------------------------------------------------------- 
    160260   rn_Cd0     =  1.e-3    !  drag coefficient [-] 
     
    166266      rn_boost =  50.         !  local boost factor  [-] 
    167267/ 
     268!----------------------------------------------------------------------- 
     269&nambbc        !   bottom temperature boundary condition                (default: OFF) 
     270!----------------------------------------------------------------------- 
     271 
     272/ 
     273!----------------------------------------------------------------------- 
     274&nambbl        !   bottom boundary layer scheme                         (default: OFF) 
     275!----------------------------------------------------------------------- 
     276/ 
     277 
    168278!!====================================================================== 
    169279!!                        Tracer (T-S) namelists                      !! 
     
    178288! 
    179289!----------------------------------------------------------------------- 
    180 &nameos        !   ocean Equation Of Seawater                           (default: OFF) 
     290&nameos        !   ocean Equation Of Seawater                           (default: NO selection) 
    181291!----------------------------------------------------------------------- 
    182292   ln_eos80    = .true.          !  = Use EOS80 equation of state 
     
    199309      rn_Ld        = 10.e+3         !  lateral diffusive length   [m] 
    200310/ 
     311!----------------------------------------------------------------------- 
     312&namtra_mle    !   mixed layer eddy parametrisation (Fox-Kemper)       (default: OFF) 
     313!----------------------------------------------------------------------- 
     314/ 
     315!----------------------------------------------------------------------- 
     316&namtra_eiv    !   eddy induced velocity param.                         (default: OFF) 
     317!----------------------------------------------------------------------- 
     318/ 
     319!----------------------------------------------------------------------- 
     320&namtra_dmp    !   tracer: T & S newtonian damping                      (default: OFF) 
     321!----------------------------------------------------------------------- 
     322/ 
     323 
    201324!!====================================================================== 
    202325!!                      ***  Dynamics namelists  ***                  !! 
     
    212335! 
    213336!----------------------------------------------------------------------- 
     337&nam_vvl       !   vertical coordinate options                          (default: z-star) 
     338!----------------------------------------------------------------------- 
     339/ 
     340!----------------------------------------------------------------------- 
    214341&namdyn_adv    !   formulation of the momentum advection                (default: NO selection) 
    215342!----------------------------------------------------------------------- 
     
    218345/ 
    219346!----------------------------------------------------------------------- 
    220 &namdyn_vor    !   Vorticity / Coriolis scheme                          (default: OFF) 
     347&namdyn_vor    !   Vorticity / Coriolis scheme                          (default: NO selection) 
    221348!----------------------------------------------------------------------- 
    222349   ln_dynvor_ene = .true.  !  energy conserving scheme 
     
    241368      rn_Lv      = 10.e+3           !  lateral viscous length   [m] 
    242369/ 
     370!----------------------------------------------------------------------- 
     371&namdta_dyn    !   offline ocean input files                            (OFF_SRC only) 
     372!----------------------------------------------------------------------- 
     373 
     374/ 
     375 
    243376!!====================================================================== 
    244377!!                     vertical physics namelists                     !! 
     
    253386! 
    254387!----------------------------------------------------------------------- 
    255 &namzdf        !   vertical physics                                     (default: NO selection) 
    256 !----------------------------------------------------------------------- 
    257    !                       ! type of vertical closure 
     388&namzdf        !   vertical physics manager                             (default: NO selection) 
     389!----------------------------------------------------------------------- 
     390   !                       ! type of vertical closure (required) 
    258391   ln_zdfcst   = .true.       !  constant mixing 
     392   ! 
     393   !                       ! convection 
    259394   ln_zdfevd   = .true.       !  enhanced vertical diffusion 
    260       nn_evdm     =    1         ! apply on tracer (=0) or on tracer and momentum (=1) 
    261       rn_evd      =   0.1        ! mixing coefficient [m2/s] 
     395   rn_evd      =  0.1         ! mixing coefficient [m2/s] 
     396   ! 
    262397   !                       ! coefficients 
    263    rn_avm0     =   1.e-3     !  vertical eddy viscosity   [m2/s] 
    264    rn_avt0     =   5.e-5     !  vertical eddy diffusivity [m2/s] 
    265 / 
     398   rn_avm0     =   1.e-3      !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst) 
     399   rn_avt0     =   5.e-5      !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst) 
     400/ 
     401!----------------------------------------------------------------------- 
     402&namzdf_ric    !   richardson number dependent vertical diffusion       (ln_zdfric =T) 
     403!----------------------------------------------------------------------- 
     404/ 
     405!----------------------------------------------------------------------- 
     406&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke =T) 
     407!----------------------------------------------------------------------- 
     408/ 
     409!----------------------------------------------------------------------- 
     410&namzdf_gls    !   GLS vertical diffusion                               (ln_zdfgls =T) 
     411!----------------------------------------------------------------------- 
     412/ 
     413!----------------------------------------------------------------------- 
     414&namzdf_osm    !   OSM vertical diffusion                               (ln_zdfosm =T) 
     415!----------------------------------------------------------------------- 
     416/ 
     417!----------------------------------------------------------------------- 
     418&namzdf_iwm    !    internal wave-driven mixing parameterization        (ln_zdfiwm =T) 
     419!----------------------------------------------------------------------- 
     420/ 
     421 
    266422!!====================================================================== 
    267423!!                  ***  Diagnostics namelists  ***                   !! 
     
    280436!!====================================================================== 
    281437! 
     438!----------------------------------------------------------------------- 
     439&namtrd        !   trend diagnostics                                    (default: OFF) 
     440!----------------------------------------------------------------------- 
     441/ 
     442!----------------------------------------------------------------------- 
     443&namptr        !   Poleward Transport Diagnostic                        (default: OFF) 
     444!----------------------------------------------------------------------- 
     445/ 
     446!----------------------------------------------------------------------- 
     447&namhsb        !  Heat and salt budgets                                 (default: OFF) 
     448!----------------------------------------------------------------------- 
     449/ 
     450!----------------------------------------------------------------------- 
     451&namdiu        !   Cool skin and warm layer models                      (default: OFF) 
     452!----------------------------------------------------------------------- 
     453/ 
     454!----------------------------------------------------------------------- 
     455&namflo        !   float parameters                                     ("key_float") 
     456!----------------------------------------------------------------------- 
     457/ 
     458!----------------------------------------------------------------------- 
     459&nam_diaharm   !   Harmonic analysis of tidal constituents              ("key_diaharm") 
     460!----------------------------------------------------------------------- 
     461/ 
     462!----------------------------------------------------------------------- 
     463&namdct        ! transports through some sections                       ("key_diadct") 
     464!----------------------------------------------------------------------- 
     465/ 
     466!----------------------------------------------------------------------- 
     467&nam_diatmb    !  Top Middle Bottom Output                              (default: OFF) 
     468!----------------------------------------------------------------------- 
     469/ 
     470!----------------------------------------------------------------------- 
     471&nam_dia25h    !  25h Mean Output                                       (default: OFF) 
     472!----------------------------------------------------------------------- 
     473/ 
     474!----------------------------------------------------------------------- 
     475&namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4") 
     476!----------------------------------------------------------------------- 
     477/ 
     478 
    282479!!====================================================================== 
    283480!!               ***  Observation & Assimilation  ***                 !! 
     
    287484!!====================================================================== 
    288485! 
     486!----------------------------------------------------------------------- 
     487&namobs        !  observation usage switch                              (default: OFF) 
     488!----------------------------------------------------------------------- 
     489/ 
     490!----------------------------------------------------------------------- 
     491&nam_asminc    !   assimilation increments                              ('key_asminc') 
     492!----------------------------------------------------------------------- 
     493/ 
     494 
    289495!!====================================================================== 
    290496!!                  ***  Miscellaneous namelists  ***                 !! 
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT2_flux_cen2_cfg

    r10075 r11835  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    ! 
    35    rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3634   ! 
    3735   rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT2_flux_ubs_cfg

    r10075 r11835  
    4646!----------------------------------------------------------------------- 
    4747   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    48    ! 
    49    rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    5048   ! 
    5149   rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT2_vect_eenH_cfg

    r10075 r11835  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    ! 
    35    rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3634   ! 
    3735   rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT2_vect_een_cfg

    r10075 r11835  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    ! 
    35    rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3634   ! 
    3735   rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT2_vect_ene_cfg

    r10075 r11835  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    ! 
    35    rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3634   ! 
    3735   rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT2_vect_ens_cfg

    r10075 r11835  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    ! 
    35    rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3634   ! 
    3735   rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT4_flux_cen2_cfg

    r10075 r11835  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    ! 
    35    rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3634   ! 
    3735   rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT4_flux_ubs_cfg

    r10075 r11835  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    ! 
    35    rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3634   ! 
    3735   rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT4_vect_eenH_cfg

    r10075 r11835  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    ! 
    35    rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3634   ! 
    3735   rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT4_vect_een_cfg

    r10075 r11835  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    ! 
    35    rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3634   ! 
    3735   rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT4_vect_ene_cfg

    r10075 r11835  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    ! 
    35    rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3634   ! 
    3735   rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT4_vect_ens_cfg

    r10075 r11835  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    ! 
    35    rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3634   ! 
    3735   rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/VORTEX/MY_SRC/domvvl.F90

    r10572 r11835  
    3636 
    3737   PUBLIC  dom_vvl_init       ! called by domain.F90 
     38   PUBLIC  dom_vvl_zgr        ! called by isfcpl.F90 
    3839   PUBLIC  dom_vvl_sf_nxt     ! called by step.F90 
    3940   PUBLIC  dom_vvl_sf_swp     ! called by step.F90 
     
    115116      !! Reference  : Leclair, M., and G. Madec, 2011, Ocean Modelling. 
    116117      !!---------------------------------------------------------------------- 
     118      ! 
     119      IF(lwp) WRITE(numout,*) 
     120      IF(lwp) WRITE(numout,*) 'dom_vvl_init : Variable volume activated' 
     121      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 
     122      ! 
     123      CALL dom_vvl_ctl     ! choose vertical coordinate (z_star, z_tilde or layer) 
     124      ! 
     125      !                    ! Allocate module arrays 
     126      IF( dom_vvl_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'dom_vvl_init : unable to allocate arrays' ) 
     127      ! 
     128      !                    ! Read or initialize e3t_(b/n), tilde_e3t_(b/n) and hdiv_lf 
     129      CALL dom_vvl_rst( nit000, 'READ' ) 
     130      e3t_a(:,:,jpk) = e3t_0(:,:,jpk)  ! last level always inside the sea floor set one for all 
     131      ! 
     132      CALL dom_vvl_zgr() ! interpolation scale factor, depth and water column 
     133      ! 
     134   END SUBROUTINE dom_vvl_init 
     135   ! 
     136   SUBROUTINE dom_vvl_zgr() 
     137      !!---------------------------------------------------------------------- 
     138      !!                ***  ROUTINE dom_vvl_init  *** 
     139      !!                    
     140      !! ** Purpose :  Interpolation of all scale factors,  
     141      !!               depths and water column heights 
     142      !! 
     143      !! ** Method  :  - interpolate scale factors 
     144      !! 
     145      !! ** Action  : - e3t_(n/b) and tilde_e3t_(n/b) 
     146      !!              - Regrid: e3(u/v)_n 
     147      !!                        e3(u/v)_b        
     148      !!                        e3w_n            
     149      !!                        e3(u/v)w_b       
     150      !!                        e3(u/v)w_n       
     151      !!                        gdept_n, gdepw_n and gde3w_n 
     152      !!              - h(t/u/v)_0 
     153      !!              - frq_rst_e3t and frq_rst_hdv 
     154      !! 
     155      !! Reference  : Leclair, M., and G. Madec, 2011, Ocean Modelling. 
     156      !!---------------------------------------------------------------------- 
    117157      INTEGER ::   ji, jj, jk 
    118158      INTEGER ::   ii0, ii1, ij0, ij1 
    119159      REAL(wp)::   zcoef 
    120160      !!---------------------------------------------------------------------- 
    121       ! 
    122       IF(lwp) WRITE(numout,*) 
    123       IF(lwp) WRITE(numout,*) 'dom_vvl_init : Variable volume activated' 
    124       IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 
    125       ! 
    126       CALL dom_vvl_ctl     ! choose vertical coordinate (z_star, z_tilde or layer) 
    127       ! 
    128       !                    ! Allocate module arrays 
    129       IF( dom_vvl_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'dom_vvl_init : unable to allocate arrays' ) 
    130       ! 
    131       !                    ! Read or initialize e3t_(b/n), tilde_e3t_(b/n) and hdiv_lf 
    132       CALL dom_vvl_rst( nit000, 'READ' ) 
    133       e3t_a(:,:,jpk) = e3t_0(:,:,jpk)  ! last level always inside the sea floor set one for all 
    134161      ! 
    135162      !                    !== Set of all other vertical scale factors  ==!  (now and before) 
     
    263290      ENDIF 
    264291      ! 
    265    END SUBROUTINE dom_vvl_init 
    266  
     292   END SUBROUTINE dom_vvl_zgr 
    267293 
    268294   SUBROUTINE dom_vvl_sf_nxt( kt, kcall )  
     
    10341060      ! 
    10351061      IF( ioptio /= 1 )   CALL ctl_stop( 'Choose ONE vertical coordinate in namelist nam_vvl' ) 
    1036       IF( .NOT. ln_vvl_zstar .AND. ln_isf ) CALL ctl_stop( 'Only vvl_zstar has been tested with ice shelf cavity' ) 
    10371062      ! 
    10381063      IF(lwp) THEN                   ! Print the choice 
Note: See TracChangeset for help on using the changeset viewer.