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 11395 for NEMO/branches/2019/ENHANCE-02_ISF_nemo/cfgs/SHARED – NEMO

Ignore:
Timestamp:
2019-08-02T16:19:00+02:00 (5 years ago)
Author:
mathiot
Message:

ENHANCE-02_ISF_nemo : Initial commit isf simplification (add ISF directory, moved isf routine in and split isf cavity and isf parametrisation, ...) (ticket #2142)

Location:
NEMO/branches/2019/ENHANCE-02_ISF_nemo
Files:
2 edited
1 copied

Legend:

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

    r10824 r11395  
    264264 
    265265          <!-- * variable related to ice shelf forcing * --> 
    266           <field id="fwfisf"       long_name="Ice shelf melting"                             unit="kg/m2/s"  /> 
    267           <field id="fwfisf3d"     long_name="Ice shelf melting"                             unit="kg/m2/s"  grid_ref="grid_T_3D" /> 
    268           <field id="qlatisf"      long_name="Ice shelf latent heat flux"                    unit="W/m2"     /> 
    269           <field id="qlatisf3d"    long_name="Ice shelf latent heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
    270           <field id="qhcisf"       long_name="Ice shelf heat content flux"                   unit="W/m2"     /> 
    271           <field id="qhcisf3d"     long_name="Ice shelf heat content flux"                   unit="W/m2"     grid_ref="grid_T_3D" /> 
    272           <field id="isfgammat"    long_name="transfert coefficient for isf (temperature) "  unit="m/s"      /> 
    273           <field id="isfgammas"    long_name="transfert coefficient for isf (salinity)    "  unit="m/s"      /> 
    274           <field id="stbl"         long_name="salinity in the Losh tbl                    "  unit="PSU"      /> 
    275           <field id="ttbl"         long_name="temperature in the Losh tbl                 "  unit="C"        /> 
    276           <field id="utbl"         long_name="zonal current in the Losh tbl at T point    "  unit="m/s"      /> 
    277           <field id="vtbl"         long_name="merid current in the Losh tbl at T point    "  unit="m/s"      /> 
    278           <field id="thermald"     long_name="thermal driving of ice shelf melting        "  unit="C"        /> 
    279           <field id="tfrz"         long_name="top freezing point (used to compute melt)   "  unit="C"        /> 
    280           <field id="tinsitu"      long_name="top insitu temperature (used to cmpt melt)  "  unit="C"        /> 
    281           <field id="ustar"        long_name="ustar at T point used in ice shelf melting  "  unit="m/s"      /> 
     266          <field id="fwfisf_cav"      long_name="Ice shelf melting"                             unit="kg/m2/s"  /> 
     267          <field id="fwfisf_par"      long_name="Ice shelf melting"                             unit="kg/m2/s"  /> 
     268          <field id="qoceisf_cav"     long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     /> 
     269          <field id="qoceisf_par"     long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     /> 
     270          <field id="qlatisf_cav"     long_name="Ice shelf latent heat flux"                    unit="W/m2"     /> 
     271          <field id="qlatisf_par"     long_name="Ice shelf latent heat flux"                    unit="W/m2"     /> 
     272          <field id="qhcisf_cav"      long_name="Ice shelf heat content flux"                   unit="W/m2"     /> 
     273          <field id="qhcisf_par"      long_name="Ice shelf heat content flux"                   unit="W/m2"     /> 
     274          <field id="fwfisf3d_cav"    long_name="Ice shelf melting"                             unit="kg/m2/s"  grid_ref="grid_T_3D" /> 
     275          <field id="fwfisf3d_par"    long_name="Ice shelf melting"                             unit="kg/m2/s"  grid_ref="grid_T_3D" /> 
     276          <field id="qoceisf3d_cav"   long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     277          <field id="qoceisf3d_par"   long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     278          <field id="qlatisf3d_cav"   long_name="Ice shelf latent heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     279          <field id="qlatisf3d_par"   long_name="Ice shelf latent heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     280          <field id="qhcisf3d_cav"    long_name="Ice shelf heat content flux"                   unit="W/m2"     grid_ref="grid_T_3D" /> 
     281          <field id="qhcisf3d_par"    long_name="Ice shelf heat content flux"                   unit="W/m2"     grid_ref="grid_T_3D" /> 
     282          <field id="ttbl_cav"        long_name="temperature in the tracer sample depth      "  unit="C"        /> 
     283          <field id="ttbl_par"        long_name="temperature in the tracer sample depth      "  unit="C"        /> 
     284          <field id="isfthermald_cav" long_name="thermal driving of ice shelf melting        "  unit="C"        /> 
     285          <field id="isfthermald_par" long_name="thermal driving of ice shelf melting        "  unit="C"        /> 
     286          <field id="isfgammat"       long_name="transfert coefficient for isf (temperature) "  unit="m/s"      /> 
     287          <field id="isfgammas"       long_name="transfert coefficient for isf (salinity)    "  unit="m/s"      /> 
     288          <field id="stbl"            long_name="salinity in the Losh tbl                    "  unit="PSU"      /> 
     289          <field id="utbl"            long_name="zonal current in the Losh tbl at T point    "  unit="m/s"      /> 
     290          <field id="vtbl"            long_name="merid current in the Losh tbl at T point    "  unit="m/s"      /> 
     291          <field id="isfustar"        long_name="ustar at T point used in ice shelf melting  "  unit="m/s"      /> 
     292          <field id="qconisf"         long_name="Conductive heat flux through the ice shelf"    unit="W/m2"     /> 
    282293 
    283294          <!-- *_oce variables available with ln_blk_clio or ln_blk_core --> 
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/cfgs/SHARED/namelist_ref

    r10808 r11395  
    5252      cn_ocerst_outdir = "."         !  directory in which to write output ocean restarts 
    5353   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    =    5840   !  frequency of creation of a restart file (modulo referenced to 1) 
     
    6868!----------------------------------------------------------------------- 
    6969   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    70    rn_isfhmin  =    1.00   !  treshold [m] to discriminate grounding ice from floating ice 
    7170   ! 
    7271   rn_rdt      = 5400.     !  time step for the dynamics and tracer 
     
    7574   ln_crs      = .false.   !  Logical switch for coarsening module      (T => fill namcrs) 
    7675   ! 
    77    ln_meshmask = .false.   !  =T create a mesh file 
     76   ln_meshmask = .true.   !  =T create a mesh file 
    7877/ 
    7978!----------------------------------------------------------------------- 
     
    180179!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T) 
    181180!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T) 
    182 !!   namsbc_isf      ice shelf melting/freezing                         (ln_isfcav  =T : read (ln_read_cfg=T) or set or usr_def_zgr ) 
    183181!!   namsbc_iscpl    coupling option between land ice model and ocean   (ln_isfcav  =T) 
    184182!!   namsbc_wave     external fields from wave model                    (ln_wave    =T) 
     
    436434/ 
    437435!----------------------------------------------------------------------- 
    438 &namsbc_isf    !  Top boundary layer (ISF)                              (ln_isfcav =T : read (ln_read_cfg=T)  
    439 !-----------------------------------------------------------------------             or set or usr_def_zgr ) 
     436&namisf       !  Top boundary layer (ISF)                              (ln_isf = T .AND. ln_isfcav: read (ln_read_cfg=T)  
     437!-----------------------------------------------------------------------                         or set or usr_def_zgr ) 
    440438   !                 ! type of top boundary layer  
    441    nn_isf      = 1         !  ice shelf melting/freezing 
    442                            !  1 = presence of ISF   ;  2 = bg03 parametrisation  
    443                            !  3 = rnf file for ISF  ;  4 = ISF specified freshwater flux 
    444                            !  options 1 and 4 need ln_isfcav = .true. (domzgr) 
    445       !              !  nn_isf = 1 or 2 cases: 
     439   ln_isfcav_mlt = .false.    ! ice shelf melting into the cavity 
     440      cn_isfcav_mlt = '3eq'   ! ice shelf melting formulation (spe/2eq/3eq/oasis) 
     441      !                       ! spe = fwfisf is read from a forcing field 
     442      !                       ! 2eq = ISOMIP  like: 2 equations formulation (Hunter et al., 2006) 
     443      !                       ! 3eq = ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) 
     444      !                       ! oasis = fwfisf is given by oasis 
     445      !              !  cn_isfcav_mlt = 2eq or 3eq cases: 
     446      cn_gammablk = 'ad15'    ! scheme to compute gammat/s (spe,ad15,hj99) 
     447      !                       ! ad15 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
     448      !                       ! hj99 = velocity and stability dependent Gamma    (Holland et al. 1999) 
    446449      rn_gammat0  = 1.e-4     ! gammat coefficient used in blk formula 
    447450      rn_gammas0  = 1.e-4     ! gammas coefficient used in blk formula 
    448       !              !  nn_isf = 1 or 4 cases: 
    449       rn_hisf_tbl =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
     451      ! 
     452      rn_htbl    =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
    450453      !                       ! 0 => thickness of the tbl = thickness of the first wet cell 
    451       !              ! nn_isf = 1 case 
    452       nn_isfblk   = 1         ! 1 ISOMIP  like: 2 equations formulation (Hunter et al., 2006) 
    453       !                       ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) 
    454       nn_gammablk = 1         ! 0 = cst Gammat (= gammat/s) 
    455       !                       ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
    456       !                       ! 2 = velocity and stability dependent Gamma    (Holland et al. 1999) 
    457  
    458    !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
    459    !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    460    !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
    461 !* nn_isf = 4 case 
    462    sn_fwfisf   = 'rnfisf'    ,         -12       ,'sowflisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    463 !* nn_isf = 3 case 
    464    sn_rnfisf   = 'rnfisf'    ,         -12       ,'sofwfisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    465 !* nn_isf = 2 and 3 cases  
    466    sn_depmax_isf ='rnfisf'    ,         -12       ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    467    sn_depmin_isf ='rnfisf'    ,         -12       ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    468 !* nn_isf = 2 case 
    469    sn_Leff_isf = 'rnfisf'    ,         -12       ,'Leff'     ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     454      ! 
     455      !* 'spe' and 'oasis' case 
     456      !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
     457      !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     458      !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
     459      sn_isfcav_fwf = 'isfmlt_cav',      -12       , 'fwflisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     460      ! 
     461   ! 
     462   ln_isfpar_mlt = .false.   ! ice shelf melting parametrised 
     463      cn_isfpar_mlt = 'spe'  ! ice shelf melting parametrisation (spe/bg03/oasis) 
     464      ! 
     465      !* all cases 
     466      !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
     467      !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     468      !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
     469      sn_isfpar_zmax = 'isfmlt_par',       0        ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     470      sn_isfpar_zmin = 'isfmlt_par',       0        ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     471      !* 'spe' and 'oasis' case 
     472      sn_isfpar_fwf = 'isfmlt_par' ,      -12       ,'fwfisf'   ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     473      !* 'bg03' case 
     474      sn_isfpar_Leff = 'isfmlt_par',       0        ,'Leff'     ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     475      ! 
     476   ! 
     477   !ln_iscpl = .false. 
     478   !   nn_drown    = 10        ! number of iteration of the extrapolation loop (fill the new wet cells) 
     479   !   ln_hsb      = .false.   ! activate conservation module (conservation exact after a time of rn_fiscpl) 
     480   !   nn_fiscpl   = 43800     ! (number of time step) conservation period (maybe should be fix to the coupling frequencey of restart frequency) 
    470481/ 
    471482!----------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.