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 12166 for NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/cfgs – NEMO

Ignore:
Timestamp:
2019-12-11T09:53:24+01:00 (4 years ago)
Author:
cetlod
Message:

dev_merge_option2 : merge in ENHANCE-02_ISF_nemo branch

Location:
NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/cfgs
Files:
4 edited
1 copied

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/cfgs/SHARED/field_def_nemo-oce.xml

    r12154 r12166  
    283283 
    284284          <!-- * variable related to ice shelf forcing * --> 
    285           <field id="fwfisf"       long_name="Ice shelf melting"                             unit="kg/m2/s"  /> 
    286           <field id="fwfisf3d"     long_name="Ice shelf melting"                             unit="kg/m2/s"  grid_ref="grid_T_3D" /> 
    287           <field id="qlatisf"      long_name="Ice shelf latent heat flux"                    unit="W/m2"     /> 
    288           <field id="qlatisf3d"    long_name="Ice shelf latent heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
    289           <field id="qhcisf"       long_name="Ice shelf heat content flux"                   unit="W/m2"     /> 
    290           <field id="qhcisf3d"     long_name="Ice shelf heat content flux"                   unit="W/m2"     grid_ref="grid_T_3D" /> 
    291           <field id="isfgammat"    long_name="transfert coefficient for isf (temperature) "  unit="m/s"      /> 
    292           <field id="isfgammas"    long_name="transfert coefficient for isf (salinity)    "  unit="m/s"      /> 
    293           <field id="stbl"         long_name="salinity in the Losh tbl                    "  unit="PSU"      /> 
    294           <field id="ttbl"         long_name="temperature in the Losh tbl                 "  unit="C"        /> 
    295           <field id="utbl"         long_name="zonal current in the Losh tbl at T point    "  unit="m/s"      /> 
    296           <field id="vtbl"         long_name="merid current in the Losh tbl at T point    "  unit="m/s"      /> 
    297           <field id="thermald"     long_name="thermal driving of ice shelf melting        "  unit="C"        /> 
    298           <field id="tfrz"         long_name="top freezing point (used to compute melt)   "  unit="C"        /> 
    299           <field id="tinsitu"      long_name="top insitu temperature (used to cmpt melt)  "  unit="C"        /> 
    300           <field id="ustar"        long_name="ustar at T point used in ice shelf melting  "  unit="m/s"      /> 
     285          <field id="isftfrz_cav"     long_name="freezing point temperature at ocean/isf interface"                unit="degC"     /> 
     286          <field id="isftfrz_par"     long_name="freezing point temperature in the parametrization boundary layer" unit="degC"     /> 
     287          <field id="fwfisf_cav"      long_name="Ice shelf melt rate"                           unit="kg/m2/s"  /> 
     288          <field id="fwfisf_par"      long_name="Ice shelf melt rate"                           unit="kg/m2/s"  /> 
     289          <field id="qoceisf_cav"     long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     /> 
     290          <field id="qoceisf_par"     long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     /> 
     291          <field id="qlatisf_cav"     long_name="Ice shelf latent heat flux"                    unit="W/m2"     /> 
     292          <field id="qlatisf_par"     long_name="Ice shelf latent heat flux"                    unit="W/m2"     /> 
     293          <field id="qhcisf_cav"      long_name="Ice shelf heat content flux of injected water" unit="W/m2"     /> 
     294          <field id="qhcisf_par"      long_name="Ice shelf heat content flux of injected water" unit="W/m2"     /> 
     295          <field id="fwfisf3d_cav"    long_name="Ice shelf melt rate"                           unit="kg/m2/s"  grid_ref="grid_T_3D" /> 
     296          <field id="fwfisf3d_par"    long_name="Ice shelf melt rate"                           unit="kg/m2/s"  grid_ref="grid_T_3D" /> 
     297          <field id="qoceisf3d_cav"   long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     298          <field id="qoceisf3d_par"   long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     299          <field id="qlatisf3d_cav"   long_name="Ice shelf latent heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     300          <field id="qlatisf3d_par"   long_name="Ice shelf latent heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     301          <field id="qhcisf3d_cav"    long_name="Ice shelf heat content flux of injected water" unit="W/m2"     grid_ref="grid_T_3D" /> 
     302          <field id="qhcisf3d_par"    long_name="Ice shelf heat content flux of injected water" unit="W/m2"     grid_ref="grid_T_3D" /> 
     303          <field id="ttbl_cav"        long_name="temperature in Losch tbl"                      unit="degC"     /> 
     304          <field id="ttbl_par"        long_name="temperature in the parametrisation boundary layer" unit="degC" /> 
     305          <field id="isfthermald_cav" long_name="thermal driving of ice shelf melting"          unit="degC"     /> 
     306          <field id="isfthermald_par" long_name="thermal driving of ice shelf melting"          unit="degC"     /> 
     307          <field id="isfgammat"       long_name="Ice shelf heat-transfert velocity"             unit="m/s"      /> 
     308          <field id="isfgammas"       long_name="Ice shelf salt-transfert velocity"             unit="m/s"      /> 
     309          <field id="stbl"            long_name="salinity in the Losh tbl"                      unit="1e-3"     /> 
     310          <field id="utbl"            long_name="zonal current in the Losh tbl at T point"      unit="m/s"      /> 
     311          <field id="vtbl"            long_name="merid current in the Losh tbl at T point"      unit="m/s"      /> 
     312          <field id="isfustar"        long_name="ustar at T point used in ice shelf melting"    unit="m/s"      /> 
     313          <field id="qconisf"         long_name="Conductive heat flux through the ice shelf"    unit="W/m2"     /> 
    301314 
    302315          <!-- *_oce variables available with ln_blk_clio or ln_blk_core --> 
  • NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/cfgs/SHARED/namelist_ice_ref

    r12112 r12166  
    6060      rn_icebfr     =  15.            !        maximum bottom stress per unit volume [N/m3] 
    6161      rn_lfrelax    =   1.e-5         !        relaxation time scale to reach static friction [s-1] 
    62       rn_tensile    =   0.2           !        isotropic tensile strength [0-0.5??] 
     62      rn_tensile    =   0.05          !        isotropic tensile strength [0-0.5??] 
    6363/ 
    6464!------------------------------------------------------------------------------ 
     
    9191!------------------------------------------------------------------------------ 
    9292   ln_rhg_EVP       = .true.          !  EVP rheology 
    93       ln_aEVP       = .false.         !     adaptive rheology (Kimmritz et al. 2016 & 2017) 
     93      ln_aEVP       = .true.          !     adaptive rheology (Kimmritz et al. 2016 & 2017) 
    9494      rn_creepl     =   2.0e-9        !     creep limit [1/s] 
    9595      rn_ecc        =   2.0           !     eccentricity of the elliptical yield curve           
  • NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/cfgs/SHARED/namelist_ref

    r12154 r12166  
    55!! namelists    2 - Surface boundary (namsbc, namsbc_flx, namsbc_blk, namsbc_cpl, 
    66!!                                    namsbc_sas, namtra_qsr, namsbc_rnf, 
    7 !!                                    namsbc_isf, namsbc_iscpl, namsbc_apr, 
     7!!                                    namsbc_isf, namsbc_iscpl, namsbc_apr,  
     8!!                                    namisf, namsbc_apr,  
    89!!                                    namsbc_ssr, namsbc_wave, namberg) 
    910!!              3 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) 
     
    5152      cn_ocerst_out   = "restart"   !  suffix of ocean restart name (output) 
    5253      cn_ocerst_outdir = "."        !  directory in which to write output ocean restarts 
    53    ln_iscpl    = .false.   !  cavity evolution forcing or coupling to ice sheet model 
    54    nn_istate   =       0   !  output the initial state (1) or not (0) 
     54   nn_istate   =       1   !  output the initial state (1) or not (0) 
    5555   ln_rst_list = .false.   !  output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F) 
    5656   nn_stock    =       0   !  used only if ln_rst_list = F: output restart freqeuncy (modulo referenced to 1) 
     
    7272!----------------------------------------------------------------------- 
    7373   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    74    rn_isfhmin  =    1.00   !  treshold [m] to discriminate grounding ice from floating ice 
    7574   ! 
    7675   rn_rdt      = 5400.     !  time step for the dynamics and tracer 
     
    7978   ln_crs      = .false.   !  Logical switch for coarsening module      (T => fill namcrs) 
    8079   ! 
    81    ln_meshmask = .false.   !  =T create a mesh file 
     80   ln_meshmask = .true.   !  =T create a mesh file 
    8281/ 
    8382!----------------------------------------------------------------------- 
     
    184183!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T) 
    185184!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T) 
    186 !!   namsbc_isf      ice shelf melting/freezing                         (ln_isfcav  =T : read (ln_read_cfg=T) or set or usr_def_zgr ) 
    187 !!   namsbc_iscpl    coupling option between land ice model and ocean   (ln_isfcav  =T) 
    188185!!   namsbc_wave     external fields from wave model                    (ln_wave    =T) 
    189186!!   namberg         iceberg floats                                     (ln_icebergs=T) 
     
    223220   ln_rnf      = .false.   !  runoffs                                   (T => fill namsbc_rnf) 
    224221   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr ) 
    225    ln_isf      = .false.   !  ice shelf                                 (T   => fill namsbc_isf & namsbc_iscpl) 
    226222   ln_wave     = .false.   !  Activate coupling with wave  (T => fill namsbc_wave) 
    227223   ln_cdgw     = .false.   !  Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave) 
     
    482478/ 
    483479!----------------------------------------------------------------------- 
    484 &namsbc_isf    !  Top boundary layer (ISF)                              (ln_isfcav =T : read (ln_read_cfg=T) 
    485 !-----------------------------------------------------------------------             or set or usr_def_zgr ) 
    486    !                 ! type of top boundary layer 
    487    nn_isf      = 1         !  ice shelf melting/freezing 
    488                            !  1 = presence of ISF   ;  2 = bg03 parametrisation 
    489                            !  3 = rnf file for ISF  ;  4 = ISF specified freshwater flux 
    490                            !  options 1 and 4 need ln_isfcav = .true. (domzgr) 
    491       !              !  nn_isf = 1 or 2 cases: 
    492       rn_gammat0  = 1.e-4     ! gammat coefficient used in blk formula 
    493       rn_gammas0  = 1.e-4     ! gammas coefficient used in blk formula 
    494       !              !  nn_isf = 1 or 4 cases: 
    495       rn_hisf_tbl =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
    496       !                       ! 0 => thickness of the tbl = thickness of the first wet cell 
    497       !              ! nn_isf = 1 case 
    498       nn_isfblk   = 1         ! 1 ISOMIP  like: 2 equations formulation (Hunter et al., 2006) 
    499       !                       ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) 
    500       nn_gammablk = 1         ! 0 = cst Gammat (= gammat/s) 
    501       !                       ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
    502       !                       ! 2 = velocity and stability dependent Gamma    (Holland et al. 1999) 
    503  
    504    !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
    505    !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    506    !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
    507 !* nn_isf = 4 case 
    508    sn_fwfisf   = 'rnfisf'    ,         -12.      ,'sowflisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    509 !* nn_isf = 3 case 
    510    sn_rnfisf   = 'rnfisf'    ,         -12.      ,'sofwfisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    511 !* nn_isf = 2 and 3 cases 
    512    sn_depmax_isf ='rnfisf'   ,         -12.      ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    513    sn_depmin_isf ='rnfisf'   ,         -12.      ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    514 !* nn_isf = 2 case 
    515    sn_Leff_isf = 'rnfisf'    ,         -12.      ,'Leff'     ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    516 / 
    517 !----------------------------------------------------------------------- 
    518 &namsbc_iscpl  !   land ice / ocean coupling option                     (ln_isfcav =T : read (ln_read_cfg=T) 
    519 !-----------------------------------------------------------------------             or set or usr_def_zgr ) 
    520    nn_drown    = 10        ! number of iteration of the extrapolation loop (fill the new wet cells) 
    521    ln_hsb      = .false.   ! activate conservation module (conservation exact after a time of rn_fiscpl) 
    522    nn_fiscpl   = 43800     ! (number of time step) conservation period (maybe should be fix to the coupling frequencey of restart frequency) 
     480&namisf       !  Top boundary layer (ISF)                               (default: OFF) 
     481!----------------------------------------------------------------------- 
     482   ! 
     483   ! ---------------- ice shelf load ------------------------------- 
     484   ! 
     485   cn_isfload = 'uniform'      ! scheme to compute ice shelf load (ln_isfcav = .true. in domain_cfg.nc) 
     486      rn_isfload_T = -1.9 
     487      rn_isfload_S = 34.4 
     488   ! 
     489   ! ---------------- ice shelf melt formulation ------------------------------- 
     490   ! 
     491   ln_isf = .false.           ! activate ice shelf module 
     492      ln_isfdebug = .false.      ! add debug print in ISF code (global min/max/sum of specific variable) 
     493      cn_isfdir   = './'         ! directory for all ice shelf input file 
     494      ! 
     495      ! ---------------- cavities opened ------------------------------- 
     496      ! 
     497      ln_isfcav_mlt = .false.    ! ice shelf melting into the cavity (need ln_isfcav = .true. in domain_cfg.nc) 
     498         cn_isfcav_mlt = '3eq'   ! ice shelf melting formulation (spe/2eq/3eq/oasis) 
     499         !                       ! spe = fwfisf is read from a forcing field 
     500         !                       ! 2eq = ISOMIP  like: 2 equations formulation (Hunter et al., 2006 for a short description) 
     501         !                       ! 3eq = ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2016 for a short description) 
     502         !                       ! oasis = fwfisf is given by oasis and pattern by file sn_isfcav_fwf 
     503         !              !  cn_isfcav_mlt = 2eq or 3eq cases: 
     504         cn_gammablk = 'vel'     ! scheme to compute gammat/s (spe,ad15,hj99) 
     505         !                       ! spe      = constant transfert velocity (rn_gammat0, rn_gammas0) 
     506         !                       ! vel      = velocity dependent transfert velocity (u* * gammat/s) (Asay-Davis et al. 2016 for a short description) 
     507         !                       ! vel_stab = velocity and stability dependent transfert coeficient (Holland et al. 1999 for a complete description) 
     508         rn_gammat0  = 1.4e-2    ! gammat coefficient used in spe, vel and vel_stab gamma computation method 
     509         rn_gammas0  = 4.0e-4    ! gammas coefficient used in spe, vel and vel_stab gamma computation method 
     510         ! 
     511         rn_htbl     =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
     512         !                       ! 0 => thickness of the tbl = thickness of the first wet cell 
     513         ! 
     514         !* 'spe' and 'oasis' case 
     515         !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
     516         !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     517         !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
     518         sn_isfcav_fwf = 'isfmlt_cav',      -12.      , 'fwflisf'  ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     519      ! 
     520      ! ---------------- cavities parametrised ------------------------------- 
     521      ! 
     522      ln_isfpar_mlt = .false.   ! ice shelf melting parametrised 
     523         cn_isfpar_mlt = 'spe'  ! ice shelf melting parametrisation (spe/bg03/oasis) 
     524         !                      ! spe   = fwfisf is read from a forcing field 
     525         !                      ! bg03  = melt computed using Beckmann and Goosse parametrisation 
     526         !                      ! oasis = fwfisf is given by oasis and pattern by file sn_isfpar_fwf 
     527         ! 
     528         !* all cases 
     529         !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
     530         !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     531         !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
     532         sn_isfpar_zmax = 'isfmlt_par',       0        ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     533         sn_isfpar_zmin = 'isfmlt_par',       0        ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     534         !* 'spe' and 'oasis' case 
     535         sn_isfpar_fwf = 'isfmlt_par' ,      -12.      ,'sofwfisf' ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
     536         !* 'bg03' case 
     537         sn_isfpar_Leff = 'isfmlt_par',       0.       ,'Leff'     ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
     538      ! 
     539      ! ---------------- ice sheet coupling ------------------------------- 
     540      ! 
     541      ln_isfcpl = .false. 
     542         nn_drown       = 10        ! number of iteration of the extrapolation loop (fill the new wet cells) 
     543         ln_isfcpl_cons = .false. 
    523544/ 
    524545!----------------------------------------------------------------------- 
     
    12381259&nam_diatmb    !  Top Middle Bottom Output                              (default: OFF) 
    12391260!----------------------------------------------------------------------- 
    1240    ln_diatmb   = .false.   !  Choose Top Middle and Bottom output or not 
     1261   ln_diatmb   = .true.   !  Choose Top Middle and Bottom output or not 
    12411262/ 
    12421263!----------------------------------------------------------------------- 
  • NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/cfgs/ref_cfgs.txt

    r12154 r12166  
    1212ORCA2_ICE OCE ICE 
    1313SPITZ12 OCE ICE 
    14 eORCA025_ICE OCE ICE 
    15 eORCA025_ICE_ABL OCE ICE ABL 
    16 eORCA025_SAS_ICE_ABL OCE SAS ICE ABL 
     14WED025 OCE ICE 
Note: See TracChangeset for help on using the changeset viewer.