Ignore:
Timestamp:
2013-09-12T11:08:28+02:00 (8 years ago)
Author:
clevy
Message:

Configuration setting/correct errors, see ticket:#1074

Location:
branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/SHARED/1_namelist_ref

    r4018 r4025  
    6464                                       !  = 5 North fold F-point pivot 
    6565                                       !  = 6 cyclic East-West AND North fold F-point pivot 
     66/ 
     67!!====================================================================== 
     68!!                      ***  Domain namelists  *** 
     69!!====================================================================== 
     70!!   namzgr       vertical coordinate 
     71!!   namzgr_sco   s-coordinate or hybrid z-s-coordinate 
     72!!   namdom       space and time domain (bathymetry, mesh, timestep) 
     73!!   namtsd       data: temperature & salinity 
     74!!====================================================================== 
     75! 
     76!----------------------------------------------------------------------- 
     77&namzgr        !   vertical coordinate 
     78!----------------------------------------------------------------------- 
     79   ln_zco      = .false.   !  z-coordinate - full    steps   (T/F)      ("key_zco" may also be defined) 
     80   ln_zps      = .true.    !  z-coordinate - partial steps   (T/F) 
     81   ln_sco      = .false.   !  s- or hybrid z-s-coordinate    (T/F) 
     82/ 
     83!----------------------------------------------------------------------- 
     84&namzgr_sco    !   s-coordinate or hybrid z-s-coordinate 
     85!----------------------------------------------------------------------- 
     86   ln_s_sh94   = .true.    !  Song & Haidvogel 1994 hybrid S-sigma   (T)| 
     87   ln_s_sf12   = .false.   !  Siddorn & Furner 2012 hybrid S-z-sigma (T)| if both are false the NEMO tanh stretching is applied 
     88   ln_sigcrit  = .false.   !  use sigma coordinates below critical depth (T) or Z coordinates (F) for Siddorn & Furner stretch 
     89                           !  stretching coefficients for all functions 
     90   rn_sbot_min =   10.0    !  minimum depth of s-bottom surface (>0) (m) 
     91   rn_sbot_max = 7000.0    !  maximum depth of s-bottom surface (= ocean depth) (>0) (m) 
     92   rn_hc       =  150.0    !  critical depth for transition to stretched coordinates 
     93                        !!!!!!!  Envelop bathymetry 
     94   rn_rmax     =    0.3    !  maximum cut-off r-value allowed (0<r_max<1) 
     95                        !!!!!!!  SH94 stretching coefficients  (ln_s_sh94 = .true.) 
     96   rn_theta    =    6.0    !  surface control parameter (0<=theta<=20) 
     97   rn_bb       =    0.8    !  stretching with SH94 s-sigma    
     98                        !!!!!!!  SF12 stretching coefficient  (ln_s_sf12 = .true.) 
     99   rn_alpha    =    4.4    !  stretching with SF12 s-sigma 
     100   rn_efold    =    0.0    !  efold length scale for transition to stretched coord 
     101   rn_zs       =    1.0    !  depth of surface grid box 
     102                           !  bottom cell depth (Zb) is a linear function of water depth Zb = H*a + b 
     103   rn_zb_a     =    0.024  !  bathymetry scaling factor for calculating Zb 
     104   rn_zb_b     =   -0.2    !  offset for calculating Zb 
     105                        !!!!!!!! Other stretching (not SH94 or SF12) [also uses rn_theta above] 
     106   rn_thetb    =    1.0    !  bottom control parameter  (0<=thetb<= 1)  
     107/ 
     108!----------------------------------------------------------------------- 
     109&namdom        !   space and time domain (bathymetry, mesh, timestep) 
     110!----------------------------------------------------------------------- 
     111   nn_bathy    =    1      !  compute (=0) or read (=1) the bathymetry file 
     112   nn_closea    =   0      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
     113   nn_msh      =    0      !  create (=1) a mesh file or not (=0) 
     114   rn_hmin     =   -3.     !  min depth of the ocean (>0) or min number of ocean level (<0) 
     115   rn_e3zps_min=   20.     !  partial step thickness is set larger than the minimum of 
     116   rn_e3zps_rat=    0.1    !  rn_e3zps_min and rn_e3zps_rat*e3t, with 0<rn_e3zps_rat<1 
     117                           ! 
     118   rn_rdt      = 5760.     !  time step for the dynamics (and tracer if nn_acc=0) 
     119   nn_baro     =   64      !  number of barotropic time step            ("key_dynspg_ts") 
     120   rn_atfp     =    0.1    !  asselin time filter parameter 
     121   nn_acc      =    0      !  acceleration of convergence : =1      used, rdt < rdttra(k) 
     122                                 !                          =0, not used, rdt = rdttra 
     123   rn_rdtmin   = 28800.          !  minimum time step on tracers (used if nn_acc=1) 
     124   rn_rdtmax   = 28800.          !  maximum time step on tracers (used if nn_acc=1) 
     125   rn_rdth     =  800.           !  depth variation of tracer time step  (used if nn_acc=1) 
    66126   jphgr_msh   =       0               !  type of horizontal mesh 
    67127                                       !  = 0 curvilinear coordinate on the sphere read in coordinate.nc 
     
    87147   ppkth2      =       48.029893720000 ! 
    88148   ppacr2      =       13.000000000000 ! 
    89 / 
    90 !!====================================================================== 
    91 !!                      ***  Domain namelists  *** 
    92 !!====================================================================== 
    93 !!   namzgr       vertical coordinate 
    94 !!   namzgr_sco   s-coordinate or hybrid z-s-coordinate 
    95 !!   namdom       space and time domain (bathymetry, mesh, timestep) 
    96 !!   namtsd       data: temperature & salinity 
    97 !!====================================================================== 
    98 ! 
    99 !----------------------------------------------------------------------- 
    100 &namzgr        !   vertical coordinate 
    101 !----------------------------------------------------------------------- 
    102    ln_zco      = .false.   !  z-coordinate - full    steps   (T/F)      ("key_zco" may also be defined) 
    103    ln_zps      = .true.    !  z-coordinate - partial steps   (T/F) 
    104    ln_sco      = .false.   !  s- or hybrid z-s-coordinate    (T/F) 
    105 / 
    106 !----------------------------------------------------------------------- 
    107 &namzgr_sco    !   s-coordinate or hybrid z-s-coordinate 
    108 !----------------------------------------------------------------------- 
    109    ln_s_sh94   = .true.    !  Song & Haidvogel 1994 hybrid S-sigma   (T)| 
    110    ln_s_sf12   = .false.   !  Siddorn & Furner 2012 hybrid S-z-sigma (T)| if both are false the NEMO tanh stretching is applied 
    111    ln_sigcrit  = .false.   !  use sigma coordinates below critical depth (T) or Z coordinates (F) for Siddorn & Furner stretch 
    112                            !  stretching coefficients for all functions 
    113    rn_sbot_min =   10.0    !  minimum depth of s-bottom surface (>0) (m) 
    114    rn_sbot_max = 7000.0    !  maximum depth of s-bottom surface (= ocean depth) (>0) (m) 
    115    rn_hc       =  150.0    !  critical depth for transition to stretched coordinates 
    116                         !!!!!!!  Envelop bathymetry 
    117    rn_rmax     =    0.3    !  maximum cut-off r-value allowed (0<r_max<1) 
    118                         !!!!!!!  SH94 stretching coefficients  (ln_s_sh94 = .true.) 
    119    rn_theta    =    6.0    !  surface control parameter (0<=theta<=20) 
    120    rn_bb       =    0.8    !  stretching with SH94 s-sigma    
    121                         !!!!!!!  SF12 stretching coefficient  (ln_s_sf12 = .true.) 
    122    rn_alpha    =    4.4    !  stretching with SF12 s-sigma 
    123    rn_efold    =    0.0    !  efold length scale for transition to stretched coord 
    124    rn_zs       =    1.0    !  depth of surface grid box 
    125                            !  bottom cell depth (Zb) is a linear function of water depth Zb = H*a + b 
    126    rn_zb_a     =    0.024  !  bathymetry scaling factor for calculating Zb 
    127    rn_zb_b     =   -0.2    !  offset for calculating Zb 
    128                         !!!!!!!! Other stretching (not SH94 or SF12) [also uses rn_theta above] 
    129    rn_thetb    =    1.0    !  bottom control parameter  (0<=thetb<= 1)  
    130 / 
    131 !----------------------------------------------------------------------- 
    132 &namdom        !   space and time domain (bathymetry, mesh, timestep) 
    133 !----------------------------------------------------------------------- 
    134    nn_bathy    =    1      !  compute (=0) or read (=1) the bathymetry file 
    135    nn_closea    =   0      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
    136    nn_msh      =    0      !  create (=1) a mesh file or not (=0) 
    137    rn_hmin     =   -3.     !  min depth of the ocean (>0) or min number of ocean level (<0) 
    138    rn_e3zps_min=   20.     !  partial step thickness is set larger than the minimum of 
    139    rn_e3zps_rat=    0.1    !  rn_e3zps_min and rn_e3zps_rat*e3t, with 0<rn_e3zps_rat<1 
    140                            ! 
    141    rn_rdt      = 5760.     !  time step for the dynamics (and tracer if nn_acc=0) 
    142    nn_baro     =   64      !  number of barotropic time step            ("key_dynspg_ts") 
    143    rn_atfp     =    0.1    !  asselin time filter parameter 
    144    nn_acc      =    0      !  acceleration of convergence : =1      used, rdt < rdttra(k) 
    145                                  !                          =0, not used, rdt = rdttra 
    146    rn_rdtmin   = 28800.          !  minimum time step on tracers (used if nn_acc=1) 
    147    rn_rdtmax   = 28800.          !  maximum time step on tracers (used if nn_acc=1) 
    148    rn_rdth     =  800.           !  depth variation of tracer time step  (used if nn_acc=1) 
    149149/ 
    150150!----------------------------------------------------------------------- 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/cfg.txt

    r3993 r4025  
    33GYRE_BFM OPA_SRC TOP_SRC 
    44GYRE OPA_SRC 
    5 ORCA2_LIM OPA_SRC LIM_SRC_2 NST_SRC 
    65AMM12 OPA_SRC 
    76ORCA2_SAS_LIM OPA_SRC SAS_SRC LIM_SRC_2 NST_SRC 
     
    109GYRE_PISCES OPA_SRC TOP_SRC 
    1110ORCA2_LIM_CFC_C14b OPA_SRC LIM_SRC_2 NST_SRC TOP_SRC 
     11ORCA2_LIM OPA_SRC LIM_SRC_2 NST_SRC 
    1212AGRIF OPA_SRC LIM_SRC_2 NST_SRC 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/NST_SRC/agrif_user.F90

    r4018 r4025  
    3131      jpi     = ( jpiglo-2*jpreci + (jpni-1+0) ) / jpni + 2*jpreci 
    3232      jpj     = ( jpjglo-2*jprecj + (jpnj-1+0) ) / jpnj + 2*jprecj 
     33      jpk     = jpkdta  
     34      jpim1   = jpi-1  
     35      jpjm1   = jpj-1  
     36      jpkm1   = jpk-1                                          
     37      jpij    = jpi*jpj  
    3338      jpidta  = jpiglo 
    3439      jpjdta  = jpjglo 
     
    6772      ENDIF 
    6873   ENDIF 
     74PRINT *, " cp_cfg = ",cp_cfg," jp_cfg= ", jp_cfg 
    6975   ! Specific fine grid Initializations 
    7076   ! no tracer damping on fine grids 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/nemogcm.F90

    r4020 r4025  
    504504         WRITE(numout,*) '      left bottom j index of the zoom (in data domain) jpizoom = ', jpjzoom 
    505505         WRITE(numout,*) '      lateral cond. type (between 0 and 6) jperio = ', jperio    
    506          WRITE(numout,*) '      type of horizontal mesh jphgr_msh           = ', jphgr_msh 
    507          WRITE(numout,*) '      longitude of first raw and column T-point ppglam0 = ', ppglam0 
    508          WRITE(numout,*) '      latitude  of first raw and column T-point ppgphi0 = ', ppgphi0 
    509          WRITE(numout,*) '      zonal      grid-spacing (degrees) ppe1_deg        = ', ppe1_deg 
    510          WRITE(numout,*) '      meridional grid-spacing (degrees) ppe2_deg        = ', ppe2_deg 
    511          WRITE(numout,*) '      zonal      grid-spacing (degrees) ppe1_m          = ', ppe1_m 
    512          WRITE(numout,*) '      meridional grid-spacing (degrees) ppe2_m          = ', ppe2_m 
    513          WRITE(numout,*) '      ORCA r4, r2 and r05 coefficients  ppsur           = ', ppsur 
    514          WRITE(numout,*) '                                        ppa0            = ', ppa0 
    515          WRITE(numout,*) '                                        ppa1            = ', ppa1 
    516          WRITE(numout,*) '                                        ppkth           = ', ppkth 
    517          WRITE(numout,*) '                                        ppacr           = ', ppacr 
    518          WRITE(numout,*) '      Minimum vertical spacing ppdzmin                  = ', ppdzmin 
    519          WRITE(numout,*) '      Maximum depth pphmax                              = ', pphmax 
    520          WRITE(numout,*) '      Use double tanf function for vertical coordinates ldbletanh = ', ldbletanh 
    521          WRITE(numout,*) '      Double tanh function parameters ppa2              = ', ppa2 
    522          WRITE(numout,*) '                                      ppkth2            = ', ppkth2 
    523          WRITE(numout,*) '                                      ppacr2            = ', ppacr2 
    524506      ENDIF 
    525507      !                             ! Parameter control 
Note: See TracChangeset for help on using the changeset viewer.