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 11586 for NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/cfgs – NEMO

Ignore:
Timestamp:
2019-09-20T17:28:02+02:00 (5 years ago)
Author:
gsamson
Message:

dev_r11265_ABL : see #2131

  • merge HPC-13_IRRMANN_BDY_optimization branch @ r11535 (last commit) with dev_r11265_ABL branch @ r11414 (except doc directory)
  • change ORCA2 results due to changes in HPC-13_IRRMANN_BDY_optimization branch
Location:
NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/cfgs
Files:
8 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/cfgs/ORCA2_ICE_PISCES/EXPREF/context_nemo.xml

    r11413 r11586  
    4242      <axis id="iax_20C" long_name="20 degC isotherm"  unit="degC"              /> 
    4343      <axis id="iax_28C" long_name="28 degC isotherm"  unit="degC"              /> 
    44       <!-- ABL vertical axis definition --> 
    45       <axis id="ght_abl" long_name="ABL Vertical T levels" unit="m" positive="up"   /> 
    46       <axis id="ghw_abl" long_name="ABL Vertical W levels" unit="m" positive="up"   /> 
    4744    </axis_definition> 
    4845  
  • NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/cfgs/SHARED/field_def_nemo-ice.xml

    r11413 r11586  
    4646      
    4747     <!-- melt ponds --> 
    48      <field id="iceapnd"      long_name="melt pond fraction"                                      standard_name="sea_ice_meltpond_fraction"                 unit="%" />  
     48     <field id="iceapnd"      long_name="melt pond concentration"                                 standard_name="sea_ice_meltpond_concentration"            unit=""  />  
     49          <field id="icehpnd"      long_name="melt pond depth"                                         standard_name="sea_ice_meltpond_depth"                    unit="m" />  
    4950          <field id="icevpnd"      long_name="melt pond volume"                                        standard_name="sea_ice_meltpond_volume"                   unit="m" />  
    5051      
     
    328329     <field field_ref="icesalt"          name="sisali" /> 
    329330     <field field_ref="iceapnd"          name="siapnd" /> 
     331     <field field_ref="icehpnd"          name="sihpnd" /> 
    330332     <field field_ref="icevpnd"          name="sivpnd" /> 
    331333          <field field_ref="iceage"           name="siage"  /> 
  • NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/cfgs/SHARED/field_def_nemo-oce.xml

    r11413 r11586  
    237237       
    238238      <!-- SBC --> 
     239 
    239240      <field_group id="SBC" > <!-- time step automaticaly defined based on nn_fsbc --> 
    240241 
     
    300301          <field id="iceshelf_cea" long_name="Iceshelf"                                  standard_name="water_flux_into_sea_water_from_iceshelf"   unit="kg/m2/s"  /> 
    301302 
    302  
    303303          <!-- available if key_oasis3 + conservative method --> 
    304304          <field id="rain"          long_name="Liquid precipitation"                                     standard_name="rainfall_flux"                                                                 unit="kg/m2/s"  /> 
     
    348348      <!-- ABL --> 
    349349      <field_group id="ABL" > <!-- time step automaticaly defined based on nn_fsbc --> 
     350 
    350351   <!-- variables available with ABL on atmospheric T grid--> 
    351352   <field_group id="grid_ABL3D" grid_ref="grid_TA_3D" > 
     
    354355          <field id="t_abl"      long_name="ABL potential temperature"     standard_name="abl_theta"      unit="K"        /> 
    355356          <field id="q_abl"      long_name="ABL specific humidity"         standard_name="abl_qspe"       unit="kg/kg"    /> 
     357          <!-- debug (to be removed) --> 
     358          <field id="u_dta"      long_name="DTA i-horizontal velocity"     standard_name="dta_x_velocity" unit="m/s"      /> 
     359          <field id="v_dta"      long_name="DTA j-horizontal velocity"     standard_name="dta_y_velocity" unit="m/s"      /> 
     360          <field id="t_dta"      long_name="DTA potential temperature"     standard_name="dta_theta"      unit="K"        /> 
     361          <field id="q_dta"      long_name="DTA specific humidity"         standard_name="dta_qspe"       unit="kg/kg"    /> 
    356362          <field id="coeft"      long_name="ABL nudging coefficient"       standard_name="coeft"          unit=""         /> 
    357363   </field_group> 
     364 
    358365   <field_group id="grid_ABL2D" grid_ref="grid_TA_2D" > 
    359366          <field id="pblh"       long_name="ABL height"                    standard_name="abl_height"     unit="m"        /> 
     367          <field id="uz1_abl"    long_name="ABL i-horizontal velocity"     standard_name="abl_x_velocity" unit="m/s"      /> 
     368          <field id="vz1_abl"    long_name="ABL j-horizontal velocity"     standard_name="abl_y_velocity" unit="m/s"      /> 
     369          <field id="uvz1_abl"   long_name="ABL wind speed module"         standard_name="abl_wind_speed" unit="m/s"       > sqrt( uz1_abl^2 + vz1_abl^2 ) </field> 
     370          <field id="tz1_abl"    long_name="ABL potential temperature"     standard_name="abl_theta"      unit="K"        /> 
     371          <field id="qz1_abl"    long_name="ABL specific humidity"         standard_name="abl_qspe"       unit="kg/kg"    /> 
     372          <field id="uz1_dta"    long_name="DTA i-horizontal velocity"     standard_name="dta_x_velocity" unit="m/s"      /> 
     373          <field id="vz1_dta"    long_name="DTA j-horizontal velocity"     standard_name="dta_y_velocity" unit="m/s"      /> 
     374          <field id="uvz1_dta"   long_name="DTA wind speed module"         standard_name="dta_wind_speed" unit="m/s"       > sqrt( uz1_dta^2 + vz1_dta^2 ) </field>  
     375          <field id="tz1_dta"    long_name="DTA potential temperature"     standard_name="dta_theta"      unit="K"        /> 
     376          <field id="qz1_dta"    long_name="DTA specific humidity"         standard_name="dta_qspe"       unit="kg/kg"    /> 
     377          <!-- debug (to be removed) --> 
     378          <field id="uz1_geo"    long_name="GEO i-horizontal velocity"     standard_name="geo_x_velocity" unit="m/s"      /> 
     379          <field id="vz1_geo"    long_name="GEO j-horizontal velocity"     standard_name="geo_y_velocity" unit="m/s"      /> 
     380          <field id="uvz1_geo"   long_name="GEO wind speed module"         standard_name="geo_wind_speed" unit="m/s"       > sqrt( uz1_geo^2 + vz1_geo^2 ) </field>  
    360381   </field_group> 
     382 
    361383      </field_group> <!-- ABL --> 
     384 
    362385       
    363386      <!-- U grid --> 
     
    375398        <field id="uocet"        long_name="ocean transport along i-axis times temperature (CRS)"                                               unit="degC*m/s"   grid_ref="grid_U_3D" /> 
    376399        <field id="uoces"        long_name="ocean transport along i-axis times salinity (CRS)"                                                  unit="1e-3*m/s"   grid_ref="grid_U_3D" /> 
     400        <field id="ssuww"        long_name="ocean surface wind work along i-axis"                   standard_name="surface_x_wind_work"         unit="N/m*s"                            > utau * ssu </field> 
    377401 
    378402        <!-- u-eddy diffusivity coefficients (available if ln_traldf_OFF=F) --> 
     
    429453        <field id="vocet"        long_name="ocean transport along j-axis times temperature (CRS)"                                               unit="degC*m/s"   grid_ref="grid_V_3D" /> 
    430454        <field id="voces"        long_name="ocean transport along j-axis times salinity (CRS)"                                                  unit="1e-3*m/s"   grid_ref="grid_V_3D" /> 
     455        <field id="ssvww"        long_name="ocean surface wind work along j-axis"                   standard_name="surface_y_wind_work"         unit="N/m*s"                            > vtau * ssv </field> 
     456 
    431457 
    432458        <!-- v-eddy diffusivity coefficients (available if ln_traldf_OFF=F) --> 
  • NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/cfgs/SHARED/namelist_ice_ref

    r11413 r11586  
    7070   ln_str_H79       = .true.          !  ice strength param.: Hibler_79   => P = pstar*<h>*exp(-c_rhg*A)                       
    7171      rn_pstar      =   2.0e+04       !     ice strength thickness parameter [N/m2] 
    72       rn_crhg       =   20.0          !     ice strength conc. parameter (-) 
     72      rn_crhg       =  20.0           !     ice strength conc. parameter (-) 
    7373                   ! -- ice_rdgrft -- ! 
    7474   rn_csrdg         =   0.5           !  fraction of shearing energy contributing to ridging 
     
    177177&namthd_pnd     !   Melt ponds 
    178178!------------------------------------------------------------------------------ 
    179    ln_pnd_H12       = .false.         !  activate evolutive melt ponds (from Holland et al 2012) 
    180    ln_pnd_CST       = .false.         !  activate constant  melt ponds 
    181       rn_apnd       =   0.2           !     prescribed pond fraction, at Tsu=0 degC 
    182       rn_hpnd       =   0.05          !     prescribed pond depth,    at Tsu=0 degC 
    183    ln_pnd_alb       = .false.         !  melt ponds affect albedo or not 
     179   ln_pnd           = .false.         !  activate melt ponds or not 
     180     ln_pnd_H12     = .false.         !  activate evolutive melt ponds (from Holland et al 2012) 
     181     ln_pnd_CST     = .false.         !  activate constant  melt ponds 
     182       rn_apnd      =   0.2           !     prescribed pond fraction, at Tsu=0 degC 
     183       rn_hpnd      =   0.05          !     prescribed pond depth,    at Tsu=0 degC 
     184     ln_pnd_alb     = .false.         !  melt ponds affect albedo or not 
    184185/ 
    185186!------------------------------------------------------------------------------ 
     
    207208   rn_hpd_ini_n     =   0.05          !  initial pond depth          (m), North 
    208209   rn_hpd_ini_s     =   0.05          !        "            "             South 
    209    !        ! if ln_iceini_file=T 
     210   ! -- for ln_iceini_file = T 
    210211   sn_hti = 'Ice_initialization'    , -12 ,'hti'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    211212   sn_hts = 'Ice_initialization'    , -12 ,'hts'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     
    215216   sn_tsu = 'Ice_initialization'    , -12 ,'tsu'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    216217   sn_tms = 'NOT USED'              , -12 ,'tms'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     218   !      melt ponds (be careful, sn_apd is the pond concentration (not fraction), so it differs from rn_apd) 
    217219   sn_apd = 'NOT USED'              , -12 ,'apd'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    218220   sn_hpd = 'NOT USED'              , -12 ,'hpd'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     
    232234&namdia         !   Diagnostics 
    233235!------------------------------------------------------------------------------ 
    234    ln_icediachk     = .false.         !  check online the heat, mass & salt budgets (T) or not (F) 
     236   ln_icediachk     = .false.         !  check online the heat, mass & salt budgets at each time step 
     237      !                               !     rate of ice spuriously gained/lost. For ex., rn_icechk=1. <=> 1mm/year, rn_icechk=0.1 <=> 1mm/10years                                    
     238      rn_icechk_cel =  1.             !     check at any gridcell           => stops the code if violated (and writes a file) 
     239      rn_icechk_glo =  0.1            !     check over the entire ice cover => only prints warnings 
    235240   ln_icediahsb     = .false.         !  output the heat, mass & salt budgets (T) or not (F) 
    236241   ln_icectl        = .false.         !  ice points output for debug (T or F) 
  • NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/cfgs/SHARED/namelist_ref

    r11413 r11586  
    199199   ln_flx      = .false.   !  flux formulation                          (T => fill namsbc_flx ) 
    200200   ln_blk      = .false.   !  Bulk formulation                          (T => fill namsbc_blk ) 
     201   ln_abl      = .false.   !  ABL  formulation                          (T => fill namsbc_abl ) 
    201202      !              ! Type of coupling (Ocean/Ice/Atmosphere) : 
    202203   ln_cpl      = .false.   !  atmosphere coupled   formulation          ( requires key_oasis3 ) 
     
    280281   sn_snow     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'SNOW'    ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    281282   sn_slp      = 'slp.15JUNE2009_fill'        ,    6.        , 'SLP'     ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     283/ 
     284!----------------------------------------------------------------------- 
     285&namsbc_abl    !   Atmospheric Boundary Layer formulation           (ln_abl = T) 
     286!----------------------------------------------------------------------- 
     287   cn_dir         = './'      !  root directory for the location of the ABL grid file 
     288   cn_dom         = 'dom_cfg_abl.nc' 
     289   ln_hpgls_frc   = .false. 
     290   ln_geos_winds  = .false. 
     291   nn_dyn_restore = 2         ! restoring option for dynamical ABL variables: = 0 no restoring 
     292                              !                                               = 1 equatorial restoring 
     293                              !                                               = 2 global restoring 
     294   rn_ldyn_min   =  4.5       !  magnitude of the nudging on ABL dynamics at the bottom of the ABL   [hour] 
     295   rn_ldyn_max   =  1.5       !  magnitude of the nudging on ABL dynamics at the top of the ABL   [hour] 
     296   rn_ltra_min   =  4.5       !  magnitude of the nudging on ABL tracers  at the bottom of the ABL   [hour] 
     297   rn_ltra_max   =  1.5       !  magnitude of the nudging on ABL tracers  at the top of the ABL   [hour] 
     298   nn_amxl        = 0         ! mixing length: = 0 Deardorff 80 length-scale 
     299                              !                = 1 length-scale based on the distance to the PBL height 
     300                              !                = 2 Bougeault & Lacarrere 89 length-scale 
     301   rn_Cm         = 0.0667     ! 0.126 in MesoNH 
     302   rn_Ct         = 0.1667     ! 0.143 in MesoNH 
     303   rn_Ce         = 0.4        ! 0.4   in MesoNH 
     304   rn_Ceps       = 0.7        ! 0.85  in MesoNH 
     305   rn_Rod        = 0.15       ! c0 in RMCA17 mixing length formulation (not yet implemented) 
     306   rn_Ric        = 0.139      !  Critical Richardson number (to compute PBL height and diffusivities) 
    282307/ 
    283308!----------------------------------------------------------------------- 
     
    600625   nn_ice_dta    =  0         !  = 0, bdy data are equal to the initial state 
    601626   !                          !  = 1, bdy data are read in 'bdydata   .nc' files 
    602    rn_ice_tem    = 270.       !  si3 only: arbitrary temperature of incoming sea ice 
    603    rn_ice_sal    = 10.        !  si3 only:      --   salinity           -- 
    604    rn_ice_age    = 30.        !  si3 only:      --   age                -- 
    605627   ! 
    606628   ln_tra_dmp    =.false.     !  open boudaries conditions for tracers 
     
    632654   bn_sal      = 'amm12_bdyT_tra'        ,         24.       , 'vosaline',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
    633655!* for si3 
    634 !   bn_a_i     = 'amm12_bdyT_ice'        ,         24.       , 'ileadfra',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
    635 !   bn_h_i     = 'amm12_bdyT_ice'        ,         24.       , 'iicethic',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
    636 !   bn_h_s     = 'amm12_bdyT_ice'        ,         24.       , 'isnowthi',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     656   bn_a_i      = 'amm12_bdyT_ice'        ,         24.       , 'siconc'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     657   bn_h_i      = 'amm12_bdyT_ice'        ,         24.       , 'sithic'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     658   bn_h_s      = 'amm12_bdyT_ice'        ,         24.       , 'snthic'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     659   bn_t_i      = 'NOT USED'              ,         24.       , 'sitemp'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     660   bn_t_s      = 'NOT USED'              ,         24.       , 'sntemp'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     661   bn_tsu      = 'NOT USED'              ,         24.       , 'sittop'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     662   bn_s_i      = 'NOT USED'              ,         24.       , 'sisalt'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     663   ! melt ponds (be careful, bn_aip is the pond concentration (not fraction), so it differs from rn_iceapnd) 
     664   bn_aip      = 'NOT USED'              ,         24.       , 'siapnd'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     665   bn_hip      = 'NOT USED'              ,         24.       , 'sihpnd'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     666   ! if bn_t_i etc are "not used", then define arbitrary temperatures and salinity and ponds 
     667   rn_ice_tem  = 270.         !  arbitrary temperature               of incoming sea ice 
     668   rn_ice_sal  = 10.          !       --   salinity                            -- 
     669   rn_ice_age  = 30.          !       --   age                                 -- 
     670   rn_ice_apnd = 0.2          !       --   pond fraction = a_ip/a_i            -- 
     671   rn_ice_hpnd = 0.05         !       --   pond depth                          -- 
    637672/ 
    638673!----------------------------------------------------------------------- 
  • NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/cfgs/SPITZ12/EXPREF/namelist_cfg

    r11413 r11586  
    170170   nn_ice_dta    =  1         !  = 0, bdy data are equal to the initial state 
    171171   !                          !  = 1, bdy data are read in 'bdydata   .nc' files 
    172    rn_ice_tem    = 267.       !  si3 only: arbitrary temperature of incoming sea ice 
    173    rn_ice_sal    =   6.       !  si3 only:      --   salinity           -- 
    174    rn_ice_age    = 365.       !  si3 only:      --   age                -- 
    175    ! 
    176172   nn_rimwidth   = 1          !  width of the relaxation zone 
    177173   ln_vol        = .false.    !  total volume correction (see nn_volctl parameter) 
  • NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/cfgs/SPITZ12/EXPREF/namelist_ice_cfg

    r11263 r11586  
    8080&namthd_pnd     !   Melt ponds 
    8181!------------------------------------------------------------------------------ 
    82    ln_pnd_H12       = .true.          !  activate evolutive melt ponds (from Holland et al 2012) 
    83    ln_pnd_alb       = .true.          !  melt ponds affect albedo or not 
     82   ln_pnd           = .true.          !  activate melt ponds or not 
     83     ln_pnd_H12     = .true.          !  activate evolutive melt ponds (from Holland et al 2012) 
     84     ln_pnd_alb     = .true.          !  melt ponds affect albedo or not 
    8485/ 
    8586 
  • NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/cfgs/ref_cfgs.txt

    r11413 r11586  
    1414eORCA025_ICE OCE ICE 
    1515eORCA025_ICE_ABL OCE ICE ABL 
     16eORCA025_SAS_ICE_ABL OCE SAS ICE ABL 
Note: See TracChangeset for help on using the changeset viewer.