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 14063 for NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/cfgs/SHARED – NEMO

Ignore:
Timestamp:
2020-12-03T17:45:05+01:00 (4 years ago)
Author:
laurent
Message:

Catch up with trunk at r14060

Location:
NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/cfgs/SHARED
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/cfgs/SHARED/axis_def_nemo.xml

    r14021 r14063  
    1414  <axis id="depthv"  long_name="Vertical V levels" unit="m" positive="down" /> 
    1515  <axis id="depthw"  long_name="Vertical W levels" unit="m" positive="down" /> 
     16  <axis id="depthf"  long_name="Vertical F levels" unit="m" positive="down" /> 
    1617  <axis id="nfloat"  long_name="Float number"      unit="-"                 /> 
    1718  <axis id="icbcla"  long_name="Iceberg class"      unit="1"               /> 
  • NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/cfgs/SHARED/domain_def_nemo.xml

    r14021 r14063  
    181181  <domain id="EqW" domain_ref="grid_W" > <zoom_domain id="EqW"/> </domain> 
    182182 
     183  <!--   F grid   --> 
     184  <domain id="grid_F" long_name="grid F"/> 
     185 
    183186  <!--   zonal mean grid   --> 
    184187  <domain_group id="gznl"> 
  • NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/cfgs/SHARED/field_def_nemo-oce.xml

    r14021 r14063  
    171171    <field id="tosmint_pot"  long_name="vertical integral of potential temperature times density"   standard_name="integral_wrt_depth_of_product_of_density_and_potential_temperature"  unit="(kg m2) degree_C" /> 
    172172 
    173  
     173    <field id="ht"           long_name="water column height at T point"                     standard_name="water_column_height_T"                      unit="m" /> 
    174174    <field id="ssh"          long_name="sea surface height"                                 standard_name="sea_surface_height_above_geoid"             unit="m" /> 
    175175    <field id="ssh2"         long_name="square of sea surface height"                       standard_name="square_of_sea_surface_height_above_geoid"   unit="m2" > ssh * ssh </field > 
     
    190190 
    191191    <!-- Energy - horizontal divergence --> 
     192    <field id="sKE"          long_name="surface kinetic energy"  standard_name="specific_kinetic_energy_of_sea_water"   unit="m2/s2"  grid_ref="grid_T_2D" /> 
    192193    <field id="hdiv"         long_name="horizontal divergence"                                                          unit="s-1"    grid_ref="grid_T_3D" /> 
    193194 
     
    270271 
    271272  <field_group id="OSMOSIS_T" grid_ref="grid_T_2D"> 
     273    <field id="hml"                 long_name="mixed layr depth"                         unit="m"       /> 
     274    <field id="hbl"                 long_name="boundary layer depth"                     unit="m"       /> 
     275    <field id="dh"                  long_name="Pycnocline thickness"                     unit=" m"      /> 
     276    <field id="ibld"                long_name="index of boundary layer depth"            unit="#"       /> 
     277    <field id="imld"                long_name="index of mixed layer depth"            unit="#"       /> 
     278    <field id="zhbl"                long_name="boundary layer depth -grid"                     unit="m"       /> 
     279    <field id="zhml"                long_name="mixed layer depth - grid"                        unit="m"       /> 
     280    <field id="zdh"                 long_name="Pycnocline  depth - grid"                 unit=" m"      /> 
     281    <field id="zustke"              long_name="magnitude of stokes drift  at T-points"   unit="m/s"     /> 
     282    <field id="us_x"        long_name="i component of active Stokes drift"                      unit="m/s"     /> 
     283    <field id="us_y"        long_name="j component of active Stokes drift"                      unit="m/s"     /> 
     284    <field id="dstokes"             long_name="stokes drift  depth scale"                unit="m"       /> 
    272285    <field id="zwth0"               long_name="surface non-local temperature flux"       unit="deg m/s" /> 
    273286    <field id="zws0"                long_name="surface non-local salinity flux"          unit="psu m/s" /> 
    274     <field id="hbl"                 long_name="boundary layer depth"                     unit="m"       /> 
    275     <field id="hbli"                long_name="initial boundary layer depth"             unit="m"       /> 
    276     <field id="dstokes"             long_name="stokes drift  depth scale"                unit="m"       /> 
    277     <field id="zustke"              long_name="magnitude of stokes drift  at T-points"   unit="m/s"     /> 
    278287    <field id="zwstrc"              long_name="convective velocity scale"                unit="m/s"     /> 
     288    <field id="zustar"              long_name="friction velocity"                        unit="m/s"     /> 
    279289    <field id="zwstrl"              long_name="langmuir velocity scale"                  unit="m/s"     /> 
    280     <field id="zustar"              long_name="friction velocity"                        unit="m/s"     /> 
    281     <field id="zhbl"                long_name="boundary layer depth"                     unit="m"       /> 
    282     <field id="zhml"                long_name="mixed layer depth"                        unit="m"       /> 
     290    <field id="zvstr"               long_name="mixed velocity scale"                     unit="m/s"     /> 
     291    <field id="zla"                 long_name="langmuir number"                          unit="m/s"     /> 
    283292    <field id="wind_wave_abs_power" long_name="\rho |U_s| x u*^2"                        unit="mW"      /> 
    284293    <field id="wind_wave_power"     long_name="U_s \dot  tau"                            unit="mW"      /> 
    285294    <field id="wind_power"          long_name="\rho  u*^3"                               unit="mW"      /> 
    286295 
    287     <!-- extra OSMOSIS diagnostics --> 
     296    <!-- interior BL OSMOSIS diagnostics --> 
    288297    <field id="zwthav"              long_name="av turb flux of T in ml"                  unit="deg m/s" /> 
    289298    <field id="zt_ml"               long_name="av T in ml"                               unit="deg"     /> 
     299    <field id="zhol"                long_name="Hoenekker number"                         unit="#"       /> 
     300    <field id="zws_ent"            long_name="entrainment turb flux of S"                unit="10^-3 m/s" /> 
    290301    <field id="zwth_ent"            long_name="entrainment turb flux of T"               unit="deg m/s" /> 
    291     <field id="zhol"                long_name="Hoenekker number"                         unit="#"       /> 
    292     <field id="zdh"                 long_name="Pycnocline  depth - grid"                 unit=" m"      /> 
    293   </field_group> 
    294  
    295   <field_group id="OSMOSIS_W" grid_ref="grid_W_3D" operation="instant" > 
     302    <field id="zwb_ent"            long_name="entrainment turb flux of buoyancy"         unit="m^2/s^-3" /> 
     303 
     304    <field id="zdt_bl"             long_name="temperature jump at base of BL"                 unit="deg"      /> 
     305    <field id="zds_bl"             long_name="salinity jump at base of BL"                 unit="10^-3"      /> 
     306    <field id="zdb_bl"             long_name="buoyancy jump at base of BL"                 unit="m/s^2"      /> 
     307    <field id="zdu_bl"             long_name="u jump at base of BL"                       unit="m/s"      /> 
     308    <field id="zdv_bl"             long_name="v jump at base of BL"                       unit="m/s"      /> 
     309 
     310    <!-- extra OSMOSIS diagnostics for debugging --> 
     311    <field id="zsc_uw_1_0"       long_name="zsc u-momentum flux on T after Stokes"                       unit="m^2/s^2" /> 
     312    <field id="zsc_uw_1_f"       long_name="zsc u-momentum flux on T after Coriolis"                       unit="m^2/s^2" /> 
     313    <field id="zsc_vw_1_f"       long_name="zsc v-momentum flux on T after Coriolis"                       unit="m^2/s^2" /> 
     314    <field id="zsc_uw_2_f"       long_name="2nd zsc u-momentum flux on T after Coriolis"                       unit="m^2/s^2" /> 
     315    <field id="zsc_vw_2_f"       long_name="2nd zsc v-momentum flux on T after Coriolis"                       unit="m^2/s^2" /> 
     316    <field id="zuw_bse"       long_name="base u-flux T-points"                          unit="m^2/s^2" /> 
     317    <field id="zvw_bse"       long_name="base v-flux T-points"                          unit="m^2/s^2" /> 
     318 
     319    <!-- FK_OSM OSMOSIS diagnostics (require also ln_osm_mle=.true.--> 
     320    <field id="hmle"          long_name="OBL FK-layer thickness"                                     unit="m"        /> 
     321    <field id="mld_prof"              long_name="FK-layer depth index"                  unit="#" /> 
     322    <field id="zmld"          long_name="target FK-layer thickness"                                     unit="m"        /> 
     323    <field id="zwb_fk"          long_name="FK b-flux"                                     unit="m^2 s^-3"        /> 
     324    <field id="zwb_fk_b"          long_name="layer averaged FK b-flux"                 unit="m^2 s^-3"       /> 
     325    <field id="zdiff_mle"          long_name="max FK diffusivity in MLE"       unit=" 10^-4 m^2 s^-1"       /> 
     326    <field id="zvel_mle"          long_name="FK velocity scale in MLE"       unit=" m s^-1"       /> 
     327  </field_group> 
     328 
     329  <field_group id="OSMOSIS_W" grid_ref="grid_W_3D" > 
     330    <field id="zviscos"       long_name="BL viscosity"   unit="m^2/s" /> 
    296331    <field id="ghamt"       long_name="non-local temperature flux"                       unit="deg m/s" /> 
    297332    <field id="ghams"       long_name="non-local salinity flux"                          unit="psu m/s" /> 
    298333    <field id="zdtdz_pyc"   long_name="Pycnocline temperature gradient"                  unit=" deg/m"  /> 
     334    <field id="zdsdz_pyc"   long_name="Pycnocline salinity gradient"                  unit=" 10^-3/m"  /> 
     335    <field id="zdbdz_pyc"   long_name="Pycnocline buoyancy gradient"                  unit=" s^-2"  /> 
     336    <field id="zdudz_pyc"   long_name="Pycnocline u gradient"                  unit=" s^-2"  /> 
     337    <field id="zdvdz_pyc"   long_name="Pycnocline v gradient"                  unit=" s^-2"  /> 
     338 
     339    <!-- extra OSMOSIS diagnostics for debugging --> 
     340    <field id="ghamu_00"       long_name="initial non-local u-momentum flux"   unit="m^2/s^2" /> 
     341    <field id="ghamv_00"       long_name="initial non-local v-momentum flux"   unit="m^2/s^2" /> 
     342    <field id="ghamu_0"       long_name="after dstokes non-local u-momentum flux"   unit="m^2/s^2" /> 
     343    <field id="ghamu_f"       long_name="after Coriolis non-local u-momentum flux"   unit="m^2/s^2" /> 
     344    <field id="ghamv_f"       long_name="after Coriolis  non-local v-momentum flux"   unit="m^2/s^2" /> 
     345    <field id="ghamu_b"       long_name="after buoyancy added non-local u-momentum flux"   unit="m^2/s^2" /> 
     346    <field id="ghamv_b"       long_name="after buoyancy added  non-local v-momentum flux"  unit="m^2/s^2" /> 
     347    <field id="ghamu_1"       long_name="after entrainment non-local u-momentum flux"   unit="m^2/s^2" /> 
     348    <field id="ghamv_1"       long_name="after entrainment  non-local v-momentum flux"  unit="m^2/s^2" /> 
    299349  </field_group> 
    300350 
    301351  <field_group id="OSMOSIS_U" grid_ref="grid_U_2D" > 
    302352    <field id="ghamu"       long_name="non-local u-momentum flux"   grid_ref="grid_U_3D" unit="m^2/s^2" /> 
    303     <field id="us_x"        long_name="i component of Stokes drift"                      unit="m/s"     /> 
     353    <!-- FK_OSM OSMOSIS diagnostics (require also ln_osm_mle=.true.--> 
     354    <field id="zdtdx"          long_name="FK  T x-gradient"                                     unit=" deg C m^-1"        /> 
     355    <field id="zdsdx"          long_name="FK  S x-gradient"                                     unit=" 10^-3 m^-1"        /> 
     356    <field id="dbdx_mle"          long_name="FK  B x-gradient"                                     unit=" s^-2"        /> 
    304357  </field_group> 
    305358 
    306359  <field_group id="OSMOSIS_V" grid_ref="grid_V_2D" > 
    307360    <field id="ghamv"       long_name="non-local v-momentum flux"   grid_ref="grid_V_3D" unit="m^2/s^2" /> 
    308     <field id="us_y"        long_name="j component of Stokes drift"                      unit="m/s"     /> 
     361    <!-- FK_OSM OSMOSIS diagnostics (require also ln_osm_mle=.true.--> 
     362    <field id="zdtdy"          long_name="FK T y-gradient"                                     unit=" deg C m^-1"        /> 
     363    <field id="zdsdy"          long_name="FK S y-gradient"                                     unit=" 10^-3 m^-1"        /> 
     364    <field id="dbdy_mle"          long_name="FK B y-gradient"                                     unit=" s^-2"        /> 
    309365  </field_group> 
    310366 
     
    503559 
    504560  <field_group id="grid_U"   grid_ref="grid_U_2D"> 
     561    <field id="hu"            long_name="water column height at U point"                         standard_name="water_column_height_U"       unit="m" /> 
    505562    <field id="e2u"           long_name="U-cell width in meridional direction"                   standard_name="cell_width"                  unit="m"                               /> 
    506563    <field id="e3u"           long_name="U-cell thickness"                                       standard_name="cell_thickness"              unit="m"          grid_ref="grid_U_3D" /> 
     
    573630    <field id="e3v"          long_name="V-cell thickness"                                       standard_name="cell_thickness"              unit="m"          grid_ref="grid_V_3D" /> 
    574631    <field id="e3v_0"        long_name="Initial V-cell thickness"                               standard_name="ref_cell_thickness"          unit="m"          grid_ref="grid_V_3D" /> 
     632    <field id="hv"            long_name="water column height at V point"                        standard_name="water_column_height_V"       unit="m" /> 
    575633    <field id="vtau"         long_name="Wind Stress along j-axis"                               standard_name="surface_downward_y_stress"   unit="N/m2"                            /> 
    576634    <field id="voce"         long_name="ocean current along j-axis"                             standard_name="sea_water_y_velocity"        unit="m/s"        grid_ref="grid_V_3D" /> 
     
    681739 
    682740  <!-- F grid --> 
     741  <field_group id="grid_F" grid_ref="grid_F_2D"> 
     742    <field id="e3f"          long_name="F-cell thickness"                    standard_name="cell_thickness"        unit="m"   grid_ref="grid_F_3D" /> 
     743    <field id="e3f_0"        long_name="F-cell thickness"                    standard_name="cell_thickness"        unit="m"   grid_ref="grid_F_3D" /> 
     744    <field id="hf"           long_name="water column height at F point"    standard_name="water_column_height_F"  unit="m"                     /> 
     745    <field id="sKEf"         long_name="surface kinetic energy at F point" standard_name="specific_kinetic_energy_of_sea_water"   unit="m2/s2" /> 
     746    <field id="relvor"       long_name="relative vorticity"                standard_name="relative_vorticity"     unit="1/s"                   /> 
     747    <field id="plavor"       long_name="planetary vorticity"               standard_name="planetary_vorticity"    unit="1/s"                   /> 
     748    <field id="relpotvor"    long_name="relative potential vorticity"      standard_name="relpot_vorticity"       unit="1/m.s"                 /> 
     749    <field id="abspotvor"    long_name="absolute potential vorticity"      standard_name="abspot_vorticity"       unit="1/m.s"                 /> 
     750    <field id="Ens"          long_name="enstrophy"                         standard_name="enstrophy"              unit="1/m2.s2"               /> 
     751  </field_group> 
     752 
    683753  <!-- AGRIF sponge --> 
    684754  <field id="agrif_spf"    long_name=" AGRIF f-sponge coefficient"   unit=" " /> 
     
    843913    <field id="strd_zdfp"     long_name="salinity   -trend: pure vert. diffusion"   unit="1e-3/s" /> 
    844914 
     915    <!-- ln_zdfosm=T only (OSMOSIS-OBL) --> 
     916    <field id="ttrd_osm"      long_name="temperature-trend: OSM-OSBL non-local forcing"                             unit="degC/s" /> 
     917    <field id="strd_osm"      long_name="salinity   -trend: OSM-OSBL non-local forcing"                             unit="1e-3/s" /> 
     918 
     919 
    845920    <!-- --> 
    846921    <field id="ttrd_dmp"      long_name="temperature-trend: interior restoring"        unit="degC/s" /> 
     
    879954    <field id="strd_zdfp_e3t"     unit="1e-3/s * m"  >  strd_zdfp * e3t </field> 
    880955 
     956    <!-- ln_zdfosm=T only (OSMOSIS-OBL) --> 
     957    <field id="ttrd_osm_e3t"      long_name="temperature-trend: OSM-OSBL non-local forcing"                             unit="degC/s * m" >  ttrd_osm * e3t </field> 
     958    <field id="strd_osm_e3t"      long_name="salinity   -trend: OSM-OSBL non-local forcing"                             unit="1e-3/s * m" >  strd_osm * e3t </field> 
     959 
    881960    <!-- --> 
    882961    <field id="ttrd_dmp_e3t"      unit="degC/s * m"  >  ttrd_dmp * e3t </field> 
     
    894973    <field id="ttrd_totad_li"    long_name="layer integrated heat-trend: total advection"         unit="W/m^2"     > ttrd_totad_e3t * 1026.0 * 3991.86795711963 </field> 
    895974    <field id="strd_totad_li"    long_name="layer integrated salt-trend: total advection"         unit="kg/(m^2 s)"    > strd_totad_e3t * 1026.0 * 0.001  </field> 
     975    <field id="ttrd_osm_li"    long_name="layer integrated heat-trend: non-local OSM"         unit="W/m^2"     > ttrd_osm_e3t * 1026.0 * 3991.86795711963 </field> 
     976    <field id="strd_osm_li"    long_name="layer integrated salt-trend: non-local OSM"         unit="kg/(m^2 s)"    > strd_osm_e3t * 1026.0 * 0.001  </field> 
    896977    <field id="ttrd_evd_li"      long_name="layer integrated heat-trend: EVD convection"          unit="W/m^2"    > ttrd_evd_e3t * 1026.0 * 3991.86795711963 </field> 
    897978    <field id="strd_evd_li"      long_name="layer integrated salt-trend: EVD convection"          unit="kg/(m^2 s)"  > strd_evd_e3t * 1026.0 * 0.001  </field> 
     
    11011182  </field_group> 
    11021183 
     1184  <!-- TMB diagnostic output --> 
     1185  <field_group  id="1h_grid_T_tmb" grid_ref="grid_T_2D" operation="instant"> 
     1186    <field id="top_temp"           name="votemper_top"  unit="degC"  /> 
     1187    <field id="mid_temp"           name="votemper_mid"  unit="degC"  /> 
     1188    <field id="bot_temp"           name="votemper_bot"  unit="degC"  /> 
     1189    <field id="top_sal"            name="vosaline_top"  unit="psu"   /> 
     1190    <field id="mid_sal"            name="vosaline_mid"  unit="psu"   /> 
     1191    <field id="bot_sal"            name="vosaline_bot"  unit="psu"   /> 
     1192    <field id="sshnmasked"         name="sossheig"      unit="m"     /> 
     1193  </field_group> 
     1194 
    11031195  <field_group  id="1h_grid_U_tmb" grid_ref="grid_U_2D" operation="instant"> 
    11041196    <field id="top_u"           name="vozocrtx_top"  unit="m/s"  /> 
  • NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/cfgs/SHARED/grid_def_nemo.xml

    r14021 r14063  
    5353    <domain domain_ref="grid_W" /> 
    5454    <axis axis_ref="depthw" /> 
     55  </grid> 
     56  <!--  --> 
     57  <grid id="grid_F_2D" > 
     58    <domain domain_ref="grid_F" /> 
     59  </grid> 
     60  <!--  --> 
     61  <grid id="grid_F_3D" > 
     62    <domain domain_ref="grid_F" /> 
     63    <axis axis_ref="depthf" /> 
    5564  </grid> 
    5665  <!--  --> 
  • NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/cfgs/SHARED/namelist_ref

    r14038 r14063  
    10131013   ln_dynvor_eeT = .false. !  energy conserving scheme (een using e3t) 
    10141014   ln_dynvor_een = .false. !  energy & enstrophy scheme 
    1015       nn_een_e3f = 0          ! =0  e3f = mi(mj(e3t))/4 
    1016       !                       ! =1  e3f = mi(mj(e3t))/mi(mj( tmask)) 
     1015   ! 
    10171016   ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T)        ==>>> PLEASE DO NOT ACTIVATE 
    1018       !                    !  (f-point vorticity schemes only) 
     1017   !                       !  (f-point vorticity schemes only) 
     1018   ! 
     1019   nn_e3f_typ = 0          !  type of e3f (EEN, ENE, ENS, MIX only)  =0  e3f = mi(mj(e3t))/4 
     1020   !                       !                                         =1  e3f = mi(mj(e3t))/mi(mj( tmask)) 
    10191021/ 
    10201022!----------------------------------------------------------------------- 
     
    10261028   ln_hpg_isf  = .false.   !  s-coordinate (sco ) adapted to isf 
    10271029   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial) 
     1030      ln_hpg_djc_vnh = .true.  !  hor.  bc type for djc scheme (T=von Neumann, F=linear extrapolation) 
     1031      ln_hpg_djc_vnv = .true.  !  vert. bc type for djc scheme (T=von Neumann, F=linear extrapolation) 
    10281032   ln_hpg_prj  = .false.   !  s-coordinate (Pressure Jacobian scheme) 
    10291033/ 
     
    10481052   !                       !  Type of the operator : 
    10491053   ln_dynldf_OFF = .false.     !  No operator (i.e. no explicit diffusion) 
     1054   nn_dynldf_typ = 0           !  =0 div-rot (default)   ;   =1 symmetric 
    10501055   ln_dynldf_lap = .false.     !    laplacian operator 
    10511056   ln_dynldf_blp = .false.     !  bilaplacian operator 
     
    11781183   ln_mxl0     = .true.    !  surface mixing length scale = F(wind stress) (T) or not (F) 
    11791184      nn_mxlice    = 2        ! type of scaling under sea-ice 
    1180                               !    = 0 no scaling under sea-ice 
    1181                               !    = 1 scaling with constant sea-ice thickness 
    1182                               !    = 2 scaling with mean sea-ice thickness ( only with SI3 sea-ice model ) 
    1183                               !    = 3 scaling with maximum sea-ice thickness 
     1185      !                       !    = 0 no scaling under sea-ice 
     1186      !                       !    = 1 scaling with constant sea-ice thickness 
     1187      !                       !    = 2 scaling with mean sea-ice thickness ( only with SI3 sea-ice model ) 
     1188      !                       !    = 3 scaling with maximum sea-ice thickness 
    11841189      rn_mxlice   = 10.       ! max constant ice thickness value when scaling under sea-ice ( nn_mxlice=1) 
    11851190   rn_mxl0     =   0.04    !  surface  buoyancy lenght scale minimum value 
     
    11881193      rn_lc       =   0.15    !  coef. associated to Langmuir cells 
    11891194   nn_etau     =   1       !  penetration of tke below the mixed layer (ML) due to NIWs 
    1190                               !        = 0 none ; = 1 add a tke source below the ML 
    1191                               !        = 2 add a tke source just at the base of the ML 
    1192                               !        = 3 as = 1 applied on HF part of the stress           (ln_cpl=T) 
     1195   !                          !        = 0 none ; = 1 add a tke source below the ML 
     1196   !                          !        = 2 add a tke source just at the base of the ML 
     1197   !                          !        = 3 as = 1 applied on HF part of the stress           (ln_cpl=T) 
    11931198      rn_efr      =   0.05    !  fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2) 
    11941199      nn_htau     =   1       !  type of exponential decrease of tke penetration below the ML 
    1195                               !        = 0  constant 10 m length scale 
    1196                               !        = 1  0.5m at the equator to 30m poleward of 40 degrees 
     1200      !                       !        = 0  constant 10 m length scale 
     1201      !                       !        = 1  0.5m at the equator to 30m poleward of 40 degrees 
    11971202   nn_eice     =   1       !  attenutaion of langmuir & surface wave breaking under ice 
    11981203   !                       !           = 0 no impact of ice cover on langmuir & surface wave breaking 
     
    12311236&namzdf_osm    !   OSM vertical diffusion                               (ln_zdfosm =T) 
    12321237!----------------------------------------------------------------------- 
    1233    ln_use_osm_la = .false.      !  Use namelist  rn_osm_la 
     1238   ln_use_osm_la = .false.     !  Use   rn_osm_la 
    12341239   rn_osm_la     = 0.3         !  Turbulent Langmuir number 
    1235    rn_osm_dstokes     = 5.     !  Depth scale of Stokes drift (m) 
     1240   rn_zdfosm_adjust_sd = 1.0   ! Stokes drift reduction factor 
     1241   rn_osm_hblfrac = 0.1        ! specify top part of hbl for nn_osm_wave = 3 or 4 
     1242   rn_osm_bl_thresh   = 5.e-5      !Threshold buoyancy for deepening of OSBL base 
    12361243   nn_ave = 0                  ! choice of horizontal averaging on avt, avmu, avmv 
    12371244   ln_dia_osm = .true.         ! output OSMOSIS-OBL variables 
     
    12411248   rn_difri  =  0.005          ! max Ri# diffusivity at Ri_g = 0 (m^2/s) 
    12421249   ln_convmix  = .true.        ! Use convective instability mixing below BL 
    1243    rn_difconv = 1.             ! diffusivity when unstable below BL  (m2/s) 
     1250   rn_difconv = 1. !0.01 !1.             ! diffusivity when unstable below BL  (m2/s) 
     1251   rn_osm_dstokes     = 5.     !  Depth scale of Stokes drift (m) 
    12441252   nn_osm_wave = 0             ! Method used to calculate Stokes drift 
    12451253      !                        !  = 2: Use ECMWF wave fields 
    12461254      !                        !  = 1: Pierson Moskowitz wave spectrum 
    12471255      !                        !  = 0: Constant La# = 0.3 
    1248 / 
     1256   nn_osm_SD_reduce = 0        ! Method used to get active Stokes drift from surface value 
     1257      !                        !  = 0: No reduction 
     1258                               !  = 1: use SD avged over top 10% hbl 
     1259                               !  = 2:use surface value of SD fit to slope at rn_osm_hblfrac*hbl below surface 
     1260   ln_zdfosm_ice_shelter = .true.  ! reduce surface SD and depth scale under ice 
     1261   ln_osm_mle = .false.        !  Use integrated FK-OSM model 
     1262/ 
     1263!----------------------------------------------------------------------- 
     1264&namosm_mle    !   mixed layer eddy parametrisation (Fox-Kemper)       (default: OFF) 
     1265!----------------------------------------------------------------------- 
     1266   rn_osm_mle_ce       = 0.06      ! magnitude of the MLE (typical value: 0.06 to 0.08) 
     1267   nn_osm_mle          = 0         ! MLE type: =0 standard Fox-Kemper ; =1 new formulation 
     1268   rn_osm_mle_lf       = 5.e+3     ! typical scale of mixed layer front (meters)                      (case rn_osm_mle=0) 
     1269   rn_osm_mle_time     = 172800.   ! time scale for mixing momentum across the mixed layer (seconds)  (case rn_osm_mle=0) 
     1270   rn_osm_mle_lat      = 20.       ! reference latitude (degrees) of MLE coef.                        (case rn_mle=1) 
     1271   rn_osm_mle_rho_c =    0.01      ! delta rho criterion used to calculate MLD for FK 
     1272   rn_osm_mle_thresh  = 0.0005     ! delta b criterion used for FK MLE criterion 
     1273   rn_osm_mle_tau     = 172800.    ! time scale for FK-OSM (seconds)  (case rn_osm_mle=0) 
     1274   ln_osm_hmle_limit   = .false.   ! limit hmle to rn_osm_hmle_limit*hbl 
     1275   rn_osm_hmle_limit   = 1.2 
     1276   / 
    12491277!----------------------------------------------------------------------- 
    12501278&namzdf_mfc     !   Mass Flux Convection 
     
    13891417   ln_sstnight = .false.             ! Logical switch for calculating night-time average for SST obs 
    13901418   ln_bound_reject  = .false.        ! Logical to remove obs near boundaries in LAMs. 
     1419   ln_default_fp_indegs = .true.     ! Logical: T=> averaging footprint is in degrees, F=> in metres 
    13911420   ln_sla_fp_indegs = .true.         ! Logical for SLA: T=> averaging footprint is in degrees, F=> in metres 
    13921421   ln_sst_fp_indegs = .true.         ! Logical for SST: T=> averaging footprint is in degrees, F=> in metres 
     
    14041433   cn_gridsearchfile ='gridsearch.nc' ! Grid search file name 
    14051434   rn_gridsearchres = 0.5            ! Grid search resolution 
     1435   rn_default_avglamscl = 0.         ! Default E/W diameter of observation footprint (metres/degrees) 
     1436   rn_default_avgphiscl = 0.         ! Default N/S diameter of observation footprint (metres/degrees) 
    14061437   rn_mdtcorr  = 1.61                ! MDT  correction 
    14071438   rn_mdtcutoff = 65.0               ! MDT cutoff for computed correction 
     
    14171448   rn_sic_avgphiscl = 0.             ! N/S diameter of SIC observation footprint (metres/degrees) 
    14181449   nn_1dint = 0                      ! Type of vertical interpolation method 
    1419    nn_2dint = 0                      ! Default horizontal interpolation method 
     1450   nn_2dint_default = 0              ! Default horizontal interpolation method 
    14201451   nn_2dint_sla = 0                  ! Horizontal interpolation method for SLA 
    14211452   nn_2dint_sst = 0                  ! Horizontal interpolation method for SST 
Note: See TracChangeset for help on using the changeset viewer.