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 12721 for NEMO/branches/UKMO/NEMO_4.0.2_GO8_package_ENHANCE-02_ISF_nemo/cfgs – NEMO

Ignore:
Timestamp:
2020-04-08T20:42:57+02:00 (4 years ago)
Author:
mathiot
Message:

NEMO_4.0.2_GO8_package_ENHANCE-02_ISF_nemo: add last year isf dev

Location:
NEMO/branches/UKMO/NEMO_4.0.2_GO8_package_ENHANCE-02_ISF_nemo/cfgs
Files:
16 added
3 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0.2_GO8_package_ENHANCE-02_ISF_nemo/cfgs/SHARED/field_def_nemo-oce.xml

    r12660 r12721  
    314314 
    315315          <!-- * variable related to ice shelf forcing * --> 
    316           <field id="fwfisf"       long_name="Ice shelf melting"                             unit="kg/m2/s"  /> 
    317           <field id="fwfisf3d"     long_name="Ice shelf melting"                             unit="kg/m2/s"  grid_ref="grid_T_3D" /> 
    318           <field id="qlatisf"      long_name="Ice shelf latent heat flux"                    unit="W/m2"     /> 
    319           <field id="qlatisf3d"    long_name="Ice shelf latent heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
    320           <field id="qhcisf"       long_name="Ice shelf heat content flux"                   unit="W/m2"     /> 
    321           <field id="qhcisf3d"     long_name="Ice shelf heat content flux"                   unit="W/m2"     grid_ref="grid_T_3D" /> 
    322           <field id="isfgammat"    long_name="transfert coefficient for isf (temperature) "  unit="m/s"      /> 
    323           <field id="isfgammas"    long_name="transfert coefficient for isf (salinity)    "  unit="m/s"      /> 
    324           <field id="stbl"         long_name="salinity in the Losh tbl                    "  unit="PSU"      /> 
    325           <field id="ttbl"         long_name="temperature in the Losh tbl                 "  unit="C"        /> 
    326           <field id="utbl"         long_name="zonal current in the Losh tbl at T point    "  unit="m/s"      /> 
    327           <field id="vtbl"         long_name="merid current in the Losh tbl at T point    "  unit="m/s"      /> 
    328           <field id="thermald"     long_name="thermal driving of ice shelf melting        "  unit="C"        /> 
    329           <field id="tfrz"         long_name="top freezing point (used to compute melt)   "  unit="C"        /> 
    330           <field id="tinsitu"      long_name="top insitu temperature (used to cmpt melt)  "  unit="C"        /> 
    331           <field id="ustar"        long_name="ustar at T point used in ice shelf melting  "  unit="m/s"      /> 
     316          <field id="isftfrz_cav"     long_name="freezing point temperature at ocean/isf interface"                unit="degC"     /> 
     317          <field id="isftfrz_par"     long_name="freezing point temperature in the parametrization boundary layer" unit="degC"     /> 
     318          <field id="fwfisf_cav"      long_name="Ice shelf melt rate"                           unit="kg/m2/s"  /> 
     319          <field id="fwfisf_par"      long_name="Ice shelf melt rate"                           unit="kg/m2/s"  /> 
     320          <field id="qoceisf_cav"     long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     /> 
     321          <field id="qoceisf_par"     long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     /> 
     322          <field id="qlatisf_cav"     long_name="Ice shelf latent heat flux"                    unit="W/m2"     /> 
     323          <field id="qlatisf_par"     long_name="Ice shelf latent heat flux"                    unit="W/m2"     /> 
     324          <field id="qhcisf_cav"      long_name="Ice shelf heat content flux of injected water" unit="W/m2"     /> 
     325          <field id="qhcisf_par"      long_name="Ice shelf heat content flux of injected water" unit="W/m2"     /> 
     326          <field id="fwfisf3d_cav"    long_name="Ice shelf melt rate"                           unit="kg/m2/s"  grid_ref="grid_T_3D" /> 
     327          <field id="fwfisf3d_par"    long_name="Ice shelf melt rate"                           unit="kg/m2/s"  grid_ref="grid_T_3D" /> 
     328          <field id="qoceisf3d_cav"   long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     329          <field id="qoceisf3d_par"   long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     330          <field id="qlatisf3d_cav"   long_name="Ice shelf latent heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     331          <field id="qlatisf3d_par"   long_name="Ice shelf latent heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     332          <field id="qhcisf3d_cav"    long_name="Ice shelf heat content flux of injected water" unit="W/m2"     grid_ref="grid_T_3D" /> 
     333          <field id="qhcisf3d_par"    long_name="Ice shelf heat content flux of injected water" unit="W/m2"     grid_ref="grid_T_3D" /> 
     334          <field id="ttbl_cav"        long_name="temperature in Losch tbl"                      unit="degC"     /> 
     335          <field id="ttbl_par"        long_name="temperature in the parametrisation boundary layer" unit="degC" /> 
     336          <field id="isfthermald_cav" long_name="thermal driving of ice shelf melting"          unit="degC"     /> 
     337          <field id="isfthermald_par" long_name="thermal driving of ice shelf melting"          unit="degC"     /> 
     338          <field id="isfgammat"       long_name="Ice shelf heat-transfert velocity"             unit="m/s"      /> 
     339          <field id="isfgammas"       long_name="Ice shelf salt-transfert velocity"             unit="m/s"      /> 
     340          <field id="stbl"            long_name="salinity in the Losh tbl"                      unit="1e-3"     /> 
     341          <field id="utbl"            long_name="zonal current in the Losh tbl at T point"      unit="m/s"      /> 
     342          <field id="vtbl"            long_name="merid current in the Losh tbl at T point"      unit="m/s"      /> 
     343          <field id="isfustar"        long_name="ustar at T point used in ice shelf melting"    unit="m/s"      /> 
     344          <field id="qconisf"         long_name="Conductive heat flux through the ice shelf"    unit="W/m2"     /> 
    332345 
    333346          <!-- *_oce variables available with ln_blk_clio or ln_blk_core --> 
  • NEMO/branches/UKMO/NEMO_4.0.2_GO8_package_ENHANCE-02_ISF_nemo/cfgs/SHARED/namelist_ref

    r12660 r12721  
    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!!                                    namisf, namsbc_apr,  
    88!!                                    namsbc_ssr, namsbc_wave, namberg) 
    99!!              3 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) 
     
    5151      cn_ocerst_out   = "restart"   !  suffix of ocean restart name (output) 
    5252      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) 
     53   nn_istate   =       1   !  output the initial state (1) or not (0) 
    5554   ln_rst_list = .false.   !  output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F) 
    5655   nn_stock    =       0   !  used only if ln_rst_list = F: output restart freqeuncy (modulo referenced to 1) 
     
    7372!----------------------------------------------------------------------- 
    7473   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    75    rn_isfhmin  =    1.00   !  treshold [m] to discriminate grounding ice from floating ice 
    7674   ! 
    7775   rn_rdt      = 5400.     !  time step for the dynamics and tracer 
     
    185183!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T) 
    186184!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T) 
    187 !!   namsbc_isf      ice shelf melting/freezing                         (ln_isfcav  =T : read (ln_read_cfg=T) or set or usr_def_zgr ) 
    188 !!   namsbc_iscpl    coupling option between land ice model and ocean   (ln_isfcav  =T) 
    189185!!   namsbc_wave     external fields from wave model                    (ln_wave    =T) 
    190186!!   namberg         iceberg floats                                     (ln_icebergs=T) 
     
    223219   ln_rnf      = .false.   !  runoffs                                   (T => fill namsbc_rnf) 
    224220   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) 
    226221   ln_wave     = .false.   !  Activate coupling with wave  (T => fill namsbc_wave) 
    227222   ln_cdgw     = .false.   !  Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave) 
     
    443438/ 
    444439!----------------------------------------------------------------------- 
    445 &namsbc_isf    !  Top boundary layer (ISF)                              (ln_isfcav =T : read (ln_read_cfg=T)  
    446 !-----------------------------------------------------------------------             or set or usr_def_zgr ) 
    447    !                 ! type of top boundary layer  
    448    nn_isf      = 1         !  ice shelf melting/freezing 
    449                            !  1 = presence of ISF   ;  2 = bg03 parametrisation  
    450                            !  3 = rnf file for ISF  ;  4 = ISF specified freshwater flux 
    451                            !  options 1 and 4 need ln_isfcav = .true. (domzgr) 
    452       !              !  nn_isf = 1 or 2 cases: 
    453       rn_gammat0  = 1.e-4     ! gammat coefficient used in blk formula 
    454       rn_gammas0  = 1.e-4     ! gammas coefficient used in blk formula 
    455       !              !  nn_isf = 1 or 4 cases: 
    456       rn_hisf_tbl =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
    457       !                       ! 0 => thickness of the tbl = thickness of the first wet cell 
    458       !              ! nn_isf = 1 case 
    459       nn_isfblk   = 1         ! 1 ISOMIP  like: 2 equations formulation (Hunter et al., 2006) 
    460       !                       ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) 
    461       nn_gammablk = 1         ! 0 = cst Gammat (= gammat/s) 
    462       !                       ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
    463       !                       ! 2 = velocity and stability dependent Gamma    (Holland et al. 1999) 
    464  
    465    !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
    466    !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    467    !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
    468 !* nn_isf = 4 case 
    469    sn_fwfisf   = 'rnfisf'    ,         -12.      ,'sowflisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    470 !* nn_isf = 3 case 
    471    sn_rnfisf   = 'rnfisf'    ,         -12.      ,'sofwfisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    472 !* nn_isf = 2 and 3 cases  
    473    sn_depmax_isf ='rnfisf'   ,         -12.      ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    474    sn_depmin_isf ='rnfisf'   ,         -12.      ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    475 !* nn_isf = 2 case 
    476    sn_Leff_isf = 'rnfisf'    ,         -12.      ,'Leff'     ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    477 / 
    478 !----------------------------------------------------------------------- 
    479 &namsbc_iscpl  !   land ice / ocean coupling option                     (ln_isfcav =T : read (ln_read_cfg=T)  
    480 !-----------------------------------------------------------------------             or set or usr_def_zgr ) 
    481    nn_drown    = 10        ! number of iteration of the extrapolation loop (fill the new wet cells) 
    482    ln_hsb      = .false.   ! activate conservation module (conservation exact after a time of rn_fiscpl) 
    483    nn_fiscpl   = 43800     ! (number of time step) conservation period (maybe should be fix to the coupling frequencey of restart frequency) 
     440&namisf       !  Top boundary layer (ISF)                               (default: OFF) 
     441!----------------------------------------------------------------------- 
     442   ! 
     443   ! ---------------- ice shelf load ------------------------------- 
     444   ! 
     445   cn_isfload = 'uniform'      ! scheme to compute ice shelf load (ln_isfcav = .true. in domain_cfg.nc) 
     446      rn_isfload_T = -1.9 
     447      rn_isfload_S = 34.4 
     448   ! 
     449   ! ---------------- ice shelf melt formulation ------------------------------- 
     450   ! 
     451   ln_isf = .false.           ! activate ice shelf module 
     452      ln_isfdebug = .false.      ! add debug print in ISF code (global min/max/sum of specific variable) 
     453      cn_isfdir   = './'         ! directory for all ice shelf input file 
     454      ! 
     455      ! ---------------- cavities opened ------------------------------- 
     456      ! 
     457      ln_isfcav_mlt = .false.    ! ice shelf melting into the cavity (need ln_isfcav = .true. in domain_cfg.nc) 
     458         cn_isfcav_mlt = '3eq'   ! ice shelf melting formulation (spe/2eq/3eq/oasis) 
     459         !                       ! spe = fwfisf is read from a forcing field 
     460         !                       ! 2eq = ISOMIP  like: 2 equations formulation (Hunter et al., 2006 for a short description) 
     461         !                       ! 3eq = ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2016 for a short description) 
     462         !                       ! oasis = fwfisf is given by oasis and pattern by file sn_isfcav_fwf 
     463         !              !  cn_isfcav_mlt = 2eq or 3eq cases: 
     464         cn_gammablk = 'vel'     ! scheme to compute gammat/s (spe,ad15,hj99) 
     465         !                       ! spe      = constant transfert velocity (rn_gammat0, rn_gammas0) 
     466         !                       ! vel      = velocity dependent transfert velocity (u* * gammat/s) (Asay-Davis et al. 2016 for a short description) 
     467         !                       ! vel_stab = velocity and stability dependent transfert coeficient (Holland et al. 1999 for a complete description) 
     468         rn_gammat0  = 1.4e-2    ! gammat coefficient used in spe, vel and vel_stab gamma computation method 
     469         rn_gammas0  = 4.0e-4    ! gammas coefficient used in spe, vel and vel_stab gamma computation method 
     470         ! 
     471         rn_htbl     =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
     472         !                       ! 0 => thickness of the tbl = thickness of the first wet cell 
     473         ! 
     474         !* 'spe' and 'oasis' case 
     475         !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
     476         !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     477         !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
     478         sn_isfcav_fwf = 'isfmlt_cav',      -12.      , 'fwflisf'  ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     479      ! 
     480      ! ---------------- cavities parametrised ------------------------------- 
     481      ! 
     482      ln_isfpar_mlt = .false.   ! ice shelf melting parametrised 
     483         cn_isfpar_mlt = 'spe'  ! ice shelf melting parametrisation (spe/bg03/oasis) 
     484         !                      ! spe   = fwfisf is read from a forcing field 
     485         !                      ! bg03  = melt computed using Beckmann and Goosse parametrisation 
     486         !                      ! oasis = fwfisf is given by oasis and pattern by file sn_isfpar_fwf 
     487         ! 
     488         !* all cases 
     489         !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
     490         !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     491         !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
     492         sn_isfpar_zmax = 'isfmlt_par',       0        ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     493         sn_isfpar_zmin = 'isfmlt_par',       0        ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     494         !* 'spe' and 'oasis' case 
     495         sn_isfpar_fwf = 'isfmlt_par' ,      -12.      ,'sofwfisf' ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
     496         !* 'bg03' case 
     497         sn_isfpar_Leff = 'isfmlt_par',       0.       ,'Leff'     ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
     498      ! 
     499      ! ---------------- ice sheet coupling ------------------------------- 
     500      ! 
     501      ln_isfcpl = .false. 
     502         nn_drown       = 10        ! number of iteration of the extrapolation loop (fill the new wet cells) 
     503         ln_isfcpl_cons = .false. 
    484504/ 
    485505!----------------------------------------------------------------------- 
  • NEMO/branches/UKMO/NEMO_4.0.2_GO8_package_ENHANCE-02_ISF_nemo/cfgs/ref_cfgs.txt

    r12658 r12721  
    99ORCA2_ICE_PISCES OCE TOP ICE NST 
    1010SPITZ12 OCE ICE 
     11WED025 OCE ICE 
Note: See TracChangeset for help on using the changeset viewer.