- Timestamp:
- 2019-07-22T14:52:29+02:00 (5 years ago)
- Location:
- NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization
- Files:
-
- 12 edited
-
cfgs/C1D_PAPA/MY_SRC/usrdef_nam.F90 (modified) (1 diff)
-
src/OCE/DOM/domain.F90 (modified) (3 diffs)
-
src/OCE/USR/usrdef_nam.F90 (modified) (1 diff)
-
tests/BENCH/MY_SRC/usrdef_nam.F90 (modified) (2 diffs)
-
tests/CANAL/MY_SRC/usrdef_nam.F90 (modified) (1 diff)
-
tests/ICE_ADV1D/MY_SRC/usrdef_nam.F90 (modified) (1 diff)
-
tests/ICE_ADV2D/MY_SRC/usrdef_nam.F90 (modified) (2 diffs)
-
tests/ICE_AGRIF/MY_SRC/usrdef_nam.F90 (modified) (2 diffs)
-
tests/ISOMIP/MY_SRC/usrdef_nam.F90 (modified) (1 diff)
-
tests/LOCK_EXCHANGE/MY_SRC/usrdef_nam.F90 (modified) (1 diff)
-
tests/VORTEX/MY_SRC/usrdef_nam.F90 (modified) (1 diff)
-
tests/WAD/MY_SRC/usrdef_nam.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/cfgs/C1D_PAPA/MY_SRC/usrdef_nam.F90
r11317 r11320 74 74 kpj = 3 75 75 kpk = 75 76 !77 ! ! control print78 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_zco84 WRITE(numout,*) ' z-partial-step coordinate flag ln_zps = ', ln_zps85 WRITE(numout,*) ' s-coordinate flag ln_sco = ', ln_sco86 WRITE(numout,*) ' C1D domain = 3 x 3 x75 grid-points '87 WRITE(numout,*) ' resulting global domain size : jpiglo = ', kpi88 WRITE(numout,*) ' jpjglo = ', kpj89 WRITE(numout,*) ' jpkglo = ', kpk90 91 !92 76 ! ! Set the lateral boundary condition of the global domain 93 77 kperio = 7 ! C1D configuration : 3x3 basin with cyclic Est-West and Norht-South condition 94 78 ! 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 97 96 ! 98 97 END SUBROUTINE usr_def_nam -
NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/OCE/DOM/domain.F90
r11317 r11320 530 530 !!---------------------------------------------------------------------- 531 531 ! 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 535 537 ! 536 538 CALL iom_open( cn_domcfg, inum ) … … 543 545 CALL iom_get( inum, 'ORCA_index', zorca_res ) ; kk_cfg = NINT( zorca_res ) 544 546 ! 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 548 552 ! 549 553 ELSE !- cd_cfg & k_cfg are not used … … 566 570 CALL iom_close( inum ) 567 571 ! 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 573 579 ! 574 580 END SUBROUTINE domain_cfg -
NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/OCE/USR/usrdef_nam.F90
r11317 r11320 80 80 #endif 81 81 kpk = jpkglo 82 !83 ! ! control print84 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_bench89 WRITE(numout,*) ' inverse resolution & implied domain size nn_GYRE = ', nn_GYRE90 #if defined key_agrif91 IF( Agrif_Root() ) THEN92 #endif93 WRITE(numout,*) ' jpiglo = 30*nn_GYRE+2 jpiglo = ', kpi94 WRITE(numout,*) ' jpjglo = 20*nn_GYRE+2 jpjglo = ', kpj95 #if defined key_agrif96 ENDIF97 #endif98 WRITE(numout,*) ' number of model levels jpkglo = ', kpk99 !100 82 ! ! Set the lateral boundary condition of the global domain 101 83 kperio = 0 ! GYRE configuration : closed domain 102 84 ! 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 105 105 ! 106 106 END SUBROUTINE usr_def_nam -
NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/tests/BENCH/MY_SRC/usrdef_nam.F90
r11317 r11320 69 69 cd_cfg = 'BENCH' ! name & resolution (not used) 70 70 kk_cfg = 0 71 71 ! 72 72 IF( nn_isize < 0 .AND. nn_jsize < 0 ) THEN 73 73 ! … … 86 86 kpj = nn_jsize 87 87 ENDIF 88 88 ! 89 89 kpk = nn_ksize 90 90 kperio = nn_perio 91 92 91 ! ! 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 103 113 ENDIF 104 IF( nn_jsize > 0 ) THEN105 WRITE(numout,*) ' global domain size-y nn_jsize = ', nn_jsize106 ELSE107 WRITE(numout,*) ' jpnj = ', jpnj108 WRITE(numout,*) ' local domain size-y -nn_jsize = ', -nn_jsize109 WRITE(numout,*) ' global domain size-y kpj = ', kpj110 ENDIF111 WRITE(numout,*) ' global domain size-z nn_ksize = ', nn_ksize112 WRITE(numout,*) ' LBC of the global domain kperio = ', kperio113 114 ! 114 115 END SUBROUTINE usr_def_nam -
NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/tests/CANAL/MY_SRC/usrdef_nam.F90
r11317 r11320 119 119 ! 120 120 zh = (kpk-1)*rn_dz 121 ! ! control print122 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_0xratio133 WRITE(numout,*) ' y-domain ratio of the 0 rn_0yratio = ', rn_0yratio134 WRITE(numout,*) ' H [m] : ', zh135 WRITE(numout,*) ' F computation nn_fcase = ', nn_fcase136 WRITE(numout,*) ' Reference latitude rn_ppgphi0 = ', rn_ppgphi0137 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_uofac141 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_botcase146 WRITE(numout,*) ' initial condition case nn_initcase= ', nn_initcase147 WRITE(numout,*) ' (0:rest, 1:zonal current, 10:shear)'148 WRITE(numout,*) ' add random noise on initial ssh ln_sshnoise= ', ln_sshnoise149 WRITE(numout,*) ' Gaussian lambda parameter rn_lambda = ', rn_lambda150 !151 121 ! ! Set the lateral boundary condition of the global domain 152 122 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 157 157 ! 158 158 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 79 79 zlx = kpi*rn_dx*1.e-3 80 80 zly = kpj*rn_dy*1.e-3 81 ! ! control print82 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]: ', zlx90 WRITE(numout,*) ' LY [km]: ', zly91 WRITE(numout,*) ' resulting global domain size : jpiglo = ', kpi92 WRITE(numout,*) ' jpjglo = ', kpj93 WRITE(numout,*) ' jpkglo = ', kpk94 WRITE(numout,*) ' Coriolis:', ln_corio95 !96 81 ! ! Set the lateral boundary condition of the global domain 97 82 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 102 102 ! 103 103 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 81 81 ! 82 82 cd_cfg = 'ICE_ADV2D' ! name & resolution (not used) 83 kk_cfg = INT( rn_dx )83 kk_cfg = NINT( rn_dx ) 84 84 ! 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 90 89 kpi = nbcellsx + 2 + 2*nbghostcells 91 90 kpj = nbcellsy + 2 + 2*nbghostcells 92 91 ENDIF 93 #endif94 92 kpk = 1 95 93 ! … … 98 96 zlx = kpi*rn_dx*1.e-3 99 97 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 100 102 ! ! 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 120 120 ENDIF 121 #endif122 !123 WRITE(numout,*) ' '124 WRITE(numout,*) ' Lateral boundary condition of the global domain'125 WRITE(numout,*) ' ICE_ADV2D : bi-periodic basin jperio = ', kperio126 121 ! 127 122 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 81 81 ! 82 82 cd_cfg = 'ICE_AGRIF' ! name & resolution (not used) 83 kk_cfg = INT( rn_dx )83 kk_cfg = NINT( rn_dx ) 84 84 ! 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 90 89 kpi = nbcellsx + 2 + 2*nbghostcells 91 90 kpj = nbcellsy + 2 + 2*nbghostcells 92 91 ENDIF 93 #endif94 92 kpk = 1 95 93 ! … … 98 96 zlx = kpi*rn_dx*1.e-3 99 97 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 100 102 ! ! 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 120 120 ENDIF 121 #endif122 !123 WRITE(numout,*) ' '124 WRITE(numout,*) ' Lateral boundary condition of the global domain'125 WRITE(numout,*) ' ICE_AGRIF : bi-periodic basin jperio = ', kperio126 121 ! 127 122 END SUBROUTINE usr_def_nam -
NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/tests/ISOMIP/MY_SRC/usrdef_nam.F90
r11317 r11320 78 78 kpk = INT( rbathy / rn_e3 ) + 1 ! add 1 for t-point in the seafloor 79 79 ! 80 ! ! control print81 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_zco87 WRITE(numout,*) ' z-partial-step coordinate flag ln_zps = ', ln_zps88 WRITE(numout,*) ' s-coordinate flag ln_sco = ', ln_sco89 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 = ', kpi95 WRITE(numout,*) ' jpjglo = ', kpj96 WRITE(numout,*) ' jpkglo = ', kpk97 !98 !99 80 ! ! Set the lateral boundary condition of the global domain 100 81 kperio = 0 ! ISOMIP configuration : close basin 101 82 ! 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 105 105 ! 106 106 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 74 74 kpj = 3 75 75 kpk = INT( 20. / rn_dz ) + 1 76 !77 ! ! control print78 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 = ', kpi86 WRITE(numout,*) ' jpjglo = ', kpj87 WRITE(numout,*) ' jpkglo = ', kpk88 !89 76 ! ! Set the lateral boundary condition of the global domain 90 77 kperio = 0 ! LOCK_EXCHANGE configuration : closed domain 91 78 ! 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 95 95 ! 96 96 END SUBROUTINE usr_def_nam -
NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/tests/VORTEX/MY_SRC/usrdef_nam.F90
r11317 r11320 80 80 ! 81 81 cd_cfg = 'VORTEX' ! name & resolution (not used) 82 kk_cfg = INT( rn_dx )82 kk_cfg = nINT( rn_dx ) 83 83 ! 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 90 88 kpi = nbcellsx + 2 + 2*nbghostcells 91 89 kpj = nbcellsy + 2 + 2*nbghostcells 92 90 ENDIF 93 #endif 91 kpk = NINT( 5000._wp / rn_dz ) + 1 94 92 ! 95 93 zlx = (kpi-2)*rn_dx*1.e-3 96 94 zly = (kpj-2)*rn_dy*1.e-3 97 95 zh = (kpk-1)*rn_dz 98 ! ! control print99 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]: ', zlx108 WRITE(numout,*) ' LY [km]: ', zly109 WRITE(numout,*) ' H [m] : ', zh110 WRITE(numout,*) ' Reference latitude rn_ppgphi0 = ', rn_ppgphi0111 !112 96 ! ! Set the lateral boundary condition of the global domain 113 97 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 118 116 ! 119 117 END SUBROUTINE usr_def_nam -
NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/tests/WAD/MY_SRC/usrdef_nam.F90
r11317 r11320 76 76 kpj = INT( 32.e3 / rn_dx ) + 2 77 77 kpk = INT( 10. / rn_dz ) + 1 78 !79 ! ! control print80 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 = ', kpi88 WRITE(numout,*) ' jpjglo = ', kpj89 WRITE(numout,*) ' jpkglo = ', kpk90 !91 78 ! ! Set the lateral boundary condition of the global domain 92 79 kperio = 0 ! WAD_TEST_CASES configuration : closed domain 93 80 IF( nn_wad_test == 8 ) kperio = 7 ! North-South cyclic test 94 81 ! 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 98 98 ! 99 99 END SUBROUTINE usr_def_nam
Note: See TracChangeset
for help on using the changeset viewer.
