Changeset 1705


Ignore:
Timestamp:
2009-11-03T17:37:00+01:00 (12 years ago)
Author:
smasson
Message:

impact of HF winds in TKE, see ticket:585

Location:
trunk
Files:
17 edited

Legend:

Unmodified
Added
Removed
  • trunk/CONFIG/GYRE/EXP00/iodef.xml

    r1704 r1705  
    4343   <field id="qns"          description="non solar Downward Heat Flux"              unit="W/m2"     /> 
    4444   <field id="qsr"          description="Shortwave Radiation"                       unit="W/m2"     /> 
     45   <field id="taum"         description="wind stress module"                        unit="N/m2"     /> 
    4546   <field id="wspd"         description="Wind speed module at 10 m "                unit="m/s"      /> 
    4647   <field id="qrp"          description="Surface Heat Flux: Damping"                unit="W/m2"     /> 
     
    114115   <field ref="mldr10_1"     name="somxl010"  /> 
    115116   <field ref="mldkz5"       name="somixhgt"  /> 
    116    <field ref="ice_cover"    name="soicecov"  /> 
    117117   <field ref="wspd"         name="sowindsp"  /> 
    118118   <field ref="qrp"          name="sohefldp"  /> 
  • trunk/CONFIG/GYRE/EXP00/namelist

    r1634 r1705  
    178178   sn_prec     = 'precip_core',       -1.         , 'precip'   ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
    179179   sn_snow     = 'snow_core'  ,       -1.         , 'snow'     ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
     180   sn_tdif     = 'taudif_core',       24.         , 'taudif'   ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
    180181! 
    181182   cn_dir      = './'      !  root directory for the location of the bulk files 
    182183   ln_2m       = .true.    !  air temperature and humidity referenced at 2m (T) instead 10m (F) 
     184   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data ? 
    183185   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    184186/ 
     
    530532   nn_etau     =   0       !  exponentially deceasing penetration of tke due to internal & intertial waves 
    531533                           !        = 0 no penetration ( O(2 km) resolution) 
    532                            !        = 1 additional tke source 
    533                            !        = 2 additional tke source applied only at the base of the mixed layer 
     534                           !        = 1 additional tke source (rn_efr * en) 
     535                           !        = 2 additional tke source (rn_efr * en) applied only at the base of the mixed layer 
     536                           !        = 3 additional tke source (HF contribution: mean of stress module - module of mean stress) 
    534537   nn_htau     =   1       !  type of exponential decrease of tke penetration 
    535538                           !        = 0  constant 10 m length scale 
    536539                           !        = 1  0.5m at the equator to 30m at high latitudes 
    537    rn_efr      =   0.05    !  fraction of surface tke value which penetrates inside the ocean 
     540                           !        = 2  30 meters constant depth penetration 
     541                           !  otion used only if nn_etau = 1 or 2:  
     542   rn_efr      =   0.05    !     fraction of surface tke value which penetrates inside the ocean 
     543                           !  otion used only if nn_etau = 3: 
     544   rn_addhft   =  -1.e-3   !     add offset   applied to the "mean of stress module - module of mean stress" (always kept > 0) 
     545   rn_sclhft   =   1.      !     scale factor applied to the "mean of stress module - module of mean stress" 
    538546   ln_lc       = .false.   !  Langmuir cell parameterisation 
    539547   rn_lc       =   0.15    !  coef. associated to Langmuir cells 
  • trunk/CONFIG/GYRE_LOBSTER/EXP00/iodef.xml

    r1704 r1705  
    4343   <field id="qns"          description="non solar Downward Heat Flux"              unit="W/m2"     /> 
    4444   <field id="qsr"          description="Shortwave Radiation"                       unit="W/m2"     /> 
     45   <field id="taum"         description="wind stress module"                        unit="N/m2"     /> 
    4546   <field id="wspd"         description="Wind speed module at 10 m "                unit="m/s"      /> 
    4647   <field id="qrp"          description="Surface Heat Flux: Damping"                unit="W/m2"     /> 
     
    153154   <field ref="mldr10_1"     name="somxl010"  /> 
    154155   <field ref="mldkz5"       name="somixhgt"  /> 
    155    <field ref="ice_cover"    name="soicecov"  /> 
    156156   <field ref="wspd"         name="sowindsp"  /> 
    157157   <field ref="qrp"          name="sohefldp"  /> 
  • trunk/CONFIG/GYRE_LOBSTER/EXP00/namelist

    r1634 r1705  
    178178   sn_prec     = 'precip_core',       -1.         , 'precip'   ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
    179179   sn_snow     = 'snow_core'  ,       -1.         , 'snow'     ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
     180   sn_tdif     = 'taudif_core',       24.         , 'taudif'   ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
    180181! 
    181182   cn_dir      = './'      !  root directory for the location of the bulk files 
    182183   ln_2m       = .true.    !  air temperature and humidity referenced at 2m (T) instead 10m (F) 
     184   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data ? 
    183185   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    184186/ 
     
    530532   nn_etau     =   0       !  exponentially deceasing penetration of tke due to internal & intertial waves 
    531533                           !        = 0 no penetration ( O(2 km) resolution) 
    532                            !        = 1 additional tke source 
    533                            !        = 2 additional tke source applied only at the base of the mixed layer 
     534                           !        = 1 additional tke source (rn_efr * en) 
     535                           !        = 2 additional tke source (rn_efr * en) applied only at the base of the mixed layer 
     536                           !        = 3 additional tke source (HF contribution: mean of stress module - module of mean stress) 
    534537   nn_htau     =   1       !  type of exponential decrease of tke penetration 
    535538                           !        = 0  constant 10 m length scale 
    536539                           !        = 1  0.5m at the equator to 30m at high latitudes 
    537    rn_efr      =   0.05    !  fraction of surface tke value which penetrates inside the ocean 
     540                           !        = 2  30 meters constant depth penetration 
     541                           !  otion used only if nn_etau = 1 or 2:  
     542   rn_efr      =   0.05    !     fraction of surface tke value which penetrates inside the ocean 
     543                           !  otion used only if nn_etau = 3: 
     544   rn_addhft   =  -1.e-3   !     add offset   applied to the "mean of stress module - module of mean stress" (always kept > 0) 
     545   rn_sclhft   =   1.      !     scale factor applied to the "mean of stress module - module of mean stress" 
    538546   ln_lc       = .false.   !  Langmuir cell parameterisation 
    539547   rn_lc       =   0.15    !  coef. associated to Langmuir cells 
  • trunk/CONFIG/ORCA2_LIM/EXP00/1_namelist

    r1634 r1705  
    178178   sn_prec     = 'precip_core',       -1.         , 'precip'   ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
    179179   sn_snow     = 'snow_core'  ,       -1.         , 'snow'     ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
     180   sn_tdif     = 'taudif_core',       24.         , 'taudif'   ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
    180181! 
    181182   cn_dir      = './'      !  root directory for the location of the bulk files 
    182183   ln_2m       = .true.    !  air temperature and humidity referenced at 2m (T) instead 10m (F) 
     184   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data ? 
    183185   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    184186/ 
     
    529531   nn_etau     =   0       !  exponentially deceasing penetration of tke due to internal & intertial waves 
    530532                           !        = 0 no penetration ( O(2 km) resolution) 
    531                            !        = 1 additional tke source 
    532                            !        = 2 additional tke source applied only at the base of the mixed layer 
     533                           !        = 1 additional tke source (rn_efr * en) 
     534                           !        = 2 additional tke source (rn_efr * en) applied only at the base of the mixed layer 
     535                           !        = 3 additional tke source (HF contribution: mean of stress module - module of mean stress) 
    533536   nn_htau     =   1       !  type of exponential decrease of tke penetration 
    534537                           !        = 0  constant 10 m length scale 
    535538                           !        = 1  0.5m at the equator to 30m at high latitudes 
    536    rn_efr      =   0.05    !  fraction of surface tke value which penetrates inside the ocean 
     539                           !        = 2  30 meters constant depth penetration 
     540                           !  otion used only if nn_etau = 1 or 2:  
     541   rn_efr      =   0.05    !     fraction of surface tke value which penetrates inside the ocean 
     542                           !  otion used only if nn_etau = 3: 
     543   rn_addhft   =  -1.e-3   !     add offset   applied to the "mean of stress module - module of mean stress" (always kept > 0) 
     544   rn_sclhft   =   1.      !     scale factor applied to the "mean of stress module - module of mean stress" 
    537545   ln_lc       = .false.   !  Langmuir cell parameterisation 
    538546   rn_lc       =   0.15    !  coef. associated to Langmuir cells 
  • trunk/CONFIG/ORCA2_LIM/EXP00/iodef.xml

    r1682 r1705  
    4949   <field id="qrp"          description="Surface Heat Flux: Damping"                                   unit="W/m2"     /> 
    5050   <field id="erp"          description="Surface Water Flux: Damping"                                  unit="Kg/m2/s"  /> 
     51   <field id="taum"         description="wind stress module"                                           unit="N/m2"     /> 
    5152   <field id="wspd"         description="Wind speed module at 10 m "                                   unit="m/s"      /> 
    5253    
     
    5657   <field id="qsb_oce"      description="Sensible Downward Heat Flux over open ocean"                  unit="W/m2"     /> 
    5758   <field id="qla_oce"      description="Latent Downward Heat Flux over open ocean"                    unit="W/m2"     /> 
     59   <field id="taum_oce"     description="wind stress module over open ocean"                           unit="N/m2"     /> 
    5860 
    5961   <field id="ice_cover"    description="Ice fraction"                                                 unit="1"        /> 
  • trunk/CONFIG/ORCA2_LIM/EXP00/namelist

    r1696 r1705  
    178178   sn_prec     = 'precip_core',       -1.         , 'precip'   ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
    179179   sn_snow     = 'snow_core'  ,       -1.         , 'snow'     ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
     180   sn_tdif     = 'taudif_core',       24.         , 'taudif'   ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
    180181! 
    181182   cn_dir      = './'      !  root directory for the location of the bulk files 
    182183   ln_2m       = .true.    !  air temperature and humidity referenced at 2m (T) instead 10m (F) 
     184   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data ? 
    183185   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    184186/ 
     
    557559   nn_etau     =   0       !  exponentially deceasing penetration of tke due to internal & intertial waves 
    558560                           !        = 0 no penetration ( O(2 km) resolution) 
    559                            !        = 1 additional tke source 
    560                            !        = 2 additional tke source applied only at the base of the mixed layer 
     561                           !        = 1 additional tke source (rn_efr * en) 
     562                           !        = 2 additional tke source (rn_efr * en) applied only at the base of the mixed layer 
     563                           !        = 3 additional tke source (HF contribution: mean of stress module - module of mean stress) 
    561564   nn_htau     =   1       !  type of exponential decrease of tke penetration 
    562565                           !        = 0  constant 10 m length scale 
    563566                           !        = 1  0.5m at the equator to 30m at high latitudes 
    564    rn_efr      =   0.05    !  fraction of surface tke value which penetrates inside the ocean 
     567                           !        = 2  30 meters constant depth penetration 
     568                           !  otion used only if nn_etau = 1 or 2:  
     569   rn_efr      =   0.05    !     fraction of surface tke value which penetrates inside the ocean 
     570                           !  otion used only if nn_etau = 3: 
     571   rn_addhft   =  -1.e-3   !     add offset   applied to the "mean of stress module - module of mean stress" (always kept > 0) 
     572   rn_sclhft   =   1.      !     scale factor applied to the "mean of stress module - module of mean stress" 
    565573   ln_lc       = .false.   !  Langmuir cell parameterisation 
    566574   rn_lc       =   0.15    !  coef. associated to Langmuir cells 
  • trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/iodef.xml

    r1682 r1705  
    4949   <field id="qrp"          description="Surface Heat Flux: Damping"                                   unit="W/m2"     /> 
    5050   <field id="erp"          description="Surface Water Flux: Damping"                                  unit="Kg/m2/s"  /> 
     51   <field id="taum"         description="wind stress module"                                           unit="N/m2"     /> 
    5152   <field id="wspd"         description="Wind speed module at 10 m "                                   unit="m/s"      /> 
    5253    
     
    5657   <field id="qsb_oce"      description="Sensible Downward Heat Flux over open ocean"                  unit="W/m2"     /> 
    5758   <field id="qla_oce"      description="Latent Downward Heat Flux over open ocean"                    unit="W/m2"     /> 
     59   <field id="taum_oce"     description="wind stress module over open ocean"                           unit="N/m2"     /> 
    5860 
    5961   <field id="ice_cover"    description="Ice fraction"                                                 unit="1"        /> 
  • trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist

    r1696 r1705  
    178178   sn_prec     = 'precip_core',       -1.         , 'precip'   ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
    179179   sn_snow     = 'snow_core'  ,       -1.         , 'snow'     ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
     180   sn_tdif     = 'taudif_core',       24.         , 'taudif'   ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
    180181! 
    181182   cn_dir      = './'      !  root directory for the location of the bulk files 
    182183   ln_2m       = .true.    !  air temperature and humidity referenced at 2m (T) instead 10m (F) 
     184   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data ? 
    183185   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    184186/ 
     
    557559   nn_etau     =   0       !  exponentially deceasing penetration of tke due to internal & intertial waves 
    558560                           !        = 0 no penetration ( O(2 km) resolution) 
    559                            !        = 1 additional tke source 
    560                            !        = 2 additional tke source applied only at the base of the mixed layer 
     561                           !        = 1 additional tke source (rn_efr * en) 
     562                           !        = 2 additional tke source (rn_efr * en) applied only at the base of the mixed layer 
     563                           !        = 3 additional tke source (HF contribution: mean of stress module - module of mean stress) 
    561564   nn_htau     =   1       !  type of exponential decrease of tke penetration 
    562565                           !        = 0  constant 10 m length scale 
    563566                           !        = 1  0.5m at the equator to 30m at high latitudes 
    564    rn_efr      =   0.05    !  fraction of surface tke value which penetrates inside the ocean 
     567                           !        = 2  30 meters constant depth penetration 
     568                           !  otion used only if nn_etau = 1 or 2:  
     569   rn_efr      =   0.05    !     fraction of surface tke value which penetrates inside the ocean 
     570                           !  otion used only if nn_etau = 3: 
     571   rn_addhft   =  -1.e-3   !     add offset   applied to the "mean of stress module - module of mean stress" (always kept > 0) 
     572   rn_sclhft   =   1.      !     scale factor applied to the "mean of stress module - module of mean stress" 
    565573   ln_lc       = .false.   !  Langmuir cell parameterisation 
    566574   rn_lc       =   0.15    !  coef. associated to Langmuir cells 
  • trunk/CONFIG/POMME/EXP00/iodef.xml

    r1682 r1705  
    1 x<?xml version="1.0"?> 
     1<?xml version="1.0"?> 
    22 
    3 <simulation> 
     3<simulation>  
    44   
    5   <context id="nemo">  
     5  <context id="nemo"> 
    66     
    77    <!-- groups of variables definition --> 
     
    1111      <!-- T grid --> 
    1212       
    13       <group id="grid_T" axis_ref="none" grid_ref="grid_T">   
     13      <group id="grid_T" axis_ref="none" grid_ref="grid_T"> 
    1414   <field id="toce"         description="temperature"                               unit="degC" axis_ref="deptht"   /> 
    15    <field id="soce"         description="salinity"                                  unit="psu"  axis_ref="deptht"   /> 
     15         <field id="soce"         description="salinity"                                  unit="psu"  axis_ref="deptht"   /> 
    1616   <field id="sst"          description="sea surface temperature"                   unit="degC"                     /> 
    1717   <field id="sss"          description="sea surface salinity"                      unit="psu"                      /> 
     
    3434   <field id="hc300"        description="Heat content 300 m"                        unit="W"                        /> 
    3535     </group> 
    36         
     36 
    3737      <!-- SBC --> 
    3838       
    39       <group id="SBC" axis_ref="none" grid_ref="grid_T"> 
    40    <field id="emp"          description="Net Upward Water Flux"                     unit="kg/m2/s"  /> 
    41    <field id="emps"         description="concentration/dilution water flux"         unit="kg/m2/s"  /> 
    42    <field id="qns+qsr"      description="Net Downward Heat Flux"                    unit="W/m2"     /> 
    43    <field id="qns"          description="non solar Downward Heat Flux"              unit="W/m2"     /> 
    44    <field id="qsr"          description="Shortwave Radiation"                       unit="W/m2"     /> 
    45    <field id="qrp"          description="Surface Heat Flux: Damping"                unit="W/m2"     /> 
    46    <field id="erp"          description="Surface Water Flux: Damping"               unit="Kg/m2/s"  /> 
     39      <group id="SBC" axis_ref="none" grid_ref="grid_T" freq_op="28800"> 
     40 
     41   <field id="emp"          description="Net Upward Water Flux"                                        unit="kg/m2/s"  /> 
     42   <field id="emps"         description="concentration/dilution water flux"                            unit="kg/m2/s"  /> 
     43   <field id="snowpre"      description="Snow precipitation"                                           unit="kg/m2/s"  /> 
     44   <field id="runoffs"      description="River Runoffs"                                                unit="Kg/m2/s"  /> 
     45 
     46   <field id="qns+qsr"      description="Net Downward Heat Flux"                                       unit="W/m2"     /> 
     47   <field id="qns"          description="non solar Downward Heat Flux"                                 unit="W/m2"     /> 
     48   <field id="qsr"          description="Shortwave Radiation"                                          unit="W/m2"     /> 
     49   <field id="qrp"          description="Surface Heat Flux: Damping"                                   unit="W/m2"     /> 
     50   <field id="erp"          description="Surface Water Flux: Damping"                                  unit="Kg/m2/s"  /> 
     51   <field id="taum"         description="wind stress module"                                           unit="N/m2"     /> 
     52   <field id="wspd"         description="Wind speed module at 10 m "                                   unit="m/s"      /> 
     53    
     54   <!-- *_oce variables available with ln_blk_clio or ln_blk_core --> 
     55   <field id="qns_oce"      description="Non solar Downward Heat Flux over open ocean"                 unit="W/m2"     /> 
     56   <field id="qlw_oce"      description="Longwave Downward Heat Flux over open ocean"                  unit="W/m2"     /> 
     57   <field id="qsb_oce"      description="Sensible Downward Heat Flux over open ocean"                  unit="W/m2"     /> 
     58   <field id="qla_oce"      description="Latent Downward Heat Flux over open ocean"                    unit="W/m2"     /> 
     59   <field id="taum_oce"     description="wind stress module over open ocean"                           unit="N/m2"     /> 
     60 
    4761      </group> 
    4862 
    4963      <!-- U grid --> 
    5064       
    51      <group id="grid_U"  axis_ref="none" grid_ref="grid_U"> 
     65      <group id="grid_U"  axis_ref="depthu" grid_ref="grid_U"> 
    5266   <field id="utau"        description="Wind Stress along i-axis"                  unit="N/m2" axis_ref="none" /> 
    5367   <field id="uoce"        description="ocean current along i-axis"                unit="m/s"                  /> 
     
    5973     </group> 
    6074       
    61      <!-- V grid --> 
    62             
    63      <group id="grid_V"  axis_ref="none" grid_ref="grid_V"> 
     75      <group id="grid_V"  axis_ref="depthv" grid_ref="grid_V"> 
    6476   <field id="vtau"        description="Wind Stress along j-axis"                  unit="N/m2" axis_ref="none" /> 
    6577   <field id="voce"        description="ocean current along j-axis"                unit="m/s"                  /> 
     
    6981   <!-- voce_eiv: available with key_trabbl_adv --> 
    7082   <field id="voce_bbl"    description="BBL ocean current along j-axis"            unit="m/s"                  /> 
    71      </group> 
    72             
    73      <!-- W grid --> 
     83      </group> 
    7484       
    75      <group id="grid_W"  axis_ref="depthw" grid_ref="grid_W"> 
     85      <!-- W grid --> 
     86       
     87      <group id="grid_W"  axis_ref="depthw" grid_ref="grid_W"> 
    7688   <field id="woce"        description="ocean vertical velocity"                   unit="m/s"                  /> 
    7789   <field id="woce_eff"    description="effective ocean vertical velocity"         unit="m/s"                  /> 
     
    113125   <field ref="mldr10_1"     name="somxl010"  /> 
    114126   <field ref="mldkz5"       name="somixhgt"  /> 
    115    <field ref="ice_cover"    name="soicecov"  /> 
    116127   <field ref="qrp"          name="sohefldp"  /> 
    117128   <field ref="erp"          name="sowafldp"  /> 
  • trunk/CONFIG/POMME/EXP00/namelist

    r1647 r1705  
    177177   sn_humi     = 'q2_1988'        ,       24.         , 'q2'       ,    .true.      , .false. ,   'yearly' , '' 
    178178   sn_prec     = 'precip_1988.nc' ,       -1.         , 'precip'   ,    .true.      , .false. ,   'yearly' , '' 
     179   sn_tdif     = 'taudif_core'    ,       24.         , 'taudif'   ,    .true.      , .true.  ,   'yearly'  ,'' 
    179180! 
    180181   cn_dir      = './'      !  root directory for the location of the bulk files 
    181182   ln_2m       = .false.   !  air temperature and humidity referenced at 2m (T) instead 10m (F) 
     183   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data ? 
    182184   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    183185/ 
     
    194196cn_snd_crt_grid   = 'T'                     ! 'T' 
    195197                                       ! receive 
    196 cn_rcv_w10m       = 'coupled'               ! 'none' 'coupled' 
     198cn_rcv_w10m       = 'none'                  ! 'none' 'coupled' 
     199cn_rcv_taumod     = 'coupled'               ! 'none' 'coupled' 
    197200cn_rcv_tau_nature = 'oce only'              ! 'oce only' 'oce and ice' 'mixed oce-ice' 
    198201cn_rcv_tau_refere = 'cartesian'             ! 'spherical' 'cartesian' 
     
    555558   nn_etau     =   0       !  exponentially deceasing penetration of tke due to internal & intertial waves 
    556559                           !        = 0 no penetration ( O(2 km) resolution) 
    557                            !        = 1 additional tke source 
    558                            !        = 2 additional tke source applied only at the base of the mixed layer 
     560                           !        = 1 additional tke source (rn_efr * en) 
     561                           !        = 2 additional tke source (rn_efr * en) applied only at the base of the mixed layer 
     562                           !        = 3 additional tke source (HF contribution: mean of stress module - module of mean stress) 
    559563   nn_htau     =   1       !  type of exponential decrease of tke penetration 
    560564                           !        = 0  constant 10 m length scale 
    561565                           !        = 1  0.5m at the equator to 30m at high latitudes 
    562    rn_efr      =   0.05    !  fraction of surface tke value which penetrates inside the ocean 
     566                           !        = 2  30 meters constant depth penetration 
     567                           !  otion used only if nn_etau = 1 or 2:  
     568   rn_efr      =   0.05    !     fraction of surface tke value which penetrates inside the ocean 
     569                           !  otion used only if nn_etau = 3: 
     570   rn_addhft   =  -1.e-3   !     add offset   applied to the "mean of stress module - module of mean stress" (always kept > 0) 
     571   rn_sclhft   =   1.      !     scale factor applied to the "mean of stress module - module of mean stress" 
    563572   ln_lc       = .false.   !  Langmuir cell parameterisation 
    564573   rn_lc       =   0.15    !  coef. associated to Langmuir cells 
  • trunk/NEMO/OPA_SRC/SBC/sbc_oce.F90

    r1695 r1705  
    2929   !                                             !:  = 1 global mean of e-p-r set to zero at each nn_fsbc time step 
    3030   !                                             !:  = 2 annual global mean of e-p-r set to zero 
    31   INTEGER , PUBLIC ::   nn_ico_cpl  = 0          !: ice-ocean coupling indicator 
     31   INTEGER , PUBLIC ::   nn_ico_cpl  = 0          !: ice-ocean coupling indicator 
    3232   !                                             !:  = 0   LIM-3 old case 
    3333   !                                             !:  = 1   stresses computed using now ocean velocity 
     
    3737   !!              Ocean Surface Boundary Condition fields 
    3838   !!---------------------------------------------------------------------- 
     39   LOGICAL , PUBLIC ::   lhftau = .FALSE.              !: HF tau contribution: mean of stress module - module of the mean stress 
    3940   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   utau      !: sea surface i-stress (ocean referential)     [N/m2] 
    4041   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   vtau      !: sea surface j-stress (ocean referential)     [N/m2] 
    4142   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   taum      !: module of sea surface stress (at T-point)    [N/m2]  
    42    REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   wndm      !: wind speed module at T-point (=|U10m-Uoce|)  [m/s] Used only in PISCES 
     43   !! wndm is used only in PISCES to compute gases exchanges at the surface of the free ocean or in the leads in sea-ice parts 
     44   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   wndm      !: wind speed module at T-point (=|U10m-Uoce|)  [m/s]  
    4345   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   qsr       !: sea heat flux:     solar                     [W/m2] 
    4446   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   qns       !: sea heat flux: non solar                     [W/m2] 
  • trunk/NEMO/OPA_SRC/SBC/sbcblk_core.F90

    r1695 r1705  
    4343   PUBLIC   blk_ice_core       ! routine called in sbc_ice_lim module 
    4444       
    45    INTEGER , PARAMETER ::   jpfld   = 8           ! maximum number of files to read  
     45   INTEGER , PARAMETER ::   jpfld   = 9           ! maximum number of files to read  
    4646   INTEGER , PARAMETER ::   jp_wndi = 1           ! index of 10m wind velocity (i-component) (m/s)    at T-point 
    4747   INTEGER , PARAMETER ::   jp_wndj = 2           ! index of 10m wind velocity (j-component) (m/s)    at T-point 
     
    5252   INTEGER , PARAMETER ::   jp_prec = 7           ! index of total precipitation (rain+snow) (Kg/m2/s) 
    5353   INTEGER , PARAMETER ::   jp_snow = 8           ! index of snow (solid prcipitation)       (kg/m2/s) 
     54   INTEGER , PARAMETER ::   jp_tdif = 9           ! index of tau diff associated to HF tau   (N/m2)   at T-point 
    5455   TYPE(FLD), ALLOCATABLE, DIMENSION(:) ::   sf   ! structure of input fields (file informations, fields read) 
    5556          
     
    6364 
    6465   !                                !!* Namelist namsbc_core : CORE bulk parameters 
    65    LOGICAL  ::   ln_2m = .FALSE.     ! logical flag for height of air temp. and hum 
    66    REAL(wp) ::   rn_pfac = 1.        ! multiplication factor for precipitation 
     66   LOGICAL  ::   ln_2m     = .FALSE.     ! logical flag for height of air temp. and hum 
     67   LOGICAL  ::   ln_taudif = .FALSE.     ! logical flag to use the "mean of stress module - module of mean stress" data 
     68   REAL(wp) ::   rn_pfac   = 1.          ! multiplication factor for precipitation 
    6769 
    6870   !! * Substitutions 
     
    9395      !!      the total precipitation (rain+snow) (Kg/m2/s) 
    9496      !!      the snow (solid prcipitation)       (kg/m2/s) 
     97      !!   OPTIONAL parameter (see ln_taudif namelist flag): 
     98      !!      the tau diff associated to HF tau   (N/m2)   at T-point  
    9599      !!              (2) CALL blk_oce_core 
    96100      !! 
     
    110114      INTEGER  ::   ierror   ! return error code 
    111115      INTEGER  ::   ifpr     ! dummy loop indice 
     116      INTEGER  ::   jfld     ! dummy loop arguments 
    112117      !! 
    113118      CHARACTER(len=100) ::  cn_dir   !   Root directory for location of core files 
     
    115120      TYPE(FLD_N) ::   sn_wndi, sn_wndj, sn_humi, sn_qsr       ! informations about the fields to be read 
    116121      TYPE(FLD_N) ::   sn_qlw , sn_tair, sn_prec, sn_snow      !   "                                 " 
    117       NAMELIST/namsbc_core/ cn_dir, ln_2m, rn_pfac, sn_wndi, sn_wndj, sn_humi, sn_qsr,   & 
    118          &                                                sn_qlw , sn_tair, sn_prec, sn_snow 
     122      TYPE(FLD_N) ::   sn_tdif                                 !   "                                 " 
     123      NAMELIST/namsbc_core/ cn_dir , ln_2m  , ln_taudif, rn_pfac,           & 
     124         &                  sn_wndi, sn_wndj, sn_humi  , sn_qsr ,           & 
     125         &                  sn_qlw , sn_tair, sn_prec  , sn_snow, sn_tdif 
    119126      !!--------------------------------------------------------------------- 
    120127 
     
    136143         sn_prec = FLD_N( 'precip'  ,    -1.    ,  'precip'  ,  .true.    , .false. ,   'yearly'  , ''       , ''         ) 
    137144         sn_snow = FLD_N( 'snow'    ,    -1.    ,  'snow'    ,  .true.    , .false. ,   'yearly'  , ''       , ''         ) 
     145         sn_tdif = FLD_N( 'taudif'  ,    24.    ,  'taudif'  ,  .true.    , .false. ,   'yearly'  , ''       , ''         ) 
    138146         ! 
    139147         REWIND( numnam )                    ! ... read in namlist namsbc_core 
     
    145153         slf_i(jp_tair) = sn_tair   ;   slf_i(jp_humi) = sn_humi 
    146154         slf_i(jp_prec) = sn_prec   ;   slf_i(jp_snow) = sn_snow 
     155         slf_i(jp_tdif) = sn_tdif 
     156         ! 
     157         ! do we use HF tau information? 
     158         lhftau = ln_taudif 
     159         jfld = jpfld - COUNT( (/.NOT. lhftau/) ) 
    147160         ! 
    148161         ! set sf structure 
    149          ALLOCATE( sf(jpfld), STAT=ierror ) 
     162         ALLOCATE( sf(jfld), STAT=ierror ) 
    150163         IF( ierror > 0 ) THEN 
    151164            CALL ctl_stop( 'sbc_blk_core: unable to allocate sf structure' )   ;   RETURN 
    152165         ENDIF 
    153          DO ifpr= 1, jpfld 
     166         DO ifpr= 1, jfld 
    154167            ALLOCATE( sf(ifpr)%fnow(jpi,jpj) ) 
    155168            ALLOCATE( sf(ifpr)%fdta(jpi,jpj,2) ) 
     
    291304         END DO 
    292305      END DO 
     306 
     307      ! ... add the HF tau contribution to the wind stress module? 
     308      IF( lhftau ) THEN  
     309!CDIR COLLAPSE 
     310         taum(:,:) = taum(:,:) + sf(jp_tdif)%fnow(:,:) 
     311      ENDIF 
     312      CALL iom_put( "taum_oce", taum )   ! output wind stress module 
     313 
    293314      ! ... utau, vtau at U- and V_points, resp. 
    294315      !     Note the use of 0.5*(2-umask) in order to unmask the stress along coastlines 
  • trunk/NEMO/OPA_SRC/SBC/sbccpl.F90

    r1698 r1705  
    419419      !                                                      ! ------------------------- ! 
    420420      srcv(jpr_taum)%clname = 'O_TauMod'   ;   IF( TRIM(cn_rcv_taumod) == 'coupled' )   srcv(jpr_taum)%laction = .TRUE. 
     421      lhftau = srcv(jpr_taum)%laction 
    421422 
    422423#if defined key_cpl_carbon_cycle 
     
    693694         taum(:,:) = frcv(:,:,jpr_taum) 
    694695         wndm(:,:) = frcv(:,:,jpr_w10m) 
     696         CALL iom_put( "taum_oce", taum )   ! output wind stress module 
    695697         !   
    696698      ENDIF 
  • trunk/NEMO/OPA_SRC/SBC/sbcmod.F90

    r1649 r1705  
    253253      CALL iom_put( "utau", utau )   ! i-wind stress   (stress can be updated at  
    254254      CALL iom_put( "vtau", vtau )   ! j-wind stress    each time step in sea-ice) 
    255       CALL iom_put( "wspd", wndm )   ! wind speed module  
     255      CALL iom_put( "taum", taum )   ! wind stress module  
     256      CALL iom_put( "wspd", wndm )   ! wind speed  module  
    256257      ! 
    257258      IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
  • trunk/NEMO/OPA_SRC/ZDF/zdftke.F90

    r1695 r1705  
    6565 
    6666   !                                       !!! ** Namelist  namzdf_tke  ** 
    67    LOGICAL  ::   ln_mxl0  = .FALSE.         ! mixing length scale surface value as function of wind stress or not 
    68    INTEGER  ::   nn_mxl   =  2              ! type of mixing length (=0/1/2/3) 
    69    REAL(wp) ::   rn_lmin0 = 0.4_wp          ! surface  min value of mixing length   [m] 
    70    REAL(wp) ::   rn_lmin  = 0.1_wp          ! interior min value of mixing length   [m] 
    71    INTEGER  ::   nn_pdl   =  1              ! Prandtl number or not (ratio avt/avm) (=0/1) 
    72    REAL(wp) ::   rn_ediff = 0.1_wp          ! coefficient for avt: avt=rn_ediff*mxl*sqrt(e) 
    73    REAL(wp) ::   rn_ediss = 0.7_wp          ! coefficient of the Kolmogoroff dissipation  
    74    REAL(wp) ::   rn_ebb   = 3.75_wp         ! coefficient of the surface input of tke 
    75    REAL(wp) ::   rn_emin  = 0.7071e-6_wp    ! minimum value of tke           [m2/s2] 
    76    REAL(wp) ::   rn_emin0 = 1.e-4_wp        ! surface minimum value of tke   [m2/s2] 
    77    REAL(wp) ::   rn_bshear= 1.e-20          ! background shear (>0) 
    78    INTEGER  ::   nn_etau  = 0               ! type of depth penetration of surface tke (=0/1/2) 
    79    INTEGER  ::   nn_htau  = 0               ! type of tke profile of penetration (=0/1) 
    80    REAL(wp) ::   rn_efr   = 1.0_wp          ! fraction of TKE surface value which penetrates in the ocean 
    81    LOGICAL  ::   ln_lc    = .FALSE.         ! Langmuir cells (LC) as a source term of TKE or not 
    82    REAL(wp) ::   rn_lc    = 0.15_wp         ! coef to compute vertical velocity of Langmuir cells 
     67   LOGICAL  ::   ln_mxl0   = .FALSE.         ! mixing length scale surface value as function of wind stress or not 
     68   INTEGER  ::   nn_mxl    =  2              ! type of mixing length (=0/1/2/3) 
     69   REAL(wp) ::   rn_lmin0  = 0.4_wp          ! surface  min value of mixing length   [m] 
     70   REAL(wp) ::   rn_lmin   = 0.1_wp          ! interior min value of mixing length   [m] 
     71   INTEGER  ::   nn_pdl    =  1              ! Prandtl number or not (ratio avt/avm) (=0/1) 
     72   REAL(wp) ::   rn_ediff  = 0.1_wp          ! coefficient for avt: avt=rn_ediff*mxl*sqrt(e) 
     73   REAL(wp) ::   rn_ediss  = 0.7_wp          ! coefficient of the Kolmogoroff dissipation  
     74   REAL(wp) ::   rn_ebb    = 3.75_wp         ! coefficient of the surface input of tke 
     75   REAL(wp) ::   rn_emin   = 0.7071e-6_wp    ! minimum value of tke           [m2/s2] 
     76   REAL(wp) ::   rn_emin0  = 1.e-4_wp        ! surface minimum value of tke   [m2/s2] 
     77   REAL(wp) ::   rn_bshear = 1.e-20          ! background shear (>0) 
     78   INTEGER  ::   nn_etau   = 0               ! type of depth penetration of surface tke (=0/1/2) 
     79   INTEGER  ::   nn_htau   = 0               ! type of tke profile of penetration (=0/1) 
     80   REAL(wp) ::   rn_efr    = 1.0_wp          ! fraction of TKE surface value which penetrates in the ocean 
     81   REAL(wp) ::   rn_addhft = 0.0_wp          ! add offset   applied to HF tau 
     82   REAL(wp) ::   rn_sclhft = 1.0_wp          ! scale factor applied to HF tau 
     83   LOGICAL  ::   ln_lc     = .FALSE.         ! Langmuir cells (LC) as a source term of TKE or not 
     84   REAL(wp) ::   rn_lc     = 0.15_wp         ! coef to compute vertical velocity of Langmuir cells 
    8385 
    8486   REAL(wp) ::   ri_cri   ! critic Richardson number (deduced from rn_ediff and rn_ediss values) 
     
    177179      USE oce,   zd_lw  =>   ta   ! use ta as workspace 
    178180      !! 
    179       INTEGER  ::   ji, jj, jk                ! dummy loop arguments 
    180       REAL(wp) ::   zbbrau, zesh2             ! temporary scalars 
    181       REAL(wp) ::   zfact1, zfact2, zfact3    !    -         - 
    182       REAL(wp) ::   ztx2  , zty2  , zcof      !    -         - 
    183       REAL(wp) ::   zus   , zwlc  , zind      !    -         - 
    184       REAL(wp) ::   zzd_up, zzd_lw            !    -         - 
     181      INTEGER  ::   ji, jj, jk                      ! dummy loop arguments 
    185182      INTEGER  ::   ikbu, ikbv, ikbum1, ikbvm1      ! temporary scalar 
    186183      INTEGER  ::   ikbt, ikbumm1, ikbvmm1          ! temporary scalar 
    187       REAL(wp) ::   zebot                           ! temporary scalars 
     184      REAL(wp) ::   zrhoa  = 1.22                   ! Air density kg/m3 
     185      REAL(wp) ::   zcdrag = 1.5e-3                 ! drag coefficient 
     186      REAL(wp) ::   zbbrau, zesh2                   ! temporary scalars 
     187      REAL(wp) ::   zfact1, zfact2, zfact3          !    -         - 
     188      REAL(wp) ::   ztx2  , zty2  , zcof            !    -         - 
     189      REAL(wp) ::   ztau  , zdif                    !    -         - 
     190      REAL(wp) ::   zus   , zwlc  , zind            !    -         - 
     191      REAL(wp) ::   zzd_up, zzd_lw                  !    -         - 
     192      REAL(wp) ::   zebot                           !    -         - 
    188193      INTEGER , DIMENSION(jpi,jpj)     ::   imlc    ! 2D workspace 
    189194      REAL(wp), DIMENSION(jpi,jpj)     ::   zhlc    !  -      - 
     
    242247         END DO 
    243248         !                        !* finite Langmuir Circulation depth 
     249         zcof = 0.5 * 0.016 * 0.016 / ( zrhoa * zcdrag ) 
    244250         imlc(:,:) = mbathy(:,:)         ! Initialization to the number of w ocean point mbathy 
    245251         DO jk = jpkm1, 2, -1 
    246252            DO jj = 1, jpj               ! Last w-level at which zpelc>=0.5*us*us  
    247253               DO ji = 1, jpi            !      with us=0.016*wind(starting from jpk-1) 
    248                   zus  = 0.000128 * wndm(ji,jj) * wndm(ji,jj) 
     254                  zus  = zcof * taum(ji,jj) 
    249255                  IF( zpelc(ji,jj,jk) > zus )   imlc(ji,jj) = jk 
    250256               END DO 
     
    265271         hlc(:,:) = zhlc(:,:) * tmask(:,:,1)      ! c1d configuration: save finite Langmuir Circulation depth 
    266272# endif 
     273         zcof = 0.016 / SQRT( zrhoa * zcdrag ) 
    267274!CDIR NOVERRCHK 
    268275         DO jk = 2, jpkm1         !* TKE Langmuir circulation source term added to en 
     
    271278!CDIR NOVERRCHK 
    272279               DO ji = fs_2, fs_jpim1   ! vector opt. 
    273 !!gm replace here wndn by a formulation with the stress module 
    274                   zus  = 0.016 * wndm(ji,jj)                  ! Stokes drift 
     280                  zus  = zcof * SQRT( taum(ji,jj) )           ! Stokes drift 
    275281                  !                                           ! vertical velocity due to LC 
    276282                  zind = 0.5 - SIGN( 0.5, fsdepw(ji,jj,jk) - zhlc(ji,jj) ) 
     
    371377      !                            !  TKE due to surface and internal wave breaking 
    372378      !                            !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
    373       IF( nn_etau == 1 ) THEN           !*  penetration throughout the water column 
     379      IF( nn_etau == 1 ) THEN           !* penetration throughout the water column 
    374380         DO jk = 2, jpkm1 
    375381            DO jj = 2, jpjm1 
     
    386392               en(ji,jj,jk) = en(ji,jj,jk) + rn_efr * en(ji,jj,1) * EXP( -fsdepw(ji,jj,jk) / htau(ji,jj) )   & 
    387393                  &                                               * ( 1.e0 - fr_i(ji,jj) )  * tmask(ji,jj,jk) 
     394            END DO 
     395         END DO 
     396      ELSEIF( nn_etau == 3 ) THEN       !* penetration throughout the water column 
     397!CDIR NOVERRCHK 
     398         DO jk = 2, jpkm1 
     399!CDIR NOVERRCHK 
     400            DO jj = 2, jpjm1 
     401!CDIR NOVERRCHK 
     402               DO ji = fs_2, fs_jpim1   ! vector opt. 
     403                  ztx2 = utau(ji-1,jj  ) + utau(ji,jj) 
     404                  zty2 = vtau(ji  ,jj-1) + vtau(ji,jj) 
     405                  ztau = 0.5 * SQRT( ztx2 * ztx2 + zty2 * zty2 )      ! module of the mean stress  
     406                  zdif = taum(ji,jj) - ztau                           ! mean of the module - module of the mean  
     407                  zdif = rn_sclhft * MAX( 0.e0, zdif + rn_addhft )    ! apply some modifications... 
     408                  en(ji,jj,jk) = en(ji,jj,jk) + zbbrau * zdif * EXP( -fsdepw(ji,jj,jk) / htau(ji,jj) )   & 
     409                     &                                        * ( 1.e0 - fr_i(ji,jj) ) * tmask(ji,jj,jk) 
     410               END DO 
    388411            END DO 
    389412         END DO 
     
    642665      INTEGER ::   ji, jj, jk   ! dummy loop indices 
    643666      !! 
    644       NAMELIST/namzdf_tke/ rn_ediff, rn_ediss , rn_ebb, rn_emin,   & 
    645          &                 rn_emin0, rn_bshear, nn_mxl, ln_mxl0,   & 
    646          &                 rn_lmin , rn_lmin0 , nn_pdl, nn_etau,   & 
    647          &                 nn_htau , rn_efr   , ln_lc , rn_lc  
     667      NAMELIST/namzdf_tke/ rn_ediff, rn_ediss , rn_ebb   , rn_emin  ,   & 
     668         &                 rn_emin0, rn_bshear, nn_mxl   , ln_mxl0  ,   & 
     669         &                 rn_lmin , rn_lmin0 , nn_pdl   , nn_etau  ,   & 
     670         &                 nn_htau , rn_efr   , rn_addhft, rn_sclhft,   & 
     671         &                 ln_lc   , rn_lc  
    648672      !!---------------------------------------------------------------------- 
    649673 
     
    658682         WRITE(numout,*) '~~~~~~~~' 
    659683         WRITE(numout,*) '   Namelist namzdf_tke : set tke mixing parameters' 
    660          WRITE(numout,*) '      coef. to compute avt                        rn_ediff = ', rn_ediff 
    661          WRITE(numout,*) '      Kolmogoroff dissipation coef.               rn_ediss = ', rn_ediss 
    662          WRITE(numout,*) '      tke surface input coef.                     rn_ebb   = ', rn_ebb 
    663          WRITE(numout,*) '      minimum value of tke                        rn_emin  = ', rn_emin 
    664          WRITE(numout,*) '      surface minimum value of tke                rn_emin0 = ', rn_emin0 
    665          WRITE(numout,*) '      background shear (>0)                       rn_bshear= ', rn_bshear 
    666          WRITE(numout,*) '      mixing length type                          nn_mxl   = ', nn_mxl 
    667          WRITE(numout,*) '      prandl number flag                          nn_pdl   = ', nn_pdl 
    668          WRITE(numout,*) '      surface mixing length = F(stress) or not    ln_mxl0  = ', ln_mxl0 
    669          WRITE(numout,*) '      surface  mixing length minimum value        rn_lmin0 = ', rn_lmin0 
    670          WRITE(numout,*) '      interior mixing length minimum value        rn_lmin0 = ', rn_lmin0 
    671          WRITE(numout,*) '      test param. to add tke induced by wind      nn_etau  = ', nn_etau 
    672          WRITE(numout,*) '      flag for computation of exp. tke profile    nn_htau  = ', nn_htau 
    673          WRITE(numout,*) '      % of rn_emin0 which pene. the thermocline   rn_efr   = ', rn_efr 
    674          WRITE(numout,*) '      flag to take into acc.  Langmuir circ.      ln_lc    = ', ln_lc 
    675          WRITE(numout,*) '      coef to compute verticla velocity of LC     rn_lc    = ', rn_lc 
     684         WRITE(numout,*) '      coef. to compute avt                        rn_ediff  = ', rn_ediff 
     685         WRITE(numout,*) '      Kolmogoroff dissipation coef.               rn_ediss  = ', rn_ediss 
     686         WRITE(numout,*) '      tke surface input coef.                     rn_ebb    = ', rn_ebb 
     687         WRITE(numout,*) '      minimum value of tke                        rn_emin   = ', rn_emin 
     688         WRITE(numout,*) '      surface minimum value of tke                rn_emin0  = ', rn_emin0 
     689         WRITE(numout,*) '      background shear (>0)                       rn_bshear = ', rn_bshear 
     690         WRITE(numout,*) '      mixing length type                          nn_mxl    = ', nn_mxl 
     691         WRITE(numout,*) '      prandl number flag                          nn_pdl    = ', nn_pdl 
     692         WRITE(numout,*) '      surface mixing length = F(stress) or not    ln_mxl0   = ', ln_mxl0 
     693         WRITE(numout,*) '      surface  mixing length minimum value        rn_lmin0  = ', rn_lmin0 
     694         WRITE(numout,*) '      interior mixing length minimum value        rn_lmin0  = ', rn_lmin0 
     695         WRITE(numout,*) '      test param. to add tke induced by wind      nn_etau   = ', nn_etau 
     696         WRITE(numout,*) '      flag for computation of exp. tke profile    nn_htau   = ', nn_htau 
     697         WRITE(numout,*) '      fraction of en which pene. the thermocline  rn_efr    = ', rn_efr 
     698         WRITE(numout,*) '      add offset   applied to HF tau              rn_addhft = ', rn_addhft 
     699         WRITE(numout,*) '      scale factor applied to HF tau              rn_sclhft = ', rn_sclhft 
     700         WRITE(numout,*) '      flag to take into acc.  Langmuir circ.      ln_lc     = ', ln_lc 
     701         WRITE(numout,*) '      coef to compute verticla velocity of LC     rn_lc     = ', rn_lc 
    676702         WRITE(numout,*) 
    677703         WRITE(numout,*) '      critical Richardson nb with your parameters  ri_cri = ', ri_cri 
     
    679705 
    680706      !                               !* Check of some namelist values 
    681       IF( nn_mxl  < 0    .OR. nn_mxl  > 3   )   CALL ctl_stop( 'bad flag: nn_mxl is  0, 1 or 2 ' ) 
    682       IF( nn_pdl  < 0    .OR. nn_pdl  > 1   )   CALL ctl_stop( 'bad flag: nn_pdl is  0 or 1    ' ) 
    683       IF( nn_htau < 0    .OR. nn_htau > 1   )   CALL ctl_stop( 'bad flag: nn_htau is 0 or 1    ' ) 
    684       IF( rn_lc   < 0.15 .OR. rn_lc   > 0.2 )   CALL ctl_stop( 'bad value: rn_lc must be between 0.15 and 0.2 ' ) 
     707      IF( nn_mxl  < 0     .OR.  nn_mxl  > 3   )   CALL ctl_stop( 'bad flag: nn_mxl is  0, 1 or 2 ' ) 
     708      IF( nn_pdl  < 0     .OR.  nn_pdl  > 1   )   CALL ctl_stop( 'bad flag: nn_pdl is  0 or 1    ' ) 
     709      IF( nn_htau < 0     .OR.  nn_htau > 2   )   CALL ctl_stop( 'bad flag: nn_htau is 0, 1 or 2 ' ) 
     710      IF( rn_lc   < 0.15  .OR.  rn_lc   > 0.2 )   CALL ctl_stop( 'bad value: rn_lc must be between 0.15 and 0.2 ' ) 
     711      IF( rn_sclhft == 0. .AND. nn_etau == 3  )   CALL ctl_stop( 'force null HF tau to penetrate the thermocline...' ) 
     712      IF( .NOT. lhftau    .AND. nn_etau == 3  )   CALL ctl_stop( 'bad flag: nn_etau == 3 must be used with HF tau' ) 
    685713 
    686714      IF( nn_etau == 2  )   CALL zdf_mxl( nit000 )      ! Initialization of nmln  
     
    697725               END DO 
    698726            END DO 
     727         CASE( 2 )                                    ! constant depth penetration (here 30 meters) 
     728            htau(:,:) = 30.e0 
    699729         END SELECT 
    700730      ENDIF 
  • trunk/UTIL/chk_iomput.sh

    r1482 r1705  
    88inxml=0 
    99insrc=0 
     10xmlfile=${1:-"../CONFIG/ORCA2_LIM/EXP00/iodef.xml"} 
     11srcdir=${2:-"../NEMO/OPA_SRC ../NEMO/LIM_SRC_2"} 
    1012while [ $# -gt 0 ]   # Until you run out of parameters . . .  
    1113do     
     
    1315        --inxml) inxml=1 ;; 
    1416        --insrc) insrc=1 ;; 
     17   *) xmlfile=${1:-"../CONFIG/ORCA2_LIM/EXP00/iodef.xml"} 
     18      srcdir=${2:-"../NEMO/OPA_SRC ../NEMO/LIM_SRC_2"} 
     19      shift 
    1520    esac 
    1621    shift       # Check next set of parameters.  
    1722done  
    18 # 
    19 xmlfile=${1:-"../CONFIG/ORCA2_LIM/EXP00/iodef.xml"} 
    20 # 
    21 srcdir=${2:-"../NEMO/OPA_SRC ../NEMO/LIM_SRC_2"} 
    2223# 
    2324#------------------------------------------------ 
Note: See TracChangeset for help on using the changeset viewer.