Changeset 15264 for NEMO/branches/2020/ticket2487/tests
- Timestamp:
- 2021-09-16T12:42:20+02:00 (3 years ago)
- Location:
- NEMO/branches/2020/ticket2487
- Files:
-
- 54 edited
- 1 copied
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/ticket2487
- Property svn:externals
-
old new 1 ^/utils/build/arch@12130 arch2 ^/utils/build/makenemo@12191 makenemo3 ^/utils/build/mk@11662 mk4 ^/utils/tools_r4.0-HEAD@1 2672tools5 ^/vendors/AGRIF/ dev@10586ext/AGRIF6 ^/vendors/FCM@10134 ext/FCM7 ^/vendors/IOIPSL@9655 ext/IOIPSL1 ^/utils/build/arch@12130 arch 2 ^/utils/build/makenemo@12191 makenemo 3 ^/utils/build/mk@11662 mk 4 ^/utils/tools_r4.0-HEAD@14974 tools 5 ^/vendors/AGRIF/stable@14105 ext/AGRIF 6 ^/vendors/FCM@10134 ext/FCM 7 ^/vendors/IOIPSL@9655 ext/IOIPSL 8 8 9 9 # SETTE mapping (inactive)
-
- Property svn:externals
-
NEMO/branches/2020/ticket2487/tests/BENCH/EXPREF/namelist_cfg_orca025_like
r11536 r15264 75 75 !! !! 76 76 !! namdrg top/bottom drag coefficient (default: NO selection) 77 !! namdrg_top top friction (ln_ OFF =F & ln_isfcav=T)78 !! namdrg_bot bottom friction (ln_ OFF =F)77 !! namdrg_top top friction (ln_drg_OFF =F & ln_isfcav=T) 78 !! namdrg_bot bottom friction (ln_drg_OFF =F) 79 79 !! nambbc bottom temperature boundary condition (default: OFF) 80 80 !! nambbl bottom boundary layer scheme (default: OFF) -
NEMO/branches/2020/ticket2487/tests/BENCH/EXPREF/namelist_cfg_orca12_like
r11536 r15264 75 75 !! !! 76 76 !! namdrg top/bottom drag coefficient (default: NO selection) 77 !! namdrg_top top friction (ln_ OFF =F & ln_isfcav=T)78 !! namdrg_bot bottom friction (ln_ OFF =F)77 !! namdrg_top top friction (ln_drg_OFF =F & ln_isfcav=T) 78 !! namdrg_bot bottom friction (ln_drg_OFF =F) 79 79 !! nambbc bottom temperature boundary condition (default: OFF) 80 80 !! nambbl bottom boundary layer scheme (default: OFF) -
NEMO/branches/2020/ticket2487/tests/BENCH/EXPREF/namelist_cfg_orca1_like
r11536 r15264 75 75 !! !! 76 76 !! namdrg top/bottom drag coefficient (default: NO selection) 77 !! namdrg_top top friction (ln_ OFF =F & ln_isfcav=T)78 !! namdrg_bot bottom friction (ln_ OFF =F)77 !! namdrg_top top friction (ln_drg_OFF =F & ln_isfcav=T) 78 !! namdrg_bot bottom friction (ln_drg_OFF =F) 79 79 !! nambbc bottom temperature boundary condition (default: OFF) 80 80 !! nambbl bottom boundary layer scheme (default: OFF) -
NEMO/branches/2020/ticket2487/tests/BENCH/MY_SRC/usrdef_sbc.F90
r10179 r15264 126 126 REAL(wp), DIMENSION(:,:,:), INTENT(in) :: phi ! ice thickness 127 127 !! 128 REAL(wp) :: zfr1, zfr2 ! local variables129 128 REAL(wp), DIMENSION(jpi,jpj) :: zsnw ! snw distribution after wind blowing 130 129 !!--------------------------------------------------------------------- … … 147 146 ! ice fields deduced from above 148 147 zsnw(:,:) = 1._wp 149 !!CALL lim_thd_snwblow( at_i_b, zsnw ) ! snow distribution over ice after150 !wind blowing151 148 emp_ice (:,:) = SUM( a_i_b(:,:,:) * evap_ice(:,:,:), dim=3 ) - sprecip(:,:) * zsnw(:,:) 152 149 emp_oce (:,:) = emp_oce(:,:) - sprecip(:,:) * (1._wp - zsnw(:,:) ) … … 161 158 qsr_tot (:,:) = at_i_b(:,:) * qsr_oce(:,:) + SUM( a_i_b(:,:,:) * qsr_ice(:,:,:), dim=3 ) 162 159 163 ! --- shortwave radiation transmitted below the surface (W/m2, see Grenfell Maykut 77) --- ! 164 zfr1 = ( 0.18 * ( 1.0 - cldf_ice ) + 0.35 * cldf_ice ) ! transmission when hi>10cm 165 zfr2 = ( 0.82 * ( 1.0 - cldf_ice ) + 0.65 * cldf_ice ) ! zfr2 such that zfr1 + zfr2 to equal 1 166 ! 167 WHERE ( phs(:,:,:) <= 0._wp .AND. phi(:,:,:) < 0.1_wp ) ! linear decrease from hi=0 to 10cm 168 qtr_ice_top(:,:,:) = qsr_ice(:,:,:) * ( zfr1 + zfr2 * ( 1._wp - phi(:,:,:) * 10._wp ) ) 169 ELSEWHERE( phs(:,:,:) <= 0._wp .AND. phi(:,:,:) >= 0.1_wp ) ! constant (zfr1) when hi>10cm 170 qtr_ice_top(:,:,:) = qsr_ice(:,:,:) * zfr1 171 ELSEWHERE ! zero when hs>0 172 qtr_ice_top(:,:,:) = 0._wp 173 END WHERE 160 ! --- shortwave radiation transmitted below the surface (W/m2) 161 qtr_ice_top(:,:,:) = 0._wp 174 162 #endif 175 163 -
NEMO/branches/2020/ticket2487/tests/CANAL/EXPREF/file_def_nemo-oce.xml
r9572 r15264 15 15 <field field_ref="soce" /> 16 16 <field field_ref="ssh" /> 17 <field field_ref="s algrad" />18 <field field_ref=" ke_zint" />17 <field field_ref="socegrad" /> 18 <field field_ref="eken_int" /> 19 19 <field field_ref="relvor" /> 20 20 <field field_ref="potvor" /> -
NEMO/branches/2020/ticket2487/tests/CANAL/EXPREF/namelist_cfg
r12206 r15264 20 20 &namusr_def ! User defined : CANAL configuration: Flat bottom, beta-plane 21 21 !----------------------------------------------------------------------- 22 rn_domszx = 3600. ! x horizontal size [km]23 rn_domszy = 1 800. ! y horizontal size [km]24 rn_domszz = 5000. ! z vertical size [m]25 rn_dx = 30. ! x horizontal resolution [km]26 rn_dy = 30. ! y horizontal resolution [km]27 rn_dz = 500. ! z vertical resolution [m]22 rn_domszx = 2000. ! x horizontal size [km] 23 rn_domszy = 1000. ! y horizontal size [km] 24 rn_domszz = 1000. ! z vertical size [m] 25 rn_dx = 10. ! x horizontal resolution [km] 26 rn_dy = 10. ! y horizontal resolution [km] 27 rn_dz = 1000. ! z vertical resolution [m] 28 28 rn_0xratio = 0.5 ! x-domain ratio of the 0 29 29 rn_0yratio = 0.5 ! y-domain ratio of the 0 … … 31 31 rn_ppgphi0 = 38.5 ! Reference latitude [degrees] 32 32 rn_u10 = 0. ! 10m wind speed [m/s] 33 rn_windszx = 4000.! longitudinal wind extension [km]34 rn_windszy = 4000.! latitudinal wind extension [km]35 rn_uofac = 0.! Uoce multiplicative factor (0.:absolute or 1.:relative winds)33 rn_windszx = 90. ! longitudinal wind extension [km] 34 rn_windszy = 90. ! latitudinal wind extension [km] 35 !!clem rn_uofac = 0. ! Uoce multiplicative factor (0.:absolute or 1.:relative winds) 36 36 rn_vtxmax = 1. ! initial vortex max current [m/s] 37 37 rn_uzonal = 1. ! initial zonal current [m/s] 38 rn_ujetszx = 4000. 39 rn_ujetszy = 400 0. ! latitudinal jet extension [km]38 rn_ujetszx = 4000. ! longitudinal jet extension [km] 39 rn_ujetszy = 400. ! latitudinal jet extension [km] 40 40 nn_botcase = 0 ! bottom definition (0:flat, 1:bump) 41 nn_initcase = 1 ! initial condition case (0:rest, 1:zonal current, 2:current shear, 3: gaussian zonal current, 42 ! ! 4: geostrophic zonal pulse, 5: vortex) 43 ln_sshnoise = .false. ! add random noise on initial ssh 44 rn_lambda = 50. ! gaussian lambda 41 nn_initcase = 1 ! initial condition case 42 ! ! -1 : stratif at rest 43 ! ! 0 : rest 44 ! ! 1 : zonal current 45 ! ! 2 : current shear 46 ! ! 3 : gaussian zonal current 47 ! ! 4 : geostrophic zonal pulse 48 ! ! 5 : baroclinic vortex 49 ln_sshnoise = .FALSE. ! add random noise on initial ssh 50 rn_lambda = 50. ! gaussian lambda 51 nn_perio = 1 45 52 / 46 53 !----------------------------------------------------------------------- … … 59 66 !----------------------------------------------------------------------- 60 67 ln_linssh = .false. ! =T linear free surface ==>> model level are fixed in time 61 rn_rdt = 1440. ! time step for the dynamics (and tracer if nn_acc=0) 62 rn_atfp = 0.05 ! asselin time filter parameter 68 rn_rdt = 1200. ! time step for the dynamics (and tracer if nn_acc=0) 69 rn_atfp = 0.0 ! asselin time filter parameter 70 / 71 !----------------------------------------------------------------------- 72 &namcfg ! parameters of the configuration (default: use namusr_def in namelist_cfg) 73 !----------------------------------------------------------------------- 74 ln_write_cfg = .false. ! (=T) create the domain configuration file 75 cn_domcfg_out = "domain_cfg" ! newly created domain configuration filename 63 76 / 64 77 !!====================================================================== … … 108 121 !! !! 109 122 !! namdrg top/bottom drag coefficient (default: NO selection) 110 !! namdrg_top top friction (ln_ OFF =F & ln_isfcav=T)111 !! namdrg_bot bottom friction (ln_ OFF =F)123 !! namdrg_top top friction (ln_drg_OFF =F & ln_isfcav=T) 124 !! namdrg_bot bottom friction (ln_drg_OFF =F) 112 125 !! nambbc bottom temperature boundary condition (default: OFF) 113 126 !! nambbl bottom boundary layer scheme (default: OFF) … … 117 130 &namdrg ! top/bottom drag coefficient (default: NO selection) 118 131 !----------------------------------------------------------------------- 119 ln_ OFF = .true. ! free-slip : Cd = 0132 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 120 133 / 121 134 !!====================================================================== … … 148 161 ln_traadv_OFF = .false. ! No tracer advection 149 162 ln_traadv_cen = .false. ! 2nd order centered scheme 150 nn_cen_h = 4! =2/4, horizontal 2nd order CEN / 4th order CEN151 nn_cen_v = 4! =2/4, vertical 2nd order CEN / 4th order COMPACT163 nn_cen_h = 2 ! =2/4, horizontal 2nd order CEN / 4th order CEN 164 nn_cen_v = 2 ! =2/4, vertical 2nd order CEN / 4th order COMPACT 152 165 ln_traadv_fct = .false. ! FCT scheme 153 nn_fct_h = 2! =2/4, horizontal 2nd / 4th order166 nn_fct_h = 4 ! =2/4, horizontal 2nd / 4th order 154 167 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 155 168 ln_traadv_mus = .false. ! MUSCL scheme … … 162 175 &namtra_ldf ! lateral diffusion scheme for tracers (default: NO selection) 163 176 !----------------------------------------------------------------------- 164 ln_traldf_OFF = .true. ! No explicit diffusion 177 ! ! Operator type: 178 ln_traldf_OFF = .true. ! No explicit diffusion 179 ln_traldf_lap = .false. ! laplacian operator 180 ln_traldf_blp = .false. ! bilaplacian operator 181 ! 182 ! ! Direction of action: 183 ln_traldf_lev = .false. ! iso-level 184 ln_traldf_hor = .true. ! horizontal (geopotential) 185 ln_traldf_iso = .false. ! iso-neutral (standard operator) 186 ln_traldf_triad = .false. ! iso-neutral (triad operator) 187 ! 188 ! ! iso-neutral options: 189 ln_traldf_msc = .false. ! Method of Stabilizing Correction (both operators) 190 rn_slpmax = 0.01 ! slope limit (both operators) 191 ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) 192 rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) 193 ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) 194 ! 195 ! ! Coefficients: 196 nn_aht_ijk_t = 31 ! space/time variation of eddy coefficient: 197 ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file 198 ! ! = 0 constant 199 ! ! = 10 F(k) =ldf_c1d 200 ! ! = 20 F(i,j) =ldf_c2d 201 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 202 ! ! = 30 F(i,j,k) =ldf_c2d * ldf_c1d 203 ! ! = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 204 ! ! time invariant coefficients: aht0 = 1/2 Ud*Ld (lap case) 205 ! ! or = 1/12 Ud*Ld^3 (blp case) 206 rn_Ud = 0.01 ! lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30) 207 rn_Ld = 200.e+3 ! lateral diffusive length [m] (nn_aht_ijk_t= 0, 10) 165 208 / 166 209 !!====================================================================== … … 183 226 nn_dynkeg = 0 ! scheme for grad(KE): =0 C2 ; =1 Hollingsworth correction 184 227 ln_dynadv_cen2 = .false. ! flux form - 2nd order centered scheme 185 ln_dynadv_ubs = .true. ! flux form - 3rd order UBS scheme228 ln_dynadv_ubs = .true. ! flux form - 3rd order UBS scheme 186 229 / 187 230 !----------------------------------------------------------------------- 188 231 &namdyn_vor ! Vorticity / Coriolis scheme (default: NO selection) 189 232 !----------------------------------------------------------------------- 190 ln_dynvor_ene = . true. ! energy conserving scheme191 ln_dynvor_ens = .false. ! enstrophy conserving scheme192 ln_dynvor_mix = .false. ! mixed scheme233 ln_dynvor_ene = .false. ! energy conserving scheme 234 ln_dynvor_ens = .false. ! enstrophy conserving scheme 235 ln_dynvor_mix = .false. ! mixed scheme 193 236 ln_dynvor_een = .false. ! energy & enstrophy scheme 237 ln_dynvor_enT = .false. ! energy conserving scheme (T-point) 238 ln_dynvor_eeT = .true. ! energy conserving scheme (een using e3t) 194 239 nn_een_e3f = 0 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 195 240 / … … 210 255 ! ! = 1 Boxcar over nn_baro sub-steps 211 256 ! ! = 2 Boxcar over 2*nn_baro " " 212 ln_bt_auto = . false. ! Number of sub-step defined from:257 ln_bt_auto = .true. ! Number of sub-step defined from: 213 258 nn_baro = 24 ! =F : the number of sub-step in rn_rdt seconds 214 259 / … … 222 267 ! ! Direction of action : 223 268 ln_dynldf_lev = .false. ! iso-level 224 ln_dynldf_hor = . true. ! horizontal (geopotential)269 ln_dynldf_hor = .false. ! horizontal (geopotential) 225 270 ln_dynldf_iso = .false. ! iso-neutral 226 271 ! ! Coefficient 227 nn_ahm_ijk_t = 20! space/time variation of eddy coef272 nn_ahm_ijk_t = 31 ! space/time variation of eddy coef 228 273 ! ! =-30 read in eddy_viscosity_3D.nc file 229 274 ! ! =-20 read in eddy_viscosity_2D.nc file … … 275 320 !! namdiu Cool skin and warm layer models (default: OFF) 276 321 !! namdiu Cool skin and warm layer models (default: OFF) 277 !! namflo float parameters (default: OFF) 278 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF) 279 !! nam_diadct transports through some sections (default: OFF) 322 !! namflo float parameters ("key_float") 323 !! nam_diaharm Harmonic analysis of tidal constituents ("key_diaharm") 324 !! namdct transports through some sections ("key_diadct") 325 !! nam_diatmb Top Middle Bottom Output (default: OFF) 280 326 !! nam_dia25h 25h Mean Output (default: OFF) 281 327 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") … … 286 332 !----------------------------------------------------------------------- 287 333 ln_glo_trd = .false. ! (T) global domain averaged diag for T, T^2, KE, and PE 288 ln_dyn_trd = .true. ! (T) 3D momentum trend output334 ln_dyn_trd = .true. ! (T) 3D momentum trend output 289 335 ln_dyn_mxl = .false. ! (T) 2D momentum trends averaged over the mixed layer (not coded yet) 290 336 ln_vor_trd = .false. ! (T) 2D barotropic vorticity trends (not coded yet) … … 313 359 &nammpp ! Massively Parallel Processing ("key_mpp_mpi") 314 360 !----------------------------------------------------------------------- 361 !! jpni = 8 ! jpni number of processors following i (set automatically if < 1) 362 !! jpnj = 1 ! jpnj number of processors following j (set automatically if < 1) 315 363 / 316 364 !----------------------------------------------------------------------- 317 365 &namctl ! Control prints (default: OFF) 318 366 !----------------------------------------------------------------------- 367 ln_timing = .true. ! timing by routine write out in timing.output file 368 !! ln_diacfl = .true. ! CFL diagnostics write out in cfl_diagnostics.ascii 319 369 / 320 370 !----------------------------------------------------------------------- -
NEMO/branches/2020/ticket2487/tests/CANAL/MY_SRC/diawri.F90
r12206 r15264 230 230 IF( iom_use('logavs') ) CALL iom_put( "logavs", LOG( MAX( 1.e-20_wp, avs(:,:,:) ) ) ) 231 231 232 IF ( iom_use("s algrad") .OR. iom_use("salgrad2") ) THEN232 IF ( iom_use("socegrad") .OR. iom_use("socegrad2") ) THEN 233 233 z3d(:,:,jpk) = 0. 234 234 DO jk = 1, jpkm1 … … 244 244 END DO 245 245 CALL lbc_lnk( 'diawri', z3d, 'T', 1. ) 246 CALL iom_put( "s algrad2", z3d ) ! square of module of sal gradient246 CALL iom_put( "socegrad2", z3d ) ! square of module of sal gradient 247 247 z3d(:,:,:) = SQRT( z3d(:,:,:) ) 248 CALL iom_put( "s algrad" , z3d ) ! module of sal gradient248 CALL iom_put( "socegrad" , z3d ) ! module of sal gradient 249 249 ENDIF 250 250 … … 299 299 END DO 300 300 END DO 301 CALL iom_put( "salt2c", rau0 * z2d ) ! vertically integrated s alt content (PSU*kg/m2)302 ENDIF 303 ! 304 IF ( iom_use("eken") ) THEN301 CALL iom_put( "salt2c", rau0 * z2d ) ! vertically integrated squared salt content (PSU*kg/m2) 302 ENDIF 303 ! 304 IF ( iom_use("eken") .OR. iom_use("eken_int") ) THEN 305 305 z3d(:,:,jpk) = 0._wp 306 306 DO jk = 1, jpkm1 307 DO jj = 2, jpj 308 DO ji = 2, jpi 307 DO jj = 2, jpjm1 308 DO ji = 2, jpim1 309 309 zztmpx = 0.5 * ( un(ji-1,jj ,jk) + un(ji,jj,jk) ) 310 310 zztmpy = 0.5 * ( vn(ji ,jj-1,jk) + vn(ji,jj,jk) ) … … 315 315 CALL lbc_lnk( 'diawri', z3d, 'T', 1. ) 316 316 CALL iom_put( "eken", z3d ) ! kinetic energy 317 ENDIF318 319 IF ( iom_use("ke") .or. iom_use("ke_zint") ) THEN320 !321 z3d(:,:,jpk) = 0._wp322 z3d(1,:, : ) = 0._wp323 z3d(:,1, : ) = 0._wp324 DO jk = 1, jpkm1325 DO jj = 2, jpj326 DO ji = 2, jpi327 z3d(ji,jj,jk) = 0.25_wp * ( un(ji ,jj,jk) * un(ji ,jj,jk) * e1e2u(ji ,jj) * e3u_n(ji ,jj,jk) &328 & + un(ji-1,jj,jk) * un(ji-1,jj,jk) * e1e2u(ji-1,jj) * e3u_n(ji-1,jj,jk) &329 & + vn(ji,jj ,jk) * vn(ji,jj ,jk) * e1e2v(ji,jj ) * e3v_n(ji,jj ,jk) &330 & + vn(ji,jj-1,jk) * vn(ji,jj-1,jk) * e1e2v(ji,jj-1) * e3v_n(ji,jj-1,jk) ) &331 & * r1_e1e2t(ji,jj) / e3t_n(ji,jj,jk) * tmask(ji,jj,jk)332 END DO333 END DO334 END DO335 336 CALL lbc_lnk( 'diawri', z3d, 'T', 1. )337 CALL iom_put( "ke", z3d ) ! kinetic energy338 317 339 318 z2d(:,:) = 0._wp … … 341 320 DO jj = 1, jpj 342 321 DO ji = 1, jpi 343 z2d(ji,jj) = z2d(ji,jj) + e3t_n(ji,jj,jk) * z3d(ji,jj,jk) * tmask(ji,jj,jk) 344 END DO 345 END DO 346 END DO 347 CALL iom_put( "ke_zint", z2d ) ! vertically integrated kinetic energy 348 322 z2d(ji,jj) = z2d(ji,jj) + e3t_n(ji,jj,jk) * z3d(ji,jj,jk) * e1e2t(ji,jj) * tmask(ji,jj,jk) 323 END DO 324 END DO 325 END DO 326 CALL iom_put( "eken_int", z2d ) ! vertically integrated kinetic energy 349 327 ENDIF 350 328 ! … … 358 336 DO ji = 1, fs_jpim1 ! vector opt. 359 337 z3d(ji,jj,jk) = ( e2v(ji+1,jj ) * vn(ji+1,jj ,jk) - e2v(ji,jj) * vn(ji,jj,jk) & 360 & 338 & - e1u(ji ,jj+1) * un(ji ,jj+1,jk) + e1u(ji,jj) * un(ji,jj,jk) ) * r1_e1e2f(ji,jj) 361 339 END DO 362 340 END DO -
NEMO/branches/2020/ticket2487/tests/CANAL/MY_SRC/usrdef_istate.F90
r10425 r15264 64 64 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~ ' 65 65 ! 66 IF (ln_sshnoise) CALL RANDOM_NUMBER(zrandom)67 66 zjetx = ABS(rn_ujetszx)/2. 68 67 zjety = ABS(rn_ujetszy)/2. 69 68 ! 69 zf0 = 2._wp * omega * SIN( rad * rn_ppgphi0 ) 70 ! 70 71 SELECT CASE(nn_initcase) 72 73 CASE(-1) ! stratif at rest 74 75 ! sea level: 76 pssh(:,:) = 0. 77 ! temperature: 78 pts(:,:,1,jp_tem) = 25. !!30._wp 79 pts(:,:,2:jpk,jp_tem) = 22. !!24._wp 80 ! salinity: 81 pts(:,:,:,jp_sal) = 35._wp 82 ! velocities: 83 pu(:,:,:) = 0. 84 pv(:,:,:) = 0. 85 71 86 CASE(0) ! rest 72 87 … … 96 111 zbeta = 2._wp * omega * COS( rad * rn_ppgphi0 ) / ra 97 112 WHERE( ABS(gphit) <= zjety ) 98 pssh(:,:) = - rn_uzonal / grav * ( ff_t(:,:)* gphit(:,:) * 1.e3 + 0.5 * zbeta * gphit(:,:) * gphit(:,:) * 1.e6 )99 ELSEWHERE 100 pssh(:,:) = - rn_uzonal / grav * ( ff_t(:,:)* SIGN(zjety, gphit(:,:)) * 1.e3 &113 pssh(:,:) = - rn_uzonal / grav * ( zf0 * gphit(:,:) * 1.e3 + 0.5 * zbeta * gphit(:,:) * gphit(:,:) * 1.e6 ) 114 ELSEWHERE 115 pssh(:,:) = - rn_uzonal / grav * ( zf0 * SIGN(zjety, gphit(:,:)) * 1.e3 & 101 116 & + 0.5 * zbeta * zjety * zjety * 1.e6 ) 102 117 END WHERE … … 107 122 pts(:,:,jpk,jp_sal) = 0. 108 123 DO jk=1, jpkm1 109 pts(:,:,jk,jp_sal) = gphit(:,:) 124 WHERE( ABS(gphit) <= zjety ) 125 !!$ WHERE( ABS(gphit) <= zjety*0.5 .AND. ABS(glamt) <= zjety*0.5 ) ! for a square of salt 126 pts(:,:,jk,jp_sal) = 35. 127 ELSEWHERE 128 pts(:,:,jk,jp_sal) = 30. 129 END WHERE 110 130 END DO 111 131 ! velocities: … … 132 152 WHERE( ABS(gphit) <= zjety ) 133 153 pssh(:,:) = - SIGN(rn_uzonal, gphit(:,:)) / grav & 134 & * ( ff_t(:,:)* gphit(:,:) * 1.e3 + 0.5 * zbeta * gphit(:,:) * gphit(:,:) * 1.e6 )154 & * ( zf0 * gphit(:,:) * 1.e3 + 0.5 * zbeta * gphit(:,:) * gphit(:,:) * 1.e6 ) 135 155 ELSEWHERE 136 156 pssh(:,:) = - SIGN(rn_uzonal, gphit(:,:)) / grav & 137 & * ( ff_t(:,:)* SIGN(zjety, gphit(:,:)) * 1.e3 + 0.5 * zbeta * zjety * zjety * 1.e6 )157 & * ( zf0 * SIGN(zjety, gphit(:,:)) * 1.e3 + 0.5 * zbeta * zjety * zjety * 1.e6 ) 138 158 END WHERE 139 159 END SELECT … … 141 161 pts(:,:,:,jp_tem) = 10._wp 142 162 ! salinity: 143 pts(:,:,:,jp_sal) = 2.144 DO jk=1, jpkm1 145 WHERE( ABS(gphiv) <= zjety ) pts(:,:,jk,jp_sal) = 2. + SIGN(1.,gphiv(:,:))163 pts(:,:,:,jp_sal) = 30. 164 DO jk=1, jpkm1 165 WHERE( ABS(gphiv) <= zjety ) pts(:,:,jk,jp_sal) = 30. + SIGN(1.,gphiv(:,:)) 146 166 END DO 147 167 ! velocities: … … 176 196 ! salinity: 177 197 DO jk=1, jpkm1 178 pts(:,:,jk,jp_sal) = gphit(:,:)198 pts(:,:,jk,jp_sal) = pssh(:,:) 179 199 END DO 180 200 ! velocities: … … 213 233 zf0 = 2._wp * omega * SIN( rad * rn_ppgphi0 ) 214 234 zumax = rn_vtxmax * SIGN(1._wp, zf0) ! Here Anticyclonic: set zumax=-1 for cyclonic 215 zlambda = SQRT(2._wp)*rn_lambda ! Horizontal scale in meters235 zlambda = SQRT(2._wp)*rn_lambda*1.e3 ! Horizontal scale in meters 216 236 zn2 = 3.e-3**2 217 237 zH = 0.5_wp * 5000._wp … … 253 273 ! velocities: 254 274 za = 2._wp * zP0 / zlambda**2 255 DO jj =1, jpj256 DO ji =1, jpim1275 DO jj = 2, jpjm1 276 DO ji = 2, jpim1 257 277 zx = glamu(ji,jj) * 1.e3 258 278 zy = gphiu(ji,jj) * 1.e3 … … 270 290 END DO 271 291 ! 272 DO jj =1, jpjm1273 DO ji =1, jpi292 DO jj = 2, jpjm1 293 DO ji = 2, jpim1 274 294 zx = glamv(ji,jj) * 1.e3 275 295 zy = gphiv(ji,jj) * 1.e3 … … 287 307 END DO 288 308 ! 309 CALL lbc_lnk_multi( 'usrdef_istate', pu, 'U', -1., pv, 'V', -1. ) 310 289 311 END SELECT 290 312 291 313 IF (ln_sshnoise) THEN 314 CALL RANDOM_SEED() 292 315 CALL RANDOM_NUMBER(zrandom) 293 316 pssh(:,:) = pssh(:,:) + ( 0.1 * zrandom(:,:) - 0.05 ) 294 317 END IF 295 CALL lbc_lnk( 'usrdef_istate', pssh, 'T', 1. ) 296 CALL lbc_lnk( 'usrdef_istate', pts, 'T', 1. ) 297 CALL lbc_lnk( 'usrdef_istate', pu, 'U', -1. ) 298 CALL lbc_lnk( 'usrdef_istate', pv, 'V', -1. ) 299 318 300 319 END SUBROUTINE usr_def_istate 301 320 -
NEMO/branches/2020/ticket2487/tests/CANAL/MY_SRC/usrdef_nam.F90
r11899 r15264 50 50 LOGICAL , PUBLIC :: ln_sshnoise=.false. ! add random noise on initial ssh 51 51 REAL(wp), PUBLIC :: rn_lambda = 50. ! gaussian lambda 52 INTEGER , PUBLIC :: nn_perio = 0 ! periodicity of the channel (0=closed, 1=E-W) 52 53 53 54 !!---------------------------------------------------------------------- … … 79 80 !! 80 81 NAMELIST/namusr_def/ rn_domszx, rn_domszy, rn_domszz, rn_dx, rn_dy, rn_dz, rn_0xratio, rn_0yratio & 81 & , nn_fcase, rn_ppgphi0, rn_ vtxmax, rn_uzonal, rn_ujetszx, rn_ujetszy&82 & , rn_ u10, rn_windszx, rn_windszy, rn_uofac&83 & , nn_botcase, nn_initcase, ln_sshnoise, rn_lambda 82 & , nn_fcase, rn_ppgphi0, rn_u10, rn_windszx, rn_windszy & !!, rn_uofac & 83 & , rn_vtxmax, rn_uzonal, rn_ujetszx, rn_ujetszy & 84 & , nn_botcase, nn_initcase, ln_sshnoise, rn_lambda, nn_perio 84 85 !!---------------------------------------------------------------------- 85 86 ! … … 151 152 WRITE(numout,*) ' add random noise on initial ssh ln_sshnoise= ', ln_sshnoise 152 153 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 154 WRITE(numout,*) ' Periodicity of the basin nn_perio = ', nn_perio 156 155 ENDIF 156 ! ! Set the lateral boundary condition of the global domain 157 kperio = nn_perio ! EW_CANAL configuration : closed basin 157 158 ! 158 159 END SUBROUTINE usr_def_nam -
NEMO/branches/2020/ticket2487/tests/CANAL/MY_SRC/usrdef_sbc.F90
r10074 r15264 17 17 USE sbc_oce ! Surface boundary condition: ocean fields 18 18 USE phycst ! physical constants 19 USE usrdef_nam, ONLY : rn_u10, rn_uofac, rn_windszy 19 USE usrdef_nam, ONLY : rn_u10, rn_uofac, rn_windszy, rn_windszx 20 20 ! 21 21 USE in_out_manager ! I/O manager … … 71 71 ! 72 72 utau(:,:) = 0._wp 73 IF( rn_u10 /= 0. .AND. rn_windszy > 0. ) THEN74 WHERE( ABS(gphit) <= rn_windszy/2. ) utau(:,:) = zrhocd * rn_u10 * rn_u1075 ENDIF76 73 vtau(:,:) = 0._wp 77 74 taum(:,:) = 0._wp … … 83 80 qsr (:,:) = 0._wp 84 81 ! 82 ENDIF 83 84 IF( rn_u10 /= 0. .AND. rn_windszy > 0. ) THEN 85 IF( nyear == 1 .AND. nmonth == 1 .AND. nday <= 10 ) THEN 86 WHERE( ABS(gphit) <= rn_windszy/2. .AND. ABS(glamt) <= rn_windszx/2. ) utau(:,:) = zrhocd * rn_u10 * rn_u10 87 ELSE 88 utau(:,:) = 0. 89 ENDIF 85 90 ENDIF 86 91 -
NEMO/branches/2020/ticket2487/tests/CANAL/MY_SRC/usrdef_zgr.F90
r10425 r15264 199 199 zmaxlam = MAXVAL(glamt) 200 200 CALL mpp_max( 'usrdef_zgr', zmaxlam ) ! max over the global domain 201 zscl = rpi / zmaxlam202 z2d(:,:) = 0.5 * ( 1. - COS( glamt(:,:) * zscl ))203 z2d(:,:) = REAL(jpkm1 - NINT( 0. 75 * REAL(jpkm1,wp) * z2d(:,:) ), wp)201 zscl = 0.5 * rpi / zmaxlam 202 z2d(:,:) = COS( glamt(:,:) * zscl ) 203 z2d(:,:) = REAL(jpkm1 - NINT( 0.5 * REAL(jpkm1,wp) * z2d(:,:) ), wp) 204 204 END SELECT 205 205 ! -
NEMO/branches/2020/ticket2487/tests/ICE_ADV1D/EXPREF/namelist_cfg
r12206 r15264 106 106 !! !! 107 107 !! namdrg top/bottom drag coefficient (default: NO selection) 108 !! namdrg_top top friction (ln_ OFF=F & ln_isfcav=T)109 !! namdrg_bot bottom friction (ln_ OFF=F)108 !! namdrg_top top friction (ln_drg_OFF=F & ln_isfcav=T) 109 !! namdrg_bot bottom friction (ln_drg_OFF=F) 110 110 !! nambbc bottom temperature boundary condition (default: OFF) 111 111 !! nambbl bottom boundary layer scheme (default: OFF) … … 115 115 &namdrg ! top/bottom drag coefficient (default: NO selection) 116 116 !----------------------------------------------------------------------- 117 ln_ OFF = .true. ! free-slip : Cd = 0117 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 118 118 / 119 119 !!====================================================================== -
NEMO/branches/2020/ticket2487/tests/ICE_ADV1D/EXPREF/namelist_cfg_120pts
r12206 r15264 106 106 !! !! 107 107 !! namdrg top/bottom drag coefficient (default: NO selection) 108 !! namdrg_top top friction (ln_ OFF=F & ln_isfcav=T)109 !! namdrg_bot bottom friction (ln_ OFF=F)108 !! namdrg_top top friction (ln_drg_OFF=F & ln_isfcav=T) 109 !! namdrg_bot bottom friction (ln_drg_OFF=F) 110 110 !! nambbc bottom temperature boundary condition (default: OFF) 111 111 !! nambbl bottom boundary layer scheme (default: OFF) … … 115 115 &namdrg ! top/bottom drag coefficient (default: NO selection) 116 116 !----------------------------------------------------------------------- 117 ln_ OFF = .true. ! free-slip : Cd = 0117 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 118 118 / 119 119 !!====================================================================== -
NEMO/branches/2020/ticket2487/tests/ICE_ADV1D/EXPREF/namelist_cfg_240pts
r12206 r15264 106 106 !! !! 107 107 !! namdrg top/bottom drag coefficient (default: NO selection) 108 !! namdrg_top top friction (ln_ OFF=F & ln_isfcav=T)109 !! namdrg_bot bottom friction (ln_ OFF=F)108 !! namdrg_top top friction (ln_drg_OFF=F & ln_isfcav=T) 109 !! namdrg_bot bottom friction (ln_drg_OFF=F) 110 110 !! nambbc bottom temperature boundary condition (default: OFF) 111 111 !! nambbl bottom boundary layer scheme (default: OFF) … … 115 115 &namdrg ! top/bottom drag coefficient (default: NO selection) 116 116 !----------------------------------------------------------------------- 117 ln_ OFF = .true. ! free-slip : Cd = 0117 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 118 118 / 119 119 !!====================================================================== -
NEMO/branches/2020/ticket2487/tests/ICE_ADV1D/EXPREF/namelist_cfg_60pts
r12206 r15264 106 106 !! !! 107 107 !! namdrg top/bottom drag coefficient (default: NO selection) 108 !! namdrg_top top friction (ln_ OFF=F & ln_isfcav=T)109 !! namdrg_bot bottom friction (ln_ OFF=F)108 !! namdrg_top top friction (ln_drg_OFF=F & ln_isfcav=T) 109 !! namdrg_bot bottom friction (ln_drg_OFF=F) 110 110 !! nambbc bottom temperature boundary condition (default: OFF) 111 111 !! nambbl bottom boundary layer scheme (default: OFF) … … 115 115 &namdrg ! top/bottom drag coefficient (default: NO selection) 116 116 !----------------------------------------------------------------------- 117 ln_ OFF = .true. ! free-slip : Cd = 0117 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 118 118 / 119 119 !!====================================================================== -
NEMO/branches/2020/ticket2487/tests/ICE_ADV1D/EXPREF/namelist_ice_cfg
r10535 r15264 88 88 !------------------------------------------------------------------------------ 89 89 ln_iceini = .true. ! activate ice initialization (T) or not (F) 90 ln_iceini_file = .true. ! netcdf file provided for initialization (T) or not (F)90 nn_iceini_file = 1 ! netcdf file provided for initialization 91 91 92 92 sn_hti = 'initice_60pts' , -12 ,'hti' , .false. , .true., 'yearly' , '' , '', '' -
NEMO/branches/2020/ticket2487/tests/ICE_ADV1D/EXPREF/namelist_ice_cfg_120pts
r10431 r15264 88 88 !------------------------------------------------------------------------------ 89 89 ln_iceini = .true. ! activate ice initialization (T) or not (F) 90 ln_iceini_file = .true. ! netcdf file provided for initialization (T) or not (F)90 nn_iceini_file = 1 ! netcdf file provided for initialization 91 91 92 92 sn_hti = 'initice_120pts' , -12 ,'hti' , .false. , .true., 'yearly' , '' , '', '' -
NEMO/branches/2020/ticket2487/tests/ICE_ADV1D/EXPREF/namelist_ice_cfg_240pts
r10431 r15264 88 88 !------------------------------------------------------------------------------ 89 89 ln_iceini = .true. ! activate ice initialization (T) or not (F) 90 ln_iceini_file = .true. ! netcdf file provided for initialization (T) or not (F)90 nn_iceini_file = 1 ! netcdf file provided for initialization 91 91 92 92 sn_hti = 'initice_240pts' , -12 ,'hti' , .false. , .true., 'yearly' , '' , '', '' -
NEMO/branches/2020/ticket2487/tests/ICE_ADV1D/EXPREF/namelist_ice_cfg_60pts
r10431 r15264 88 88 !------------------------------------------------------------------------------ 89 89 ln_iceini = .true. ! activate ice initialization (T) or not (F) 90 ln_iceini_file = .true. ! netcdf file provided for initialization (T) or not (F)90 nn_iceini_file = 1 ! netcdf file provided for initialization 91 91 92 92 sn_hti = 'initice_60pts' , -12 ,'hti' , .false. , .true., 'yearly' , '' , '', '' -
NEMO/branches/2020/ticket2487/tests/ICE_ADV1D/MY_SRC/usrdef_sbc.F90
r10513 r15264 108 108 REAL(wp), DIMENSION(:,:,:), INTENT(in) :: phi ! ice thickness 109 109 !! 110 INTEGER :: jl 110 111 REAL(wp) :: zfr1, zfr2 ! local variables 111 112 REAL(wp), DIMENSION(jpi,jpj) :: zsnw ! snw distribution after wind blowing 113 REAL(wp), DIMENSION(jpi,jpj) :: ztri 112 114 !!--------------------------------------------------------------------- 113 115 ! … … 142 144 143 145 ! --- shortwave radiation transmitted below the surface (W/m2, see Grenfell Maykut 77) --- ! 144 zfr1 = ( 0.18 * ( 1.0 - cldf_ice ) + 0.35 * cldf_ice ) ! transmission when hi>10cm145 z fr2 = ( 0.82 * ( 1.0 - cldf_ice ) + 0.65 * cldf_ice ) ! zfr2 such that zfr1 + zfr2 to equal 1146 cloud_fra(:,:) = pp_cldf 147 ztri(:,:) = 0.18 * ( 1.0 - cloud_fra(:,:) ) + 0.35 * cloud_fra(:,:) ! surface transmission when hi>10cm 146 148 ! 147 WHERE ( phs(:,:,:) <= 0._wp .AND. phi(:,:,:) < 0.1_wp ) ! linear decrease from hi=0 to 10cm 148 qtr_ice_top(:,:,:) = qsr_ice(:,:,:) * ( zfr1 + zfr2 * ( 1._wp - phi(:,:,:) * 10._wp ) ) 149 ELSEWHERE( phs(:,:,:) <= 0._wp .AND. phi(:,:,:) >= 0.1_wp ) ! constant (zfr1) when hi>10cm 150 qtr_ice_top(:,:,:) = qsr_ice(:,:,:) * zfr1 151 ELSEWHERE ! zero when hs>0 152 qtr_ice_top(:,:,:) = 0._wp 153 END WHERE 154 149 DO jl = 1, jpl 150 WHERE ( phs(:,:,jl) <= 0._wp .AND. phi(:,:,jl) < 0.1_wp ) ! linear decrease from hi=0 to 10cm 151 qtr_ice_top(:,:,jl) = qsr_ice(:,:,jl) * ( ztri(:,:) + ( 1._wp - ztri(:,:) ) * ( 1._wp - phi(:,:,jl) * 10._wp ) ) 152 ELSEWHERE( phs(:,:,jl) <= 0._wp .AND. phi(:,:,jl) >= 0.1_wp ) ! constant (ztri) when hi>10cm 153 qtr_ice_top(:,:,jl) = qsr_ice(:,:,jl) * ztri(:,:) 154 ELSEWHERE ! zero when hs>0 155 qtr_ice_top(:,:,jl) = 0._wp 156 END WHERE 157 ENDDO 158 159 155 160 END SUBROUTINE usrdef_sbc_ice_flx 156 161 -
NEMO/branches/2020/ticket2487/tests/ICE_ADV2D/EXPREF/file_def_nemo-ice.xml
r10516 r15264 55 55 <field field_ref="normstr" name="normstr" /> 56 56 <field field_ref="sheastr" name="sheastr" /> 57 <field field_ref="isig1" name="isig1" />58 <field field_ref="isig2" name="isig2" />59 <field field_ref="isig3" name="isig3" />60 57 61 58 <!-- heat fluxes --> -
NEMO/branches/2020/ticket2487/tests/ICE_ADV2D/EXPREF/namelist_cfg
r12206 r15264 106 106 !! !! 107 107 !! namdrg top/bottom drag coefficient (default: NO selection) 108 !! namdrg_top top friction (ln_ OFF=F & ln_isfcav=T)109 !! namdrg_bot bottom friction (ln_ OFF=F)108 !! namdrg_top top friction (ln_drg_OFF=F & ln_isfcav=T) 109 !! namdrg_bot bottom friction (ln_drg_OFF=F) 110 110 !! nambbc bottom temperature boundary condition (default: OFF) 111 111 !! nambbl bottom boundary layer scheme (default: OFF) … … 115 115 &namdrg ! top/bottom drag coefficient (default: NO selection) 116 116 !----------------------------------------------------------------------- 117 ln_ OFF = .true. ! free-slip : Cd = 0117 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 118 118 / 119 119 !!====================================================================== -
NEMO/branches/2020/ticket2487/tests/ICE_ADV2D/EXPREF/namelist_ice_cfg
r10535 r15264 86 86 !------------------------------------------------------------------------------ 87 87 ln_iceini = .true. ! activate ice initialization (T) or not (F) 88 ln_iceini_file = .true. ! netcdf file provided for initialization (T) or not (F)88 nn_iceini_file = 1 ! netcdf file provided for initialization 89 89 90 90 sn_hti = 'initice' , -12 ,'hti' , .false. , .true., 'yearly' , '' , '', '' -
NEMO/branches/2020/ticket2487/tests/ICE_ADV2D/MY_SRC/usrdef_sbc.F90
r10515 r15264 108 108 REAL(wp), DIMENSION(:,:,:), INTENT(in) :: phi ! ice thickness 109 109 !! 110 INTEGER :: jl 110 111 REAL(wp) :: zfr1, zfr2 ! local variables 111 112 REAL(wp), DIMENSION(jpi,jpj) :: zsnw ! snw distribution after wind blowing 113 REAL(wp), DIMENSION(jpi,jpj) :: ztri 112 114 !!--------------------------------------------------------------------- 113 115 ! … … 142 144 143 145 ! --- shortwave radiation transmitted below the surface (W/m2, see Grenfell Maykut 77) --- ! 144 zfr1 = ( 0.18 * ( 1.0 - cldf_ice ) + 0.35 * cldf_ice ) ! transmission when hi>10cm145 z fr2 = ( 0.82 * ( 1.0 - cldf_ice ) + 0.65 * cldf_ice ) ! zfr2 such that zfr1 + zfr2 to equal 1146 cloud_fra(:,:) = pp_cldf 147 ztri(:,:) = 0.18 * ( 1.0 - cloud_fra(:,:) ) + 0.35 * cloud_fra(:,:) ! surface transmission when hi>10cm 146 148 ! 147 WHERE ( phs(:,:,:) <= 0._wp .AND. phi(:,:,:) < 0.1_wp ) ! linear decrease from hi=0 to 10cm 148 qtr_ice_top(:,:,:) = qsr_ice(:,:,:) * ( zfr1 + zfr2 * ( 1._wp - phi(:,:,:) * 10._wp ) ) 149 ELSEWHERE( phs(:,:,:) <= 0._wp .AND. phi(:,:,:) >= 0.1_wp ) ! constant (zfr1) when hi>10cm 150 qtr_ice_top(:,:,:) = qsr_ice(:,:,:) * zfr1 151 ELSEWHERE ! zero when hs>0 152 qtr_ice_top(:,:,:) = 0._wp 153 END WHERE 154 149 DO jl = 1, jpl 150 WHERE ( phs(:,:,jl) <= 0._wp .AND. phi(:,:,jl) < 0.1_wp ) ! linear decrease from hi=0 to 10cm 151 qtr_ice_top(:,:,jl) = qsr_ice(:,:,jl) * ( ztri(:,:) + ( 1._wp - ztri(:,:) ) * ( 1._wp - phi(:,:,jl) * 10._wp ) ) 152 ELSEWHERE( phs(:,:,jl) <= 0._wp .AND. phi(:,:,jl) >= 0.1_wp ) ! constant (ztri) when hi>10cm 153 qtr_ice_top(:,:,jl) = qsr_ice(:,:,jl) * ztri(:,:) 154 ELSEWHERE ! zero when hs>0 155 qtr_ice_top(:,:,jl) = 0._wp 156 END WHERE 157 ENDDO 158 159 155 160 END SUBROUTINE usrdef_sbc_ice_flx 156 161 -
NEMO/branches/2020/ticket2487/tests/ICE_AGRIF/EXPREF/1_namelist_cfg
r12206 r15264 106 106 !! !! 107 107 !! namdrg top/bottom drag coefficient (default: NO selection) 108 !! namdrg_top top friction (ln_ OFF=F & ln_isfcav=T)109 !! namdrg_bot bottom friction (ln_ OFF=F)108 !! namdrg_top top friction (ln_drg_OFF=F & ln_isfcav=T) 109 !! namdrg_bot bottom friction (ln_drg_OFF=F) 110 110 !! nambbc bottom temperature boundary condition (default: OFF) 111 111 !! nambbl bottom boundary layer scheme (default: OFF) … … 115 115 &namdrg ! top/bottom drag coefficient (default: NO selection) 116 116 !----------------------------------------------------------------------- 117 ln_ OFF = .true. ! free-slip : Cd = 0117 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 118 118 / 119 119 !!====================================================================== -
NEMO/branches/2020/ticket2487/tests/ICE_AGRIF/EXPREF/file_def_nemo-ice.xml
r11159 r15264 53 53 <field field_ref="normstr" name="normstr" /> 54 54 <field field_ref="sheastr" name="sheastr" /> 55 <field field_ref="isig1" name="isig1" />56 <field field_ref="isig2" name="isig2" />57 <field field_ref="isig3" name="isig3" />58 55 59 56 <!-- heat fluxes --> -
NEMO/branches/2020/ticket2487/tests/ICE_AGRIF/EXPREF/namelist_cfg
r12206 r15264 106 106 !! !! 107 107 !! namdrg top/bottom drag coefficient (default: NO selection) 108 !! namdrg_top top friction (ln_ OFF=F & ln_isfcav=T)109 !! namdrg_bot bottom friction (ln_ OFF=F)108 !! namdrg_top top friction (ln_drg_OFF=F & ln_isfcav=T) 109 !! namdrg_bot bottom friction (ln_drg_OFF=F) 110 110 !! nambbc bottom temperature boundary condition (default: OFF) 111 111 !! nambbl bottom boundary layer scheme (default: OFF) … … 115 115 &namdrg ! top/bottom drag coefficient (default: NO selection) 116 116 !----------------------------------------------------------------------- 117 ln_ OFF = .true. ! free-slip : Cd = 0117 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 118 118 / 119 119 !!====================================================================== -
NEMO/branches/2020/ticket2487/tests/ICE_AGRIF/EXPREF/namelist_ice_cfg
r10535 r15264 86 86 !------------------------------------------------------------------------------ 87 87 ln_iceini = .true. ! activate ice initialization (T) or not (F) 88 ln_iceini_file = .true. ! netcdf file provided for initialization (T) or not (F)88 nn_iceini_file = 1 ! netcdf file provided for initialization 89 89 90 90 sn_hti = 'initice' , -12 ,'hti' , .false. , .true., 'yearly' , '' , '', '' -
NEMO/branches/2020/ticket2487/tests/ICE_AGRIF/MY_SRC/usrdef_sbc.F90
r10516 r15264 108 108 REAL(wp), DIMENSION(:,:,:), INTENT(in) :: phi ! ice thickness 109 109 !! 110 INTEGER :: jl 110 111 REAL(wp) :: zfr1, zfr2 ! local variables 111 112 REAL(wp), DIMENSION(jpi,jpj) :: zsnw ! snw distribution after wind blowing 113 REAL(wp), DIMENSION(jpi,jpj) :: ztri 112 114 !!--------------------------------------------------------------------- 113 115 ! … … 142 144 143 145 ! --- shortwave radiation transmitted below the surface (W/m2, see Grenfell Maykut 77) --- ! 144 zfr1 = ( 0.18 * ( 1.0 - cldf_ice ) + 0.35 * cldf_ice ) ! transmission when hi>10cm145 z fr2 = ( 0.82 * ( 1.0 - cldf_ice ) + 0.65 * cldf_ice ) ! zfr2 such that zfr1 + zfr2 to equal 1146 cloud_fra(:,:) = pp_cldf 147 ztri(:,:) = 0.18 * ( 1.0 - cloud_fra(:,:) ) + 0.35 * cloud_fra(:,:) ! surface transmission when hi>10cm 146 148 ! 147 WHERE ( phs(:,:,:) <= 0._wp .AND. phi(:,:,:) < 0.1_wp ) ! linear decrease from hi=0 to 10cm 148 qtr_ice_top(:,:,:) = qsr_ice(:,:,:) * ( zfr1 + zfr2 * ( 1._wp - phi(:,:,:) * 10._wp ) ) 149 ELSEWHERE( phs(:,:,:) <= 0._wp .AND. phi(:,:,:) >= 0.1_wp ) ! constant (zfr1) when hi>10cm 150 qtr_ice_top(:,:,:) = qsr_ice(:,:,:) * zfr1 151 ELSEWHERE ! zero when hs>0 152 qtr_ice_top(:,:,:) = 0._wp 153 END WHERE 149 DO jl = 1, jpl 150 WHERE ( phs(:,:,jl) <= 0._wp .AND. phi(:,:,jl) < 0.1_wp ) ! linear decrease from hi=0 to 10cm 151 qtr_ice_top(:,:,jl) = qsr_ice(:,:,jl) * ( ztri(:,:) + ( 1._wp - ztri(:,:) ) * ( 1._wp - phi(:,:,jl) * 10._wp ) ) 152 ELSEWHERE( phs(:,:,jl) <= 0._wp .AND. phi(:,:,jl) >= 0.1_wp ) ! constant (ztri) when hi>10cm 153 qtr_ice_top(:,:,jl) = qsr_ice(:,:,jl) * ztri(:,:) 154 ELSEWHERE ! zero when hs>0 155 qtr_ice_top(:,:,jl) = 0._wp 156 END WHERE 157 ENDDO 154 158 155 159 END SUBROUTINE usrdef_sbc_ice_flx -
NEMO/branches/2020/ticket2487/tests/ISOMIP/EXPREF/namelist_cfg
r12206 r15264 128 128 !! !! 129 129 !! namdrg top/bottom drag coefficient (default: NO selection) 130 !! namdrg_top top friction (ln_ OFF=F & ln_isfcav=T)131 !! namdrg_bot bottom friction (ln_ OFF=F)130 !! namdrg_top top friction (ln_drg_OFF=F & ln_isfcav=T) 131 !! namdrg_bot bottom friction (ln_drg_OFF=F) 132 132 !! nambbc bottom temperature boundary condition (default: OFF) 133 133 !! nambbl bottom boundary layer scheme (default: OFF) … … 137 137 &namdrg ! top/bottom drag coefficient (default: NO selection) 138 138 !----------------------------------------------------------------------- 139 ln_ OFF = .false.! free-slip : Cd = 0 (F => fill namdrg_bot139 ln_drg_OFF = .false. ! free-slip : Cd = 0 (F => fill namdrg_bot 140 140 ln_lin = .false. ! linear drag: Cd = Cd0 Uc0 & namdrg_top) 141 141 ln_non_lin = .true. ! non-linear drag: Cd = Cd0 |U| -
NEMO/branches/2020/ticket2487/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT2_flux_cen2_cfg
r10075 r15264 67 67 &namdrg ! top/bottom drag coefficient (default: NO selection) 68 68 !----------------------------------------------------------------------- 69 ln_ OFF = .true. ! free-slip : Cd = 069 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 70 70 / 71 71 !----------------------------------------------------------------------- -
NEMO/branches/2020/ticket2487/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT2_flux_ubs_cfg
r12206 r15264 112 112 !! !! 113 113 !! namdrg top/bottom drag coefficient (default: NO selection) 114 !! namdrg_top top friction (ln_ OFF=F & ln_isfcav=T)115 !! namdrg_bot bottom friction (ln_ OFF=F)114 !! namdrg_top top friction (ln_drg_OFF=F & ln_isfcav=T) 115 !! namdrg_bot bottom friction (ln_drg_OFF=F) 116 116 !! nambbc bottom temperature boundary condition (default: OFF) 117 117 !! nambbl bottom boundary layer scheme (default: OFF) … … 121 121 &namdrg ! top/bottom drag coefficient (default: NO selection) 122 122 !----------------------------------------------------------------------- 123 ln_ OFF = .true. ! free-slip : Cd = 0123 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 124 124 / 125 125 !!====================================================================== -
NEMO/branches/2020/ticket2487/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT2_vect_eenH_cfg
r10075 r15264 67 67 &namdrg ! top/bottom drag coefficient (default: NO selection) 68 68 !----------------------------------------------------------------------- 69 ln_ OFF = .true. ! free-slip : Cd = 069 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 70 70 / 71 71 !----------------------------------------------------------------------- -
NEMO/branches/2020/ticket2487/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT2_vect_een_cfg
r10075 r15264 67 67 &namdrg ! top/bottom drag coefficient (default: NO selection) 68 68 !----------------------------------------------------------------------- 69 ln_ OFF = .true. ! free-slip : Cd = 069 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 70 70 / 71 71 !----------------------------------------------------------------------- -
NEMO/branches/2020/ticket2487/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT2_vect_ene_cfg
r10075 r15264 67 67 &namdrg ! top/bottom drag coefficient (default: NO selection) 68 68 !----------------------------------------------------------------------- 69 ln_ OFF = .true. ! free-slip : Cd = 069 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 70 70 / 71 71 !----------------------------------------------------------------------- -
NEMO/branches/2020/ticket2487/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT2_vect_ens_cfg
r10075 r15264 67 67 &namdrg ! top/bottom drag coefficient (default: NO selection) 68 68 !----------------------------------------------------------------------- 69 ln_ OFF = .true. ! free-slip : Cd = 069 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 70 70 / 71 71 !----------------------------------------------------------------------- -
NEMO/branches/2020/ticket2487/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT4_flux_cen2_cfg
r10075 r15264 67 67 &namdrg ! top/bottom drag coefficient (default: NO selection) 68 68 !----------------------------------------------------------------------- 69 ln_ OFF = .true. ! free-slip : Cd = 069 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 70 70 / 71 71 !----------------------------------------------------------------------- -
NEMO/branches/2020/ticket2487/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT4_flux_ubs_cfg
r10075 r15264 67 67 &namdrg ! top/bottom drag coefficient (default: NO selection) 68 68 !----------------------------------------------------------------------- 69 ln_ OFF = .true. ! free-slip : Cd = 069 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 70 70 / 71 71 !----------------------------------------------------------------------- -
NEMO/branches/2020/ticket2487/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT4_vect_eenH_cfg
r10075 r15264 67 67 &namdrg ! top/bottom drag coefficient (default: NO selection) 68 68 !----------------------------------------------------------------------- 69 ln_ OFF = .true. ! free-slip : Cd = 069 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 70 70 / 71 71 !----------------------------------------------------------------------- -
NEMO/branches/2020/ticket2487/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT4_vect_een_cfg
r10075 r15264 67 67 &namdrg ! top/bottom drag coefficient (default: NO selection) 68 68 !----------------------------------------------------------------------- 69 ln_ OFF = .true. ! free-slip : Cd = 069 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 70 70 / 71 71 !----------------------------------------------------------------------- -
NEMO/branches/2020/ticket2487/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT4_vect_ene_cfg
r10075 r15264 67 67 &namdrg ! top/bottom drag coefficient (default: NO selection) 68 68 !----------------------------------------------------------------------- 69 ln_ OFF = .true. ! free-slip : Cd = 069 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 70 70 / 71 71 !----------------------------------------------------------------------- -
NEMO/branches/2020/ticket2487/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT4_vect_ens_cfg
r10075 r15264 67 67 &namdrg ! top/bottom drag coefficient (default: NO selection) 68 68 !----------------------------------------------------------------------- 69 ln_ OFF = .true. ! free-slip : Cd = 069 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 70 70 / 71 71 !----------------------------------------------------------------------- -
NEMO/branches/2020/ticket2487/tests/OVERFLOW/EXPREF/namelist_sco_FCT2_flux_cen-ahm1000_cfg
r11848 r15264 71 71 &namdrg ! top/bottom drag coefficient (default: NO selection) 72 72 !----------------------------------------------------------------------- 73 ln_ OFF = .true.! free-slip : Cd = 0 (F => fill namdrg_bot73 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 74 74 ln_lin = .false. ! linear drag: Cd = Cd0 Uc0 & namdrg_top) 75 75 ln_non_lin = .false. ! non-linear drag: Cd = Cd0 |U| -
NEMO/branches/2020/ticket2487/tests/OVERFLOW/EXPREF/namelist_sco_FCT2_flux_ubs_cfg
r11848 r15264 71 71 &namdrg ! top/bottom drag coefficient (default: NO selection) 72 72 !----------------------------------------------------------------------- 73 ln_ OFF= .true. ! free-slip : Cd = 0 (F => fill namdrg_bot73 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 74 74 ln_lin = .false. ! linear drag: Cd = Cd0 Uc0 & namdrg_top) 75 75 ln_non_lin = .false. ! non-linear drag: Cd = Cd0 |U| -
NEMO/branches/2020/ticket2487/tests/OVERFLOW/EXPREF/namelist_sco_FCT4_flux_cen-ahm1000_cfg
r10075 r15264 71 71 &namdrg ! top/bottom drag coefficient (default: NO selection) 72 72 !----------------------------------------------------------------------- 73 ln_ OFF = .true.! free-slip : Cd = 0 (F => fill namdrg_bot73 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 74 74 ln_lin = .false. ! linear drag: Cd = Cd0 Uc0 & namdrg_top) 75 75 ln_non_lin = .false. ! non-linear drag: Cd = Cd0 |U| -
NEMO/branches/2020/ticket2487/tests/OVERFLOW/EXPREF/namelist_sco_FCT4_flux_ubs_cfg
r11848 r15264 71 71 &namdrg ! top/bottom drag coefficient (default: NO selection) 72 72 !----------------------------------------------------------------------- 73 ln_ OFF = .true.! free-slip : Cd = 0 (F => fill namdrg_bot73 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 74 74 ln_lin = .false. ! linear drag: Cd = Cd0 Uc0 & namdrg_top) 75 75 ln_non_lin = .false. ! non-linear drag: Cd = Cd0 |U| -
NEMO/branches/2020/ticket2487/tests/OVERFLOW/EXPREF/namelist_zps_FCT2_flux_ubs_cfg
r11848 r15264 71 71 &namdrg ! top/bottom drag coefficient (default: NO selection) 72 72 !----------------------------------------------------------------------- 73 ln_ OFF = .true.! free-slip : Cd = 0 (F => fill namdrg_bot73 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 74 74 ln_lin = .false. ! linear drag: Cd = Cd0 Uc0 & namdrg_top) 75 75 ln_non_lin = .false. ! non-linear drag: Cd = Cd0 |U| -
NEMO/branches/2020/ticket2487/tests/OVERFLOW/EXPREF/namelist_zps_FCT4_flux_ubs_cfg
r12206 r15264 105 105 !! !! 106 106 !! namdrg top/bottom drag coefficient (default: NO selection) 107 !! namdrg_top top friction (ln_ OFF=F & ln_isfcav=T)108 !! namdrg_bot bottom friction (ln_ OFF=F)107 !! namdrg_top top friction (ln_drg_OFF=F & ln_isfcav=T) 108 !! namdrg_bot bottom friction (ln_drg_OFF=F) 109 109 !! nambbc bottom temperature boundary condition (default: OFF) 110 110 !! nambbl bottom boundary layer scheme (default: OFF) … … 114 114 &namdrg ! top/bottom drag coefficient (default: NO selection) 115 115 !----------------------------------------------------------------------- 116 ln_ OFF= .true. ! free-slip : Cd = 0 (F => fill namdrg_bot116 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 117 117 ln_lin = .false. ! linear drag: Cd = Cd0 Uc0 & namdrg_top) 118 118 ln_non_lin = .false. ! non-linear drag: Cd = Cd0 |U| -
NEMO/branches/2020/ticket2487/tests/OVERFLOW/EXPREF/namelist_zps_FCT4_vect_een_cfg
r11848 r15264 71 71 &namdrg ! top/bottom drag coefficient (default: NO selection) 72 72 !----------------------------------------------------------------------- 73 ln_ OFF = .true.! free-slip : Cd = 0 (F => fill namdrg_bot73 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 74 74 ln_lin = .false. ! linear drag: Cd = Cd0 Uc0 & namdrg_top) 75 75 ln_non_lin = .false. ! non-linear drag: Cd = Cd0 |U| -
NEMO/branches/2020/ticket2487/tests/README.rst
r11743 r15264 53 53 List of test cases 54 54 ================== 55 56 BENCH 57 ----- 58 | Benchmark configuration. Allow to run any configuration (including ORCA type or BDY) with idealized grid 59 and initial state so it does not need any input file other than the namelists. 60 As usual, all configuration changes can be done through the namelist. 61 We provide 3 example of namelist_cfg to mimic ORCA1, OR025 or ORCA12 configurations. 62 By default do not produce any output file. An extensive description of BENCH will be abailable in 63 Irrmann et al. 2021. 55 64 56 65 ICE_AGRIF -
NEMO/branches/2020/ticket2487/tests/VORTEX/EXPREF/1_namelist_cfg
r12206 r15264 107 107 !! !! 108 108 !! namdrg top/bottom drag coefficient (default: NO selection) 109 !! namdrg_top top friction (ln_ OFF=F & ln_isfcav=T)110 !! namdrg_bot bottom friction (ln_ OFF=F)109 !! namdrg_top top friction (ln_drg_OFF=F & ln_isfcav=T) 110 !! namdrg_bot bottom friction (ln_drg_OFF=F) 111 111 !! nambbc bottom temperature boundary condition (default: OFF) 112 112 !! nambbl bottom boundary layer scheme (default: OFF) … … 116 116 &namdrg ! top/bottom drag coefficient (default: NO selection) 117 117 !----------------------------------------------------------------------- 118 ln_ OFF = .true. ! free-slip : Cd = 0118 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 119 119 / 120 120 !!====================================================================== -
NEMO/branches/2020/ticket2487/tests/VORTEX/EXPREF/namelist_cfg
r12206 r15264 99 99 !! !! 100 100 !! namdrg top/bottom drag coefficient (default: NO selection) 101 !! namdrg_top top friction (ln_ OFF=F & ln_isfcav=T)102 !! namdrg_bot bottom friction (ln_ OFF=F)101 !! namdrg_top top friction (ln_drg_OFF=F & ln_isfcav=T) 102 !! namdrg_bot bottom friction (ln_drg_OFF=F) 103 103 !! nambbc bottom temperature boundary condition (default: OFF) 104 104 !! nambbl bottom boundary layer scheme (default: OFF) … … 108 108 &namdrg ! top/bottom drag coefficient (default: NO selection) 109 109 !----------------------------------------------------------------------- 110 ln_ OFF = .true. ! free-slip : Cd = 0110 ln_drg_OFF = .true. ! free-slip : Cd = 0 (F => fill namdrg_bot 111 111 / 112 112 !!====================================================================== -
NEMO/branches/2020/ticket2487/tests/WAD/EXPREF/namelist_cfg
r12206 r15264 200 200 !! !! 201 201 !! namdrg top/bottom drag coefficient (default: NO selection) 202 !! namdrg_top top friction (ln_ OFF=F & ln_isfcav=T)203 !! namdrg_bot bottom friction (ln_ OFF=F)202 !! namdrg_top top friction (ln_drg_OFF=F & ln_isfcav=T) 203 !! namdrg_bot bottom friction (ln_drg_OFF=F) 204 204 !! nambbc bottom temperature boundary condition (default: OFF) 205 205 !! nambbl bottom boundary layer scheme (default: OFF)
Note: See TracChangeset
for help on using the changeset viewer.