Changeset 11628


Ignore:
Timestamp:
2019-10-01T18:12:05+02:00 (13 months ago)
Author:
mathiot
Message:

ENHANCE-03_domcfg: add cpp_ENHANCE-03_domcfg.fcm and remove make_closea_masks.py (as now managed in the F90 code and remove compilation file (ticket #2143)

Location:
NEMO/branches/2019/ENHANCE-03_domcfg
Files:
1 added
5 deleted
3 edited
1 moved

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/ENHANCE-03_domcfg/README.rst

    r11599 r11628  
    1919NOTA: it can be run in multiproc mode, but in output there will be domain_cfg_00xx.nc files 
    2020 
    21  === Closed seas (closea module) === 
    22  
    23 If you want to define closed seas in the bathymetry either to suppress them at runtime or  
    24 redistribute freshwater fluxes, then you need to run make_closea_masks.py after you have 
    25 created the basic domain_cfg file. This utility will add "closea_mask*" fields to the 
    26 domain_cfg file to define the closed seas in the configuration. (If you have closed seas 
    27 but don't want to treat them in a special way then you can ignore this step).  
    28  
    2921================================ 
    3022= HOW TO USE 
     
    3325IMPORTANT : keep the namelist_ref committed inchanged. !!! 
    3426 
    35 NEW OPTION ln_e3_dep in the namelist_ref: 
     27NEW OPTION: 
     28* ln_e3_dep 
    3629till nemo_v3.6_stable e3 were done like an analytical derivative of depth function 
    3730now  e3=dk[depth] in discret sens 
     
    4437   !                       !      only there for backward compatibility test with v3.6 
    4538   ! 
     39 
     40* ln_domclo 
     41The definition of closed sea are now integrated into ./make_domain_cfg.exe. The need on an external python script and hard coded indices has been removed. 
     42ln_domclo=T will generate the masks needed to mask or correct  
     43the fwf unbalance over closed sea. If non defined closed sea are detected, a mask containing this will also be generated. 
     44The detection of the lake is done using a lat/lon seed and a flood filling algorithm (the default namelist should work for any resolution). 
     45The definition of the target area is done by defining a lat/lon position center, then all wet points at a defined distance of it will be selected. 
     46Options are available to select only coastal points if needed. 
    4647 
    47482) copy in DOMAINcfg directory same input files (of related configuration) required in v3.6_stable. 
  • NEMO/branches/2019/ENHANCE-03_domcfg/namelist_cfg

    r10727 r11628  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! NEMO/OPA  Configuration namelist : used to overwrite defaults values defined in SHARED/namelist_ref 
     2!! NEMO/OCE :   Reference namelist_ref                                !! 
    33!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    4 ! 
     4!! NEMO/OCE  :  1 - Domain & run manager (namrun, namcfg, namdom, namzgr, namzgr_sco ) 
     5!!              2 - diagnostics      (namnc4) 
     6!!              3 - miscellaneous    (nammpp, namctl) 
     7!! 
     8!! namelist skeleton : egrep -E '(^/ *$|^! *$|^ *$|&nam.*|!---.*|!! .*|!!==.*|!!>>>.*)' namelist_ref > namelist_skl 
     9!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    510!----------------------------------------------------------------------- 
    611&namrun        !   parameters of the run 
    712!----------------------------------------------------------------------- 
    8    nn_no       =       0   !  job number (no more used...) 
    9    cn_exp      =  "domaincfg"  !  experience name 
    10    nn_it000    =       1   !  first time step 
    11    nn_itend    =      75   !  last  time step (std 5475) 
     13/ 
     14!----------------------------------------------------------------------- 
     15&namdom        !   space and time domain (bathymetry, mesh, timestep) 
     16!----------------------------------------------------------------------- 
    1217/ 
    1318!----------------------------------------------------------------------- 
    1419&namcfg        !   parameters of the configuration 
    1520!----------------------------------------------------------------------- 
    16    ! 
    17    ln_e3_dep   = .true.    ! =T : e3=dk[depth] in discret sens.  
    18    !                       !      ===>>> will become the only possibility in v4.0 
    19    !                       ! =F : e3 analytical derivative of depth function 
    20    !                       !      only there for backward compatibility test with v3.6 
    21    !                       !       
    22    cp_cfg      =  "orca"   !  name of the configuration 
    23    jp_cfg      =       2   !  resolution of the configuration 
    24    jpidta      =     182   !  1st lateral dimension ( >= jpi ) 
    25    jpjdta      =     149   !  2nd    "         "    ( >= jpj ) 
    26    jpkdta      =      31   !  number of levels      ( >= jpk ) 
    27    jpiglo      =     182   !  1st dimension of global domain --> i =jpidta 
    28    jpjglo      =     149   !  2nd    -                  -    --> j  =jpjdta 
    29    jpizoom     =       1   !  left bottom (i,j) indices of the zoom 
    30    jpjzoom     =       1   !  in data domain indices 
    31    jperio      =       4   !  lateral cond. type (between 0 and 6) 
    3221/ 
    3322!----------------------------------------------------------------------- 
    34 &namzgr        !   vertical coordinate 
    35 !----------------------------------------------------------------------- 
    36    ln_zps      = .true.    !  z-coordinate - partial steps 
    37    ln_linssh   = .true.    !  linear free surface 
    38 / 
    39 &namzgr_sco 
    40 / 
    41 &namlbc 
    42 / 
    43 &namagrif 
    44 / 
    45 &nambdy 
    46 / 
    47 &nam_vvl 
    48 / 
    49  
    50 !----------------------------------------------------------------------- 
    51 &namdom        !   space and time domain (bathymetry, mesh, timestep) 
    52 !----------------------------------------------------------------------- 
    53    jphgr_msh   =       0               !  type of horizontal mesh 
    54    ppglam0     =  999999.0             !  longitude of first raw and column T-point (jphgr_msh = 1) 
    55    ppgphi0     =  999999.0             ! latitude  of first raw and column T-point (jphgr_msh = 1) 
    56    ppe1_deg    =  999999.0             !  zonal      grid-spacing (degrees) 
    57    ppe2_deg    =  999999.0             !  meridional grid-spacing (degrees) 
    58    ppe1_m      =  999999.0             !  zonal      grid-spacing (degrees) 
    59    ppe2_m      =  999999.0             !  meridional grid-spacing (degrees) 
    60    ppsur       =   -4762.96143546300   !  ORCA r4, r2 and r05 coefficients 
    61    ppa0        =     255.58049070440   ! (default coefficients) 
    62    ppa1        =     245.58132232490   ! 
    63    ppkth       =      21.43336197938   ! 
    64    ppacr       =       3.0             ! 
    65    ppdzmin     =  999999.              !  Minimum vertical spacing 
    66    pphmax      =  999999.              !  Maximum depth 
    67    ldbletanh   =  .FALSE.              !  Use/do not use double tanf function for vertical coordinates 
    68    ppa2        =  999999.              !  Double tanh function parameters 
    69    ppkth2      =  999999.              ! 
    70    ppacr2      =  999999.              ! 
    71 / 
    72 !----------------------------------------------------------------------- 
    73 &nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
     23&namzgr        !   vertical coordinate                                  (default: NO selection) 
    7424!----------------------------------------------------------------------- 
    7525/ 
    7626!----------------------------------------------------------------------- 
    77 &namctl        !   Control prints & Benchmark 
     27&namzgr_sco    !   s-coordinate or hybrid z-s-coordinate                (default F) 
     28!----------------------------------------------------------------------- 
     29                        !!!!!!!  Envelop bathymetry 
     30                        !!!!!!!  SH94 stretching coefficients  (ln_s_sh94 = .true.) 
     31                        !!!!!!!  SF12 stretching coefficient  (ln_s_sf12 = .true.) 
     32                        !!!!!!!! Other stretching (not SH94 or SF12) [also uses rn_theta above] 
     33/ 
     34!----------------------------------------------------------------------- 
     35&namclo ! (closed sea : need ln_domclo = .true. in namcfg) 
    7836!----------------------------------------------------------------------- 
    7937/ 
     38!----------------------------------------------------------------------- 
     39&namlbc        !   lateral momentum boundary condition                  (default: NO selection) 
     40!----------------------------------------------------------------------- 
     41/ 
     42!----------------------------------------------------------------------- 
     43&namagrif      !  AGRIF zoom                                            ("key_agrif") 
     44!----------------------------------------------------------------------- 
     45/ 
     46!----------------------------------------------------------------------- 
     47&nambdy        !  unstructured open boundaries                          (default: OFF) 
     48!----------------------------------------------------------------------- 
     49/ 
     50!----------------------------------------------------------------------- 
     51&namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4") 
     52!----------------------------------------------------------------------- 
     53/ 
     54!----------------------------------------------------------------------- 
     55&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi") 
     56!----------------------------------------------------------------------- 
     57/ 
     58!----------------------------------------------------------------------- 
     59&namctl        !   Control prints                                       (default: OFF) 
     60!----------------------------------------------------------------------- 
     61/ 
  • NEMO/branches/2019/ENHANCE-03_domcfg/namelist_ref

    r11201 r11628  
    33!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    44!! NEMO/OCE  :  1 - Domain & run manager (namrun, namcfg, namdom, namzgr, namzgr_sco ) 
    5 !!              8 - diagnostics      (namnc4) 
    6 !!             10 - miscellaneous    (nammpp, namctl) 
     5!!              2 - diagnostics      (namnc4) 
     6!!              3 - miscellaneous    (nammpp, namctl) 
    77!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    88!----------------------------------------------------------------------- 
     
    139139&namclo ! (closed sea : need ln_domclo = .true. in namcfg) 
    140140!----------------------------------------------------------------------- 
    141    rn_lon_opnsea = -2.0     ! longitude seed of open ocean 
    142    rn_lat_opnsea = -2.0     ! latitude  seed of open ocean 
     141   rn_lon_opnsea = -0.0    ! longitude seed of open ocean 
     142   rn_lat_opnsea = -0.0    ! latitude  seed of open ocean 
    143143   nn_closea = 8           ! number of closed seas ( = 0; only the open_sea mask will be computed) 
    144    !                name   ! lon_src ! lat_src ! lon_trg ! lat_trg ! river mouth area   ! net evap/precip correction scheme ! radius tgt   ! id trg 
    145    !                       ! (degree)! (degree)! (degree)! (degree)! local/coast/global ! (glo/rnf/emp)                     !     (m)      ! 
     144   ! 
     145   !                name   ! lon_src ! lat_src ! lon_trg ! lat_trg ! river mouth area   ! correction scheme ! radius trg   ! id trg 
     146   !                       ! (degree)! (degree)! (degree)! (degree)! local/coast/global ! (glo/rnf/emp)     !     (m)      ! 
    146147   ! North American lakes 
    147    sn_lake(1) = 'superior' ,  -86.57 ,  47.30  , -66.49  , 50.45   , 'local'            , 'rnf'                             ,   550000.0 , 2     
    148    sn_lake(2) = 'michigan' ,  -87.06 ,  42.74  , -66.49  , 50.45   , 'local'            , 'rnf'                             ,   550000.0 , 2     
    149    sn_lake(3) = 'huron'    ,  -82.51 ,  44.74  , -66.49  , 50.45   , 'local'            , 'rnf'                             ,   550000.0 , 2     
    150    sn_lake(4) = 'erie'     ,  -81.13 ,  42.25  , -66.49  , 50.45   , 'local'            , 'rnf'                             ,   550000.0 , 2     
    151    sn_lake(5) = 'ontario'  ,  -77.72 ,  43.62  , -66.49  , 50.45   , 'local'            , 'rnf'                             ,   550000.0 , 2     
     148   sn_lake(1) = 'superior' ,  -86.57 ,  47.30  , -66.49  , 50.45   , 'local'            , 'rnf'             ,   550000.0  , 2     
     149   sn_lake(2) = 'michigan' ,  -87.06 ,  42.74  , -66.49  , 50.45   , 'local'            , 'rnf'             ,   550000.0  , 2     
     150   sn_lake(3) = 'huron'    ,  -82.51 ,  44.74  , -66.49  , 50.45   , 'local'            , 'rnf'             ,   550000.0  , 2     
     151   sn_lake(4) = 'erie'     ,  -81.13 ,  42.25  , -66.49  , 50.45   , 'local'            , 'rnf'             ,   550000.0  , 2     
     152   sn_lake(5) = 'ontario'  ,  -77.72 ,  43.62  , -66.49  , 50.45   , 'local'            , 'rnf'             ,   550000.0  , 2     
    152153   ! African Lake 
    153    sn_lake(6) = 'victoria' ,   32.93 ,  -1.08  ,  30.44  , 31.37   , 'coast'            , 'emp'                             ,   100000.0 , 3     
     154   sn_lake(6) = 'victoria' ,   32.93 ,  -1.08  ,  30.44  , 31.37   , 'coast'            , 'emp'             ,   100000.0  , 3     
    154155   ! Asian Lakes 
    155    sn_lake(7) = 'caspian'  ,   50.0  ,  44.0   ,   0.0   ,  0.0    , 'global'           , 'glo'                             ,        0.0 , 1      
    156    sn_lake(8) = 'aral'     ,   60.0  ,  45.0   ,   0.0   ,  0.0    , 'global'           , 'glo'                             ,        0.0 , 1     
     156   sn_lake(7) = 'caspian'  ,   50.0  ,  44.0   ,   0.0   ,  0.0    , 'global'           , 'glo'             ,        0.0  , 1      
     157   sn_lake(8) = 'aral'     ,   60.0  ,  45.0   ,   0.0   ,  0.0    , 'global'           , 'glo'             ,        0.0  , 1     
    157158/ 
    158159!----------------------------------------------------------------------- 
  • NEMO/branches/2019/ENHANCE-03_domcfg/src/domclo.F90

    r11604 r11628  
    238238 
    239239            !! define estuary 
    240             !! deal with local/coast cases 
     240            !! deal with global/local/coastal cases 
    241241            SELECT CASE (cloc) 
    242242            CASE ('global') 
Note: See TracChangeset for help on using the changeset viewer.