Changeset 11320


Ignore:
Timestamp:
2019-07-22T14:52:29+02:00 (13 months ago)
Author:
smasson
Message:

dev_r10984_HPC-13 : add missing if(lwp) in [11317], see #2307 and #2285

Location:
NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization
Files:
12 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/cfgs/C1D_PAPA/MY_SRC/usrdef_nam.F90

    r11317 r11320  
    7474      kpj = 3 
    7575      kpk = 75  
    76       ! 
    77       !                             ! control print 
    78       WRITE(numout,*) '   ' 
    79       WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
    80       WRITE(numout,*) '~~~~~~~~~~~ ' 
    81       WRITE(numout,*) '   Namelist namusr_def : C1 case' 
    82       WRITE(numout,*) '      type of vertical coordinate : ' 
    83       WRITE(numout,*) '         z-coordinate flag                     ln_zco = ', ln_zco 
    84       WRITE(numout,*) '         z-partial-step coordinate flag        ln_zps = ', ln_zps 
    85       WRITE(numout,*) '         s-coordinate flag                     ln_sco = ', ln_sco 
    86       WRITE(numout,*) '      C1D domain = 3 x 3 x75 grid-points                ' 
    87       WRITE(numout,*) '         resulting global domain size :        jpiglo = ', kpi 
    88       WRITE(numout,*) '                                               jpjglo = ', kpj 
    89       WRITE(numout,*) '                                               jpkglo = ', kpk 
    90  
    91       ! 
    9276      !                             ! Set the lateral boundary condition of the global domain 
    9377      kperio =  7                   ! C1D configuration : 3x3 basin with cyclic Est-West and Norht-South condition 
    9478      ! 
    95       WRITE(numout,*) '   Lateral boundary condition of the global domain' 
    96       WRITE(numout,*) '      C1D : closed basin                 jperio = ', kperio 
     79      !                             ! control print 
     80      IF(lwp) THEN 
     81         WRITE(numout,*) '   ' 
     82         WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
     83         WRITE(numout,*) '~~~~~~~~~~~ ' 
     84         WRITE(numout,*) '   Namelist namusr_def : C1 case' 
     85         WRITE(numout,*) '      type of vertical coordinate : ' 
     86         WRITE(numout,*) '         z-coordinate flag                     ln_zco = ', ln_zco 
     87         WRITE(numout,*) '         z-partial-step coordinate flag        ln_zps = ', ln_zps 
     88         WRITE(numout,*) '         s-coordinate flag                     ln_sco = ', ln_sco 
     89         WRITE(numout,*) '      C1D domain = 3 x 3 x75 grid-points                ' 
     90         WRITE(numout,*) '         resulting global domain size :        jpiglo = ', kpi 
     91         WRITE(numout,*) '                                               jpjglo = ', kpj 
     92         WRITE(numout,*) '                                               jpkglo = ', kpk 
     93         WRITE(numout,*) '   Lateral boundary condition of the global domain' 
     94         WRITE(numout,*) '      C1D : closed basin                       jperio = ', kperio 
     95      ENDIF 
    9796      ! 
    9897   END SUBROUTINE usr_def_nam 
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/OCE/DOM/domain.F90

    r11317 r11320  
    530530      !!---------------------------------------------------------------------- 
    531531      ! 
    532       WRITE(numout,*) '           ' 
    533       WRITE(numout,*) 'domain_cfg : domain size read in ', TRIM( cn_domcfg ), ' file' 
    534       WRITE(numout,*) '~~~~~~~~~~ ' 
     532      IF(lwp) THEN 
     533         WRITE(numout,*) '           ' 
     534         WRITE(numout,*) 'domain_cfg : domain size read in ', TRIM( cn_domcfg ), ' file' 
     535         WRITE(numout,*) '~~~~~~~~~~ ' 
     536      ENDIF 
    535537      ! 
    536538      CALL iom_open( cn_domcfg, inum ) 
     
    543545         CALL iom_get( inum, 'ORCA_index', zorca_res )   ;   kk_cfg = NINT( zorca_res ) 
    544546         ! 
    545          WRITE(numout,*) '   .' 
    546          WRITE(numout,*) '   ==>>>   ORCA configuration ' 
    547          WRITE(numout,*) '   .' 
     547         IF(lwp) THEN 
     548            WRITE(numout,*) '   .' 
     549            WRITE(numout,*) '   ==>>>   ORCA configuration ' 
     550            WRITE(numout,*) '   .' 
     551         ENDIF 
    548552         ! 
    549553      ELSE                                !- cd_cfg & k_cfg are not used 
     
    566570      CALL iom_close( inum ) 
    567571      ! 
    568       WRITE(numout,*) '      cn_cfg = ', TRIM(cd_cfg), '   nn_cfg = ', kk_cfg 
    569       WRITE(numout,*) '      jpiglo = ', kpi 
    570       WRITE(numout,*) '      jpjglo = ', kpj 
    571       WRITE(numout,*) '      jpkglo = ', kpk 
    572       WRITE(numout,*) '      type of global domain lateral boundary   jperio = ', kperio 
     572      IF(lwp) THEN 
     573         WRITE(numout,*) '      cn_cfg = ', TRIM(cd_cfg), '   nn_cfg = ', kk_cfg 
     574         WRITE(numout,*) '      jpiglo = ', kpi 
     575         WRITE(numout,*) '      jpjglo = ', kpj 
     576         WRITE(numout,*) '      jpkglo = ', kpk 
     577         WRITE(numout,*) '      type of global domain lateral boundary   jperio = ', kperio 
     578      ENDIF 
    573579      !         
    574580   END SUBROUTINE domain_cfg 
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/OCE/USR/usrdef_nam.F90

    r11317 r11320  
    8080#endif 
    8181      kpk = jpkglo 
    82       ! 
    83       !                             ! control print 
    84       WRITE(numout,*) '   ' 
    85       WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
    86       WRITE(numout,*) '~~~~~~~~~~~ ' 
    87       WRITE(numout,*) '   Namelist namusr_def : GYRE case' 
    88       WRITE(numout,*) '      GYRE used as Benchmark (=T)                      ln_bench  = ', ln_bench 
    89       WRITE(numout,*) '      inverse resolution & implied domain size         nn_GYRE   = ', nn_GYRE 
    90 #if defined key_agrif 
    91       IF( Agrif_Root() ) THEN 
    92 #endif 
    93       WRITE(numout,*) '         jpiglo = 30*nn_GYRE+2                            jpiglo = ', kpi 
    94       WRITE(numout,*) '         jpjglo = 20*nn_GYRE+2                            jpjglo = ', kpj 
    95 #if defined key_agrif 
    96       ENDIF 
    97 #endif 
    98       WRITE(numout,*) '      number of model levels                              jpkglo = ', kpk 
    99       ! 
    10082      !                             ! Set the lateral boundary condition of the global domain 
    10183      kperio = 0                    ! GYRE configuration : closed domain 
    10284      ! 
    103       WRITE(numout,*) '   ' 
    104       WRITE(numout,*) '   Lateral b.c. of the global domain set to closed     jperio = ', kperio 
     85      !                             ! control print 
     86      IF(lwp) THEN 
     87         WRITE(numout,*) '   ' 
     88         WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
     89         WRITE(numout,*) '~~~~~~~~~~~ ' 
     90         WRITE(numout,*) '   Namelist namusr_def : GYRE case' 
     91         WRITE(numout,*) '      GYRE used as Benchmark (=T)                      ln_bench  = ', ln_bench 
     92         WRITE(numout,*) '      inverse resolution & implied domain size         nn_GYRE   = ', nn_GYRE 
     93#if defined key_agrif 
     94         IF( Agrif_Root() ) THEN 
     95#endif 
     96         WRITE(numout,*) '         jpiglo = 30*nn_GYRE+2                            jpiglo = ', kpi 
     97         WRITE(numout,*) '         jpjglo = 20*nn_GYRE+2                            jpjglo = ', kpj 
     98#if defined key_agrif 
     99         ENDIF 
     100#endif 
     101         WRITE(numout,*) '      number of model levels                              jpkglo = ', kpk 
     102         WRITE(numout,*) '   ' 
     103         WRITE(numout,*) '   Lateral b.c. of the global domain set to closed        jperio = ', kperio 
     104      ENDIF 
    105105      ! 
    106106   END SUBROUTINE usr_def_nam 
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/tests/BENCH/MY_SRC/usrdef_nam.F90

    r11317 r11320  
    6969      cd_cfg = 'BENCH'             ! name & resolution (not used) 
    7070      kk_cfg = 0 
    71  
     71      ! 
    7272      IF( nn_isize < 0 .AND. nn_jsize < 0 ) THEN 
    7373      ! 
     
    8686         kpj = nn_jsize 
    8787      ENDIF 
    88   
     88      ! 
    8989      kpk = nn_ksize 
    9090      kperio = nn_perio 
    91  
    9291      !                             ! control print 
    93       WRITE(numout,*) '   ' 
    94       WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
    95       WRITE(numout,*) '~~~~~~~~~~~ ' 
    96       WRITE(numout,*) '   Namelist namusr_def : BENCH test case' 
    97       IF( nn_isize > 0 ) THEN 
    98          WRITE(numout,*) '      global domain size-x            nn_isize = ',  nn_isize 
    99       ELSE 
    100          WRITE(numout,*) '                                          jpni = ', jpni 
    101          WRITE(numout,*) '       local domain size-x           -nn_isize = ', -nn_isize 
    102          WRITE(numout,*) '      global domain size-x                 kpi = ', kpi 
     92      IF(lwp) THEN 
     93         WRITE(numout,*) '   ' 
     94         WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
     95         WRITE(numout,*) '~~~~~~~~~~~ ' 
     96         WRITE(numout,*) '   Namelist namusr_def : BENCH test case' 
     97         IF( nn_isize > 0 ) THEN 
     98            WRITE(numout,*) '      global domain size-x            nn_isize = ',  nn_isize 
     99         ELSE 
     100            WRITE(numout,*) '                                          jpni = ', jpni 
     101            WRITE(numout,*) '       local domain size-x           -nn_isize = ', -nn_isize 
     102            WRITE(numout,*) '      global domain size-x                 kpi = ', kpi 
     103         ENDIF 
     104         IF( nn_jsize > 0 ) THEN 
     105            WRITE(numout,*) '      global domain size-y            nn_jsize = ', nn_jsize 
     106         ELSE 
     107            WRITE(numout,*) '                                          jpnj = ', jpnj 
     108            WRITE(numout,*) '       local domain size-y           -nn_jsize = ', -nn_jsize 
     109            WRITE(numout,*) '      global domain size-y                 kpj = ', kpj 
     110         ENDIF 
     111         WRITE(numout,*) '      global domain size-z            nn_ksize = ', nn_ksize 
     112         WRITE(numout,*) '      LBC of the global domain          kperio = ', kperio 
    103113      ENDIF 
    104       IF( nn_jsize > 0 ) THEN 
    105          WRITE(numout,*) '      global domain size-y            nn_jsize = ', nn_jsize 
    106       ELSE 
    107          WRITE(numout,*) '                                          jpnj = ', jpnj 
    108          WRITE(numout,*) '       local domain size-y           -nn_jsize = ', -nn_jsize 
    109          WRITE(numout,*) '      global domain size-y                 kpj = ', kpj 
    110       ENDIF 
    111       WRITE(numout,*) '      global domain size-z            nn_ksize = ', nn_ksize 
    112       WRITE(numout,*) '      LBC of the global domain          kperio = ', kperio 
    113114      ! 
    114115   END SUBROUTINE usr_def_nam 
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/tests/CANAL/MY_SRC/usrdef_nam.F90

    r11317 r11320  
    119119      ! 
    120120      zh  = (kpk-1)*rn_dz 
    121       !                             ! control print 
    122       WRITE(numout,*) '   ' 
    123       WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
    124       WRITE(numout,*) '~~~~~~~~~~~ ' 
    125       WRITE(numout,*) '   Namelist namusr_def : EW_CANAL test case' 
    126       WRITE(numout,*) '      horizontal domain size-x          rn_domszx  = ', rn_domszx, ' km' 
    127       WRITE(numout,*) '      horizontal domain size-y          rn_domszy  = ', rn_domszy, ' km' 
    128       WRITE(numout,*) '      vertical   domain size-z          rn_domszz  = ', rn_domszz, '  m' 
    129       WRITE(numout,*) '      horizontal x-resolution           rn_dx      = ',     rn_dx, ' km' 
    130       WRITE(numout,*) '      horizontal y-resolution           rn_dy      = ',     rn_dy, ' km' 
    131       WRITE(numout,*) '      vertical resolution               rn_dz      = ',     rn_dz, '  m' 
    132       WRITE(numout,*) '      x-domain ratio of the 0           rn_0xratio = ', rn_0xratio 
    133       WRITE(numout,*) '      y-domain ratio of the 0           rn_0yratio = ', rn_0yratio 
    134       WRITE(numout,*) '          H [m] : ', zh 
    135       WRITE(numout,*) '      F computation                     nn_fcase   = ',   nn_fcase 
    136       WRITE(numout,*) '      Reference latitude                rn_ppgphi0 = ', rn_ppgphi0 
    137       WRITE(numout,*) '      10m wind speed                    rn_u10     = ',     rn_u10, ' m/s' 
    138       WRITE(numout,*) '         wind latitudinal extension     rn_windszy = ', rn_windszy, ' km' 
    139       WRITE(numout,*) '         wind longitudinal extension    rn_windszx = ', rn_windszx, ' km' 
    140       WRITE(numout,*) '         Uoce multiplicative factor     rn_uofac   = ',   rn_uofac 
    141       WRITE(numout,*) '      initial Canal max current        rn_vtxmax  = ',  rn_vtxmax, ' m/s' 
    142       WRITE(numout,*) '      initial zonal current             rn_uzonal  = ',  rn_uzonal, ' m/s' 
    143       WRITE(numout,*) '         Jet latitudinal extension      rn_ujetszy = ', rn_ujetszy, ' km' 
    144       WRITE(numout,*) '         Jet longitudinal extension     rn_ujetszx = ', rn_ujetszx, ' km' 
    145       WRITE(numout,*) '      bottom definition (0:flat)        nn_botcase = ', nn_botcase 
    146       WRITE(numout,*) '      initial condition case            nn_initcase= ', nn_initcase 
    147       WRITE(numout,*) '                   (0:rest, 1:zonal current, 10:shear)' 
    148       WRITE(numout,*) '      add random noise on initial ssh   ln_sshnoise= ', ln_sshnoise 
    149       WRITE(numout,*) '      Gaussian lambda parameter          rn_lambda = ', rn_lambda 
    150       ! 
    151121      !                             ! Set the lateral boundary condition of the global domain 
    152122      kperio = 1                    ! EW_CANAL configuration : closed basin 
    153       ! 
    154       WRITE(numout,*) '   ' 
    155       WRITE(numout,*) '   Lateral boundary condition of the global domain' 
    156       WRITE(numout,*) '      EW_CANAL : closed basin            jperio = ', kperio 
     123      !                             ! control print 
     124      IF(lwp) THEN 
     125         WRITE(numout,*) '   ' 
     126         WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
     127         WRITE(numout,*) '~~~~~~~~~~~ ' 
     128         WRITE(numout,*) '   Namelist namusr_def : EW_CANAL test case' 
     129         WRITE(numout,*) '      horizontal domain size-x          rn_domszx  = ', rn_domszx, ' km' 
     130         WRITE(numout,*) '      horizontal domain size-y          rn_domszy  = ', rn_domszy, ' km' 
     131         WRITE(numout,*) '      vertical   domain size-z          rn_domszz  = ', rn_domszz, '  m' 
     132         WRITE(numout,*) '      horizontal x-resolution           rn_dx      = ',     rn_dx, ' km' 
     133         WRITE(numout,*) '      horizontal y-resolution           rn_dy      = ',     rn_dy, ' km' 
     134         WRITE(numout,*) '      vertical resolution               rn_dz      = ',     rn_dz, '  m' 
     135         WRITE(numout,*) '      x-domain ratio of the 0           rn_0xratio = ', rn_0xratio 
     136         WRITE(numout,*) '      y-domain ratio of the 0           rn_0yratio = ', rn_0yratio 
     137         WRITE(numout,*) '          H [m] : ', zh 
     138         WRITE(numout,*) '      F computation                     nn_fcase   = ',   nn_fcase 
     139         WRITE(numout,*) '      Reference latitude                rn_ppgphi0 = ', rn_ppgphi0 
     140         WRITE(numout,*) '      10m wind speed                    rn_u10     = ',     rn_u10, ' m/s' 
     141         WRITE(numout,*) '         wind latitudinal extension     rn_windszy = ', rn_windszy, ' km' 
     142         WRITE(numout,*) '         wind longitudinal extension    rn_windszx = ', rn_windszx, ' km' 
     143         WRITE(numout,*) '         Uoce multiplicative factor     rn_uofac   = ',   rn_uofac 
     144         WRITE(numout,*) '      initial Canal max current         rn_vtxmax  = ',  rn_vtxmax, ' m/s' 
     145         WRITE(numout,*) '      initial zonal current             rn_uzonal  = ',  rn_uzonal, ' m/s' 
     146         WRITE(numout,*) '         Jet latitudinal extension      rn_ujetszy = ', rn_ujetszy, ' km' 
     147         WRITE(numout,*) '         Jet longitudinal extension     rn_ujetszx = ', rn_ujetszx, ' km' 
     148         WRITE(numout,*) '      bottom definition (0:flat)        nn_botcase = ', nn_botcase 
     149         WRITE(numout,*) '      initial condition case            nn_initcase= ', nn_initcase 
     150         WRITE(numout,*) '                   (0:rest, 1:zonal current, 10:shear)' 
     151         WRITE(numout,*) '      add random noise on initial ssh   ln_sshnoise= ', ln_sshnoise 
     152         WRITE(numout,*) '      Gaussian lambda parameter          rn_lambda = ', rn_lambda 
     153         WRITE(numout,*) '   ' 
     154         WRITE(numout,*) '   Lateral boundary condition of the global domain' 
     155         WRITE(numout,*) '      EW_CANAL : closed basin               jperio = ', kperio 
     156      ENDIF 
    157157      ! 
    158158   END SUBROUTINE usr_def_nam 
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/tests/ICE_ADV1D/MY_SRC/usrdef_nam.F90

    r11317 r11320  
    7979      zlx = kpi*rn_dx*1.e-3 
    8080      zly = kpj*rn_dy*1.e-3 
    81       !                             ! control print 
    82       WRITE(numout,*) '   ' 
    83       WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
    84       WRITE(numout,*) '~~~~~~~~~~~ ' 
    85       WRITE(numout,*) '   Namelist namusr_def : ICE_ADV1D test case' 
    86       WRITE(numout,*) '      horizontal resolution                    rn_dx  = ', rn_dx, ' meters' 
    87       WRITE(numout,*) '      horizontal resolution                    rn_dy  = ', rn_dy, ' meters' 
    88       WRITE(numout,*) '      ICE_ADV1D domain  ' 
    89       WRITE(numout,*) '         LX [km]: ', zlx 
    90       WRITE(numout,*) '         LY [km]: ', zly 
    91       WRITE(numout,*) '         resulting global domain size :        jpiglo = ', kpi 
    92       WRITE(numout,*) '                                               jpjglo = ', kpj 
    93       WRITE(numout,*) '                                               jpkglo = ', kpk 
    94       WRITE(numout,*) '         Coriolis:', ln_corio 
    95       ! 
    9681      !                             ! Set the lateral boundary condition of the global domain 
    9782      kperio = 0                    ! ICE_ADV1D configuration : bi-periodic basin 
    98       ! 
    99       WRITE(numout,*) '   ' 
    100       WRITE(numout,*) '   Lateral boundary condition of the global domain' 
    101       WRITE(numout,*) '      ICE_ADV1D : closed basin                    jperio = ', kperio 
     83      !                             ! control print 
     84      IF(lwp) THEN 
     85         WRITE(numout,*) '   ' 
     86         WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
     87         WRITE(numout,*) '~~~~~~~~~~~ ' 
     88         WRITE(numout,*) '   Namelist namusr_def : ICE_ADV1D test case' 
     89         WRITE(numout,*) '      horizontal resolution                    rn_dx  = ', rn_dx, ' meters' 
     90         WRITE(numout,*) '      horizontal resolution                    rn_dy  = ', rn_dy, ' meters' 
     91         WRITE(numout,*) '      ICE_ADV1D domain  ' 
     92         WRITE(numout,*) '         LX [km]: ', zlx 
     93         WRITE(numout,*) '         LY [km]: ', zly 
     94         WRITE(numout,*) '         resulting global domain size :        jpiglo = ', kpi 
     95         WRITE(numout,*) '                                               jpjglo = ', kpj 
     96         WRITE(numout,*) '                                               jpkglo = ', kpk 
     97         WRITE(numout,*) '         Coriolis:', ln_corio 
     98         WRITE(numout,*) '   ' 
     99         WRITE(numout,*) '   Lateral boundary condition of the global domain' 
     100         WRITE(numout,*) '      ICE_ADV1D : closed basin                 jperio = ', kperio 
     101      ENDIF 
    102102      ! 
    103103   END SUBROUTINE usr_def_nam 
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/tests/ICE_ADV2D/MY_SRC/usrdef_nam.F90

    r11317 r11320  
    8181      ! 
    8282      cd_cfg = 'ICE_ADV2D'           ! name & resolution (not used) 
    83       kk_cfg = INT( rn_dx ) 
     83      kk_cfg = NINT( rn_dx ) 
    8484      ! 
    85       ! Global Domain size:  ICE_ADV2D domain is  300 km x 300 Km x 10 m 
    86       kpi = INT( 300.e3 / rn_dx ) -1 
    87       kpj = INT( 300.e3 / rn_dy ) -1 
    88 #if defined key_agrif 
    89       IF( .NOT. Agrif_Root() ) THEN 
     85      IF( Agrif_Root() ) THEN        ! Global Domain size:  ICE_AGRIF domain is  300 km x 300 Km x 10 m 
     86         kpi = NINT( 300.e3 / rn_dx ) - 1 
     87         kpj = NINT( 300.e3 / rn_dy ) - 1 
     88      ELSE 
    9089         kpi = nbcellsx + 2 + 2*nbghostcells 
    9190         kpj = nbcellsy + 2 + 2*nbghostcells 
    9291      ENDIF 
    93 #endif 
    9492      kpk = 1 
    9593      ! 
     
    9896      zlx = kpi*rn_dx*1.e-3 
    9997      zly = kpj*rn_dy*1.e-3 
     98      ! 
     99      IF( Agrif_Root() ) THEN   ;   kperio = 7     ! ICE_AGRIF configuration : bi-periodic basin 
     100      ELSE                      ;   kperio = 0     ! closed periodicity for the zoom 
     101      ENDIF 
    100102      !                             ! control print 
    101       WRITE(numout,*) '   ' 
    102       WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
    103       WRITE(numout,*) '~~~~~~~~~~~ ' 
    104       WRITE(numout,*) '   Namelist namusr_def : ICE_ADV2D test case' 
    105       WRITE(numout,*) '      horizontal resolution                    rn_dx  = ', rn_dx, ' meters' 
    106       WRITE(numout,*) '      horizontal resolution                    rn_dy  = ', rn_dy, ' meters' 
    107       WRITE(numout,*) '      ICE_ADV2D domain = 300 km x 300Km x 1 grid-point ' 
    108       WRITE(numout,*) '         LX [km]: ', zlx 
    109       WRITE(numout,*) '         LY [km]: ', zly 
    110       WRITE(numout,*) '         resulting global domain size :        jpiglo = ', kpi 
    111       WRITE(numout,*) '                                               jpjglo = ', kpj 
    112       WRITE(numout,*) '                                               jpkglo = ', kpk 
    113       WRITE(numout,*) '         Coriolis:', ln_corio 
    114       ! 
    115       !                             ! Set the lateral boundary condition of the global domain 
    116       kperio = 7                    ! ICE_ADV2D configuration : bi-periodic basin 
    117 #if defined key_agrif 
    118       IF( .NOT. Agrif_Root() ) THEN 
    119       kperio = 0 
     103      IF(lwp) THEN 
     104         WRITE(numout,*) '   ' 
     105         WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
     106         WRITE(numout,*) '~~~~~~~~~~~ ' 
     107         WRITE(numout,*) '   Namelist namusr_def : ICE_ADV2D test case' 
     108         WRITE(numout,*) '      horizontal resolution                    rn_dx  = ', rn_dx, ' meters' 
     109         WRITE(numout,*) '      horizontal resolution                    rn_dy  = ', rn_dy, ' meters' 
     110         WRITE(numout,*) '      ICE_ADV2D domain = 300 km x 300Km x 1 grid-point ' 
     111         WRITE(numout,*) '         LX [km]: ', zlx 
     112         WRITE(numout,*) '         LY [km]: ', zly 
     113         WRITE(numout,*) '         resulting global domain size :        jpiglo = ', kpi 
     114         WRITE(numout,*) '                                               jpjglo = ', kpj 
     115         WRITE(numout,*) '                                               jpkglo = ', kpk 
     116         WRITE(numout,*) '         Coriolis:', ln_corio 
     117         WRITE(numout,*) '   ' 
     118         WRITE(numout,*) '   Lateral boundary condition of the global domain' 
     119         WRITE(numout,*) '      ICE_ADV2D : bi-periodic basin            jperio = ', kperio 
    120120      ENDIF 
    121 #endif 
    122       ! 
    123       WRITE(numout,*) '   ' 
    124       WRITE(numout,*) '   Lateral boundary condition of the global domain' 
    125       WRITE(numout,*) '      ICE_ADV2D : bi-periodic basin               jperio = ', kperio 
    126121      ! 
    127122   END SUBROUTINE usr_def_nam 
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/tests/ICE_AGRIF/MY_SRC/usrdef_nam.F90

    r11317 r11320  
    8181      ! 
    8282      cd_cfg = 'ICE_AGRIF'           ! name & resolution (not used) 
    83       kk_cfg = INT( rn_dx ) 
     83      kk_cfg = NINT( rn_dx ) 
    8484      ! 
    85       ! Global Domain size:  ICE_AGRIF domain is  300 km x 300 Km x 10 m 
    86       kpi = INT( 300.e3 / rn_dx ) -1 
    87       kpj = INT( 300.e3 / rn_dy ) -1 
    88 #if defined key_agrif 
    89       IF( .NOT. Agrif_Root() ) THEN 
     85      IF( Agrif_Root() ) THEN        ! Global Domain size:  ICE_AGRIF domain is  300 km x 300 Km x 10 m 
     86         kpi = NINT( 300.e3 / rn_dx ) - 1 
     87         kpj = NINT( 300.e3 / rn_dy ) - 1 
     88      ELSE 
    9089         kpi = nbcellsx + 2 + 2*nbghostcells 
    9190         kpj = nbcellsy + 2 + 2*nbghostcells 
    9291      ENDIF 
    93 #endif 
    9492      kpk = 1 
    9593      ! 
     
    9896      zlx = kpi*rn_dx*1.e-3 
    9997      zly = kpj*rn_dy*1.e-3 
     98      ! 
     99      IF( Agrif_Root() ) THEN   ;   kperio = 7     ! ICE_AGRIF configuration : bi-periodic basin 
     100      ELSE                      ;   kperio = 0     ! closed periodicity for the zoom 
     101      ENDIF 
    100102      !                             ! control print 
    101       WRITE(numout,*) '   ' 
    102       WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
    103       WRITE(numout,*) '~~~~~~~~~~~ ' 
    104       WRITE(numout,*) '   Namelist namusr_def : ICE_AGRIF test case' 
    105       WRITE(numout,*) '      horizontal resolution                    rn_dx  = ', rn_dx, ' meters' 
    106       WRITE(numout,*) '      horizontal resolution                    rn_dy  = ', rn_dy, ' meters' 
    107       WRITE(numout,*) '      ICE_AGRIF domain = 300 km x 300Km x 1 grid-point ' 
    108       WRITE(numout,*) '         LX [km]: ', zlx 
    109       WRITE(numout,*) '         LY [km]: ', zly 
    110       WRITE(numout,*) '         resulting global domain size :        jpiglo = ', kpi 
    111       WRITE(numout,*) '                                               jpjglo = ', kpj 
    112       WRITE(numout,*) '                                               jpkglo = ', kpk 
    113       WRITE(numout,*) '         Coriolis:', ln_corio 
    114       ! 
    115       !                             ! Set the lateral boundary condition of the global domain 
    116       kperio = 7                    ! ICE_AGRIF configuration : bi-periodic basin 
    117 #if defined key_agrif 
    118       IF( .NOT. Agrif_Root() ) THEN 
    119       kperio = 0 
     103      IF(lwp) THEN 
     104         WRITE(numout,*) '   ' 
     105         WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
     106         WRITE(numout,*) '~~~~~~~~~~~ ' 
     107         WRITE(numout,*) '   Namelist namusr_def : ICE_AGRIF test case' 
     108         WRITE(numout,*) '      horizontal resolution                    rn_dx  = ', rn_dx, ' meters' 
     109         WRITE(numout,*) '      horizontal resolution                    rn_dy  = ', rn_dy, ' meters' 
     110         WRITE(numout,*) '      ICE_AGRIF domain = 300 km x 300Km x 1 grid-point ' 
     111         WRITE(numout,*) '         LX [km]: ', zlx 
     112         WRITE(numout,*) '         LY [km]: ', zly 
     113         WRITE(numout,*) '         resulting global domain size :        jpiglo = ', kpi 
     114         WRITE(numout,*) '                                               jpjglo = ', kpj 
     115         WRITE(numout,*) '                                               jpkglo = ', kpk 
     116         WRITE(numout,*) '         Coriolis:', ln_corio 
     117         WRITE(numout,*) '   ' 
     118         WRITE(numout,*) '   Lateral boundary condition of the global domain' 
     119         WRITE(numout,*) '      ICE_AGRIF : bi-periodic basin            jperio = ', kperio 
    120120      ENDIF 
    121 #endif 
    122       ! 
    123       WRITE(numout,*) '   ' 
    124       WRITE(numout,*) '   Lateral boundary condition of the global domain' 
    125       WRITE(numout,*) '      ICE_AGRIF : bi-periodic basin               jperio = ', kperio 
    126121      ! 
    127122   END SUBROUTINE usr_def_nam 
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/tests/ISOMIP/MY_SRC/usrdef_nam.F90

    r11317 r11320  
    7878      kpk = INT( rbathy / rn_e3    ) + 1     ! add 1 for t-point in the seafloor 
    7979      ! 
    80       !                             ! control print 
    81       WRITE(numout,*) '   ' 
    82       WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
    83       WRITE(numout,*) '~~~~~~~~~~~ ' 
    84       WRITE(numout,*) '   Namelist namusr_def : ISOMIP test case' 
    85       WRITE(numout,*) '      type of vertical coordinate : ' 
    86       WRITE(numout,*) '         z-coordinate flag                     ln_zco   = ', ln_zco 
    87       WRITE(numout,*) '         z-partial-step coordinate flag        ln_zps   = ', ln_zps 
    88       WRITE(numout,*) '         s-coordinate flag                     ln_sco   = ', ln_sco 
    89       WRITE(numout,*) '      resolution' 
    90       WRITE(numout,*) '         zonal      resolution                 rn_e1deg = ', rn_e1deg, ' degrees' 
    91       WRITE(numout,*) '         meridional resolution                 rn_e1deg = ', rn_e1deg, ' degrees' 
    92       WRITE(numout,*) '         vertical   resolution                 rn_e3    = ', rn_e3   , ' meters' 
    93       WRITE(numout,*) '      ISOMIP domain = 15° x 10° x 900 m' 
    94       WRITE(numout,*) '         resulting global domain size :        jpiglo   = ', kpi 
    95       WRITE(numout,*) '                                               jpjglo   = ', kpj 
    96       WRITE(numout,*) '                                               jpkglo   = ', kpk 
    97       ! 
    98       ! 
    9980      !                             ! Set the lateral boundary condition of the global domain 
    10081      kperio = 0                    ! ISOMIP configuration : close basin 
    10182      ! 
    102       WRITE(numout,*) '   ' 
    103       WRITE(numout,*) '   Lateral boundary condition of the global domain' 
    104       WRITE(numout,*) '      ISOMIP : closed basin                    jperio   = ', kperio 
     83      !                             ! control print 
     84      IF(lwp) THEN 
     85         WRITE(numout,*) '   ' 
     86         WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
     87         WRITE(numout,*) '~~~~~~~~~~~ ' 
     88         WRITE(numout,*) '   Namelist namusr_def : ISOMIP test case' 
     89         WRITE(numout,*) '      type of vertical coordinate : ' 
     90         WRITE(numout,*) '         z-coordinate flag                     ln_zco   = ', ln_zco 
     91         WRITE(numout,*) '         z-partial-step coordinate flag        ln_zps   = ', ln_zps 
     92         WRITE(numout,*) '         s-coordinate flag                     ln_sco   = ', ln_sco 
     93         WRITE(numout,*) '      resolution' 
     94         WRITE(numout,*) '         zonal      resolution                 rn_e1deg = ', rn_e1deg, ' degrees' 
     95         WRITE(numout,*) '         meridional resolution                 rn_e1deg = ', rn_e1deg, ' degrees' 
     96         WRITE(numout,*) '         vertical   resolution                 rn_e3    = ', rn_e3   , ' meters' 
     97         WRITE(numout,*) '      ISOMIP domain = 15° x 10° x 900 m' 
     98         WRITE(numout,*) '         resulting global domain size :        jpiglo   = ', kpi 
     99         WRITE(numout,*) '                                               jpjglo   = ', kpj 
     100         WRITE(numout,*) '                                               jpkglo   = ', kpk 
     101         WRITE(numout,*) '   ' 
     102         WRITE(numout,*) '   Lateral boundary condition of the global domain' 
     103         WRITE(numout,*) '      ISOMIP : closed basin                    jperio   = ', kperio 
     104      ENDIF 
    105105      ! 
    106106   END SUBROUTINE usr_def_nam 
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/tests/LOCK_EXCHANGE/MY_SRC/usrdef_nam.F90

    r11317 r11320  
    7474      kpj = 3 
    7575      kpk = INT(  20.  / rn_dz ) + 1 
    76       ! 
    77       !                             ! control print 
    78       WRITE(numout,*) '   ' 
    79       WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
    80       WRITE(numout,*) '~~~~~~~~~~~ ' 
    81       WRITE(numout,*) '   Namelist namusr_def : LOCK_EXCHANGE test case' 
    82       WRITE(numout,*) '      horizontal resolution                    rn_dx  = ', rn_dx, ' meters' 
    83       WRITE(numout,*) '      vertical   resolution                    rn_dz  = ', rn_dz, ' meters' 
    84       WRITE(numout,*) '      LOCK_EXCHANGE domain = 64 km  x  3 grid-points  x  20 m' 
    85       WRITE(numout,*) '         resulting global domain size :        jpiglo = ', kpi 
    86       WRITE(numout,*) '                                               jpjglo = ', kpj 
    87       WRITE(numout,*) '                                               jpkglo = ', kpk 
    88       ! 
    8976      !                             ! Set the lateral boundary condition of the global domain 
    9077      kperio = 0                    ! LOCK_EXCHANGE configuration : closed domain 
    9178      ! 
    92       WRITE(numout,*) '   ' 
    93       WRITE(numout,*) '   Lateral boundary condition of the global domain' 
    94       WRITE(numout,*) '      closed                                   jperio = ', kperio 
     79      !                             ! control print 
     80      IF(lwp) THEN 
     81         WRITE(numout,*) '   ' 
     82         WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
     83         WRITE(numout,*) '~~~~~~~~~~~ ' 
     84         WRITE(numout,*) '   Namelist namusr_def : LOCK_EXCHANGE test case' 
     85         WRITE(numout,*) '      horizontal resolution                    rn_dx  = ', rn_dx, ' meters' 
     86         WRITE(numout,*) '      vertical   resolution                    rn_dz  = ', rn_dz, ' meters' 
     87         WRITE(numout,*) '      LOCK_EXCHANGE domain = 64 km  x  3 grid-points  x  20 m' 
     88         WRITE(numout,*) '         resulting global domain size :        jpiglo = ', kpi 
     89         WRITE(numout,*) '                                               jpjglo = ', kpj 
     90         WRITE(numout,*) '                                               jpkglo = ', kpk 
     91         WRITE(numout,*) '   ' 
     92         WRITE(numout,*) '   Lateral boundary condition of the global domain' 
     93         WRITE(numout,*) '      closed                                   jperio = ', kperio 
     94      ENDIF 
    9595      ! 
    9696   END SUBROUTINE usr_def_nam 
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/tests/VORTEX/MY_SRC/usrdef_nam.F90

    r11317 r11320  
    8080      ! 
    8181      cd_cfg = 'VORTEX'             ! name & resolution (not used) 
    82       kk_cfg = INT( rn_dx ) 
     82      kk_cfg = nINT( rn_dx ) 
    8383      ! 
    84       ! Global Domain size:  VORTEX global domain is  1800 km x 1800 Km x 5000 m 
    85       kpi = INT( 1800.e3  / rn_dx ) + 3   
    86       kpj = INT( 1800.e3  / rn_dy ) + 3  
    87       kpk = INT( 5000._wp / rn_dz ) + 1 
    88 #if defined key_agrif 
    89       IF( .NOT. Agrif_Root() ) THEN 
     84      IF( Agrif_Root() ) THEN       ! Global Domain size:  VORTEX global domain is  1800 km x 1800 Km x 5000 m 
     85         kpi = NINT( 1800.e3  / rn_dx ) + 3   
     86         kpj = NINT( 1800.e3  / rn_dy ) + 3  
     87      ELSE 
    9088         kpi  = nbcellsx + 2 + 2*nbghostcells 
    9189         kpj  = nbcellsy + 2 + 2*nbghostcells 
    9290      ENDIF 
    93 #endif 
     91      kpk = NINT( 5000._wp / rn_dz ) + 1 
    9492      ! 
    9593      zlx = (kpi-2)*rn_dx*1.e-3 
    9694      zly = (kpj-2)*rn_dy*1.e-3 
    9795      zh  = (kpk-1)*rn_dz 
    98       !                             ! control print 
    99       WRITE(numout,*) '   ' 
    100       WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
    101       WRITE(numout,*) '~~~~~~~~~~~ ' 
    102       WRITE(numout,*) '   Namelist namusr_def : VORTEX test case' 
    103       WRITE(numout,*) '      horizontal resolution             rn_dx  = ', rn_dx, ' m' 
    104       WRITE(numout,*) '      horizontal resolution             rn_dy  = ', rn_dy, ' m' 
    105       WRITE(numout,*) '      vertical resolution               rn_dz  = ', rn_dz, ' m' 
    106       WRITE(numout,*) '      VORTEX domain: ' 
    107       WRITE(numout,*) '         LX [km]: ', zlx 
    108       WRITE(numout,*) '         LY [km]: ', zly 
    109       WRITE(numout,*) '          H [m] : ', zh 
    110       WRITE(numout,*) '      Reference latitude            rn_ppgphi0 = ', rn_ppgphi0 
    111       ! 
    11296      !                             ! Set the lateral boundary condition of the global domain 
    11397      kperio = 0                    ! VORTEX configuration : closed basin 
    114       ! 
    115       WRITE(numout,*) '   ' 
    116       WRITE(numout,*) '   Lateral boundary condition of the global domain' 
    117       WRITE(numout,*) '      VORTEX : closed basin            jperio = ', kperio 
     98      !                             ! control print 
     99      IF(lwp) THEN 
     100         WRITE(numout,*) '   ' 
     101         WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
     102         WRITE(numout,*) '~~~~~~~~~~~ ' 
     103         WRITE(numout,*) '   Namelist namusr_def : VORTEX test case' 
     104         WRITE(numout,*) '      horizontal resolution             rn_dx  = ', rn_dx, ' m' 
     105         WRITE(numout,*) '      horizontal resolution             rn_dy  = ', rn_dy, ' m' 
     106         WRITE(numout,*) '      vertical resolution               rn_dz  = ', rn_dz, ' m' 
     107         WRITE(numout,*) '      VORTEX domain: ' 
     108         WRITE(numout,*) '         LX [km]: ', zlx 
     109         WRITE(numout,*) '         LY [km]: ', zly 
     110         WRITE(numout,*) '          H [m] : ', zh 
     111         WRITE(numout,*) '      Reference latitude            rn_ppgphi0 = ', rn_ppgphi0 
     112         WRITE(numout,*) '   ' 
     113         WRITE(numout,*) '   Lateral boundary condition of the global domain' 
     114         WRITE(numout,*) '      VORTEX : closed basin            jperio = ', kperio 
     115      ENDIF 
    118116      ! 
    119117   END SUBROUTINE usr_def_nam 
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/tests/WAD/MY_SRC/usrdef_nam.F90

    r11317 r11320  
    7676      kpj = INT(  32.e3 / rn_dx ) + 2 
    7777      kpk = INT(  10.  / rn_dz ) + 1 
    78       ! 
    79       !                             ! control print 
    80       WRITE(numout,*) '   ' 
    81       WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
    82       WRITE(numout,*) '~~~~~~~~~~~ ' 
    83       WRITE(numout,*) '   Namelist namusr_def : WAD_TEST_CASES test case' 
    84       WRITE(numout,*) '      horizontal resolution                    rn_dx  = ', rn_dx, ' meters' 
    85       WRITE(numout,*) '      vertical   resolution                    rn_dz  = ', rn_dz, ' meters' 
    86       WRITE(numout,*) '      WAD_TEST_CASES domain = 52 km  x  34 km x 10 m' 
    87       WRITE(numout,*) '         resulting global domain size :        jpiglo = ', kpi 
    88       WRITE(numout,*) '                                               jpjglo = ', kpj 
    89       WRITE(numout,*) '                                               jpkglo = ', kpk 
    90       ! 
    9178      !                             ! Set the lateral boundary condition of the global domain 
    9279      kperio = 0                    ! WAD_TEST_CASES configuration : closed domain 
    9380      IF( nn_wad_test == 8 ) kperio = 7 ! North-South cyclic test 
    9481      ! 
    95       WRITE(numout,*) '   ' 
    96       WRITE(numout,*) '   Lateral boundary condition of the global domain' 
    97       WRITE(numout,*) '      closed                                   jperio = ', kperio 
     82      !                             ! control print 
     83      IF(lwp) THEN 
     84         WRITE(numout,*) '   ' 
     85         WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
     86         WRITE(numout,*) '~~~~~~~~~~~ ' 
     87         WRITE(numout,*) '   Namelist namusr_def : WAD_TEST_CASES test case' 
     88         WRITE(numout,*) '      horizontal resolution                    rn_dx  = ', rn_dx, ' meters' 
     89         WRITE(numout,*) '      vertical   resolution                    rn_dz  = ', rn_dz, ' meters' 
     90         WRITE(numout,*) '      WAD_TEST_CASES domain = 52 km  x  34 km x 10 m' 
     91         WRITE(numout,*) '         resulting global domain size :        jpiglo = ', kpi 
     92         WRITE(numout,*) '                                               jpjglo = ', kpj 
     93         WRITE(numout,*) '                                               jpkglo = ', kpk 
     94         WRITE(numout,*) '   ' 
     95         WRITE(numout,*) '   Lateral boundary condition of the global domain' 
     96         WRITE(numout,*) '      closed                                   jperio = ', kperio 
     97      ENDIF 
    9898      ! 
    9999   END SUBROUTINE usr_def_nam 
Note: See TracChangeset for help on using the changeset viewer.