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 9019 for branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER – NEMO

Ignore:
Timestamp:
2017-12-13T15:58:53+01:00 (7 years ago)
Author:
timgraham
Message:

Merge of dev_CNRS_2017 into branch

Location:
branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER
Files:
8 edited
6 copied

Legend:

Unmodified
Added
Removed

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER/EXP00/1_namelist_cfg

    r7821 r9019  
    1414   cn_exp      =  "SASBIPER"  !  experience name 
    1515   nn_it000    =       1      !  first time step 
    16    nn_itend    =    30300     !  last  time step (std 5475) 
     16   nn_itend    =    1500      !  last  time step (std 5475) 
     17   nn_istate   =       1      !  output the initial state (1) or not (0) 
    1718/ 
    1819!----------------------------------------------------------------------- 
    … …  
    2627&namdom        !   space and time domain (bathymetry, mesh, timestep) 
    2728!----------------------------------------------------------------------- 
    28    ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    29    ! 
    30    nn_msh      =    1      !  create (>0) a mesh file or not (=0) 
     29   ln_linssh   = .true.   !  =T  linear free surface  ==>>  model level are fixed in time 
     30   ! 
     31   nn_msh      =    1     !  create (>0) a mesh file or not (=0) 
    3132   ! 
    3233   rn_rdt      = 400.     !  time step for the dynamics (and tracer if nn_acc=0) 
    3334/ 
    3435!----------------------------------------------------------------------- 
    35 &namcrs        !   Grid coarsening for dynamics output and/or 
    36                !   passive tracer coarsened online simulations 
     36&namcrs        !   coarsened grid (for outputs and/or TOP)              (ln_crs =T) 
    3737!----------------------------------------------------------------------- 
    3838/ 
    … …  
    4545!----------------------------------------------------------------------- 
    4646   ln_usr      = .true.   !  user defined formulation                  (T => check usrdef_sbc) 
     47   ln_blk      = .false.  !  Bulk formulation                          (T => fill namsbc_blk ) 
    4748   nn_fsbc     =  1       !  frequency of surface boundary condition computation 
     49   nn_ice      =  2       !  sea-ice model 
    4850/ 
    4951!----------------------------------------------------------------------- 
    … …  
    5254/ 
    5355!----------------------------------------------------------------------- 
     56&namsbc_sas    !   Stand Alone Surface boundary condition 
     57!----------------------------------------------------------------------- 
     58   l_sasread   = .false.   ! Read fields in a file if .TRUE. , or initialize to 0. in sbcssm.F90 if .FALSE. 
     59/ 
     60!----------------------------------------------------------------------- 
    5461&namtra_qsr    !   penetrative solar radiation 
    5562!----------------------------------------------------------------------- 
    … …  
    7077&namberg       !   iceberg parameters 
    7178!----------------------------------------------------------------------- 
     79   ln_bergdia               = .false.               ! Calculate budgets 
    7280/ 
    7381!----------------------------------------------------------------------- 
    … …  
    7886/ 
    7987!----------------------------------------------------------------------- 
    80 &nambfr        !   bottom friction 
    81 !----------------------------------------------------------------------- 
     88&namagrif      !  AGRIF zoom                                            ("key_agrif") 
     89!----------------------------------------------------------------------- 
     90   nn_cln_update =    1    !  baroclinic update frequency 
     91/ 
     92!----------------------------------------------------------------------- 
     93&namdrg            !   top/bottom drag coefficient                      (default: NO selection) 
     94!----------------------------------------------------------------------- 
     95   ln_NONE    = .true.     !  free-slip       : Cd = 0                   
    8296/ 
    8397!----------------------------------------------------------------------- 
    … …  
    95109/ 
    96110!----------------------------------------------------------------------- 
    97 &namtra_adv    !   advection scheme for tracer 
     111&namtra_adv    !   advection scheme for tracer                          (default: NO selection) 
    98112!----------------------------------------------------------------------- 
    99113   ln_traadv_fct =  .true.    !  FCT scheme 
    100114      nn_fct_h   =  4               !  =2/4, horizontal 2nd / 4th order  
    101115      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
    102       nn_fct_zts =  0               !  > 1 , 2nd order FCT scheme with vertical sub-timestepping 
    103       !                             !        (number of sub-timestep = nn_fct_zts) 
    104116/ 
    105117!----------------------------------------------------------------------- 
    … …  
    155167&namtra_dmp    !   tracer: T & S newtonian damping                      (default: NO) 
    156168!----------------------------------------------------------------------- 
    157 !----------------------------------------------------------------------- 
    158 &namdyn_adv    !   formulation of the momentum advection 
    159 !----------------------------------------------------------------------- 
     169/ 
     170!----------------------------------------------------------------------- 
     171&namdyn_adv    !   formulation of the momentum advection                (default: NO selection) 
     172!----------------------------------------------------------------------- 
     173   ln_dynadv_NONE= .false. !  linear dynamics (no momentum advection) 
     174   ln_dynadv_vec = .true. !  vector form - 2nd centered scheme 
     175     nn_dynkeg     = 0        ! grad(KE) scheme: =0   C2  ;  =1   Hollingsworth correction 
     176   ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme 
     177   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme 
    160178/ 
    161179!----------------------------------------------------------------------- 
    … …  
    205223/ 
    206224!----------------------------------------------------------------------- 
    207 &namzdf        !   vertical physics 
    208 !----------------------------------------------------------------------- 
    209 / 
    210 !----------------------------------------------------------------------- 
    211 &namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke") 
    212 !----------------------------------------------------------------------- 
    213 / 
    214 !----------------------------------------------------------------------- 
    215 &namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm") 
    216 !----------------------------------------------------------------------- 
    217 / 
    218 !----------------------------------------------------------------------- 
    219 &namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx") 
    220 !----------------------------------------------------------------------- 
    221 / 
    222 !----------------------------------------------------------------------- 
    223 &nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
     225&namzdf        !   vertical physics                                     (default: NO selection) 
     226!----------------------------------------------------------------------- 
     227   !                       ! type of vertical closure 
     228   ln_zdfcst   = .true.       !  constant mixing 
     229/ 
     230!----------------------------------------------------------------------- 
     231&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke =T) 
     232!----------------------------------------------------------------------- 
     233/ 
     234!----------------------------------------------------------------------- 
     235&namzdf_iwm    !   tidal mixing parameterization                        (ln_zdfiwm =T) 
     236!----------------------------------------------------------------------- 
     237/ 
     238!----------------------------------------------------------------------- 
     239&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi") 
    224240!----------------------------------------------------------------------- 
    225241/ 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER/EXP00/1_namelist_ice_cfg

    r7821 r9019  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! LIM3 namelist:   
    3 !!              1 - Generic parameters                 (namicerun) 
    4 !!              2 - Diagnostics                        (namicediag) 
    5 !!              3 - Ice initialization                 (namiceini) 
    6 !!              4 - Ice discretization                 (namiceitd) 
    7 !!              5 - Ice dynamics and transport         (namicedyn) 
    8 !!              6 - Ice diffusion                      (namicehdf) 
    9 !!              7 - Ice thermodynamics                 (namicethd) 
    10 !!              8 - Ice salinity                       (namicesal) 
    11 !!              9 - Ice mechanical redistribution      (namiceitdme) 
     2!! ESIM namelist:   
     3!!              1 - Generic parameters                 (nampar) 
     4!!              2 - Ice thickness discretization       (namitd) 
     5!!              3 - Ice dynamics                       (namdyn) 
     6!!              4 - Ice ridging/rafting                (namdyn_rdgrft) 
     7!!              5 - Ice rheology                       (namdyn_rhg) 
     8!!              6 - Ice advection                      (namdyn_adv) 
     9!!              7 - Ice surface forcing                (namforcing) 
     10!!              8 - Ice thermodynamics                 (namthd) 
     11!!              9 - Ice heat diffusion                 (namthd_zdf) 
     12!!             10 - Ice lateral melting                (namthd_da) 
     13!!             11 - Ice growth in open water           (namthd_do) 
     14!!             12 - Ice salinity                       (namthd_sal) 
     15!!             13 - Ice melt ponds                     (namthd_pnd) 
     16!!             14 - Ice initialization                 (namini) 
     17!!             15 - Ice/snow albedos                   (namalb) 
     18!!             16 - Ice diagnostics                    (namdia) 
    1219!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    1320! 
    1421!------------------------------------------------------------------------------ 
    15 &namicerun     !   Generic parameters 
     22&nampar         !   Generic parameters 
    1623!------------------------------------------------------------------------------ 
    17    jpl              =    1          !  number of ice  categories 
    18    nlay_i           =    1          !  number of ice  layers 
    19    ln_limthd        =  .false.       !  ice thermo   (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    20    ln_limdyn        =  .true.       !  ice dynamics (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    21    nn_limdyn        =   0           !     (ln_limdyn=T) switch for ice dynamics 
    22                                     !      2: total 
    23                                     !      1: advection only (no diffusion, no ridging/rafting) 
    24                                     !      0: advection only (as 1 but with prescribed velocity, bypass rheology) 
    25    rn_uice          =   0.5    !     (nn_limdyn=0) ice u-velocity 
    26    rn_vice          =   0.0    !     (nn_limdyn=0) ice v-velocity 
     24   ln_icedyn        = .true.          !  ice dynamics (T) or not (F) 
     25   ln_icethd        = .false.         !  ice thermo   (T) or not (F) 
    2726/ 
    2827!------------------------------------------------------------------------------ 
    29 &namicediag    !   Diagnostics 
     28&namitd         !   Ice discretization 
     29!------------------------------------------------------------------------------ 
     30   rn_himin         =   0.1           !  minimum ice thickness (m) used in remapping 
     31/ 
     32!------------------------------------------------------------------------------ 
     33&namdyn         !   Ice dynamics 
     34!------------------------------------------------------------------------------ 
     35   ln_dynFULL       = .false.         !  dyn.: full ice dynamics               (rheology + advection + ridging/rafting + correction) 
     36   ln_dynRHGADV     = .true.          !  dyn.: no ridge/raft & no corrections  (rheology + advection) 
     37   ln_dynADV        = .false.         !  dyn.: only advection w prescribed vel.(rn_uvice + advection) 
     38      rn_uice       =   0.5           !        prescribed ice u-velocity 
     39      rn_vice       =   0.            !        prescribed ice v-velocity 
     40/ 
     41!------------------------------------------------------------------------------ 
     42&namdyn_rdgrft  !   Ice ridging/rafting 
    3043!------------------------------------------------------------------------------ 
    3144/ 
    3245!------------------------------------------------------------------------------ 
    33 &namiceini     !   Ice initialization 
    34 !------------------------------------------------------------------------------ 
    35                   ! -- limistate -- ! 
    36    ln_limini      = .false.         !  activate ice initialization (T) or not (F) 
    37    ln_limini_file = .false.         !  netcdf file provided for initialization (T) or not (F) 
    38    cn_dir="./" 
    39    sn_hti = 'initice'    , -12 ,'hti'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    40    sn_hts = 'initice'    , -12 ,'hts'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    41    sn_ati = 'initice'    , -12 ,'ati'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    42    sn_tsu = 'initice'    , -12 ,'tsu'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    43    sn_tmi = 'initice'    , -12 ,'tmi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    44    sn_smi = 'initice'    , -12 ,'smi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    45 / 
    46 !------------------------------------------------------------------------------ 
    47 &namiceitd     !   Ice discretization 
     46&namdyn_rhg     !   Ice rheology 
    4847!------------------------------------------------------------------------------ 
    4948/ 
    5049!------------------------------------------------------------------------------ 
    51 &namicedyn     !   Ice dynamics and transport 
     50&namdyn_adv     !   Ice advection 
    5251!------------------------------------------------------------------------------ 
    53                   ! -- limtrp & limadv -- ! 
    54    nn_limadv      =    0            !  choose the advection scheme (-1=Prather ; 0=Ultimate-Macho) 
    55    nn_limadv_ord  =    5            !  choose the order of the advection scheme (if nn_limadv=0) 
     52   ln_adv_Pra       = .false.         !  Advection scheme (Prather) 
     53   ln_adv_UMx       = .true.          !  Advection scheme (Ultimate-Macho) 
     54      nn_UMx        =   5             !     order of the scheme for UMx (1-5 ; 20=centered 2nd order) 
    5655/ 
    5756!------------------------------------------------------------------------------ 
    58 &namicehdf     !   Ice horizontal diffusion 
     57&namforcing     !   Ice surface forcing 
    5958!------------------------------------------------------------------------------ 
    6059/ 
    6160!------------------------------------------------------------------------------ 
    62 &namicethd     !   Ice thermodynamics 
     61&namthd         !   Ice thermodynamics 
    6362!------------------------------------------------------------------------------ 
    64                   ! -- limthd_dh -- ! 
    65    ln_limdH       = .true.          !  activate ice thickness change from growing/melting (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    66                   ! -- limthd_da -- ! 
    67    ln_limdA       = .true.          !  activate lateral melting param. (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    68                  ! -- limthd_lac -- ! 
    69    ln_limdO       = .true.          !  activate ice growth in open-water (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    70    rn_hnewice     = 0.02            !  thickness for new ice formation in open water (m) 
    71                   ! -- limitd_th -- ! 
    72    rn_himin       = 0.01            !  minimum ice thickness (m) used in remapping, must be smaller than rn_hnewice 
    7363/ 
    7464!------------------------------------------------------------------------------ 
    75 &namicesal     !   Ice salinity 
     65&namthd_zdf     !   Ice heat diffusion 
    7666!------------------------------------------------------------------------------ 
    77                  ! -- limthd_sal -- ! 
    78    ln_limdS       = .true.          !  activate gravity drainage and flushing (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    7967/ 
    8068!------------------------------------------------------------------------------ 
    81 &namiceitdme   !   Ice mechanical redistribution (ridging and rafting) 
     69&namthd_da      !   Ice lateral melting 
    8270!------------------------------------------------------------------------------ 
    83                   ! -- limitd_me -- ! 
    84    ln_ridging     =   .true.        !  ridging activated (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    85    ln_rafting     =   .true.        !  rafting activated (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    8671/ 
     72!------------------------------------------------------------------------------ 
     73&namthd_do      !   Ice growth in open water 
     74!------------------------------------------------------------------------------ 
     75/ 
     76!------------------------------------------------------------------------------ 
     77&namthd_sal     !   Ice salinity 
     78!------------------------------------------------------------------------------ 
     79/ 
     80!------------------------------------------------------------------------------ 
     81&namthd_pnd     !   Melt ponds 
     82!------------------------------------------------------------------------------ 
     83/ 
     84!------------------------------------------------------------------------------ 
     85&namini         !   Ice initialization 
     86!------------------------------------------------------------------------------ 
     87   ln_iceini        = .false.          !  activate ice initialization (T) or not (F) 
     88/ 
     89!------------------------------------------------------------------------------ 
     90&namalb         !   albedo parameters 
     91!------------------------------------------------------------------------------ 
     92/ 
     93!------------------------------------------------------------------------------ 
     94&namdia         !   Diagnostics 
     95!------------------------------------------------------------------------------ 
     96/ 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER/EXP00/file_def_nemo-lim.xml

    r7821 r9019  
    1717   <file id="file1" name_suffix="_icemod" description="ice variables" enabled=".true." > 
    1818 
    19         <field field_ref="snowthic_cea"     name="snthic"     long_name="surface_snow_thickness"   /> 
    20         <field field_ref="snowvol"          name="snvolu" /> 
    21         <field field_ref="isnowhc"          name="snheco" /> 
     19        <field field_ref="snwvolu"          name="snvolu" /> 
     20        <field field_ref="icethic"          name="sithic" /> 
     21        <field field_ref="icevolu"          name="sivolu" /> 
     22        <field field_ref="iceconc"          name="siconc" /> 
     23        <field field_ref="icesalt"          name="sisali" /> 
     24        <field field_ref="icetemp"          name="sitemp" /> 
     25        <field field_ref="icettop"          name="sittop" /> 
     26        <field field_ref="uice"             name="sivelu" /> 
     27        <field field_ref="vice"             name="sivelv" /> 
     28        <field field_ref="icevel"           name="sivelo" /> 
     29        <field field_ref="icediv"           name="sidive" /> 
     30        <field field_ref="iceshe"           name="sishea" /> 
     31        <field field_ref="icestr"           name="sistre" /> 
     32        <field field_ref="sst_m"            name="sst_m" /> 
    2233 
    23         <field field_ref="icethic_cea"      name="sithic"     long_name="sea_ice_thickness"        /> 
    24         <field field_ref="icevolu"          name="sivolu" /> 
    25         <field field_ref="iceconc"          name="siconc"  /> 
    26         <field field_ref="micesalt"         name="sisali" /> 
    27         <field field_ref="micet"            name="sitemp" /> 
    28         <field field_ref="icest"            name="sistem" /> 
    29         <field field_ref="icehc"            name="siheco" /> 
    30         <field field_ref="uice_ipa"         name="sivelu" /> 
    31         <field field_ref="vice_ipa"         name="sivelv" /> 
    32         <field field_ref="icevel"           name="sivelo" /> 
    33         <field field_ref="idive"            name="sidive" /> 
    34         <field field_ref="ishear"           name="sishea" /> 
    35         <field field_ref="icestr"           name="sistre" /> 
     34        <!-- fluxes --> 
     35        <field field_ref="qt_oce_ai"        name="qt_oce_ai" /> 
     36        <field field_ref="qt_atm_oi"        name="qt_atm_oi" /> 
     37        <field field_ref="hfxcndtop"        name="hfxcndtop" /> 
     38        <field field_ref="hfxcndbot"        name="hfxcndbot" /> 
     39        <field field_ref="hfxsensib"        name="hfxsensib" /> 
     40        <field field_ref="hfxbom"           name="hfxbom" /> 
     41        <field field_ref="hfxbog"           name="hfxbog" /> 
     42        <field field_ref="hfxsum"           name="hfxsum" /> 
     43        <field field_ref="hfxopw"           name="hfxopw" /> 
     44        <field field_ref="hfxdif"           name="hfxdif" /> 
     45        <field field_ref="hfxsnw"           name="hfxsnw" /> 
     46        <field field_ref="hfxerr"           name="hfxerr" /> 
     47 
    3648        <!-- categories --> 
    37         <field field_ref="snowthic_cat"     name="snthicat"/> 
     49        <field field_ref="snwthic_cat"      name="snthicat"/> 
    3850        <field field_ref="iceconc_cat"      name="siconcat"/> 
    3951        <field field_ref="icethic_cat"      name="sithicat"/> 
    40         <field field_ref="salinity_cat"     name="salincat"/> 
    41         <field field_ref="iceage_cat"       name="siagecat"/> 
    42         <field field_ref="brinevol_cat"     name="sibricat"/> 
    43         <field field_ref="icetemp_cat"      name="sitemcat"/> 
    44         <field field_ref="snwtemp_cat"      name="sntemcat"/> 
    4552 
    4653        </file> 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER/EXP00/make_INITICE.py

    r7830 r9019  
    136136# ---------------------------------------------- 
    137137# for basin=99x99km with dx=1km ; dy=1km + AGRIF 
    138 sigx=-0.04 
    139 sigy=-0.04 
     138sigx=-0.012 
     139sigy=-0.012 
    140140xshift=20.-1. 
    141141yshift=50.-1. 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER/EXP00/namelist_cfg

    r8599 r9019  
    1414   cn_exp      =  "SASBIPER"  !  experience name 
    1515   nn_it000    =       1      !  first time step 
    16    nn_itend    =    10100     !  last  time step (std 5475) 
     16   nn_itend    =     500      !  last  time step (std 5475) 
     17   nn_istate   =       1      !  output the initial state (1) or not (0) 
    1718/ 
    1819!----------------------------------------------------------------------- 
    … …  
    2627&namdom        !   space and time domain (bathymetry, mesh, timestep) 
    2728!----------------------------------------------------------------------- 
    28    ln_linssh   = .true.   !  =T  linear free surface  ==>>  model level are fixed in time 
     29   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time 
    2930   ! 
    3031   nn_msh      =    1      !  create (>0) a mesh file or not (=0) 
    … …  
    3334/ 
    3435!----------------------------------------------------------------------- 
    35 &namcrs        !   Grid coarsening for dynamics output and/or 
    36                !   passive tracer coarsened online simulations 
     36&namcrs        !   coarsened grid (for outputs and/or TOP)              (ln_crs =T) 
    3737!----------------------------------------------------------------------- 
    3838/ 
    … …  
    4747   ln_blk      = .false.  !  Bulk formulation                          (T => fill namsbc_blk ) 
    4848   nn_fsbc     =  1       !  frequency of surface boundary condition computation 
     49   nn_ice      =  2       !  sea-ice model 
    4950/ 
    5051!----------------------------------------------------------------------- 
    … …  
    5354/ 
    5455!----------------------------------------------------------------------- 
     56&namsbc_sas    !   Stand Alone Surface boundary condition 
     57!----------------------------------------------------------------------- 
     58   l_sasread   = .false.   ! Read fields in a file if .TRUE. , or initialize to 0. in sbcssm.F90 if .FALSE. 
     59/ 
     60!----------------------------------------------------------------------- 
    5561&namtra_qsr    !   penetrative solar radiation 
    5662!----------------------------------------------------------------------- 
    … …  
    6975/ 
    7076!----------------------------------------------------------------------- 
    71 &namsbc_sas    !   Stand Alone Surface boundary condition 
    72 !----------------------------------------------------------------------- 
    73    l_sasread   = .false.   ! Read fields in a file if .TRUE. , or initialize to 0. in sbcssm.F90 if .FALSE. 
    74 !----------------------------------------------------------------------- 
    7577&namberg       !   iceberg parameters 
    7678!----------------------------------------------------------------------- 
     79   ln_bergdia               = .false.               ! Calculate budgets 
    7780/ 
    7881!----------------------------------------------------------------------- 
    … …  
    8386/ 
    8487!----------------------------------------------------------------------- 
    85 &nambfr        !   bottom friction 
    86 !----------------------------------------------------------------------- 
     88&namagrif      !  AGRIF zoom                                            ("key_agrif") 
     89!----------------------------------------------------------------------- 
     90   nn_cln_update =    1    !  baroclinic update frequency 
     91/ 
     92!----------------------------------------------------------------------- 
     93&namdrg            !   top/bottom drag coefficient                      (default: NO selection) 
     94!----------------------------------------------------------------------- 
     95   ln_NONE    = .true.     !  free-slip       : Cd = 0                   
    8796/ 
    8897!----------------------------------------------------------------------- 
    … …  
    100109/ 
    101110!----------------------------------------------------------------------- 
    102 &namtra_adv    !   advection scheme for tracer 
     111&namtra_adv    !   advection scheme for tracer                          (default: NO selection) 
    103112!----------------------------------------------------------------------- 
    104113   ln_traadv_fct =  .true.    !  FCT scheme 
    105114      nn_fct_h   =  4               !  =2/4, horizontal 2nd / 4th order  
    106115      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
    107       nn_fct_zts =  0               !  > 1 , 2nd order FCT scheme with vertical sub-timestepping 
    108       !                             !        (number of sub-timestep = nn_fct_zts) 
    109116/ 
    110117!----------------------------------------------------------------------- 
    … …  
    160167&namtra_dmp    !   tracer: T & S newtonian damping                      (default: NO) 
    161168!----------------------------------------------------------------------- 
    162 !----------------------------------------------------------------------- 
    163 &namdyn_adv    !   formulation of the momentum advection 
    164 !----------------------------------------------------------------------- 
     169/ 
     170!----------------------------------------------------------------------- 
     171&namdyn_adv    !   formulation of the momentum advection                (default: NO selection) 
     172!----------------------------------------------------------------------- 
     173   ln_dynadv_NONE= .false. !  linear dynamics (no momentum advection) 
     174   ln_dynadv_vec = .true. !  vector form - 2nd centered scheme 
     175     nn_dynkeg     = 0        ! grad(KE) scheme: =0   C2  ;  =1   Hollingsworth correction 
     176   ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme 
     177   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme 
    165178/ 
    166179!----------------------------------------------------------------------- 
    … …  
    210223/ 
    211224!----------------------------------------------------------------------- 
    212 &namzdf        !   vertical physics 
    213 !----------------------------------------------------------------------- 
    214 / 
    215 !----------------------------------------------------------------------- 
    216 &namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke") 
    217 !----------------------------------------------------------------------- 
    218 / 
    219 !----------------------------------------------------------------------- 
    220 &namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm") 
    221 !----------------------------------------------------------------------- 
    222 / 
    223 !----------------------------------------------------------------------- 
    224 &namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx") 
    225 !----------------------------------------------------------------------- 
    226 / 
    227 !----------------------------------------------------------------------- 
    228 &nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
     225&namzdf        !   vertical physics                                     (default: NO selection) 
     226!----------------------------------------------------------------------- 
     227   !                       ! type of vertical closure 
     228   ln_zdfcst   = .true.       !  constant mixing 
     229/ 
     230!----------------------------------------------------------------------- 
     231&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke =T) 
     232!----------------------------------------------------------------------- 
     233/ 
     234!----------------------------------------------------------------------- 
     235&namzdf_iwm    !   tidal mixing parameterization                        (ln_zdfiwm =T) 
     236!----------------------------------------------------------------------- 
     237/ 
     238!----------------------------------------------------------------------- 
     239&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi") 
    229240!----------------------------------------------------------------------- 
    230241/ 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER/EXP00/namelist_ice_cfg

    r8599 r9019  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! LIM3 namelist:   
    3 !!              1 - Generic parameters                 (namicerun) 
    4 !!              2 - Diagnostics                        (namicediag) 
    5 !!              3 - Ice initialization                 (namiceini) 
    6 !!              4 - Ice discretization                 (namiceitd) 
    7 !!              5 - Ice dynamics and transport         (namicedyn) 
    8 !!              6 - Ice diffusion                      (namicehdf) 
    9 !!              7 - Ice thermodynamics                 (namicethd) 
    10 !!              8 - Ice salinity                       (namicesal) 
    11 !!              9 - Ice mechanical redistribution      (namiceitdme) 
     2!! ESIM namelist:   
     3!!              1 - Generic parameters                 (nampar) 
     4!!              2 - Ice thickness discretization       (namitd) 
     5!!              3 - Ice dynamics                       (namdyn) 
     6!!              4 - Ice ridging/rafting                (namdyn_rdgrft) 
     7!!              5 - Ice rheology                       (namdyn_rhg) 
     8!!              6 - Ice advection                      (namdyn_adv) 
     9!!              7 - Ice surface forcing                (namforcing) 
     10!!              8 - Ice thermodynamics                 (namthd) 
     11!!              9 - Ice heat diffusion                 (namthd_zdf) 
     12!!             10 - Ice lateral melting                (namthd_da) 
     13!!             11 - Ice growth in open water           (namthd_do) 
     14!!             12 - Ice salinity                       (namthd_sal) 
     15!!             13 - Ice melt ponds                     (namthd_pnd) 
     16!!             14 - Ice initialization                 (namini) 
     17!!             15 - Ice/snow albedos                   (namalb) 
     18!!             16 - Ice diagnostics                    (namdia) 
    1219!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    1320! 
    1421!------------------------------------------------------------------------------ 
    15 &namicerun     !   Generic parameters 
     22&nampar         !   Generic parameters 
    1623!------------------------------------------------------------------------------ 
    17    jpl              =    1          !  number of ice  categories 
    18    nlay_i           =    1          !  number of ice  layers 
    19    ln_limthd        =  .false.       !  ice thermo   (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    20    ln_limdyn        =  .true.       !  ice dynamics (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    21    nn_limdyn        =   0           !     (ln_limdyn=T) switch for ice dynamics 
    22                                     !      2: total 
    23                                     !      1: advection only (no diffusion, no ridging/rafting) 
    24                                     !      0: advection only (as 1 but with prescribed velocity, bypass rheology) 
    25    rn_uice          =   0.5    !     (nn_limdyn=0) ice u-velocity 
    26    rn_vice          =   0.0    !     (nn_limdyn=0) ice v-velocity 
     24   ln_icedyn        = .true.          !  ice dynamics (T) or not (F) 
     25   ln_icethd        = .false.         !  ice thermo   (T) or not (F) 
    2726/ 
    2827!------------------------------------------------------------------------------ 
    29 &namicediag    !   Diagnostics 
     28&namitd         !   Ice discretization 
     29!------------------------------------------------------------------------------ 
     30   rn_himin         =   0.1           !  minimum ice thickness (m) used in remapping 
     31/ 
     32!------------------------------------------------------------------------------ 
     33&namdyn         !   Ice dynamics 
     34!------------------------------------------------------------------------------ 
     35   ln_dynFULL       = .false.         !  dyn.: full ice dynamics               (rheology + advection + ridging/rafting + correction) 
     36   ln_dynRHGADV     = .true.          !  dyn.: no ridge/raft & no corrections  (rheology + advection) 
     37   ln_dynADV        = .false.         !  dyn.: only advection w prescribed vel.(rn_uvice + advection) 
     38      rn_uice       =   0.5           !        prescribed ice u-velocity 
     39      rn_vice       =   0.            !        prescribed ice v-velocity 
     40/ 
     41!------------------------------------------------------------------------------ 
     42&namdyn_rdgrft  !   Ice ridging/rafting 
    3043!------------------------------------------------------------------------------ 
    3144/ 
    3245!------------------------------------------------------------------------------ 
    33 &namiceini     !   Ice initialization 
    34 !------------------------------------------------------------------------------ 
    35                   ! -- limistate -- ! 
    36    ln_limini      = .true.         !  activate ice initialization (T) or not (F) 
    37    ln_limini_file = .true.         !  netcdf file provided for initialization (T) or not (F) 
    38    cn_dir="./" 
    39    sn_hti = 'initice'    , -12 ,'hti'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    40    sn_hts = 'initice'    , -12 ,'hts'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    41    sn_ati = 'initice'    , -12 ,'ati'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    42    sn_tsu = 'initice'    , -12 ,'tsu'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    43    sn_tmi = 'initice'    , -12 ,'tmi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    44    sn_smi = 'initice'    , -12 ,'smi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    45 / 
    46 !------------------------------------------------------------------------------ 
    47 &namiceitd     !   Ice discretization 
     46&namdyn_rhg     !   Ice rheology 
    4847!------------------------------------------------------------------------------ 
    4948/ 
    5049!------------------------------------------------------------------------------ 
    51 &namicedyn     !   Ice dynamics and transport 
     50&namdyn_adv     !   Ice advection 
    5251!------------------------------------------------------------------------------ 
    53                   ! -- limtrp & limadv -- ! 
    54    nn_limadv      =    0            !  choose the advection scheme (-1=Prather ; 0=Ultimate-Macho) 
    55    nn_limadv_ord  =    5            !  choose the order of the advection scheme (if nn_limadv=0) 
     52   ln_adv_Pra       = .false.         !  Advection scheme (Prather) 
     53   ln_adv_UMx       = .true.          !  Advection scheme (Ultimate-Macho) 
     54      nn_UMx        =   5             !     order of the scheme for UMx (1-5 ; 20=centered 2nd order) 
    5655/ 
    5756!------------------------------------------------------------------------------ 
    58 &namicehdf     !   Ice horizontal diffusion 
     57&namforcing     !   Ice surface forcing 
    5958!------------------------------------------------------------------------------ 
    6059/ 
    6160!------------------------------------------------------------------------------ 
    62 &namicethd     !   Ice thermodynamics 
     61&namthd         !   Ice thermodynamics 
    6362!------------------------------------------------------------------------------ 
    64                   ! -- limthd_dh -- ! 
    65    ln_limdH       = .true.          !  activate ice thickness change from growing/melting (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    66                   ! -- limthd_da -- ! 
    67    ln_limdA       = .true.          !  activate lateral melting param. (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    68                  ! -- limthd_lac -- ! 
    69    ln_limdO       = .true.          !  activate ice growth in open-water (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    70    rn_hnewice     = 0.02            !  thickness for new ice formation in open water (m) 
    71                   ! -- limitd_th -- ! 
    72    rn_himin       = 0.01            !  minimum ice thickness (m) used in remapping, must be smaller than rn_hnewice 
    7363/ 
    7464!------------------------------------------------------------------------------ 
    75 &namicesal     !   Ice salinity 
     65&namthd_zdf     !   Ice heat diffusion 
    7666!------------------------------------------------------------------------------ 
    77                  ! -- limthd_sal -- ! 
    78    ln_limdS       = .true.          !  activate gravity drainage and flushing (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    7967/ 
    8068!------------------------------------------------------------------------------ 
    81 &namiceitdme   !   Ice mechanical redistribution (ridging and rafting) 
     69&namthd_da      !   Ice lateral melting 
    8270!------------------------------------------------------------------------------ 
    83                   ! -- limitd_me -- ! 
    84    ln_ridging     =   .true.        !  ridging activated (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    85    ln_rafting     =   .true.        !  rafting activated (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    8671/ 
     72!------------------------------------------------------------------------------ 
     73&namthd_do      !   Ice growth in open water 
     74!------------------------------------------------------------------------------ 
     75/ 
     76!------------------------------------------------------------------------------ 
     77&namthd_sal     !   Ice salinity 
     78!------------------------------------------------------------------------------ 
     79/ 
     80!------------------------------------------------------------------------------ 
     81&namthd_pnd     !   Melt ponds 
     82!------------------------------------------------------------------------------ 
     83/ 
     84!------------------------------------------------------------------------------ 
     85&namini         !   Ice initialization 
     86!------------------------------------------------------------------------------ 
     87   ln_iceini        = .true.          !  activate ice initialization (T) or not (F) 
     88   ln_iceini_file   = .true.          !  netcdf file provided for initialization (T) or not (F) 
     89 
     90   sn_hti = 'initice'                 , -12 ,'hti'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     91   sn_hts = 'initice'                 , -12 ,'hts'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     92   sn_ati = 'initice'                 , -12 ,'ati'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     93   sn_tsu = 'initice'                 , -12 ,'tsu'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     94   sn_tmi = 'initice'                 , -12 ,'tmi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     95   sn_smi = 'initice'                 , -12 ,'smi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     96   cn_dir='./' 
     97/ 
     98!------------------------------------------------------------------------------ 
     99&namalb         !   albedo parameters 
     100!------------------------------------------------------------------------------ 
     101/ 
     102!------------------------------------------------------------------------------ 
     103&namdia         !   Diagnostics 
     104!------------------------------------------------------------------------------ 
     105/ 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER/MY_SRC/usrdef_sbc.F90

    r7377 r9019  
    1818   USE sbc_ice         ! Surface boundary condition: ice fields 
    1919   USE phycst          ! physical constants 
    20    USE ice, ONLY       : pfrld, a_i_b 
    21    USE limthd_dh       ! for CALL lim_thd_snwblow 
     20   USE ice, ONLY       : at_i_b, a_i_b 
     21   USE icethd_dh       ! for CALL ice_thd_snwblow 
    2222   ! 
    2323   USE in_out_manager  ! I/O manager 
    … …  
    2525   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
    2626   USE lib_fortran     ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined)  
    27    USE wrk_nemo 
    2827 
    2928   IMPLICIT NONE 
    … …  
    8685      INTEGER, INTENT(in) ::   kt   ! ocean time step 
    8786      !!--------------------------------------------------------------------- 
    88       IF( kt==nit000 .AND. lwp)   WRITE(numout,*)' usrdef_sbc_ice : SAS_BIPER case: NO stress forcing' 
     87      IF( kt==nit000 .AND. lwp)   WRITE(numout,*)' usrdef_sbc_ice : SAS_BIPER case: constant stress forcing' 
    8988      ! 
    90       utau_ice(:,:) = 0._wp 
     89      utau_ice(:,:) = 1.3_wp   ! <=> 0.5 m/s 
    9190      vtau_ice(:,:) = 0._wp 
    9291      ! 
    9392   END SUBROUTINE usrdef_sbc_ice_tau 
    9493 
    95    SUBROUTINE usrdef_sbc_ice_flx( kt ) 
     94   SUBROUTINE usrdef_sbc_ice_flx( kt, phs, phi ) 
    9695      !!--------------------------------------------------------------------- 
    9796      !!                     ***  ROUTINE usrdef_sbc_ice_flx  *** 
    … …  
    9998      !! ** Purpose :   provide the surface boundary (flux) condition over sea-ice 
    10099      !!--------------------------------------------------------------------- 
    101       REAL(wp), DIMENSION(:,:), POINTER ::   zsnw       ! snw distribution after wind blowing 
    102100      INTEGER, INTENT(in) ::   kt   ! ocean time step 
     101      REAL(wp), DIMENSION(:,:,:), INTENT(in)  ::   phs    ! snow thickness 
     102      REAL(wp), DIMENSION(:,:,:), INTENT(in)  ::   phi    ! ice thickness 
     103      !! 
     104      REAL(wp) ::   zfr1, zfr2                 ! local variables 
     105      REAL(wp), DIMENSION(jpi,jpj) ::   zsnw   ! snw distribution after wind blowing 
    103106      !!--------------------------------------------------------------------- 
    104       CALL wrk_alloc( jpi,jpj, zsnw ) 
    105107      ! 
    106108      IF( kt==nit000 .AND. lwp)   WRITE(numout,*)' usrdef_sbc_ice : SAS_BIPER case: NO flux forcing' 
    … …  
    120122      ! ice fields deduced from above 
    121123      zsnw(:,:) = 1._wp 
    122       !!CALL lim_thd_snwblow( pfrld, zsnw )  ! snow distribution over ice after wind blowing  
     124      !!CALL lim_thd_snwblow( at_i_b, zsnw )  ! snow distribution over ice after wind blowing  
    123125      emp_ice  (:,:)   = SUM( a_i_b(:,:,:) * evap_ice(:,:,:), dim=3 ) - sprecip(:,:) * zsnw(:,:) 
    124126      emp_oce  (:,:)   = emp_oce(:,:) - sprecip(:,:) * (1._wp - zsnw(:,:) ) 
    … …  
    130132      ! total fluxes 
    131133      emp_tot (:,:) = emp_ice  + emp_oce 
    132       qns_tot (:,:) = pfrld(:,:) * qns_oce(:,:) + SUM( a_i_b(:,:,:) * qns_ice(:,:,:), dim=3 ) + qemp_ice(:,:) + qemp_oce(:,:) 
    133       qsr_tot (:,:) = pfrld(:,:) * qsr_oce(:,:) + SUM( a_i_b(:,:,:) * qsr_ice(:,:,:), dim=3 ) 
     134      qns_tot (:,:) = at_i_b(:,:) * qns_oce(:,:) + SUM( a_i_b(:,:,:) * qns_ice(:,:,:), dim=3 ) + qemp_ice(:,:) + qemp_oce(:,:) 
     135      qsr_tot (:,:) = at_i_b(:,:) * qsr_oce(:,:) + SUM( a_i_b(:,:,:) * qsr_ice(:,:,:), dim=3 ) 
    134136 
    135       !-------------------------------------------------------------------- 
    136       ! FRACTIONs of net shortwave radiation which is not absorbed in the 
    137       ! thin surface layer and penetrates inside the ice cover 
    138       ! ( Maykut and Untersteiner, 1971 ; Ebert and Curry, 1993 ) 
    139       fr1_i0(:,:) = ( 0.18 * ( 1.0 - cldf_ice ) + 0.35 * cldf_ice ) 
    140       fr2_i0(:,:) = ( 0.82 * ( 1.0 - cldf_ice ) + 0.65 * cldf_ice ) 
    141  
    142       CALL wrk_dealloc( jpi,jpj, zsnw ) 
    143  
     137      ! --- shortwave radiation transmitted below the surface (W/m2, see Grenfell Maykut 77) --- ! 
     138      zfr1 = ( 0.18 * ( 1.0 - cldf_ice ) + 0.35 * cldf_ice )            ! transmission when hi>10cm 
     139      zfr2 = ( 0.82 * ( 1.0 - cldf_ice ) + 0.65 * cldf_ice )            ! zfr2 such that zfr1 + zfr2 to equal 1 
     140      ! 
     141      WHERE    ( phs(:,:,:) <= 0._wp .AND. phi(:,:,:) <  0.1_wp )       ! linear decrease from hi=0 to 10cm   
     142         qsr_ice_tr(:,:,:) = qsr_ice(:,:,:) * ( zfr1 + zfr2 * ( 1._wp - phi(:,:,:) * 10._wp ) ) 
     143      ELSEWHERE( phs(:,:,:) <= 0._wp .AND. phi(:,:,:) >= 0.1_wp )       ! constant (zfr1) when hi>10cm 
     144         qsr_ice_tr(:,:,:) = qsr_ice(:,:,:) * zfr1 
     145      ELSEWHERE                                                         ! zero when hs>0 
     146         qsr_ice_tr(:,:,:) = 0._wp  
     147      END WHERE 
     148           
    144149   END SUBROUTINE usrdef_sbc_ice_flx 
    145150 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER/MY_SRC/usrdef_zgr.F90

    r7355 r9019  
    6262      ! 
    6363      IF(lwp) WRITE(numout,*) 
    64       IF(lwp) WRITE(numout,*) 'usr_def_zgr : LOCK_EXCHANGE configuration (z-coordinate closed box ocean without cavities)' 
     64      IF(lwp) WRITE(numout,*) 'usr_def_zgr : SAS_BIPER configuration (slab ocean - advection of an ice patch in a biperiodic square box domain)' 
    6565      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
    6666      ! 
    … …  
    6868      ! type of vertical coordinate 
    6969      ! --------------------------- 
    70       ! set in usrdef_nam.F90 by reading the namusr_def namelist only ln_zco 
    71       ln_zco    = .TRUE.       ! z-partial-step coordinate 
    72       ln_zps    = .FALSE.      ! z-partial-step coordinate 
    73       ln_sco    = .FALSE.      ! s-coordinate 
     70      ld_zco    = .TRUE.       ! z-full-step coordinate 
     71      ld_zps    = .FALSE.      ! z-partial-step coordinate 
     72      ld_sco    = .FALSE.      ! s-coordinate 
    7473      ld_isfcav = .FALSE.      ! ISF Ice Shelves Flag 
    7574      ! 
Note: See TracChangeset for help on using the changeset viewer.