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 9169 for branches/2017/dev_merge_2017/NEMOGCM/CONFIG – NEMO

Ignore:
Timestamp:
2017-12-26T17:32:56+01:00 (6 years ago)
Author:
gm
Message:

dev_merge_2017: all SRC: finalize the removal of useless warning when reading namelist_cfg + remove all nn_closea + nn_msh replaced by a logical

Location:
branches/2017/dev_merge_2017/NEMOGCM/CONFIG
Files:
25 edited

Legend:

Unmodified
Added
Removed

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/C1D_PAPA/EXP00/namelist_cfg

    r9019 r9169  
    3939   nn_bathy    =    0      !  compute (=0) or read (=1) the bathymetry file 
    4040   rn_bathy    = 4198.     !  value of the bathymetry. if (=0) bottom flat at jpkm1 
    41    nn_msh      =    0      !  create (=1) a mesh file or not (=0) 
     41   ! 
    4242   rn_rdt      =  360.     !  time step for the dynamics  
    4343   ppglam0     =    -150.0               !  longitude of first raw and column T-point (jphgr_msh = 1) 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/GYRE_BFM/EXP00/namelist_cfg

    r9019 r9169  
    3737!----------------------------------------------------------------------- 
    3838   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time 
    39    ! 
    40    nn_msh      =    0      !  create (>0) a mesh file or not (=0) 
    4139   ! 
    4240   rn_rdt      = 7200.     !  time step for the dynamics (and tracer if nn_acc=0) 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/GYRE_PISCES/EXP00/namelist_cfg

    r9019 r9169  
    3737!----------------------------------------------------------------------- 
    3838   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time 
    39    nn_msh      =    0      !  create (>0) a mesh file or not (=0) 
     39   ! 
    4040   rn_rdt      = 7200.     !  time step for the dynamics (and tracer if nn_acc=0) 
    4141/ 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/ORCA2_LIM3_PISCES/EXP00/1_namelist_cfg

    r9063 r9169  
    2323!-----------------------------------------------------------------------   
    2424   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    25    nn_closea   =    0      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
    2625   ! 
    2726   rn_rdt      = 2880.     !  time step for the dynamics (and tracer if nn_acc=0) 
    28    ! 
    2927/ 
    3028!----------------------------------------------------------------------- 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/ORCA2_LIM3_PISCES/EXP00/namelist_cfg

    r9061 r9169  
    2121!----------------------------------------------------------------------- 
    2222   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    23    ! 
    24    nn_msh      =    0      !  create (>0) a mesh file or not (=0) 
    25    ! 
    2623/ 
    2724!----------------------------------------------------------------------- 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/ORCA2_LIM3_PISCES/IDL_scripts/README

    r7732 r9169  
    1919wget --recursive -l2 --no-directories --no-parent -A.nc -erobots=off http://dodsp.idris.fr/reee512/NEMO_OUT/ORCA2_LIM/ 
    2020 
    21 You will aslo need the meshmask file (set nn_msh = 1 or nn_msh = 7 in your namelist and run the model for at least 1 time step). 
     21You will aslo need the meshmask file (set ln_meshmask = TRUE in your namelist and run the model for at least 1 time step). 
    2222 
    2323#---------------------------------------------------------------------- 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/ORCA2_SAS_LIM3/EXP00/namelist_cfg

    r9019 r9169  
    2222!----------------------------------------------------------------------- 
    2323   ln_linssh   = .true.   !  =T  linear free surface  ==>>  model level are fixed in time 
    24    ! 
    25    nn_msh      =    0      !  create (>0) a mesh file or not (=0) 
    26    ! 
    2724/ 
    2825!----------------------------------------------------------------------- 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/SHARED/namelist_ref

    r9168 r9169  
    9292!----------------------------------------------------------------------- 
    9393   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    94    nn_msh      =    0      !  create (>0) a mesh file or not (=0) 
    95    rn_isfhmin  =    1.00   !  treshold (m) to discriminate grounding ice to floating ice 
     94   rn_isfhmin  =    1.00   !  treshold [m] to discriminate grounding ice from floating ice 
    9695   ! 
    9796   rn_rdt      = 5760.     !  time step for the dynamics and tracer 
    … …  
    9998   ! 
    10099   ln_crs      = .false.   !  Logical switch for coarsening module      (T => fill namcrs) 
     100   ! 
     101   ln_meshmask = .false.   !  =T create a mesh file 
    101102/ 
    102103!----------------------------------------------------------------------- 
    … …  
    134135                           !  1, coarse grid is binned with centering at the equator 
    135136                           !    Symmetry with nn_facty being odd-numbered. Asymmetry with even-numbered nn_facty. 
    136    nn_msh_crs  = 1         !  create (=1) a mesh file or not (=0) 
     137   ln_msh_crs  = .false.   ! =T create a mesh & mask file 
    137138   nn_crs_kz   = 0         ! 0, MEAN of volume boxes 
    138139                           ! 1, MAX of boxes 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/ISOMIP/EXP00/namelist_cfg

    r9135 r9169  
    4040   ln_linssh   = .false.    !  =T  linear free surface  ==>>  model level are fixed in time 
    4141   ! 
    42    nn_msh      =    0      !  create (>0) a mesh file or not (=0) 
    43    ! 
    4442   rn_rdt      = 1800.     !  time step for the dynamics (and tracer if nn_acc=0) 
     43   ! 
     44   ln_meshmask = .false.   !  =T create a mesh file 
    4545/ 
    4646!----------------------------------------------------------------------- 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT2_flux_cen2_cfg

    r9019 r9169  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    nn_closea   =    1      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
    35    ! 
    36    nn_msh      =    1      !  create (>0) a mesh file or not (=0) 
     34   ! 
    3735   rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3836   ! 
    … …  
    4139   ! 
    4240   ln_crs      = .false.   !  Logical switch for coarsening module 
     41   ! 
     42   ln_meshmask = .false.   !  =T create a mesh file 
    4343/ 
    4444!----------------------------------------------------------------------- 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT2_flux_ubs_cfg

    r9019 r9169  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    nn_closea   =    1      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
    35    ! 
    36    nn_msh      =    1      !  create (>0) a mesh file or not (=0) 
     34   ! 
    3735   rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3836   ! 
    … …  
    4139   ! 
    4240   ln_crs      = .false.   !  Logical switch for coarsening module 
     41   ! 
     42   ln_meshmask = .false.   !  =T create a mesh file 
    4343/ 
    4444!----------------------------------------------------------------------- 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT2_vect_eenH_cfg

    r9019 r9169  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    nn_closea   =    1      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
    35    ! 
    36    nn_msh      =    1      !  create (>0) a mesh file or not (=0) 
     34   ! 
    3735   rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3836   ! 
    … …  
    4139   ! 
    4240   ln_crs      = .false.   !  Logical switch for coarsening module 
     41   ! 
     42   ln_meshmask = .false.   !  =T create a mesh file 
    4343/ 
    4444!----------------------------------------------------------------------- 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT2_vect_een_cfg

    r9019 r9169  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    nn_closea   =    1      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
    35    ! 
    36    nn_msh      =    1      !  create (>0) a mesh file or not (=0) 
     34   ! 
    3735   rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3836   ! 
    … …  
    4139   ! 
    4240   ln_crs      = .false.   !  Logical switch for coarsening module 
     41   ! 
     42   ln_meshmask = .false.   !  =T create a mesh file 
    4343/ 
    4444!----------------------------------------------------------------------- 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT2_vect_ene_cfg

    r9019 r9169  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    nn_closea   =    1      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
    35    ! 
    36    nn_msh      =    1      !  create (>0) a mesh file or not (=0) 
     34   ! 
    3735   rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3836   ! 
    … …  
    4139   ! 
    4240   ln_crs      = .false.   !  Logical switch for coarsening module 
     41   ! 
     42   ln_meshmask = .false.   !  =T create a mesh file 
    4343/ 
    4444!----------------------------------------------------------------------- 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT2_vect_ens_cfg

    r9019 r9169  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    nn_closea   =    1      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
    35    ! 
    36    nn_msh      =    1      !  create (>0) a mesh file or not (=0) 
     34   ! 
    3735   rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3836   ! 
    … …  
    4139   ! 
    4240   ln_crs      = .false.   !  Logical switch for coarsening module 
     41   ! 
     42   ln_meshmask = .false.   !  =T create a mesh file 
    4343/ 
    4444!----------------------------------------------------------------------- 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT4_flux_cen2_cfg

    r9019 r9169  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    nn_closea   =    1      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
    35    ! 
    36    nn_msh      =    1      !  create (>0) a mesh file or not (=0) 
     34   ! 
    3735   rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3836   ! 
    … …  
    4139   ! 
    4240   ln_crs      = .false.   !  Logical switch for coarsening module 
     41   ! 
     42   ln_meshmask = .false.   !  =T create a mesh file 
    4343/ 
    4444!----------------------------------------------------------------------- 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT4_flux_ubs_cfg

    r9019 r9169  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    nn_closea   =    1      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
    35    ! 
    36    nn_msh      =    1      !  create (>0) a mesh file or not (=0) 
     34   ! 
    3735   rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3836   ! 
    … …  
    4139   ! 
    4240   ln_crs      = .false.   !  Logical switch for coarsening module 
     41   ! 
     42   ln_meshmask = .false.   !  =T create a mesh file 
    4343/ 
    4444!----------------------------------------------------------------------- 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT4_vect_eenH_cfg

    r9019 r9169  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    nn_closea   =    1      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
    35    ! 
    36    nn_msh      =    1      !  create (>0) a mesh file or not (=0) 
     34   ! 
    3735   rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3836   ! 
    … …  
    4139   ! 
    4240   ln_crs      = .false.   !  Logical switch for coarsening module 
     41   ! 
     42   ln_meshmask = .false.   !  =T create a mesh file 
    4343/ 
    4444!----------------------------------------------------------------------- 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT4_vect_een_cfg

    r9019 r9169  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    nn_closea   =    1      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
    35    ! 
    36    nn_msh      =    1      !  create (>0) a mesh file or not (=0) 
     34   ! 
    3735   rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3836   ! 
    … …  
    4139   ! 
    4240   ln_crs      = .false.   !  Logical switch for coarsening module 
     41   ! 
     42   ln_meshmask = .false.   !  =T create a mesh file 
    4343/ 
    4444!----------------------------------------------------------------------- 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT4_vect_ene_cfg

    r9019 r9169  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    nn_closea   =    1      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
    35    ! 
    36    nn_msh      =    1      !  create (>0) a mesh file or not (=0) 
     34   ! 
    3735   rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3836   ! 
    … …  
    4139   ! 
    4240   ln_crs      = .false.   !  Logical switch for coarsening module 
     41   ! 
     42   ln_meshmask = .false.   !  =T create a mesh file 
    4343/ 
    4444!----------------------------------------------------------------------- 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT4_vect_ens_cfg

    r9019 r9169  
    3232!----------------------------------------------------------------------- 
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    34    nn_closea   =    1      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
    35    ! 
    36    nn_msh      =    1      !  create (>0) a mesh file or not (=0) 
     34   ! 
    3735   rn_isfhmin  =    0.00   !  treshold (m) to discriminate grounding ice to floating ice 
    3836   ! 
    … …  
    4139   ! 
    4240   ln_crs      = .false.   !  Logical switch for coarsening module 
     41   ! 
     42   ln_meshmask = .false.   !  =T create a mesh file 
    4343/ 
    4444!----------------------------------------------------------------------- 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER/EXP00/1_namelist_cfg

    r9060 r9169  
    2929   ln_linssh   = .true.   !  =T  linear free surface  ==>>  model level are fixed in time 
    3030   ! 
    31    nn_msh      =    1     !  create (>0) a mesh file or not (=0) 
    32    ! 
    3331   rn_rdt      = 400.     !  time step for the dynamics (and tracer if nn_acc=0) 
     32   ! 
     33   ln_meshmask = .true.   !  =T  create a mesh file 
    3434/ 
    3535!----------------------------------------------------------------------- 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER/EXP00/namelist_cfg

    r9019 r9169  
    2929   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time 
    3030   ! 
    31    nn_msh      =    1      !  create (>0) a mesh file or not (=0) 
    32    ! 
    3331   rn_rdt      = 1200.     !  time step for the dynamics (and tracer if nn_acc=0) 
     32   ! 
     33   ln_meshmask = .true.   !  =T  create a mesh file 
    3434/ 
    3535!----------------------------------------------------------------------- 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/WAD/EXP00/namelist_cfg

    r9024 r9169  
    1515   cn_exp      =  "WAD"    !  experience name 
    1616   nn_it000    =       1   !  first time step 
    17    nn_itend    =      3840  !  last  time step 
    18    !nn_itend    =         6  !  last  time step 
     17   nn_itend    =     3840  !  last  time step 
     18   !nn_itend   =       6   !  last  time step 
    1919   nn_leapy    =      30   !  Leap year calendar (1) or not (0) 
    2020   nn_stock    =    48000  !  frequency of creation of a restart file (modulo referenced to 1) 
    … …  
    3636   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    3737   ! 
    38    nn_msh      =    1      !  create (>0) a mesh file or not (=0) 
     38   ln_meshmask = .true.    !  create (>0) a mesh file or not (=0) 
    3939   rn_rdt      =    18.    !  time step for the dynamics  
    4040/ 

  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/WAD/MY_SRC/domain.F90

    r9168 r9169  
    3030   USE trc_oce        ! shared ocean & passive tracers variab 
    3131   USE phycst         ! physical constants 
    32    USE usrdef_closea  ! closed seas 
     32   USE closea         ! closed seas 
    3333   USE domhgr         ! domain: set the horizontal mesh 
    3434   USE domzgr         ! domain: set the vertical mesh 
    … …  
    3737   USE domvvl         ! variable volume 
    3838   USE c1d            ! 1D configuration 
    39    USE domc1d         ! 1D configuration: column location 
    4039   USE dyncor_c1d     ! 1D configuration: Coriolis term    (cor_c1d routine) 
    4140   USE wet_dry        ! wetting and drying 
    … …  
    7170      !!                         and scale factors, and the coriolis factor 
    7271      !!              - dom_zgr: define the vertical coordinate and the bathymetry 
    73       !!              - dom_wri: create the meshmask file if nn_msh=1 
     72      !!              - dom_wri: create the meshmask file (ln_meshmask=T) 
    7473      !!              - 1D configuration, move Coriolis, u and v at T-point 
    7574      !!---------------------------------------------------------------------- 
    … …  
    9493         WRITE(numout,cform) '       ' ,'   jpij    : ', jpij 
    9594         WRITE(numout,*)     '      mpp local domain info (mpp):' 
    96          WRITE(numout,*)     '              jpni    : ', jpni, '   jpreci  : ', jpreci 
    97          WRITE(numout,*)     '              jpnj    : ', jpnj, '   jprecj  : ', jprecj 
     95         WRITE(numout,*)     '              jpni    : ', jpni, '   nn_hls  : ', nn_hls 
     96         WRITE(numout,*)     '              jpnj    : ', jpnj, '   nn_hls  : ', nn_hls 
    9897         WRITE(numout,*)     '              jpnij   : ', jpnij 
    9998         WRITE(numout,*)     '      lateral boundary of the Global domain : jperio  = ', jperio 
    … …  
    106105         CASE( 5 )   ;   WRITE(numout,*) '         (i.e. north fold with F-point pivot)' 
    107106         CASE( 6 )   ;   WRITE(numout,*) '         (i.e. cyclic east-west and north fold with F-point pivot)' 
    108          CASE( 7 )   ;   WRITE(numout,*) '         (i.e. cyclic north-south (jpnij=1 only))' 
     107         CASE( 7 )   ;   WRITE(numout,*) '         (i.e. cyclic east-west and north-south)' 
    109108         CASE DEFAULT 
    110109            CALL ctl_stop( 'jperio is out of range' ) 
    111110         END SELECT 
    112111         WRITE(numout,*)     '      Ocean model configuration used:' 
    113          WRITE(numout,*)     '              cn_cfg = ', cn_cfg 
    114          WRITE(numout,*)     '              nn_cfg = ', nn_cfg 
    115       ENDIF 
    116       ! 
    117       !       
    118 !!gm  This should be removed with the new configuration interface 
    119       IF( lk_c1d .AND. ln_c1d_locpt )  CALL dom_c1d( rn_lat1d, rn_lon1d ) 
    120 !!gm end 
     112         WRITE(numout,*)     '         cn_cfg = ', TRIM( cn_cfg ), '   nn_cfg = ', nn_cfg 
     113      ENDIF 
    121114      ! 
    122115      !           !==  Reference coordinate system  ==! 
    … …  
    124117      CALL dom_glo                     ! global domain versus local domain 
    125118      CALL dom_nam                     ! read namelist ( namrun, namdom ) 
    126       CALL dom_clo( cn_cfg, nn_cfg )   ! Closed seas and lake 
    127119      CALL dom_hgr                     ! Horizontal mesh 
    128120      CALL dom_zgr( ik_top, ik_bot )   ! Vertical mesh and bathymetry 
    129       IF( nn_closea == 0 )   CALL clo_bat( ik_top, ik_bot )    !==  remove closed seas or lakes  ==! 
    130121      CALL dom_msk( ik_top, ik_bot )   ! Masks 
     122      IF( ln_closea ) CALL dom_clo     ! ln_closea=T : closed seas included in the simulation 
     123                                       ! Read in masks to define closed seas and lakes  
    131124      ! 
    132125      DO jj = 1, jpj                   ! depth of the iceshelves 
    … …  
    182175      IF( lk_c1d         )   CALL cor_c1d       ! 1D configuration: Coriolis set at T-point 
    183176      ! 
    184       IF( nn_msh > 0 .AND. .NOT. ln_iscpl )                         CALL dom_wri      ! Create a domain file 
    185       IF( nn_msh > 0 .AND.       ln_iscpl .AND. .NOT. ln_rstart )   CALL dom_wri      ! Create a domain file 
    186       IF( .NOT.ln_rstart )   CALL dom_ctl       ! Domain control 
    187       ! 
    188        
     177      IF( ln_meshmask .AND. .NOT.ln_iscpl )                        CALL dom_wri     ! Create a domain file 
     178      IF( ln_meshmask .AND.      ln_iscpl .AND. .NOT.ln_rstart )   CALL dom_wri     ! Create a domain file 
     179      IF(                                       .NOT.ln_rstart )   CALL dom_ctl     ! Domain control 
     180      ! 
     181      IF( ln_write_cfg )   CALL cfg_write         ! create the configuration file 
     182      ! 
    189183      IF(lwp) THEN 
    190184         WRITE(numout,*) 
    191          WRITE(numout,*) 'dom_init : end of domain initialization nn_msh=', nn_msh 
     185         WRITE(numout,*) 'dom_init :   ==>>>   END of domain initialization' 
     186         WRITE(numout,*) '~~~~~~~~' 
    192187         WRITE(numout,*)  
    193188      ENDIF 
    194       ! 
    195       IF( ln_write_cfg )   CALL cfg_write         ! create the configuration file 
    196189      ! 
    197190   END SUBROUTINE dom_init 
    … …  
    240233         IF( nn_print >= 1 ) THEN 
    241234            WRITE(numout,*) 
    242             WRITE(numout,*) '          conversion local  ==> global i-index domain' 
     235            WRITE(numout,*) '          conversion local  ==> global i-index domain (mig)' 
    243236            WRITE(numout,25)              (mig(ji),ji = 1,jpi) 
    244237            WRITE(numout,*) 
    245238            WRITE(numout,*) '          conversion global ==> local  i-index domain' 
    246             WRITE(numout,*) '             starting index' 
     239            WRITE(numout,*) '             starting index (mi0)' 
    247240            WRITE(numout,25)              (mi0(ji),ji = 1,jpiglo) 
    248             WRITE(numout,*) '             ending index' 
     241            WRITE(numout,*) '             ending index (mi1)' 
    249242            WRITE(numout,25)              (mi1(ji),ji = 1,jpiglo) 
    250243            WRITE(numout,*) 
    251             WRITE(numout,*) '          conversion local  ==> global j-index domain' 
     244            WRITE(numout,*) '          conversion local  ==> global j-index domain (mjg)' 
    252245            WRITE(numout,25)              (mjg(jj),jj = 1,jpj) 
    253246            WRITE(numout,*) 
    254247            WRITE(numout,*) '          conversion global ==> local  j-index domain' 
    255             WRITE(numout,*) '             starting index' 
     248            WRITE(numout,*) '             starting index (mj0)' 
    256249            WRITE(numout,25)              (mj0(jj),jj = 1,jpjglo) 
    257             WRITE(numout,*) '             ending index' 
     250            WRITE(numout,*) '             ending index (mj1)' 
    258251            WRITE(numout,25)              (mj1(jj),jj = 1,jpjglo) 
    259252         ENDIF 
    … …  
    275268      !!---------------------------------------------------------------------- 
    276269      USE ioipsl 
     270      !! 
     271      INTEGER  ::   ios   ! Local integer 
     272      ! 
    277273      NAMELIST/namrun/ cn_ocerst_indir, cn_ocerst_outdir, nn_stocklist, ln_rst_list,                 & 
    278274         &             nn_no   , cn_exp   , cn_ocerst_in, cn_ocerst_out, ln_rstart , nn_rstctl ,     & 
    … …  
    280276         &             nn_stock, nn_write , ln_mskland  , ln_clobber   , nn_chunksz, nn_euler  ,     & 
    281277         &             ln_cfmeta, ln_iscpl 
    282       NAMELIST/namdom/ ln_linssh, nn_closea, nn_msh, rn_isfhmin, rn_rdt, rn_atfp, ln_crs 
     278      NAMELIST/namdom/ ln_linssh, rn_isfhmin, rn_rdt, rn_atfp, ln_crs, ln_meshmask 
    283279#if defined key_netcdf4 
    284280      NAMELIST/namnc4/ nn_nchunks_i, nn_nchunks_j, nn_nchunks_k, ln_nc4zip 
    285281#endif 
    286       INTEGER  ::   ios                 ! Local integer output status for namelist read 
    287       !!---------------------------------------------------------------------- 
     282      !!---------------------------------------------------------------------- 
     283      ! 
     284      IF(lwp) THEN 
     285         WRITE(numout,*) 
     286         WRITE(numout,*) 'dom_nam  : domain initialization through namelist read' 
     287         WRITE(numout,*) '~~~~~~~ ' 
     288      ENDIF 
    288289      ! 
    289290      REWIND( numnam_ref )              ! Namelist namrun in reference namelist : Parameters of the run 
    290291      READ  ( numnam_ref, namrun, IOSTAT = ios, ERR = 901) 
    291292901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namrun in reference namelist', lwp ) 
    292       ! 
    293293      REWIND( numnam_cfg )              ! Namelist namrun in configuration namelist : Parameters of the run 
    294294      READ  ( numnam_cfg, namrun, IOSTAT = ios, ERR = 902 ) 
    … …  
    297297      ! 
    298298      IF(lwp) THEN                  ! control print 
    299          WRITE(numout,*) 
    300          WRITE(numout,*) 'dom_nam  : domain initialization through namelist read' 
    301          WRITE(numout,*) '~~~~~~~ ' 
    302          WRITE(numout,*) '   Namelist namrun' 
    303          WRITE(numout,*) '      job number                      nn_no      = ', nn_no 
    304          WRITE(numout,*) '      experiment name for output      cn_exp     = ', cn_exp 
    305          WRITE(numout,*) '      file prefix restart input       cn_ocerst_in= ', cn_ocerst_in 
    306          WRITE(numout,*) '      restart input directory         cn_ocerst_indir= ', cn_ocerst_indir 
    307          WRITE(numout,*) '      file prefix restart output      cn_ocerst_out= ', cn_ocerst_out 
    308          WRITE(numout,*) '      restart output directory        cn_ocerst_outdir= ', cn_ocerst_outdir 
    309          WRITE(numout,*) '      restart logical                 ln_rstart  = ', ln_rstart 
    310          WRITE(numout,*) '      start with forward time step    nn_euler   = ', nn_euler 
    311          WRITE(numout,*) '      control of time step            nn_rstctl  = ', nn_rstctl 
    312          WRITE(numout,*) '      number of the first time step   nn_it000   = ', nn_it000 
    313          WRITE(numout,*) '      number of the last time step    nn_itend   = ', nn_itend 
    314          WRITE(numout,*) '      initial calendar date aammjj    nn_date0   = ', nn_date0 
    315          WRITE(numout,*) '      initial time of day in hhmm     nn_time0   = ', nn_time0 
    316          WRITE(numout,*) '      leap year calendar (0/1)        nn_leapy   = ', nn_leapy 
    317          WRITE(numout,*) '      initial state output            nn_istate  = ', nn_istate 
     299         WRITE(numout,*) '   Namelist : namrun' 
     300         WRITE(numout,*) '      job number                      nn_no           = ', nn_no 
     301         WRITE(numout,*) '      experiment name for output      cn_exp          = ', TRIM( cn_exp           ) 
     302         WRITE(numout,*) '      file prefix restart input       cn_ocerst_in    = ', TRIM( cn_ocerst_in     ) 
     303         WRITE(numout,*) '      restart input directory         cn_ocerst_indir = ', TRIM( cn_ocerst_indir  ) 
     304         WRITE(numout,*) '      file prefix restart output      cn_ocerst_out   = ', TRIM( cn_ocerst_out    ) 
     305         WRITE(numout,*) '      restart output directory        cn_ocerst_outdir= ', TRIM( cn_ocerst_outdir ) 
     306         WRITE(numout,*) '      restart logical                 ln_rstart       = ', ln_rstart 
     307         WRITE(numout,*) '      start with forward time step    nn_euler        = ', nn_euler 
     308         WRITE(numout,*) '      control of time step            nn_rstctl       = ', nn_rstctl 
     309         WRITE(numout,*) '      number of the first time step   nn_it000        = ', nn_it000 
     310         WRITE(numout,*) '      number of the last time step    nn_itend        = ', nn_itend 
     311         WRITE(numout,*) '      initial calendar date aammjj    nn_date0        = ', nn_date0 
     312         WRITE(numout,*) '      initial time of day in hhmm     nn_time0        = ', nn_time0 
     313         WRITE(numout,*) '      leap year calendar (0/1)        nn_leapy        = ', nn_leapy 
     314         WRITE(numout,*) '      initial state output            nn_istate       = ', nn_istate 
    318315         IF( ln_rst_list ) THEN 
    319             WRITE(numout,*) '      list of restart dump times      nn_stocklist   =', nn_stocklist 
     316            WRITE(numout,*) '      list of restart dump times      nn_stocklist    =', nn_stocklist 
    320317         ELSE 
    321             WRITE(numout,*) '      frequency of restart file       nn_stock   = ', nn_stock 
     318            WRITE(numout,*) '      frequency of restart file       nn_stock        = ', nn_stock 
    322319         ENDIF 
    323          WRITE(numout,*) '      frequency of output file        nn_write   = ', nn_write 
    324          WRITE(numout,*) '      mask land points                ln_mskland = ', ln_mskland 
    325          WRITE(numout,*) '      additional CF standard metadata ln_cfmeta  = ', ln_cfmeta 
    326          WRITE(numout,*) '      overwrite an existing file      ln_clobber = ', ln_clobber 
    327          WRITE(numout,*) '      NetCDF chunksize (bytes)        nn_chunksz = ', nn_chunksz 
    328          WRITE(numout,*) '      IS coupling at the restart step ln_iscpl   = ', ln_iscpl 
     320         WRITE(numout,*) '      frequency of output file        nn_write        = ', nn_write 
     321         WRITE(numout,*) '      mask land points                ln_mskland      = ', ln_mskland 
     322         WRITE(numout,*) '      additional CF standard metadata ln_cfmeta       = ', ln_cfmeta 
     323         WRITE(numout,*) '      overwrite an existing file      ln_clobber      = ', ln_clobber 
     324         WRITE(numout,*) '      NetCDF chunksize (bytes)        nn_chunksz      = ', nn_chunksz 
     325         WRITE(numout,*) '      IS coupling at the restart step ln_iscpl        = ', ln_iscpl 
    329326      ENDIF 
    330327 
    … …  
    341338      nwrite = nn_write 
    342339      neuler = nn_euler 
    343       IF ( neuler == 1 .AND. .NOT. ln_rstart ) THEN 
    344          WRITE(ctmp1,*) 'ln_rstart =.FALSE., nn_euler is forced to 0 ' 
    345          CALL ctl_warn( ctmp1 ) 
     340      IF( neuler == 1 .AND. .NOT. ln_rstart ) THEN 
     341         IF(lwp) WRITE(numout,*)   
     342         IF(lwp) WRITE(numout,*)'   ==>>>   Start from rest (ln_rstart=F)' 
     343         IF(lwp) WRITE(numout,*)'           an Euler initial time step is used : nn_euler is forced to 0 '    
    346344         neuler = 0 
    347345      ENDIF 
    348346      !                             ! control of output frequency 
    349       IF ( nstock == 0 .OR. nstock > nitend ) THEN 
     347      IF( nstock == 0 .OR. nstock > nitend ) THEN 
    350348         WRITE(ctmp1,*) 'nstock = ', nstock, ' it is forced to ', nitend 
    351349         CALL ctl_warn( ctmp1 ) 
    … …  
    379377      READ  ( numnam_ref, namdom, IOSTAT = ios, ERR = 903) 
    380378903   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namdom in reference namelist', lwp ) 
    381       ! 
    382379      REWIND( numnam_cfg )              ! Namelist namdom in configuration namelist : space & time domain (bathymetry, mesh, timestep) 
    383380      READ  ( numnam_cfg, namdom, IOSTAT = ios, ERR = 904 ) 
    384381904   IF( ios >  0 )   CALL ctl_nam ( ios , 'namdom in configuration namelist', lwp ) 
    385       IF(lwm) WRITE ( numond, namdom ) 
     382      IF(lwm) WRITE( numond, namdom ) 
    386383      ! 
    387384      IF(lwp) THEN 
    388385         WRITE(numout,*) 
    389          WRITE(numout,*) '   Namelist namdom : space & time domain' 
    390          WRITE(numout,*) '      linear free surface (=T)              ln_linssh  = ', ln_linssh 
    391          WRITE(numout,*) '      suppression of closed seas (=0)       nn_closea  = ', nn_closea 
    392          WRITE(numout,*) '      create mesh/mask file(s)              nn_msh     = ', nn_msh 
    393          WRITE(numout,*) '           = 0   no file created           ' 
    394          WRITE(numout,*) '           = 1   mesh_mask                 ' 
    395          WRITE(numout,*) '           = 2   mesh and mask             ' 
    396          WRITE(numout,*) '           = 3   mesh_hgr, msh_zgr and mask' 
    397          WRITE(numout,*) '      treshold to open the isf cavity       rn_isfhmin = ', rn_isfhmin, ' (m)' 
    398          WRITE(numout,*) '      ocean time step                       rn_rdt     = ', rn_rdt 
    399          WRITE(numout,*) '      asselin time filter parameter         rn_atfp    = ', rn_atfp 
    400          WRITE(numout,*) '      online coarsening of dynamical fields ln_crs     = ', ln_crs 
    401       ENDIF 
    402        
    403       call flush( numout ) 
    404       ! 
    405 !     !          ! conversion DOCTOR names into model names (this should disappear soon) 
    406       atfp      = rn_atfp 
    407       rdt       = rn_rdt 
     386         WRITE(numout,*) '   Namelist : namdom   ---   space & time domain' 
     387         WRITE(numout,*) '      linear free surface (=T)                ln_linssh   = ', ln_linssh 
     388         WRITE(numout,*) '      create mesh/mask file                   ln_meshmask = ', ln_meshmask 
     389         WRITE(numout,*) '      treshold to open the isf cavity         rn_isfhmin  = ', rn_isfhmin, ' [m]' 
     390         WRITE(numout,*) '      ocean time step                         rn_rdt      = ', rn_rdt 
     391         WRITE(numout,*) '      asselin time filter parameter           rn_atfp     = ', rn_atfp 
     392         WRITE(numout,*) '      online coarsening of dynamical fields   ln_crs      = ', ln_crs 
     393      ENDIF 
     394      ! 
     395      !          ! conversion DOCTOR names into model names (this should disappear soon) 
     396      atfp = rn_atfp 
     397      rdt  = rn_rdt 
    408398 
    409399#if defined key_netcdf4 
    … …  
    412402      READ  ( numnam_ref, namnc4, IOSTAT = ios, ERR = 907) 
    413403907   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namnc4 in reference namelist', lwp ) 
    414       ! 
    415404      REWIND( numnam_cfg )              ! Namelist namnc4 in configuration namelist : NETCDF 
    416405      READ  ( numnam_cfg, namnc4, IOSTAT = ios, ERR = 908 ) 
    … …  
    421410         WRITE(numout,*) 
    422411         WRITE(numout,*) '   Namelist namnc4 - Netcdf4 chunking parameters' 
    423          WRITE(numout,*) '      number of chunks in i-dimension      nn_nchunks_i  = ', nn_nchunks_i 
    424          WRITE(numout,*) '      number of chunks in j-dimension      nn_nchunks_j  = ', nn_nchunks_j 
    425          WRITE(numout,*) '      number of chunks in k-dimension      nn_nchunks_k  = ', nn_nchunks_k 
    426          WRITE(numout,*) '      apply netcdf4/hdf5 chunking & compression ln_nc4zip = ', ln_nc4zip 
     412         WRITE(numout,*) '      number of chunks in i-dimension             nn_nchunks_i = ', nn_nchunks_i 
     413         WRITE(numout,*) '      number of chunks in j-dimension             nn_nchunks_j = ', nn_nchunks_j 
     414         WRITE(numout,*) '      number of chunks in k-dimension             nn_nchunks_k = ', nn_nchunks_k 
     415         WRITE(numout,*) '      apply netcdf4/hdf5 chunking & compression   ln_nc4zip    = ', ln_nc4zip 
    427416      ENDIF 
    428417 
    … …  
    586575      !                       ! ============================= ! 
    587576      !          
    588       clnam = 'domcfg_out'  ! filename (configuration information) 
     577      clnam = cn_domcfg_out  ! filename (configuration information) 
    589578      CALL iom_open( TRIM(clnam), inum, ldwrt = .TRUE., kiolib = jprstlib ) 
    590579       
Note: See TracChangeset for help on using the changeset viewer.