Changeset 6909


Ignore:
Timestamp:
08/07/24 10:06:13 (6 months ago)
Author:
cetlod
Message:

NEMOv6.5 : minor improvments

Location:
CONFIG/UNIFORM/v6/NEMO_v6.5
Files:
3 added
1 deleted
20 edited

Legend:

Unmodified
Added
Removed
  • CONFIG/UNIFORM/v6/NEMO_v6.5/EXPERIMENTS/ORCA1_ICE_PISCES/core/clim/COMP/opa9.card

    r6905 r6909  
    116116GatherWithInternal = (nav_lon, nav_lat, deptht, time_counter, time_centered, time_centered_bounds) 
    117117TimeSeriesVars2D = (zos, tos, sos, mldr10_1, nshfls, rsntds, rsds, friver, hc300, wfo) 
    118 ChunckJob2D = NONE 
     118ChunckJob2D = 100Y 
    119119TimeSeriesVars3D = (thetao, so) 
    120 ChunckJob3D = NONE 
     120ChunckJob3D = 100Y 
    121121Seasonal=ON 
    122122 
  • CONFIG/UNIFORM/v6/NEMO_v6.5/EXPERIMENTS/ORCA1_ICE_PISCES/core/clim/COMP/si3.card

    r6142 r6909  
    3838GatherWithInternal = (nav_lon, nav_lat, time_counter, time_centered, time_centered_bounds) 
    3939TimeSeriesVars2D = (snthic, sithic, sivolu, snvolu, siconc, sisali, sitemp, sistem, siheco, snheco, sivelo, sishea, sfxice) 
    40 ChunckJob2D = NONE 
     40ChunckJob2D = 100Y 
    4141TimeSeriesVars3D = () 
    4242ChunckJob3D = NONE 
  • CONFIG/UNIFORM/v6/NEMO_v6.5/EXPERIMENTS/ORCA1_ICE_PISCES/core/ia/COMP/opa9.card

    r6905 r6909  
    111111GatherWithInternal = (nav_lon, nav_lat, deptht, time_counter, time_centered, time_centered_bounds) 
    112112TimeSeriesVars2D = (zos, tos, sos, mldr10_1, nshfls, rsntds, rsds, friver, hc300, wfo) 
    113 ChunckJob2D = NONE 
     113ChunckJob2D = 100Y 
    114114TimeSeriesVars3D = (thetao, so) 
    115 ChunckJob3D = NONE 
     115ChunckJob3D = 100Y 
    116116Seasonal=ON 
    117117 
  • CONFIG/UNIFORM/v6/NEMO_v6.5/EXPERIMENTS/ORCA1_ICE_PISCES/core/ia/COMP/si3.card

    r6142 r6909  
    3838GatherWithInternal = (nav_lon, nav_lat, time_counter, time_centered, time_centered_bounds) 
    3939TimeSeriesVars2D = (snthic, sithic, sivolu, snvolu, siconc, sisali, sitemp, sistem, siheco, snheco, sivelo, sishea, sfxice) 
    40 ChunckJob2D = NONE 
     40ChunckJob2D = 100Y 
    4141TimeSeriesVars3D = () 
    4242ChunckJob3D = NONE 
  • CONFIG/UNIFORM/v6/NEMO_v6.5/EXPERIMENTS/ORCA1_ICE_PISCES/jra/ia/COMP/si3.card

    r6142 r6909  
    3838GatherWithInternal = (nav_lon, nav_lat, time_counter, time_centered, time_centered_bounds) 
    3939TimeSeriesVars2D = (snthic, sithic, sivolu, snvolu, siconc, sisali, sitemp, sistem, siheco, snheco, sivelo, sishea, sfxice) 
    40 ChunckJob2D = NONE 
     40ChunckJob2D = 100Y 
    4141TimeSeriesVars3D = () 
    4242ChunckJob3D = NONE 
  • CONFIG/UNIFORM/v6/NEMO_v6.5/EXPERIMENTS/ORCA1_ICE_TRC/core/clim/COMP/opa9.card

    r6905 r6909  
    116116GatherWithInternal = (nav_lon, nav_lat, deptht, time_counter, time_centered, time_centered_bounds) 
    117117TimeSeriesVars2D = (zos, tos, sos, mldr10_1, nshfls, rsntds, rsds, friver, hc300, wfo) 
    118 ChunckJob2D = NONE 
     118ChunckJob2D = 100Y 
    119119TimeSeriesVars3D = (thetao, so) 
    120 ChunckJob3D = NONE 
     120ChunckJob3D = 100Y 
    121121Seasonal=ON 
    122122 
  • CONFIG/UNIFORM/v6/NEMO_v6.5/EXPERIMENTS/ORCA1_ICE_TRC/core/clim/COMP/si3.card

    r6142 r6909  
    3838GatherWithInternal = (nav_lon, nav_lat, time_counter, time_centered, time_centered_bounds) 
    3939TimeSeriesVars2D = (snthic, sithic, sivolu, snvolu, siconc, sisali, sitemp, sistem, siheco, snheco, sivelo, sishea, sfxice) 
    40 ChunckJob2D = NONE 
     40ChunckJob2D = 100Y 
    4141TimeSeriesVars3D = () 
    4242ChunckJob3D = NONE 
  • CONFIG/UNIFORM/v6/NEMO_v6.5/EXPERIMENTS/ORCA1_ICE_TRC/core/ia/COMP/opa9.card

    r6905 r6909  
    111111GatherWithInternal = (nav_lon, nav_lat, deptht, time_counter, time_centered, time_centered_bounds) 
    112112TimeSeriesVars2D = (zos, tos, sos, mldr10_1, nshfls, rsntds, rsds, friver, hc300, wfo) 
    113 ChunckJob2D = NONE 
     113ChunckJob2D = 100Y 
    114114TimeSeriesVars3D = (thetao, so) 
    115 ChunckJob3D = NONE 
     115ChunckJob3D = 100Y 
    116116Seasonal=ON 
    117117 
  • CONFIG/UNIFORM/v6/NEMO_v6.5/EXPERIMENTS/ORCA1_ICE_TRC/core/ia/COMP/si3.card

    r6142 r6909  
    3838GatherWithInternal = (nav_lon, nav_lat, time_counter, time_centered, time_centered_bounds) 
    3939TimeSeriesVars2D = (snthic, sithic, sivolu, snvolu, siconc, sisali, sitemp, sistem, siheco, snheco, sivelo, sishea, sfxice) 
    40 ChunckJob2D = NONE 
     40ChunckJob2D = 100Y 
    4141TimeSeriesVars3D = () 
    4242ChunckJob3D = NONE 
  • CONFIG/UNIFORM/v6/NEMO_v6.5/EXPERIMENTS/ORCA1_ICE_TRC/jra/ia/COMP/opa9.card

    r6905 r6909  
    111111GatherWithInternal = (nav_lon, nav_lat, deptht, time_counter, time_centered, time_centered_bounds) 
    112112TimeSeriesVars2D = (zos, tos, sos, mldr10_1, nshfls, rsntds, rsds, friver, hc300, wfo) 
    113 ChunckJob2D = NONE 
     113ChunckJob2D = 100Y 
    114114TimeSeriesVars3D = (thetao, so) 
    115 ChunckJob3D = NONE 
     115ChunckJob3D = 100Y 
    116116Seasonal=ON 
    117117 
  • CONFIG/UNIFORM/v6/NEMO_v6.5/EXPERIMENTS/ORCA1_ICE_TRC/jra/ia/COMP/si3.card

    r6142 r6909  
    3838GatherWithInternal = (nav_lon, nav_lat, time_counter, time_centered, time_centered_bounds) 
    3939TimeSeriesVars2D = (snthic, sithic, sivolu, snvolu, siconc, sisali, sitemp, sistem, siheco, snheco, sivelo, sishea, sfxice) 
    40 ChunckJob2D = NONE 
     40ChunckJob2D = 100Y 
    4141TimeSeriesVars3D = () 
    4242ChunckJob3D = NONE 
  • CONFIG/UNIFORM/v6/NEMO_v6.5/GENERAL/DRIVER/pisces.driver

    r6618 r6909  
    258258      if [ -f atmn2o.txt ] ; then 
    259259         line_year=`awk '/'${year}'/{print NR}' atmn2o.txt` 
    260          PisNo2=`awk 'NR == '$line_year' {print $2 }' atmn2o.txt` 
     260         PisN2o=`awk 'NR == '$line_year' {print $2 }' atmn2o.txt` 
    261261      else 
    262262         PisN2o=${pisces_UserChoices_AtmN2o} 
  • CONFIG/UNIFORM/v6/NEMO_v6.5/GENERAL/PARAM/NAMELIST/ORCA1/namelist_core_clim_cfg

    r6601 r6909  
    2525   ln_rstart   =  _AUTO_ 
    2626   nn_rstctl   =  _AUTO_ 
     27   cn_ocerst_in    = "restart_oce_in"   !  suffix of ocean restart name (input) 
    2728   nn_it000    =  _AUTO_ 
    2829   nn_itend    =  _AUTO_ 
     
    3637!----------------------------------------------------------------------- 
    3738   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    38    rn_Dt      = 2700.     !  time step for the dynamics and tracer 
     39   rn_Dt      = 3600.     !  time step for the dynamics and tracer 
    3940   ln_meshmask =  _AUTO_ 
    4041/ 
     
    8485&namsbc        !   Surface Boundary Condition manager                   (default: NO selection) 
    8586!----------------------------------------------------------------------- 
    86    nn_fsbc     = 2         !  frequency of SBC module call 
     87   nn_fsbc     = 3         !  frequency of SBC module call 
    8788                           !     (also = the frequency of sea-ice & iceberg model call) 
    8889                     ! Type of air-sea fluxes  
     
    9495   ln_dm2dc    = .true.   !  daily mean to diurnal cycle on short wave 
    9596   ln_ssr      = .true.   !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
    96    nn_fwb      = 2         !  FreshWater Budget: =0 unchecked 
     97   nn_fwb      = 1         !  FreshWater Budget: =0 unchecked 
    9798      !                    !     =1 global mean of e-p-r set to zero at each time step 
    9899      !                    !     =2 annual global mean of e-p-r set to zero 
     
    105106   !                    !  bulk algorithm : 
    106107   ln_NCAR    = .true.     ! "NCAR"      algorithm   (Large and Yeager 2008) 
     108   !! Bulk transfer coefficients over sea-ice: (relevant IF: nn_ice >=1 ) 
     109   ln_Cx_ice_cst = .true.     ! use constant ice-air bulk transfer coefficients (value given below) 
     110      rn_Cd_i  = 1.e-3       ! sea-ice drag coefficient 
     111      rn_Ce_i  = 1.e-3       !    "    sublimation coefficient 
     112      rn_Ch_i  = 1.e-3       !    "    sensible heat flux coefficient 
    107113 
    108114   cn_dir = './'  !  root directory for the bulk data location 
     
    313319!----------------------------------------------------------------------- 
    314320   ln_trabbl   = .true.    !  Bottom Boundary Layer parameterisation flag 
    315       nn_bbl_ldf  =  1        !  diffusive bbl (=1)   or not (=0) 
    316       nn_bbl_adv  =  0        !  advective bbl (=1/2) or not (=0) 
     321      nn_bbl_ldf  =  0        !  diffusive bbl (=1)   or not (=0) 
     322      nn_bbl_adv  =  2        !  advective bbl (=1/2) or not (=0) 
    317323      rn_ahtbbl   =  1000.    !  lateral mixing coefficient in the bbl  [m2/s] 
    318324      rn_gambbl   =  10.      !  advective bbl coefficient                 [s] 
     
    347353   ln_traldf_iso   = .true.    !  iso-neutral (Standard operator) 
    348354   ln_traldf_msc   = .true.    !  Method of Stabilizing Correction      (both operators) 
     355   ln_ldfCdL       = .true.   !  revisited magnitude of aht & eiv coef by CdL 
    349356   !                       !  Coefficients: 
    350357   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     
    437444      nn_evdm  =    0            !  evd apply on tracer (=0) or on tracer and momentum (=1) 
    438445      rn_evd   =  100.           !  evd mixing coefficient [m2/s] 
    439    ln_zdfddm   = .true.    ! double diffusive mixing 
     446   ln_zdfddm   = .false.    ! double diffusive mixing 
    440447      rn_avts  =    1.e-4     !  maximum avs (vertical mixing on salinity) 
    441448      rn_hsbfr =    1.6       !  heat/salt buoyancy flux ratio 
     
    454461      nn_mxlice   = 2         ! type of scaling under sea-ice 
    455462                              !    = 2  scaling with mean sea-ice thickness 
    456    nn_etau     =   0       !  penetration of tke below the mixed layer (ML) due to NIWs 
     463   nn_etau     =   1       !  penetration of tke below the mixed layer (ML) due to NIWs 
    457464                               !        = 0 none ; = 1 add a tke source below the ML 
     465   rn_efr     = 0.08       ! fraction of surface TKE value which penetrates below the mixed layer 
    458466   ln_lc      = .true.     !  Langmuir cell parameterisation (Axell 2002) 
    459    rn_lc      =  0.20      !  coef. associated to Langmuir cells 
     467   rn_lc      =  0.25      !  coef. associated to Langmuir cells 
    460468   nn_eice     =   3       !  attenutaion of langmuir & surface wave breaking under ice 
    461469   !                       !           = 3 weighted by 1-MIN(1,4*fr_i) 
     
    465473!----------------------------------------------------------------------- 
    466474   ln_mevar    = .false.    !  variable (T) or constant (F) mixing efficiency 
    467    ln_tsdiff   = .true.     !  account for differential T/S mixing (T) or not (F) 
     475   ln_tsdiff   = .false.     !  account for differential T/S mixing (T) or not (F) 
    468476 
    469477   cn_dir      = './'      !  root directory for the iwm data location 
  • CONFIG/UNIFORM/v6/NEMO_v6.5/GENERAL/PARAM/NAMELIST/ORCA1/namelist_core_ia_cfg

    r6876 r6909  
    2525   ln_rstart   =  _AUTO_ 
    2626   nn_rstctl   =  _AUTO_ 
     27   cn_ocerst_in    = "restart_oce_in"   !  suffix of ocean restart name (input) 
    2728   nn_it000    =  _AUTO_ 
    2829   nn_itend    =  _AUTO_ 
     
    3637!----------------------------------------------------------------------- 
    3738   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    38    rn_Dt      = 2700.     !  time step for the dynamics and tracer 
     39   rn_Dt      = 3600.     !  time step for the dynamics and tracer 
    3940   ln_meshmask =  _AUTO_ 
    4041/ 
     
    8485&namsbc        !   Surface Boundary Condition manager                   (default: NO selection) 
    8586!----------------------------------------------------------------------- 
    86    nn_fsbc     = 2         !  frequency of SBC module call 
     87   nn_fsbc     = 3         !  frequency of SBC module call 
    8788                           !     (also = the frequency of sea-ice & iceberg model call) 
    8889                     ! Type of air-sea fluxes  
     
    9495   ln_dm2dc    = .true.   !  daily mean to diurnal cycle on short wave 
    9596   ln_ssr      = .true.   !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
    96    nn_fwb      = 2         !  FreshWater Budget: =0 unchecked 
     97   nn_fwb      = 1         !  FreshWater Budget: =0 unchecked 
    9798      !                    !     =1 global mean of e-p-r set to zero at each time step 
    9899      !                    !     =2 annual global mean of e-p-r set to zero 
     
    105106   !                    !  bulk algorithm : 
    106107   ln_NCAR    = .true.     ! "NCAR"      algorithm   (Large and Yeager 2008) 
     108   !! Bulk transfer coefficients over sea-ice: (relevant IF: nn_ice >=1 ) 
     109   ln_Cx_ice_cst = .true.     ! use constant ice-air bulk transfer coefficients (value given below) 
     110      rn_Cd_i  = 1.e-3       ! sea-ice drag coefficient 
     111      rn_Ce_i  = 1.e-3       !    "    sublimation coefficient 
     112      rn_Ch_i  = 1.e-3       !    "    sensible heat flux coefficient 
    107113 
    108114   cn_dir = './'  !  root directory for the bulk data location 
     
    231237         sn_isfpar_zmin = 'runoff-icb_DaiTrenberth_Depoorter.nc' ,   -1      ,'sodepmin_isf' ,   .false.  , .true.  , 'yearly'  ,  ''      ,   '' 
    232238         !* 'spe' and 'oasis' case 
    233          sn_isfpar_fwf = 'runoff-icb_DaiTrenberth_Depoorter.nc' ,   -1      ,'sornfisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   '' 
     239         sn_isfpar_fwf = 'runoff-icb_DaiTrenberth_Depoorter.nc' ,    -1      ,'sornfisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   '' 
    234240         !* 'bg03' case 
    235241         sn_isfpar_Leff = 'isfmlt_par',       0.       ,'Leff'     ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     
    313319!----------------------------------------------------------------------- 
    314320   ln_trabbl   = .true.    !  Bottom Boundary Layer parameterisation flag 
    315       nn_bbl_ldf  =  1        !  diffusive bbl (=1)   or not (=0) 
    316       nn_bbl_adv  =  0        !  advective bbl (=1/2) or not (=0) 
     321      nn_bbl_ldf  =  0        !  diffusive bbl (=1)   or not (=0) 
     322      nn_bbl_adv  =  2        !  advective bbl (=1/2) or not (=0) 
    317323      rn_ahtbbl   =  1000.    !  lateral mixing coefficient in the bbl  [m2/s] 
    318324      rn_gambbl   =  10.      !  advective bbl coefficient                 [s] 
     
    347353   ln_traldf_iso   = .true.    !  iso-neutral (Standard operator) 
    348354   ln_traldf_msc   = .true.    !  Method of Stabilizing Correction      (both operators) 
     355   ln_ldfCdL       = .true.   !  revisited magnitude of aht & eiv coef by CdL 
    349356   !                       !  Coefficients: 
    350357   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     
    437444      nn_evdm  =    0            !  evd apply on tracer (=0) or on tracer and momentum (=1) 
    438445      rn_evd   =  100.           !  evd mixing coefficient [m2/s] 
    439    ln_zdfddm   = .true.    ! double diffusive mixing 
     446   ln_zdfddm   = .false.    ! double diffusive mixing 
    440447      rn_avts  =    1.e-4     !  maximum avs (vertical mixing on salinity) 
    441448      rn_hsbfr =    1.6       !  heat/salt buoyancy flux ratio 
     
    454461      nn_mxlice   = 2         ! type of scaling under sea-ice 
    455462                              !    = 2  scaling with mean sea-ice thickness 
    456    nn_etau     =   0       !  penetration of tke below the mixed layer (ML) due to NIWs 
     463   nn_etau     =   1       !  penetration of tke below the mixed layer (ML) due to NIWs 
    457464                               !        = 0 none ; = 1 add a tke source below the ML 
     465   rn_efr     = 0.08       ! fraction of surface TKE value which penetrates below the mixed layer 
    458466   ln_lc      = .true.     !  Langmuir cell parameterisation (Axell 2002) 
    459    rn_lc      =  0.20      !  coef. associated to Langmuir cells 
     467   rn_lc      =  0.25      !  coef. associated to Langmuir cells 
    460468   nn_eice     =   3       !  attenutaion of langmuir & surface wave breaking under ice 
    461469   !                       !           = 3 weighted by 1-MIN(1,4*fr_i) 
     
    465473!----------------------------------------------------------------------- 
    466474   ln_mevar    = .false.    !  variable (T) or constant (F) mixing efficiency 
    467    ln_tsdiff   = .true.     !  account for differential T/S mixing (T) or not (F) 
     475   ln_tsdiff   = .false.     !  account for differential T/S mixing (T) or not (F) 
    468476 
    469477   cn_dir      = './'      !  root directory for the iwm data location 
  • CONFIG/UNIFORM/v6/NEMO_v6.5/GENERAL/PARAM/NAMELIST/ORCA1/namelist_jra_ia_cfg

    r6876 r6909  
    2525   ln_rstart   =  _AUTO_ 
    2626   nn_rstctl   =  _AUTO_ 
     27   cn_ocerst_in    = "restart_oce_in"   !  suffix of ocean restart name (input) 
    2728   nn_it000    =  _AUTO_ 
    2829   nn_itend    =  _AUTO_ 
     
    3637!----------------------------------------------------------------------- 
    3738   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    38    rn_Dt      = 2700.     !  time step for the dynamics and tracer 
     39   rn_Dt      = 3600.     !  time step for the dynamics and tracer 
    3940   ln_meshmask =  _AUTO_ 
    4041/ 
     
    8485&namsbc        !   Surface Boundary Condition manager                   (default: NO selection) 
    8586!----------------------------------------------------------------------- 
    86    nn_fsbc     = 2         !  frequency of SBC module call 
     87   nn_fsbc     = 3         !  frequency of SBC module call 
    8788                           !     (also = the frequency of sea-ice & iceberg model call) 
    8889                     ! Type of air-sea fluxes  
     
    9293                     ! Misc. options of sbc :  
    9394   ln_traqsr   = .true.   !  Light penetration in the ocean            (T => fill namtra_qsr) 
    94    ln_dm2dc    = .false.  !  daily mean to diurnal cycle on short wave 
     95   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave 
    9596   ln_ssr      = .true.   !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
    96    nn_fwb      = 2        !  FreshWater Budget: =0 unchecked 
    97       !                   !     =1 global mean of e-p-r set to zero at each time step 
    98       !                   !     =2 annual global mean of e-p-r set to zero 
     97   nn_fwb      = 1         !  FreshWater Budget: =0 unchecked 
     98      !                    !     =1 global mean of e-p-r set to zero at each time step 
     99      !                    !     =2 annual global mean of e-p-r set to zero 
    99100   ln_rnf      = .true.   !  runoffs                                   (T => fill namsbc_rnf) 
    100101/ 
     
    102103&namsbc_blk    !   namsbc_blk  generic Bulk formula                     (ln_blk =T) 
    103104!----------------------------------------------------------------------- 
     105   !                    !  bulk algorithm : 
    104106   !                    !  bulk algorithm : 
    105107   ln_NCAR    = .true.     ! "NCAR"      algorithm   (Large and Yeager 2008) 
     108   !! Bulk transfer coefficients over sea-ice: (relevant IF: nn_ice >=1 ) 
     109   ln_Cx_ice_cst = .true.     ! use constant ice-air bulk transfer coefficients (value given below) 
     110      rn_Cd_i  = 1.e-3       ! sea-ice drag coefficient 
     111      rn_Ce_i  = 1.e-3       !    "    sublimation coefficient 
     112      rn_Ch_i  = 1.e-3       !    "    sensible heat flux coefficient 
    106113 
    107114   cn_dir = './'  !  root directory for the bulk data location 
     
    230237         sn_isfpar_zmin = 'runoff-icb_DaiTrenberth_Depoorter.nc' ,   -1      ,'sodepmin_isf' ,   .false.  , .true.  , 'yearly'  ,  ''      ,   '' 
    231238         !* 'spe' and 'oasis' case 
    232          sn_isfpar_fwf = 'runoff-icb_DaiTrenberth_Depoorter.nc' ,   -1      ,'sornfisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   '' 
     239         sn_isfpar_fwf = 'runoff-icb_DaiTrenberth_Depoorter.nc' ,    -1      ,'sornfisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   '' 
    233240         !* 'bg03' case 
    234241         sn_isfpar_Leff = 'isfmlt_par',       0.       ,'Leff'     ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     
    312319!----------------------------------------------------------------------- 
    313320   ln_trabbl   = .true.    !  Bottom Boundary Layer parameterisation flag 
    314       nn_bbl_ldf  =  1        !  diffusive bbl (=1)   or not (=0) 
    315       nn_bbl_adv  =  0        !  advective bbl (=1/2) or not (=0) 
     321      nn_bbl_ldf  =  0        !  diffusive bbl (=1)   or not (=0) 
     322      nn_bbl_adv  =  2        !  advective bbl (=1/2) or not (=0) 
    316323      rn_ahtbbl   =  1000.    !  lateral mixing coefficient in the bbl  [m2/s] 
    317324      rn_gambbl   =  10.      !  advective bbl coefficient                 [s] 
     
    346353   ln_traldf_iso   = .true.    !  iso-neutral (Standard operator) 
    347354   ln_traldf_msc   = .true.    !  Method of Stabilizing Correction      (both operators) 
     355   ln_ldfCdL       = .true.   !  revisited magnitude of aht & eiv coef by CdL 
    348356   !                       !  Coefficients: 
    349357   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     
    436444      nn_evdm  =    0            !  evd apply on tracer (=0) or on tracer and momentum (=1) 
    437445      rn_evd   =  100.           !  evd mixing coefficient [m2/s] 
    438    ln_zdfddm   = .true.    ! double diffusive mixing 
     446   ln_zdfddm   = .false.    ! double diffusive mixing 
    439447      rn_avts  =    1.e-4     !  maximum avs (vertical mixing on salinity) 
    440448      rn_hsbfr =    1.6       !  heat/salt buoyancy flux ratio 
     
    453461      nn_mxlice   = 2         ! type of scaling under sea-ice 
    454462                              !    = 2  scaling with mean sea-ice thickness 
    455    nn_etau     =   0       !  penetration of tke below the mixed layer (ML) due to NIWs 
     463   nn_etau     =   1       !  penetration of tke below the mixed layer (ML) due to NIWs 
    456464                               !        = 0 none ; = 1 add a tke source below the ML 
     465   rn_efr     = 0.08       ! fraction of surface TKE value which penetrates below the mixed layer 
    457466   ln_lc      = .true.     !  Langmuir cell parameterisation (Axell 2002) 
    458    rn_lc      =  0.20      !  coef. associated to Langmuir cells 
     467   rn_lc      =  0.25      !  coef. associated to Langmuir cells 
    459468   nn_eice     =   3       !  attenutaion of langmuir & surface wave breaking under ice 
    460469   !                       !           = 3 weighted by 1-MIN(1,4*fr_i) 
     
    464473!----------------------------------------------------------------------- 
    465474   ln_mevar    = .false.    !  variable (T) or constant (F) mixing efficiency 
    466    ln_tsdiff   = .true.     !  account for differential T/S mixing (T) or not (F) 
     475   ln_tsdiff   = .false.     !  account for differential T/S mixing (T) or not (F) 
    467476 
    468477   cn_dir      = './'      !  root directory for the iwm data location 
  • CONFIG/UNIFORM/v6/NEMO_v6.5/GENERAL/PARAM/NAMELIST/ORCA1/namelist_ncep_ia_cfg

    r6876 r6909  
    2525   ln_rstart   =  _AUTO_ 
    2626   nn_rstctl   =  _AUTO_ 
     27   cn_ocerst_in    = "restart_oce_in"   !  suffix of ocean restart name (input) 
    2728   nn_it000    =  _AUTO_ 
    2829   nn_itend    =  _AUTO_ 
     
    3637!----------------------------------------------------------------------- 
    3738   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    38    rn_Dt      = 2700.     !  time step for the dynamics and tracer 
     39   rn_Dt      = 3600.     !  time step for the dynamics and tracer 
    3940   ln_meshmask =  _AUTO_ 
    4041/ 
     
    8485&namsbc        !   Surface Boundary Condition manager                   (default: NO selection) 
    8586!----------------------------------------------------------------------- 
    86    nn_fsbc     = 2         !  frequency of SBC module call 
     87   nn_fsbc     = 3         !  frequency of SBC module call 
    8788                           !     (also = the frequency of sea-ice & iceberg model call) 
    8889                     ! Type of air-sea fluxes  
     
    9293                     ! Misc. options of sbc :  
    9394   ln_traqsr   = .true.   !  Light penetration in the ocean            (T => fill namtra_qsr) 
    94    ln_dm2dc    = .false.  !  daily mean to diurnal cycle on short wave 
     95   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave 
    9596   ln_ssr      = .true.   !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
    96    nn_fwb      = 2        !  FreshWater Budget: =0 unchecked 
    97       !                   !     =1 global mean of e-p-r set to zero at each time step 
    98       !                   !     =2 annual global mean of e-p-r set to zero 
     97   nn_fwb      = 1         !  FreshWater Budget: =0 unchecked 
     98      !                    !     =1 global mean of e-p-r set to zero at each time step 
     99      !                    !     =2 annual global mean of e-p-r set to zero 
    99100   ln_rnf      = .true.   !  runoffs                                   (T => fill namsbc_rnf) 
    100101/ 
     
    102103&namsbc_blk    !   namsbc_blk  generic Bulk formula                     (ln_blk =T) 
    103104!----------------------------------------------------------------------- 
     105   !                    !  bulk algorithm : 
    104106   !                    !  bulk algorithm : 
    105107   ln_NCAR    = .true.     ! "NCAR"      algorithm   (Large and Yeager 2008) 
    106       rn_zqt       = 2.     !  Air temperature & humidity reference height (m) 
    107       rn_zu        = 10.     !  Wind vector reference height (m) 
     108   !! Bulk transfer coefficients over sea-ice: (relevant IF: nn_ice >=1 ) 
     109   ln_Cx_ice_cst = .true.     ! use constant ice-air bulk transfer coefficients (value given below) 
     110      rn_Cd_i  = 1.e-3       ! sea-ice drag coefficient 
     111      rn_Ce_i  = 1.e-3       !    "    sublimation coefficient 
     112      rn_Ch_i  = 1.e-3       !    "    sensible heat flux coefficient 
    108113 
    109114   cn_dir = './'  !  root directory for the bulk data location 
     
    232237         sn_isfpar_zmin = 'runoff-icb_DaiTrenberth_Depoorter.nc' ,   -1      ,'sodepmin_isf' ,   .false.  , .true.  , 'yearly'  ,  ''      ,   '' 
    233238         !* 'spe' and 'oasis' case 
    234          sn_isfpar_fwf = 'runoff-icb_DaiTrenberth_Depoorter.nc' ,   -1      ,'sornfisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   '' 
     239         sn_isfpar_fwf = 'runoff-icb_DaiTrenberth_Depoorter.nc' ,    -1      ,'sornfisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   '' 
    235240         !* 'bg03' case 
    236241         sn_isfpar_Leff = 'isfmlt_par',       0.       ,'Leff'     ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     
    314319!----------------------------------------------------------------------- 
    315320   ln_trabbl   = .true.    !  Bottom Boundary Layer parameterisation flag 
    316       nn_bbl_ldf  =  1        !  diffusive bbl (=1)   or not (=0) 
    317       nn_bbl_adv  =  0        !  advective bbl (=1/2) or not (=0) 
     321      nn_bbl_ldf  =  0        !  diffusive bbl (=1)   or not (=0) 
     322      nn_bbl_adv  =  2        !  advective bbl (=1/2) or not (=0) 
    318323      rn_ahtbbl   =  1000.    !  lateral mixing coefficient in the bbl  [m2/s] 
    319324      rn_gambbl   =  10.      !  advective bbl coefficient                 [s] 
     
    348353   ln_traldf_iso   = .true.    !  iso-neutral (Standard operator) 
    349354   ln_traldf_msc   = .true.    !  Method of Stabilizing Correction      (both operators) 
     355   ln_ldfCdL       = .true.   !  revisited magnitude of aht & eiv coef by CdL 
    350356   !                       !  Coefficients: 
    351357   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     
    438444      nn_evdm  =    0            !  evd apply on tracer (=0) or on tracer and momentum (=1) 
    439445      rn_evd   =  100.           !  evd mixing coefficient [m2/s] 
    440    ln_zdfddm   = .true.    ! double diffusive mixing 
     446   ln_zdfddm   = .false.    ! double diffusive mixing 
    441447      rn_avts  =    1.e-4     !  maximum avs (vertical mixing on salinity) 
    442448      rn_hsbfr =    1.6       !  heat/salt buoyancy flux ratio 
     
    455461      nn_mxlice   = 2         ! type of scaling under sea-ice 
    456462                              !    = 2  scaling with mean sea-ice thickness 
    457    nn_etau     =   0       !  penetration of tke below the mixed layer (ML) due to NIWs 
     463   nn_etau     =   1       !  penetration of tke below the mixed layer (ML) due to NIWs 
    458464                               !        = 0 none ; = 1 add a tke source below the ML 
     465   rn_efr     = 0.08       ! fraction of surface TKE value which penetrates below the mixed layer 
    459466   ln_lc      = .true.     !  Langmuir cell parameterisation (Axell 2002) 
    460    rn_lc      =  0.20      !  coef. associated to Langmuir cells 
     467   rn_lc      =  0.25      !  coef. associated to Langmuir cells 
    461468   nn_eice     =   3       !  attenutaion of langmuir & surface wave breaking under ice 
    462469   !                       !           = 3 weighted by 1-MIN(1,4*fr_i) 
     
    466473!----------------------------------------------------------------------- 
    467474   ln_mevar    = .false.    !  variable (T) or constant (F) mixing efficiency 
    468    ln_tsdiff   = .true.     !  account for differential T/S mixing (T) or not (F) 
     475   ln_tsdiff   = .false.     !  account for differential T/S mixing (T) or not (F) 
    469476 
    470477   cn_dir      = './'      !  root directory for the iwm data location 
  • CONFIG/UNIFORM/v6/NEMO_v6.5/GENERAL/PARAM/NAMELIST/namelist_pisces_cfg

    r6904 r6909  
    182182   bureffmin   =  0.02     ! Minimum burial efficiency 
    183183   bureffvar   =  0.08     ! Variable coef. for burial efficiency 
    184    sedsilfrac  = 0.03      ! percentage of silica loss in the sediments  sedsilfrac 
    185    sedcalfrac  = 0.99      ! percentage of calcite loss in the sediments  sedsilfrac 
     184   sedsilfrac  = 0.03      ! percentage of silica loss in the sediments 
     185   sedcalfrac  = 0.78      ! percentage of calcite loss in the sediments ( 0.78 is for river x1.5 ) 
     186   sedfactcalmin = 0.3     !: Minimum value for dissolving calcite at the bottom 
     187   sedfactcalvar = 0.7     !: Variable  value for dissolving calcite at the bottom 
    186188/ 
    187189!----------------------------------------------------------------------- 
     
    192194&nampisdmp     !   Damping 
    193195!----------------------------------------------------------------------- 
    194   nn_pisdmp   =  _AUTO_ 
     196  ln_pisdmp     = .true. 
     197  ln_pisdmp_alk = .false. 
     198  nn_pisdmp     =  _AUTO_ 
    195199/ 
    196200!----------------------------------------------------------------------- 
  • CONFIG/UNIFORM/v6/NEMO_v6.5/GENERAL/PARAM/NAMELIST/namelist_top_p4z_cfg

    r6812 r6909  
    77   ln_rsttr      = _AUTO_ 
    88   nn_rsttr      = _AUTO_ 
     9   cn_trcrst_in  = "restart_trc_in"   !  suffix of pass. sn_tracer restart name (input) 
    910   ln_top_euler  = .true. 
    1011/ 
     
    137138   sn_trccbc(14) = 'river.orca' ,    -1            , 'riverdic'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    138139   sn_trccbc(23) = 'river.orca' ,    -1            , 'riverdin'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    139    rn_trcfac(1)  = 8.333333e+01      !  ( data in Mg/m2/yr : 1e3/12/ryyss) 
    140    rn_trcfac(2)  = 8.333333e+01      !  ( 1e3 /12 ) 
    141    rn_trcfac(5)  = 3.774193e+03   !  ( 1e3 / 31. * 117 ) 
    142    rn_trcfac(7)  = 3.558719e+01   !  ( 1e3 / 28.1 ) 
    143    rn_trcfac(10) = 8.333333e+01   !  ( 1e3 / 12 
    144    rn_trcfac(14) = 4.166667e-03   !  ( 1e3 / 12 * 5e-5 ) 
    145    rn_trcfac(23) = 5.223214e+02   !  (  1e3 / 14 * 7.3125 ) 
     140!  Inputs river are multiply by 1.5 
     141   rn_trcfac(1)  = 1.25e+02 !   8.333333e+01      !  ( data in Mg/yr : 1e3/12/ryyss) 
     142   rn_trcfac(2)  = 1.25e+02 !   8.333333e+01      !  ( 1e3 /12 ) 
     143   rn_trcfac(5)  = 5.66129e+03   !   3.774193e+03   !  ( 1e3 / 31. * 117 ) 
     144   rn_trcfac(7)  = 5.33807e+01  !   3.558719e+01   !  ( 1e3 / 28.1 ) 
     145   rn_trcfac(10) = 1.25e+02 !   8.333333e+01   !  ( 1e3 / 12 
     146   rn_trcfac(14) = 6.25e-03 !   4.166667e-03   !  ( 1e3 / 12 * 5e-5 ) 
     147   rn_trcfac(23) = 7.834821e+02 !    5.223214e+02   !  (  1e3 / 14 * 7.3125 ) 
    146148   rn_cbc_time   = 3.1536e+7      !  Time scaling factor for CBC data (seconds in a year) 
    147149/ 
  • CONFIG/UNIFORM/v6/NEMO_v6.5/GENERAL/PARAM/NAMELIST/namelist_top_p4z_gas_cfg

    r6812 r6909  
    77   ln_rsttr      = _AUTO_ 
    88   nn_rsttr      = _AUTO_ 
     9   cn_trcrst_in  = "restart_trc_in"   !  suffix of pass. sn_tracer restart name (input) 
    910   ln_top_euler  = .true. 
    1011/ 
     
    151152   sn_trccbc(14) = 'river.orca' ,    -1            , 'riverdic'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    152153   sn_trccbc(23) = 'river.orca' ,    -1            , 'riverdin'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    153    rn_trcfac(1)  = 8.333333e+01      !  ( data in Mg/m2/yr : 1e3/12/ryyss) 
    154    rn_trcfac(2)  = 8.333333e+01      !  ( 1e3 /12 ) 
    155    rn_trcfac(5)  = 3.774193e+03   !  ( 1e3 / 31. * 117 ) 
    156    rn_trcfac(7)  = 3.558719e+01   !  ( 1e3 / 28.1 ) 
    157    rn_trcfac(10) = 8.333333e+01   !  ( 1e3 / 12 
    158    rn_trcfac(14) = 4.166667e-03   !  ( 1e3 / 12 * 5e-5 ) 
    159    rn_trcfac(23) = 5.223214e+02   !  (  1e3 / 14 * 7.3125 ) 
     154!  Inputs river are multiply by 1.5 
     155   rn_trcfac(1)  = 1.25e+02 !   8.333333e+01      !  ( data in Mg/yr : 1e3/12/ryyss) 
     156   rn_trcfac(2)  = 1.25e+02 !   8.333333e+01      !  ( 1e3 /12 ) 
     157   rn_trcfac(5)  = 5.66129e+03   !   3.774193e+03   !  ( 1e3 / 31. * 117 ) 
     158   rn_trcfac(7)  = 5.33807e+01  !   3.558719e+01   !  ( 1e3 / 28.1 ) 
     159   rn_trcfac(10) = 1.25e+02 !   8.333333e+01   !  ( 1e3 / 12 
     160   rn_trcfac(14) = 6.25e-03 !   4.166667e-03   !  ( 1e3 / 12 * 5e-5 ) 
     161   rn_trcfac(23) = 7.834821e+02 !    5.223214e+02   !  (  1e3 / 14 * 7.3125 ) 
    160162   rn_cbc_time   = 3.1536e+7      !  Time scaling factor for CBC data (seconds in a year) 
    161163/ 
  • CONFIG/UNIFORM/v6/NEMO_v6.5/SOURCES/ldftra.F90

    r6601 r6909  
    7373   REAL(wp), PUBLIC ::      rn_Le               !: lateral diffusive length    [m] 
    7474 
     75   LOGICAL , PUBLIC ::   ln_ldfCdL  = .FALSE.          !: revisited magnitude of aht & eiv coef by CdL 
     76 
    7577   !                                  ! Flag to control the type of lateral diffusive operator 
    7678   INTEGER, PARAMETER, PUBLIC ::   np_ERROR  =-10   ! error in specification of lateral diffusion 
     
    140142         &                 ln_traldf_iso, ln_traldf_msc  ,  rn_slpmax     ,   &   ! option for iso-neutral operator 
    141143         &                 ln_triad_iso , ln_botmix_triad, rn_sw_triad    ,   &   ! option for triad operator 
    142          &                 nn_aht_ijk_t , rn_Ud          , rn_Ld                  ! lateral eddy coefficient 
     144         &                 nn_aht_ijk_t , rn_Ud          , rn_Ld,  ln_ldfCdL      ! lateral eddy coefficient 
    143145      !!---------------------------------------------------------------------- 
    144146      ! 
     
    175177         WRITE(numout,*) '            switching triad or not               rn_sw_triad     = ', rn_sw_triad 
    176178         WRITE(numout,*) '            lateral mixing on bottom             ln_botmix_triad = ', ln_botmix_triad 
     179         WRITE(numout,*) '        magnitude of lateral mixing revisiting   ln_ldf_CdL      = ', ln_ldf_CdL  
    177180         WRITE(numout,*) '      coefficients :' 
    178181         WRITE(numout,*) '         type of time-space variation            nn_aht_ijk_t    = ', nn_aht_ijk_t 
     
    418421         !                                             !   increase to aht0 within 20N-20S 
    419422         IF( ln_ldfeiv .AND. nn_aei_ijk_t == 21 ) THEN   ! use the already computed aei. 
    420             ahtu(:,:,1) = 3._wp * aeiu(:,:,1)  ! /!\ Coeffs modified by CdL 
    421             ahtv(:,:,1) = 3._wp * aeiv(:,:,1)  ! /!\ Coeffs modified by CdL 
     423            IF( ln_ldf_CdL ) THEN 
     424               ahtu(:,:,1) = 3._wp * aeiu(:,:,1)  ! /!\ Coeffs modified by CdL 
     425               ahtv(:,:,1) = 3._wp * aeiv(:,:,1)  ! /!\ Coeffs modified by CdL 
     426            ELSE 
     427               ahtu(:,:,1) = aeiu(:,:,1)   
     428               ahtv(:,:,1) = aeiv(:,:,1)   
     429            ENDIF         
    422430         ELSE                                            ! compute aht. 
    423431            CALL ldf_eiv( kt, aht0, ahtu, ahtv, Kmm ) 
     
    693701      IF( nn_hls == 1 )   CALL lbc_lnk( 'ldftra', zaeiw(:,:), 'W', 1.0_wp )   ! lateral boundary condition 
    694702      ! 
    695       DO_2D( 0, 0, 0, 0 ) ! /!\ Coeffs halved by CdL 
    696          paeiu(ji,jj,1) = r1_4 * ( zaeiw(ji,jj) + zaeiw(ji+1,jj  ) ) * umask(ji,jj,1) 
    697          paeiv(ji,jj,1) = r1_4 * ( zaeiw(ji,jj) + zaeiw(ji  ,jj+1) ) * vmask(ji,jj,1) 
    698       END_2D 
     703      IF( ln_ldf_CdL ) THEN 
     704         DO_2D( 0, 0, 0, 0 ) ! /!\ Coeffs halved by CdL 
     705            paeiu(ji,jj,1) = r1_4 * ( zaeiw(ji,jj) + zaeiw(ji+1,jj  ) ) * umask(ji,jj,1) 
     706            paeiv(ji,jj,1) = r1_4 * ( zaeiw(ji,jj) + zaeiw(ji  ,jj+1) ) * vmask(ji,jj,1) 
     707         END_2D 
     708      ELSE 
     709         DO_2D( 0, 0, 0, 0 )  
     710            paeiu(ji,jj,1) = r1_2 * ( zaeiw(ji,jj) + zaeiw(ji+1,jj  ) ) * umask(ji,jj,1) 
     711            paeiv(ji,jj,1) = r1_2 * ( zaeiw(ji,jj) + zaeiw(ji  ,jj+1) ) * vmask(ji,jj,1) 
     712         END_2D 
     713      ENDIF 
    699714      CALL lbc_lnk( 'ldftra', paeiu(:,:,1), 'U', 1.0_wp , paeiv(:,:,1), 'V', 1.0_wp )      ! lateral boundary condition 
    700715 
Note: See TracChangeset for help on using the changeset viewer.