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 4147 for branches/2013/dev_LOCEAN_2013/NEMOGCM/NEMO/OPA_SRC/DOM/dom_oce.F90 – NEMO

Ignore:
Timestamp:
2013-11-04T12:51:55+01:00 (10 years ago)
Author:
cetlod
Message:

merge in dev_LOCEAN_2013, the 1st development branch dev_r3853_CNRS9_Confsetting, from its starting point ( r3853 ) on the trunk: see ticket #1169

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_LOCEAN_2013/NEMOGCM/NEMO/OPA_SRC/DOM/dom_oce.F90

    r3851 r4147  
    2828   !! time & space domain namelist 
    2929   !! ---------------------------- 
    30    !                                              !!* Namelist namdom : time & space domain * 
    31    INTEGER , PUBLIC ::   nn_bathy     =    0       !: = 0/1 ,compute/read the bathymetry file 
    32    REAL(wp), PUBLIC ::   rn_hmin      =   -3.0_wp  !: minimum ocean depth (>0) or minimum number of ocean levels (<0) 
    33    REAL(wp), PUBLIC ::   rn_e3zps_min =    5.0_wp  !: miminum thickness for partial steps (meters) 
    34    REAL(wp), PUBLIC ::   rn_e3zps_rat =    0.1_wp  !: minimum thickness ration for partial steps 
    35    INTEGER , PUBLIC ::   nn_msh       =    0       !: = 1 create a mesh-mask file 
    36    INTEGER , PUBLIC ::   nn_acc       =    0       !: = 0/1 use of the acceleration of convergence technique 
    37    REAL(wp), PUBLIC ::   rn_atfp      =    0.1_wp  !: asselin time filter parameter 
    38    REAL(wp), PUBLIC ::   rn_rdt       = 3600._wp   !: time step for the dynamics (and tracer if nacc=0) 
    39    REAL(wp), PUBLIC ::   rn_rdtmin    = 3600._wp   !: minimum time step on tracers 
    40    REAL(wp), PUBLIC ::   rn_rdtmax    = 3600._wp   !: maximum time step on tracers 
    41    REAL(wp), PUBLIC ::   rn_rdth      =  800._wp   !: depth variation of tracer step 
    42    INTEGER , PUBLIC ::   nn_baro      =   64       !: number of barotropic time steps (key_dynspg_ts) 
    43    INTEGER , PUBLIC ::   nn_closea    =    0       !: =0 suppress closed sea/lake from the ORCA domain or not (=1) 
     30   !                                    !!* Namelist namdom : time & space domain * 
     31   INTEGER , PUBLIC ::   nn_bathy        !: = 0/1 ,compute/read the bathymetry file 
     32   REAL(wp), PUBLIC ::   rn_hmin         !: minimum ocean depth (>0) or minimum number of ocean levels (<0) 
     33   REAL(wp), PUBLIC ::   rn_e3zps_min    !: miminum thickness for partial steps (meters) 
     34   REAL(wp), PUBLIC ::   rn_e3zps_rat    !: minimum thickness ration for partial steps 
     35   INTEGER , PUBLIC ::   nn_msh          !: = 1 create a mesh-mask file 
     36   INTEGER , PUBLIC ::   nn_acc          !: = 0/1 use of the acceleration of convergence technique 
     37   REAL(wp), PUBLIC ::   rn_atfp         !: asselin time filter parameter 
     38   REAL(wp), PUBLIC ::   rn_rdt          !: time step for the dynamics (and tracer if nacc=0) 
     39   REAL(wp), PUBLIC ::   rn_rdtmin       !: minimum time step on tracers 
     40   REAL(wp), PUBLIC ::   rn_rdtmax       !: maximum time step on tracers 
     41   REAL(wp), PUBLIC ::   rn_rdth         !: depth variation of tracer step 
     42   INTEGER , PUBLIC ::   nn_baro         !: number of barotropic time steps (key_dynspg_ts) 
     43   INTEGER , PUBLIC ::   nn_closea       !: =0 suppress closed sea/lake from the ORCA domain or not (=1) 
     44   !! Horizontal grid parameters for domhgr 
     45   !! ===================================== 
     46   INTEGER       ::   jphgr_msh        !: type of horizontal mesh 
     47   !                                       !  = 0 curvilinear coordinate on the sphere read in coordinate.nc 
     48   !                                       !  = 1 geographical mesh on the sphere with regular grid-spacing 
     49   !                                       !  = 2 f-plane with regular grid-spacing 
     50   !                                       !  = 3 beta-plane with regular grid-spacing 
     51   !                                       !  = 4 Mercator grid with T/U point at the equator 
     52 
     53   REAL(wp)      ::   ppglam0              !: longitude of first raw and column T-point (jphgr_msh = 1) 
     54   REAL(wp)      ::   ppgphi0              !: latitude  of first raw and column T-point (jphgr_msh = 1) 
     55   !                                                        !  used for Coriolis & Beta parameters (jphgr_msh = 2 or 3) 
     56   REAL(wp)      ::   ppe1_deg             !: zonal      grid-spacing (degrees) 
     57   REAL(wp)      ::   ppe2_deg             !: meridional grid-spacing (degrees) 
     58   REAL(wp)      ::   ppe1_m               !: zonal      grid-spacing (degrees) 
     59   REAL(wp)      ::   ppe2_m               !: meridional grid-spacing (degrees) 
     60 
     61   !! Vertical grid parameter for domzgr 
     62   !! ================================== 
     63   REAL(wp)      ::   ppsur                !: ORCA r4, r2 and r05 coefficients 
     64   REAL(wp)      ::   ppa0                 !: (default coefficients) 
     65   REAL(wp)      ::   ppa1                 !: 
     66   REAL(wp)      ::   ppkth                !: 
     67   REAL(wp)      ::   ppacr                !: 
     68   ! 
     69   !  If both ppa0 ppa1 and ppsur are specified to 0, then 
     70   !  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 
     71   REAL(wp)      ::   ppdzmin              !: Minimum vertical spacing 
     72   REAL(wp)      ::   pphmax               !: Maximum depth 
     73   ! 
     74   LOGICAL       ::   ldbletanh            !: Use/do not use double tanf function for vertical coordinates 
     75   REAL(wp)      ::   ppa2                 !: Double tanh function parameters 
     76   REAL(wp)      ::   ppkth2               !: 
     77   REAL(wp)      ::   ppacr2               !: 
    4478 
    4579   !                                    !! old non-DOCTOR names still used in the model 
     
    6296 
    6397   !                                         !!* Namelist namcla : cross land advection 
    64    INTEGER, PUBLIC ::   nn_cla = 0            !: =1 cross land advection for exchanges through some straits (ORCA2) 
     98   INTEGER, PUBLIC ::   nn_cla               !: =1 cross land advection for exchanges through some straits (ORCA2) 
    6599 
    66100   !!---------------------------------------------------------------------- 
     
    72106   LOGICAL, PUBLIC ::   lzoom_s    =  .FALSE.   !: South zoom type flag 
    73107   LOGICAL, PUBLIC ::   lzoom_n    =  .FALSE.   !: North zoom type flag 
    74    LOGICAL, PUBLIC ::   lzoom_arct =  .FALSE.   !: ORCA    arctic zoom flag 
    75    LOGICAL, PUBLIC ::   lzoom_anta =  .FALSE.   !: ORCA antarctic zoom flag 
    76108 
    77109   !                                     !!! domain parameters linked to mpp 
     
    100132   INTEGER, PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   mig        !: local  ==> global domain i-index 
    101133   INTEGER, PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   mjg        !: local  ==> global domain j-index 
    102    INTEGER, PUBLIC,               DIMENSION(jpidta) ::   mi0, mi1   !: global ==> local  domain i-index    !!bug ==> other solution? 
     134   INTEGER, PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   mi0, mi1   !: global ==> local  domain i-index    !!bug ==> other solution? 
    103135   !                                                  ! (mi0=1 and mi1=0 if the global index is not in the local domain) 
    104    INTEGER, PUBLIC,               DIMENSION(jpjdta) ::   mj0, mj1   !: global ==> local  domain j-index     !!bug ==> other solution? 
     136   INTEGER, PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   mj0, mj1   !: global ==> local  domain j-index     !!bug ==> other solution? 
    105137   !                                                  ! (mi0=1 and mi1=0 if the global index is not in the local domain) 
    106138   INTEGER, PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   nimppt, njmppt   !: i-, j-indexes for each processor 
     
    127159   !! vertical coordinate and scale factors 
    128160   !! --------------------------------------------------------------------- 
    129    !                                           !!* Namelist namzgr : vertical coordinate * 
    130    LOGICAL, PUBLIC ::   ln_zco     =  .TRUE.    !: z-coordinate - full step 
    131    LOGICAL, PUBLIC ::   ln_zps     =  .FALSE.   !: z-coordinate - partial step 
    132    LOGICAL, PUBLIC ::   ln_sco     =  .FALSE.   !: s-coordinate or hybrid z-s coordinate 
     161   !                                 !!* Namelist namzgr : vertical coordinate * 
     162   LOGICAL, PUBLIC ::   ln_zco        !: z-coordinate - full step 
     163   LOGICAL, PUBLIC ::   ln_zps        !: z-coordinate - partial step 
     164   LOGICAL, PUBLIC ::   ln_sco        !: s-coordinate or hybrid z-s coordinate 
    133165 
    134166   !! All coordinates 
     
    195227   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: tmask, umask, vmask, fmask   !: land/ocean mask at T-, U-, V- and F-pts 
    196228 
    197    REAL(wp), PUBLIC, DIMENSION(jpiglo) ::   tpol, fpol          !: north fold mask (jperio= 3 or 4) 
     229   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   tpol, fpol          !: north fold mask (jperio= 3 or 4) 
    198230 
    199231#if defined key_noslip_accurate 
     
    272304      ALLOCATE( nimppt(jpnij) , ibonit(jpnij) , nlcit(jpnij) , nlcjt(jpnij) ,     & 
    273305         &      njmppt(jpnij) , ibonjt(jpnij) , nldit(jpnij) , nldjt(jpnij) ,     & 
    274          &                                      nleit(jpnij) , nlejt(jpnij) , STAT=ierr(2) ) 
     306         &                                      nleit(jpnij) , nlejt(jpnij) ,     & 
     307         &      mi0(jpidta)   , mi1 (jpidta),  mj0(jpjdta)   , mj1 (jpjdta),      & 
     308         &      tpol(jpiglo)  , fpol(jpiglo)                               , STAT=ierr(2) ) 
    275309         ! 
    276310      ALLOCATE( glamt(jpi,jpj) , gphit(jpi,jpj) , e1t(jpi,jpj) , e2t(jpi,jpj) ,                      &  
Note: See TracChangeset for help on using the changeset viewer.