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 778 for branches/dev_001_GM/NEMO/OPA_SRC/par_eel_r2.F90 – NEMO

Ignore:
Timestamp:
2007-12-22T11:30:56+01:00 (16 years ago)
Author:
gm
Message:

dev_001_GM - par_oce evolution: suppress the par_...h90 and create new module for all ocean configuration -compilation OK

File:
1 moved

Legend:

Unmodified
Added
Removed
  • branches/dev_001_GM/NEMO/OPA_SRC/par_eel_r2.F90

    r772 r778  
     1MODULE par_eel_r2 
     2   !!====================================================================== 
     3   !!                        ***  par_eel_r2  *** 
     4   !! Ocean Domain parameter : 2 km resolution cyclic Channel (EEL_R2 configuration) 
     5   !!====================================================================== 
     6   !! History :   9.0  !  2007-11  (G. Madec)  Original code from par_EEL_R2.h90 
     7   !!---------------------------------------------------------------------- 
     8   !! NEMO/OPA 2.4 , LOCEAN-IPSL (2007)  
     9   !! $Header: $  
     10   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 
     11   !!---------------------------------------------------------------------- 
     12#if defined key_eel_r2 
    113   !!--------------------------------------------------------------------- 
    2    !!                     ***  par_EEL_R2.h90  *** 
    3    !!   Ocean Domain : 2 km resolution Channel (EEL_R2 configuration) 
     14   !!   'key_eel_r2'   :                            EEL configuration : R2 
    415   !!--------------------------------------------------------------------- 
    5    !!---------------------------------------------------------------------- 
    6    !!  OPA 9.0 , LOCEAN-IPSL (2005)  
    7    !! $Header$  
    8    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
    9    !!---------------------------------------------------------------------- 
    10    CHARACTER (len=16)      & 
    11 #if !defined key_agrif 
    12       , PARAMETER  & 
     16   USE par_kind          ! kind parameters 
     17 
     18   IMPLICIT NONE 
     19   PUBLIC 
     20 
     21# if defined key_agrif 
     22   !------!!-------------------------------------------------------------- 
     23   !      !!   'key_agrif'   :                                       AGRIF  
     24   !------!!-------------------------------------------------------------- 
     25 
     26   ! Configuration, and domain size are not assigned to a parameter value 
     27   !    as they are different in mother and children models 
     28   INTEGER ::   cp_cfg, jp_cfg                          !: configuration id       (AGRIF case) 
     29   INTEGER ::   jpidta, jpjdta, jpkdta                  !: input data size        (AGRIF case) 
     30   INTEGER ::   jpiglo, jpjglo, jpk, jpizoom, jpjzoom   !: domain size and zoom   (AGRIF case) 
     31   INTEGER ::   jperio, jpisl, jpnisl                   !: domain characteristics (AGRIF case) 
     32 
     33# else 
     34   !------!!-------------------------------------------------------------- 
     35   !      !!   Default       :                  no AGRIF: model parameters 
     36   !------!!-------------------------------------------------------------- 
     37 
     38   ! Configuration name and resolution      !!! *  GYRE R2  * 
     39   CHARACTER (len=16) ::   cp_cfg = "eel"    !: name of the configuration  
     40   INTEGER            ::   jp_cfg = 2        !: resolution of the configuration [km] 
     41 
     42   ! input data size                        !!! *  size of all input files  * 
     43   INTEGER, PARAMETER ::   jpidta  =  83                !: 1st lateral dimension ( >= jpiglo ) 
     44   INTEGER, PARAMETER ::   jpjdta  = 242                !: 2nd    "       "      ( >= jpjglo ) 
     45   INTEGER, PARAMETER ::   jpkdta  =  30                !: number of levels      ( >= jpk    )  
     46    
     47 
     48#  if defined key_cfg_1d 
     49   !------------!!-------------------------------------------------------- 
     50   !            !!   'key_cfg_1d' :            1D configuration of GYRE_R2 
     51   !------------!!-------------------------------------------------------- 
     52   ! global domain size                                 !!! *  1D configuration  * 
     53   INTEGER, PARAMETER ::   jpiglo  = 3                   !: 1st dimension of global domain --> i 
     54   INTEGER, PARAMETER ::   jpjglo  = 3                   !: 2nd    "                  "    --> j 
     55   INTEGER, PARAMETER ::   jpk     = jpkdta              !: number of vertical levels 
     56   ! starting position of the zoom  
     57   INTEGER, PARAMETER ::   jpizoom =   16                !: left bottom (i,j) indices of the zoom 
     58   INTEGER, PARAMETER ::   jpjzoom =   11                !: in data domain indices 
     59   ! Domain characteristics 
     60   INTEGER, PARAMETER ::   jperio  =   0                 !: lateral cond. type (between 0 and 6) 
     61   INTEGER, PARAMETER ::   jpisl   =   0                 !: number of islands 
     62   INTEGER, PARAMETER ::   jpnisl  =   0                 !: maximum number of points per island 
     63 
     64#  else 
     65   !------------!!-------------------------------------------------------- 
     66   !            !!   Default :               Full global domain of GYRE_R2 
     67   !------------!!-------------------------------------------------------- 
     68   ! global domain size                                 !!! *  global domain  * 
     69   INTEGER, PARAMETER ::   jpiglo  = jpidta              !: 1st dimension of global domain --> i 
     70   INTEGER, PARAMETER ::   jpjglo  = jpjdta              !: 2nd    "                  "    --> j 
     71   INTEGER, PARAMETER ::   jpk     = jpkdta              !: number of vertical levels 
     72   ! starting position of the zoom  
     73   INTEGER, PARAMETER ::   jpizoom =   1                 !: left bottom (i,j) indices of the zoom 
     74   INTEGER, PARAMETER ::   jpjzoom =   1                 !: in data domain indices 
     75   ! Domain characteristics 
     76   INTEGER, PARAMETER ::   jperio  =   1                 !: lateral cond. type (between 0 and 6) 
     77   INTEGER, PARAMETER ::   jpisl   =   1                 !: number of islands 
     78   INTEGER, PARAMETER ::   jpnisl  = jpiglo              !: maximum number of points per island 
     79#  endif 
     80 
     81# endif 
     82 
     83   !!--------------------------------------------------------------------- 
     84   !! grid parameters for domhgr and domzgr 
     85   !!--------------------------------------------------------------------- 
     86   ! Values set to pp_not_used indicates that this parameter is not used in THIS config. 
     87   ! Values set to pp_to_be_computed  indicates that variables will be computed in domzgr 
     88   REAL(wp), PARAMETER ::   pp_not_used       = 999999_wp   !: default value for "not used" 
     89   REAL(wp), PARAMETER ::   pp_to_be_computed = 0._wp       !: default value for "to be computed" 
     90 
     91   ! horizontal mesh (domhgr) 
     92   !    Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 ) 
     93   INTEGER, PARAMETER ::   jphgr_msh = 3   !: type of horizontal mesh 
     94   !                                       ! = 3 beta-plane with regular grid-spacing 
     95 
     96   !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid. 
     97   !   The mercator grid starts only approximately at gphi0 because 
     98   !   of the constraint that the equator be a T point. 
     99   REAL(wp) ,PARAMETER ::   ppglam0  =    0.0_wp     !: longitude of first raw and column T-point (jphgr_msh=1) 
     100   REAL(wp) ,PARAMETER ::   ppgphi0  =   35.0_wp     !: latitude  of first raw and column T-point (jphgr_msh=1) 
     101   !                                                 !  latitude  of the Coriolis/Beta parameter  (jphgr_msh=2 or 3) 
     102   REAL(wp) ,PARAMETER ::   ppe1_deg = pp_not_used   !: zonal      grid-spacing (degrees) 
     103   REAL(wp) ,PARAMETER ::   ppe2_deg = pp_not_used   !: meridional grid-spacing (degrees) 
     104   ! 
     105   REAL(wp) ,PARAMETER ::   ppe1_m   = 2000.0_wp     !: zonal      grid-spacing (meters ) 
     106   REAL(wp) ,PARAMETER ::   ppe2_m   = 2000.0_wp     !: meridional grid-spacing (meters ) 
     107 
     108   ! Vertical mesh (domzgr) 
     109   !    parameter used in the analytical function defining the repartition of the level and their thickness 
     110   REAL(wp), PARAMETER ::   ppsur = -2033.194295283385_wp  !: ORCA r4, r2 and r05 coefficients 
     111   REAL(wp), PARAMETER ::   ppa0  =  155.8325369664153_wp  !: (default coefficients) 
     112   REAL(wp), PARAMETER ::   ppa1  =  146.3615918601890_wp  !: 
     113   REAL(wp), PARAMETER ::   ppkth =  17.28520372419791_wp  !: (non dimensional): gives the approximate 
     114   !                                                       !: layer number above which  stretching will 
     115   !                                                       !: be maximum. Usually of order jpk/2. 
     116   REAL(wp), PARAMETER ::   ppacr =     5.00000000000_wp   !: (non dimensional): stretching factor 
     117   !                                                       !: for the grid. The highest zacr, the smallest 
     118   !                                                       !: the stretching. 
     119   !  If both ppa0 ppa1 and ppsur are specified to 0, then 
     120   !  they are computed from ppdzmin, pphmax, ppkth, ppacr in dom_zgr 
     121   REAL(wp), PARAMETER ::   ppdzmin = pp_not_used   !: (meters) vertical thickness of the top layer 
     122   REAL(wp), PARAMETER ::   pphmax  = pp_not_used   !: (meters) Maximum depth of the ocean gdepw(jpk) 
     123 
     124#else 
     125   !!--------------------------------------------------------------------- 
     126   !!   default option  :                                     Empty module 
     127   !!--------------------------------------------------------------------- 
    13128#endif 
    14       ::    &   
    15       cp_cfg = "eel"            !: name of the configuration 
    16    INTEGER     & 
    17 #if !defined key_agrif 
    18       , PARAMETER  & 
    19 #endif 
    20       :: & 
    21       jp_cfg = 2   ,         &  !: resolution of the configuration (km) 
    22129 
    23       ! data size              !!! * size of all the input files * 
    24       jpidta  = 83,          &  !: 1st horizontal dimension ( >= jpi ) 
    25       jpjdta  = 242,         &  !: 2nd    "            "    ( >= jpj ) 
    26       jpkdta  = 30,          &  !: number of levels         ( >= jpk ) 
    27  
    28       ! global domain size     !!! * full domain * 
    29       jpiglo  = jpidta,      &  !: 1st dimension of global domain --> i 
    30       jpjglo  = jpjdta,      &  !: 2nd    "                  "    --> j 
    31       jpk     = jpkdta,      &  !: number of vertical levels 
    32       ! zoom starting position 
    33       jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom 
    34       jpjzoom =   1   ,      &  !: in data indices 
    35  
    36       ! Domain characteristics 
    37       jperio  =     1 ,      &  !: lateral cond. type (between 0 and 6) 
    38       jpisl   =     1 ,      &  !: number of islands 
    39       jpnisl  = jpiglo          !: maximum number of points per island 
    40  
    41    !!  Values set to pp_not_used indicates that this parameter is not used in THIS config. 
    42    !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr 
    43    REAL(wp), PARAMETER ::   &  !: 
    44       pp_not_used       = 999999._wp  , & !: ??? 
    45       pp_to_be_computed =      0._wp      !: ??? 
    46    !! 
    47    !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 ) 
    48    !! 
    49    INTEGER,PARAMETER   ::    & !: 
    50       jphgr_msh = 3            !: type of horizontal mesh 
    51       !                        ! = 0 curvilinear coordinate on the sphere 
    52       !                        !     read in coordinate.nc file 
    53       !                        ! = 1 geographical mesh on the sphere 
    54       !                        !     with regular grid-spacing 
    55       !                        ! = 2 f-plane with regular grid-spacing 
    56       !                        ! = 3 beta-plane with regular grid-spacing 
    57       !                        ! = 4 Mercator grid with T/U point at the equator  with 
    58       !                        !     isotropic resolution (e1_deg) 
    59  
    60       !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid. 
    61       !   The mercator grid starts only approximately at gphi0 because 
    62       !   of the constraint that the equator be a T point. 
    63    REAL(wp) ,PARAMETER ::     &  !: 
    64       ppglam0  =    0.0_wp,   &  !: longitude of first raw and column T-point (jphgr_msh = 1) 
    65       ppgphi0  =   35.0_wp,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1) 
    66       !                          ! latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3) 
    67       ppe1_deg = pp_not_used ,   &  !: zonal      grid-spacing (degrees) 
    68       ppe2_deg = pp_not_used ,   &  !: meridional grid-spacing (degrees) 
    69       ! 
    70       ppe1_m   = 2000.0_wp,   &  !: zonal      grid-spacing (meters ) 
    71       ppe2_m   = 2000.0_wp       !: meridional grid-spacing (meters ) 
    72    !! 
    73    !!  Coefficients associated with the vertical coordinate system 
    74    !! 
    75    REAL(wp) & 
    76 #if !defined key_agrif 
    77       , PARAMETER  & 
    78 #endif 
    79       ::     &  !: 
    80       &     ppsur = -2033.194295283385_wp   ,  &  !: Computed in domzgr 
    81       &     ppa0  =  155.8325369664153_wp   ,  &  !: 
    82       &     ppa1  =  146.3615918601890_wp   ,  &  !: 
    83       ! 
    84       &     ppkth =  17.28520372419791_wp   ,  &  !: (non dimensional): gives the approximate 
    85       !                                           !: layer number above which  stretching will 
    86       !                                           !: be maximum. Usually of order jpk/2. 
    87       &     ppacr =  5.000000000000000_wp         !: (non dimensional): stretching factor 
    88       !                                           !: for the grid. The highest zacr, the smallest 
    89       !                                           !: the stretching. 
    90  
    91    !! 
    92    !!  If both ppa0 ppa1 and ppsur are specified to 0, then 
    93    !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 
    94    !! 
    95    REAL(wp), PARAMETER ::        & 
    96       &     ppdzmin = pp_not_used           ,  &  !: (meters) vertical thickness of the top layer 
    97       &     pphmax  = pp_not_used                 !: (meters) Maximum depth of the ocean gdepw(jpk) 
    98    !!--------------------------------------------------------------------- 
     130   !!====================================================================== 
     131END MODULE par_eel_r2 
Note: See TracChangeset for help on using the changeset viewer.