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 12749 for NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs – NEMO

Ignore:
Timestamp:
2020-04-15T15:49:42+02:00 (4 years ago)
Author:
gsamson
Message:

ABL code cleaning and small bugfixes; modify shared and ORCA2_ICE_ABL cfgs files accordingly (ticket #2419)

Location:
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/ORCA2_ICE_ABL/EXPREF/file_def_nemo-oce.xml

    r12063 r12749  
    5656     <field field_ref="t_abl" /> 
    5757     <field field_ref="q_abl" /> 
     58     <field field_ref="uvz1_abl" /> 
     59     <field field_ref="tz1_abl" /> 
     60     <field field_ref="qz1_abl" /> 
     61     <field field_ref="uvz1_dta" /> 
     62     <field field_ref="tz1_dta" /> 
     63     <field field_ref="qz1_dta" /> 
    5864     <field field_ref="pblh" /> 
    5965   </file> 

  • NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/ORCA2_ICE_ABL/EXPREF/namelist_cfg

    r12489 r12749  
    2222   cn_exp      =  "ORCA2"  !  experience name 
    2323   nn_it000    =       1   !  first time step 
    24    nn_itend    =      16   !  last  time step (std 5475) 
     24   nn_itend    =    5840   !  last  time step (std 5840) 
    2525   nn_date0    =  20130101 !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) 
    2626   nn_write    =       4   !  frequency of write in the output file   (modulo referenced to nn_it000) 
    … …  
    3131&namdom        !   time and space domain 
    3232!----------------------------------------------------------------------- 
    33    rn_Dt      = 5400.     !  time step for the dynamics and tracer 
     33   rn_Dt       = 5400.     !  time step for the dynamics and tracer 
    3434/ 
    3535!----------------------------------------------------------------------- 
    … …  
    5656   sn_sal = 'data_1m_salinity_nomask'             ,     -1.     ,'vosaline',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    5757/ 
     58 
    5859!!====================================================================== 
    5960!!            ***  Surface Boundary Condition namelists  ***          !! 
    … …  
    6970!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T) 
    7071!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T) 
    71 !!   namsbc_isf      ice shelf melting/freezing                         (ln_isfcav  =T : read (ln_read_cfg=T) or set or usr_def_zgr ) 
    72 !!   namsbc_iscpl    coupling option between land ice model and ocean   (ln_isfcav  =T) 
    7372!!   namsbc_wave     external fields from wave model                    (ln_wave    =T) 
    7473!!   namberg         iceberg floats                                     (ln_icebergs=T) 
    … …  
    7978!----------------------------------------------------------------------- 
    8079   nn_fsbc     = 1         !  frequency of SBC module call 
    81                            !     (also = the frequency of sea-ice & iceberg model call) 
    82                      ! Type of air-sea fluxes  
     80      !                    !  (control sea-ice & iceberg model call) 
     81                     ! Type of air-sea fluxes 
    8382   ln_blk      = .false.   !  Bulk formulation                          (T => fill namsbc_blk ) 
    8483   ln_abl      = .true.    !  ABL  formulation                          (T => fill namsbc_abl ) 
    8584                     ! Sea-ice : 
    86    nn_ice      = 2         !  =2 or 3 automatically for SI3 or CICE    ("key_si3" or "key_cice") 
    87                            !          except in AGRIF zoom where it has to be specified 
    88                      ! Misc. options of sbc :  
     85   nn_ice      = 2         !  =0 no ice boundary condition 
     86      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif ) 
     87      !                    !  =2 or 3 for SI3 and CICE, respectively 
     88   ln_ice_embd = .false.   !  =T embedded sea-ice (pressure + mass and salt exchanges) 
     89      !                    !  =F levitating ice (no pressure, mass and salt exchanges) 
     90                     ! Misc. options of sbc : 
    8991   ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr) 
     92   ln_dm2dc    = .true.    !  daily mean to diurnal cycle on short wave 
    9093   ln_ssr      = .true.    !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
    91    ln_dm2dc    = .true.    !  daily mean to diurnal cycle on short wave 
     94   nn_fwb      = 2         !  FreshWater Budget: =0 unchecked 
     95      !                    !     =1 global mean of e-p-r set to zero at each time step 
     96      !                    !     =2 annual global mean of e-p-r set to zero 
    9297   ln_rnf      = .true.    !  runoffs                                   (T => fill namsbc_rnf) 
    93    nn_fwb      = 2         !  FreshWater Budget:  
    94    !                       !    =2 annual global mean of e-p-r set to zero 
    9598   ln_wave     = .false.   !  Activate coupling with wave  (T => fill namsbc_wave) 
    9699   ln_cdgw     = .false.   !  Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave) 
    97    ln_sdw      = .false.   !  Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave)  
     100   ln_sdw      = .false.   !  Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave) 
    98101   nn_sdrift   =  0        !  Parameterization for the calculation of 3D-Stokes drift from the surface Stokes drift 
    99102      !                    !   = 0 Breivik 2015 parameterization: v_z=v_0*[exp(2*k*z)/(1-8*k*z)] 
    … …  
    111114   ln_COARE_3p0 = .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003) 
    112115   ln_COARE_3p6 = .false.   ! "COARE 3.6" algorithm   (Edson et al. 2013) 
    113    ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 31) 
    114       rn_zqt      = 10.     !  Air temperature & humidity reference height (m) 
    115       rn_zu       = 10.     !  Wind vector reference height (m) 
     116   ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 45r1) 
    116117      ! 
    117       ! Skin is ONLY available in ECMWF and COARE algorithms: 
    118       ln_skin_cs = .false.  !  use the cool-skin parameterization => set nn_fsbc=1 and ln_dm2dc=.true.! 
    119       ln_skin_wl = .false.  !  use the warm-layer        "        => set nn_fsbc=1 and ln_dm2dc=.true.! 
    120       ! 
    121       ln_humi_sph = .true.  !  humidity specified below in "sn_humi" is specific humidity     [kg/kg] if .true. 
    122       ln_humi_dpt = .false. !  humidity specified below in "sn_humi" is dew-point temperature   [K]   if .true. 
    123       ln_humi_rlh = .false. !  humidity specified below in "sn_humi" is relative humidity       [%]   if .true. 
     118      ln_humi_sph = .true.  !  humidity "sn_humi" is specific humidity  [kg/kg] 
     119      ln_tpot     = .false. !!GS: compute potential temperature or not 
    124120   ! 
    125    cn_dir = './'  !  root directory for the bulk data location 
     121   cn_dir      = './'      !  root directory for the bulk data location 
    126122   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________! 
    127123   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !          weights filename            ! rotation ! land/sea mask ! 
    … …  
    146142   cn_dir         = './'      !  root directory for the location of the ABL grid file 
    147143   cn_dom         = 'dom_cfg_abl_L25Z10.nc' 
     144 
     145   ln_rstart_abl  = .false. 
    148146   ln_hpgls_frc   = .true. 
    149147   ln_geos_winds  = .false. 
    150148   nn_dyn_restore = 1 
    151    rn_ldyn_min   =  7.5       !  magnitude of the nudging on ABL dynamics at the bottom of the ABL   [hour] 
     149 
     150   rn_ldyn_min   =  4.5       !  magnitude of the nudging on ABL dynamics at the bottom of the ABL   [hour] 
    152151   rn_ldyn_max   =  1.5       !  magnitude of the nudging on ABL dynamics at the top of the ABL   [hour] 
    153    rn_ltra_min   =  7.5       !  magnitude of the nudging on ABL tracers  at the bottom of the ABL   [hour] 
     152   rn_ltra_min   =  4.5       !  magnitude of the nudging on ABL tracers  at the bottom of the ABL   [hour] 
    154153   rn_ltra_max   =  1.5       !  magnitude of the nudging on ABL tracers  at the top of the ABL   [hour] 
    155154/ 
    … …  
    201200&namberg       !   iceberg parameters                                   (default: OFF) 
    202201!----------------------------------------------------------------------- 
     202   ln_icebergs = .true.    ! activate iceberg floats (force =F with "key_agrif") 
     203 
     204   cn_dir = './'  !  root directory for the location of drag coefficient files 
     205   !______!___________!___________________!______________!______________!_________!___________!__________!__________!_______________! 
     206   !      ! file name ! frequency (hours) !   variable   ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     207   !      !           !  (if <0  months)  !     name     !   (logical)  !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
     208   sn_icb =  'calving',       -1.         , 'calving'    ,   .true.     , .true.  , 'yearly'  , ''       , ''       , '' 
    203209/ 
    204210!!====================================================================== 

  • NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/SHARED/field_def_nemo-oce.xml

    r12377 r12749  
    440440          <field id="t_dta"      long_name="DTA potential temperature"     standard_name="dta_theta"      unit="K"        /> 
    441441          <field id="q_dta"      long_name="DTA specific humidity"         standard_name="dta_qspe"       unit="kg/kg"    /> 
    442           <field id="coeft"      long_name="ABL nudging coefficient"       standard_name="coeft"          unit=""         /> 
     442          <field id="u_geo"      long_name="GEO i-horizontal velocity"     standard_name="geo_x_velocity" unit="m/s"      /> 
     443          <field id="v_geo"      long_name="GEO j-horizontal velocity"     standard_name="geo_y_velocity" unit="m/s"      /> 
    443444          <field id="tke_abl"    long_name="ABL turbulent kinetic energy"  standard_name="abl_tke"        unit="m2/s2"    /> 
    444445          <field id="avm_abl"    long_name="ABL turbulent viscosity"       standard_name="abl_avm"        unit="m2/s"     /> 
    445446          <field id="avt_abl"    long_name="ABL turbulent diffusivity"     standard_name="abl_avt"        unit="m2/s"     /> 
    446           <field id="mxl_abl"    long_name="ABL mixing length"             standard_name="abl_mxl"        unit="m"        /> 
     447          <field id="mxlm_abl"   long_name="ABL master mixing length"      standard_name="abl_mxlm"       unit="m"        /> 
     448          <field id="mxld_abl"   long_name="ABL dissipative mixing length" standard_name="abl_mxld"       unit="m"        /> 
    447449   </field_group> 
    448450 

  • NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/SHARED/namelist_ref

    r12589 r12749  
    259259!----------------------------------------------------------------------- 
    260260   !                    !  bulk algorithm : 
    261    ln_NCAR      = .true.    ! "NCAR"      algorithm   (Large and Yeager 2008) 
     261   ln_NCAR      = .false.   ! "NCAR"      algorithm   (Large and Yeager 2008) 
    262262   ln_COARE_3p0 = .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003) 
    263263   ln_COARE_3p6 = .false.   ! "COARE 3.6" algorithm   (Edson et al. 2013) 
    … …  
    271271      rn_pfac    = 1.       !  multipl. factor for precipitation (total & snow) 
    272272      rn_efac    = 1.       !  multipl. factor for evaporation (0. or 1.) 
    273       rn_vfac    = 0.       !  multipl. factor for ocean & ice velocity 
     273      rn_vfac    = 1.       !  multipl. factor for ocean & ice velocity 
    274274      !                     !  used to calculate the wind stress 
    275275      !                     ! (0. => absolute or 1. => relative winds) 
    … …  
    277277      ln_skin_wl = .false.  !  use the warm-layer parameterization 
    278278      !                     !   ==> only available in ECMWF and COARE algorithms 
    279       ln_humi_sph = .true. !  humidity "sn_humi" is specific humidity  [kg/kg] 
     279      ln_humi_sph = .false. !  humidity "sn_humi" is specific humidity  [kg/kg] 
    280280      ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K] 
    281281      ln_humi_rlh = .false. !  humidity "sn_humi" is relative humidity     [%] 
     282      ln_tpot     = .true.  !!GS: compute potential temperature or not 
    282283   ! 
    283284   cn_dir      = './'      !  root directory for the bulk data location 
    … …  
    315316                              !                                               = 1 equatorial restoring 
    316317                              !                                               = 2 global restoring 
    317    rn_ldyn_min   =  4.5       !  magnitude of the nudging on ABL dynamics at the bottom of the ABL   [hour] 
    318    rn_ldyn_max   =  1.5       !  magnitude of the nudging on ABL dynamics at the top of the ABL   [hour] 
    319    rn_ltra_min   =  4.5       !  magnitude of the nudging on ABL tracers  at the bottom of the ABL   [hour] 
    320    rn_ltra_max   =  1.5       !  magnitude of the nudging on ABL tracers  at the top of the ABL   [hour] 
     318   rn_ldyn_min   =  4.5       ! dynamics nudging magnitude inside the ABL [hour] (~3 rn_Dt) 
     319   rn_ldyn_max   =  1.5       ! dynamics nudging magnitude above  the ABL [hour] (~1 rn_Dt) 
     320   rn_ltra_min   =  4.5       ! tracers  nudging magnitude inside the ABL [hour] (~3 rn_Dt) 
     321   rn_ltra_max   =  1.5       ! tracers  nudging magnitude above  the ABL [hour] (~1 rn_Dt) 
    321322   nn_amxl       =  0         ! mixing length: = 0 Deardorff 80 length-scale 
    322323                              !                = 1 length-scale based on the distance to the PBL height 

  • NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/ref_cfgs.txt

    r12377 r12749  
    77ORCA2_OFF_TRC OCE TOP OFF 
    88ORCA2_SAS_ICE OCE ICE NST SAS 
    9 ORCA2_ICE_PISCES OCE TOP ICE NST 
     9ORCA2_ICE_PISCES OCE TOP ICE NST ABL 
    1010ORCA2_ICE_ABL OCE ICE ABL 
    11 ORCA2_SAS_ICE_ABL OCE SAS ICE ABL 
    12 ORCA2_ICE OCE ICE 
    1311SPITZ12 OCE ICE 
    1412WED025 OCE ICE 
    15 eORCA025_ICE OCE ICE 
    16 eORCA025_ICE_ABL OCE ICE ABL 
    17 eORCA025_SAS_ICE_ABL OCE SAS ICE ABL 
Note: See TracChangeset for help on using the changeset viewer.