Changeset 5836
- Timestamp:
- 2015-10-26T15:49:40+01:00 (7 years ago)
- Location:
- trunk/NEMOGCM
- Files:
-
- 31 deleted
- 217 edited
- 8 copied
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/CONFIG/AMM12/EXP00/namelist_cfg
r5501 r5836 196 196 / 197 197 !----------------------------------------------------------------------- 198 &namcla ! cross land advection199 !-----------------------------------------------------------------------200 /201 !-----------------------------------------------------------------------202 &namobc ! open boundaries parameters ("key_obc")203 !-----------------------------------------------------------------------204 /205 !-----------------------------------------------------------------------206 198 &namagrif ! AGRIF zoom ("key_agrif") 207 199 !----------------------------------------------------------------------- … … 281 273 &namtra_adv ! advection scheme for tracer 282 274 !----------------------------------------------------------------------- 275 ln_traadv_fct = .true. ! FCT scheme 276 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 277 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 278 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping 279 ! ! (number of sub-timestep = nn_fct_zts) 283 280 / 284 281 !----------------------------------------------------------------------- … … 289 286 &namtra_ldf ! lateral diffusion scheme for tracers 290 287 !---------------------------------------------------------------------------------- 291 ln_traldf_hor = .true. ! horizontal (geopotential) (needs "key_ldfslp" when ln_sco=T) 292 ln_traldf_iso = .false. ! iso-neutral (needs "key_ldfslp") 293 rn_aeiv_0 = 0. ! eddy induced velocity coefficient [m2/s] 294 rn_aht_0 = 50. ! horizontal eddy diffusivity for tracers [m2/s] 288 ! ! Operator type: 289 ln_traldf_lap = .true. ! laplacian operator 290 ln_traldf_blp = .false. ! bilaplacian operator 291 ! ! Direction of action: 292 ln_traldf_lev = .false. ! iso-level 293 ln_traldf_hor = .true. ! horizontal (geopotential) 294 ln_traldf_iso = .false. ! iso-neutral 295 ln_traldf_triad = .false. ! iso-neutral using Griffies triads 296 ! 297 ! ! iso-neutral options: 298 ln_traldf_msc = .true. ! Method of Stabilizing Correction (both operators) 299 rn_slpmax = 0.01 ! slope limit (both operators) 300 ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) 301 rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) 302 ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) 303 ! 304 ! ! Coefficients: 305 nn_aht_ijk_t = 0 ! space/time variation of eddy coef 306 ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file 307 ! ! = 0 constant 308 ! ! = 10 F(k) =ldf_c1d 309 ! ! = 20 F(i,j) =ldf_c2d 310 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 311 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d 312 ! ! = 31 F(i,j,k,t)=F(local velocity) 313 rn_aht_0 = 50. ! lateral eddy diffusivity (lap. operator) [m2/s] 314 rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] 315 / 316 !---------------------------------------------------------------------------------- 317 &namtra_ldfeiv ! eddy induced velocity param. 318 !---------------------------------------------------------------------------------- 319 ln_ldfeiv =.false. ! use eddy induced velocity parameterization 295 320 / 296 321 !----------------------------------------------------------------------- … … 306 331 &namdyn_vor ! option of physics/algorithm (not control by CPP keys) 307 332 !----------------------------------------------------------------------- 333 ln_dynvor_ene = .false. ! enstrophy conserving scheme 334 ln_dynvor_ens = .false. ! energy conserving scheme 335 ln_dynvor_mix = .false. ! mixed scheme 336 ln_dynvor_een = .true. ! energy & enstrophy scheme 337 nn_een_e3f = 1 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 308 338 / 309 339 !----------------------------------------------------------------------- … … 324 354 !----------------------------------------------------------------------- 325 355 ! ! Type of the operator : 326 ln_dynldf_bilap = .true. ! bilaplacian operator 327 ln_dynldf_lap = .false. ! bilaplacian operator 356 ! ! no diffusion: set ln_dynldf_lap=..._blp=F 357 ln_dynldf_lap = .false. ! laplacian operator 358 ln_dynldf_blp = .true. ! bilaplacian operator 328 359 ! ! Direction of action : 329 ln_dynldf_level = .true. ! iso-level 330 ln_dynldf_hor = .false. ! horizontal (geopotential) (require "key_ldfslp" in s-coord.) 331 ! Coefficient 332 rn_ahm_0_lap = 60.0 ! horizontal laplacian eddy viscosity [m2/s] 333 rn_ahm_0_blp = -1.0e+10 ! horizontal bilaplacian eddy viscosity [m4/s] 360 ln_dynldf_lev = .true. ! iso-level 361 ln_dynldf_hor = .false. ! horizontal (geopotential) 362 ln_dynldf_iso = .false. ! iso-neutral 363 ! ! Coefficient 364 nn_ahm_ijk_t = 0 ! space/time variation of eddy coef 365 ! ! =-30 read in eddy_viscosity_3D.nc file 366 ! ! =-20 read in eddy_viscosity_2D.nc file 367 ! ! = 0 constant 368 ! ! = 10 F(k)=c1d 369 ! ! = 20 F(i,j)=F(grid spacing)=c2d 370 ! ! = 30 F(i,j,k)=c2d*c1d 371 ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) 372 rn_ahm_0 = 60. ! horizontal laplacian eddy viscosity [m2/s] 373 rn_ahm_b = 0. ! background eddy viscosity for ldf_iso [m2/s] 374 rn_bhm_0 = 1.0e+10 ! horizontal bilaplacian eddy viscosity [m4/s] 375 ! 376 ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 334 377 / 335 378 !----------------------------------------------------------------------- … … 349 392 !----------------------------------------------------------------------- 350 393 / 351 !------------------------------------------------------------------------352 &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionally:353 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb")354 /355 394 !----------------------------------------------------------------------- 356 395 &namzdf_gls ! GLS vertical diffusion ("key_zdfgls") -
trunk/NEMOGCM/CONFIG/AMM12/cpp_AMM12.fcm
r4245 r5836 1 bld::tool::fppkeys key_bdy key_tide key_dynspg_ts key_ldfslpkey_zdfgls key_vvl key_diainstant key_mpp_mpi key_iomput1 bld::tool::fppkeys key_bdy key_tide key_dynspg_ts key_zdfgls key_vvl key_diainstant key_mpp_mpi key_iomput -
trunk/NEMOGCM/CONFIG/C1D_PAPA/EXP00/namelist_cfg
r5407 r5836 172 172 / 173 173 !----------------------------------------------------------------------- 174 &namcla ! cross land advection175 !-----------------------------------------------------------------------176 /177 !-----------------------------------------------------------------------178 &namobc ! open boundaries parameters ("key_obc")179 !-----------------------------------------------------------------------180 /181 !-----------------------------------------------------------------------182 174 &namagrif ! AGRIF zoom ("key_agrif") 183 175 !----------------------------------------------------------------------- … … 226 218 &namtra_adv ! advection scheme for tracer 227 219 !----------------------------------------------------------------------- 220 ! C1D : no advection scheme 228 221 / 229 222 !----------------------------------------------------------------------- … … 231 224 !----------------------------------------------------------------------- 232 225 / 233 !----------------------------------------------------------------------- -----------226 !----------------------------------------------------------------------- 234 227 &namtra_ldf ! lateral diffusion scheme for tracers 235 !---------------------------------------------------------------------------------- 236 !---------------------------------------------------------------------------------- 237 ln_traldf_hor = .true. ! horizontal (geopotential) (needs "key_ldfslp" when ln_sco=T) 238 ln_traldf_iso = .false. ! iso-neutral (needs "key_ldfslp") 239 rn_aeiv_0 = 0. ! eddy induced velocity coefficient [m2/s] 240 rn_aht_0 = 0. ! horizontal eddy diffusivity for tracers [m2/s] 228 !----------------------------------------------------------------------- 229 ! C1D : no lateral diffusion 230 / 231 !----------------------------------------------------------------------- 232 &namtra_ldfeiv ! eddy induced velocity param. 233 !----------------------------------------------------------------------- 234 ! C1D : no eiv 241 235 / 242 236 !----------------------------------------------------------------------- … … 248 242 &namdyn_adv ! formulation of the momentum advection 249 243 !----------------------------------------------------------------------- 244 ! C1D : no advection scheme 250 245 / 251 246 !----------------------------------------------------------------------- … … 285 280 &namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion ("key_zdftke") 286 281 !----------------------------------------------------------------------- 287 /288 !------------------------------------------------------------------------289 &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionally:290 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb")291 282 / 292 283 !----------------------------------------------------------------------- -
trunk/NEMOGCM/CONFIG/GYRE/EXP00/namelist_cfg
r5407 r5836 156 156 / 157 157 !----------------------------------------------------------------------- 158 &namcla ! cross land advection159 !-----------------------------------------------------------------------160 /161 !-----------------------------------------------------------------------162 &namobc ! open boundaries parameters ("key_obc")163 !-----------------------------------------------------------------------164 /165 !-----------------------------------------------------------------------166 158 &namagrif ! AGRIF zoom ("key_agrif") 167 159 !----------------------------------------------------------------------- … … 223 215 &namtra_adv ! advection scheme for tracer 224 216 !----------------------------------------------------------------------- 225 ln_traadv_cen2 = .false. ! 2nd order centered scheme 226 ln_traadv_tvd = .true. ! TVD scheme 227 ln_traadv_muscl = .false. ! MUSCL scheme 228 ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 229 ln_traadv_ubs = .false. ! UBS scheme 230 ln_traadv_qck = .false. ! QUICKEST scheme 231 ln_traadv_msc_ups= .false. ! use upstream scheme within muscl 217 ln_traadv_fct = .true. ! FCT scheme 218 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 219 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 220 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping 221 ! ! (number of sub-timestep = nn_fct_zts) 232 222 / 233 223 !----------------------------------------------------------------------- … … 238 228 &namtra_ldf ! lateral diffusion scheme for tracers 239 229 !---------------------------------------------------------------------------------- 240 rn_aeiv_0 = 0. ! eddy induced velocity coefficient [m2/s] 241 rn_aht_0 = 1000. ! horizontal eddy diffusivity for tracers [m2/s] 230 ! ! Operator type: 231 ln_traldf_lap = .true. ! laplacian operator 232 ln_traldf_blp = .false. ! bilaplacian operator 233 ! ! Direction of action: 234 ln_traldf_lev = .false. ! iso-level 235 ln_traldf_hor = .false. ! horizontal (geopotential) 236 ln_traldf_iso = .true. ! iso-neutral 237 ln_traldf_triad = .false. ! iso-neutral using Griffies triads 238 ! 239 ! ! iso-neutral options: 240 ln_traldf_msc = .false. ! Method of Stabilizing Correction (both operators) 241 rn_slpmax = 0.01 ! slope limit (both operators) 242 ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) 243 rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) 244 ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) 245 ! 246 ! ! Coefficients: 247 nn_aht_ijk_t = 0 ! space/time variation of eddy coef 248 ! ! =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 249 ! ! = 0 constant 250 ! ! = 10 F(k) =ldf_c1d 251 ! ! = 20 F(i,j) =ldf_c2d 252 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 253 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d 254 ! ! = 31 F(i,j,k,t)=F(local velocity) 255 rn_aht_0 = 1000. ! lateral eddy diffusivity (lap. operator) [m2/s] 256 rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] 257 / 258 !---------------------------------------------------------------------------------- 259 &namtra_ldfeiv ! eddy induced velocity param. 260 !---------------------------------------------------------------------------------- 261 ln_ldfeiv =.false. ! use eddy induced velocity parameterization 242 262 / 243 263 !----------------------------------------------------------------------- … … 253 273 &namdyn_vor ! option of physics/algorithm (not control by CPP keys) 254 274 !----------------------------------------------------------------------- 255 ln_dynvor_ene = .true. ! energy conserving scheme 256 ln_dynvor_ens = .false. ! enstrophy conserving scheme 275 ln_dynvor_ene = .true. ! enstrophy conserving scheme 276 ln_dynvor_ens = .false. ! energy conserving scheme 277 ln_dynvor_mix = .false. ! mixed scheme 257 278 ln_dynvor_een = .false. ! energy & enstrophy scheme 279 nn_een_e3f = 1 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 258 280 / 259 281 !----------------------------------------------------------------------- … … 270 292 &namdyn_ldf ! lateral diffusion on momentum 271 293 !----------------------------------------------------------------------- 272 rn_ahm_0_lap = 100000. ! horizontal laplacian eddy viscosity [m2/s] 294 ! ! Type of the operator : 295 ! ! no diffusion: set ln_dynldf_lap=..._blp=F 296 ln_dynldf_lap = .true. ! laplacian operator 297 ln_dynldf_blp = .false. ! bilaplacian operator 298 ! ! Direction of action : 299 ln_dynldf_lev = .true. ! iso-level 300 ln_dynldf_hor = .false. ! horizontal (geopotential) 301 ln_dynldf_iso = .false. ! iso-neutral 302 ! ! Coefficient 303 nn_ahm_ijk_t = 0 ! space/time variation of eddy coef 304 ! ! =-30 read in eddy_viscosity_3D.nc file 305 ! ! =-20 read in eddy_viscosity_2D.nc file 306 ! ! = 0 constant 307 ! ! = 10 F(k)=c1d 308 ! ! = 20 F(i,j)=F(grid spacing)=c2d 309 ! ! = 30 F(i,j,k)=c2d*c1d 310 ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) 311 rn_ahm_0 = 100000. ! horizontal laplacian eddy viscosity [m2/s] 312 rn_ahm_b = 0. ! background eddy viscosity for ldf_iso [m2/s] 313 rn_bhm_0 = 0. ! horizontal bilaplacian eddy viscosity [m4/s] 273 314 / 274 315 !----------------------------------------------------------------------- … … 285 326 !----------------------------------------------------------------------- 286 327 nn_etau = 0 ! penetration of tke below the mixed layer (ML) due to internal & intertial waves 287 /288 !------------------------------------------------------------------------289 &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionally:290 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb")291 328 / 292 329 !----------------------------------------------------------------------- -
trunk/NEMOGCM/CONFIG/GYRE/cpp_GYRE.fcm
r4990 r5836 1 bld::tool::fppkeys key_dynspg_flt key_ ldfslp key_zdftke key_iomput key_mpp_mpi key_nosignedzero1 bld::tool::fppkeys key_dynspg_flt key_zdftke key_iomput key_mpp_mpi -
trunk/NEMOGCM/CONFIG/GYRE_BFM/EXP00/namelist_cfg
r5407 r5836 161 161 / 162 162 !----------------------------------------------------------------------- 163 &namcla ! cross land advection164 !-----------------------------------------------------------------------165 /166 !-----------------------------------------------------------------------167 &namobc ! open boundaries parameters ("key_obc")168 !-----------------------------------------------------------------------169 /170 !-----------------------------------------------------------------------171 163 &namagrif ! AGRIF zoom ("key_agrif") 172 164 !----------------------------------------------------------------------- … … 228 220 &namtra_adv ! advection scheme for tracer 229 221 !----------------------------------------------------------------------- 230 ln_traadv_msc_ups= .false. ! use upstream scheme within muscl 222 ln_traadv_fct = .true. ! FCT scheme 223 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 224 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 225 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping 226 ! ! (number of sub-timestep = nn_fct_zts) 231 227 / 232 228 !----------------------------------------------------------------------- … … 237 233 &namtra_ldf ! lateral diffusion scheme for tracers 238 234 !---------------------------------------------------------------------------------- 239 rn_aeiv_0 = 0. ! eddy induced velocity coefficient [m2/s] 240 rn_aht_0 = 1000. ! horizontal eddy diffusivity for tracers [m2/s] 235 ! ! Operator type: 236 ln_traldf_lap = .true. ! laplacian operator 237 ln_traldf_blp = .false. ! bilaplacian operator 238 ! ! Direction of action: 239 ln_traldf_lev = .false. ! iso-level 240 ln_traldf_hor = .false. ! horizontal (geopotential) 241 ln_traldf_iso = .true. ! iso-neutral (standard operator) 242 ln_traldf_triad = .false. ! iso-neutral (triad operator) 243 ! 244 ! ! iso-neutral options: 245 ln_traldf_msc = .true. ! Method of Stabilizing Correction (both operators) 246 rn_slpmax = 0.01 ! slope limit (both operators) 247 ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) 248 rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) 249 ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) 250 ! 251 ! ! Coefficients: 252 nn_aht_ijk_t = 0 ! space/time variation of eddy coef 253 ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file 254 ! ! = 0 constant 255 ! ! = 10 F(k) =ldf_c1d 256 ! ! = 20 F(i,j) =ldf_c2d 257 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 258 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d 259 ! ! = 31 F(i,j,k,t)=F(local velocity) 260 rn_aht_0 = 1000. ! lateral eddy diffusivity (lap. operator) [m2/s] 261 rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] 262 / 263 !---------------------------------------------------------------------------------- 264 &namtra_ldfeiv ! eddy induced velocity param. 265 !---------------------------------------------------------------------------------- 266 ln_ldfeiv =.false. ! use eddy induced velocity parameterization 241 267 / 242 268 !----------------------------------------------------------------------- … … 252 278 &namdyn_vor ! option of physics/algorithm (not control by CPP keys) 253 279 !----------------------------------------------------------------------- 254 ln_dynvor_ene = .true. ! energy conserving scheme 255 ln_dynvor_ens = .false. ! enstrophy conserving scheme 256 ln_dynvor_een = .false. ! energy & enstrophy scheme 280 ln_dynvor_ene = .true. ! enstrophy conserving scheme 257 281 / 258 282 !----------------------------------------------------------------------- … … 268 292 &namdyn_ldf ! lateral diffusion on momentum 269 293 !----------------------------------------------------------------------- 294 ! ! Type of the operator : 295 ! ! no diffusion: set ln_dynldf_lap=..._blp=F 296 ln_dynldf_lap = .true. ! laplacian operator 297 ln_dynldf_blp = .false. ! bilaplacian operator 298 ! ! Direction of action : 299 ln_dynldf_lev = .true. ! iso-level 300 ln_dynldf_hor = .false. ! horizontal (geopotential) 301 ln_dynldf_iso = .false. ! iso-neutral 302 ! ! Coefficient 303 nn_ahm_ijk_t = 0 ! space/time variation of eddy coef 304 ! ! =-30 read in eddy_viscosity_3D.nc file 305 ! ! =-20 read in eddy_viscosity_2D.nc file 306 ! ! = 0 constant 307 ! ! = 10 F(k)=c1d 308 ! ! = 20 F(i,j)=F(grid spacing)=c2d 309 ! ! = 30 F(i,j,k)=c2d*c1d 310 ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) 311 rn_ahm_0 = 100000. ! horizontal laplacian eddy viscosity [m2/s] 312 rn_ahm_b = 0. ! background eddy viscosity for ldf_iso [m2/s] 313 rn_bhm_0 = 0. ! horizontal bilaplacian eddy viscosity [m4/s] 314 ! 315 ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 316 / 270 317 rn_ahm_0_lap = 100000. ! horizontal laplacian eddy viscosity [m2/s] 271 318 / … … 283 330 !----------------------------------------------------------------------- 284 331 nn_etau = 0 ! penetration of tke below the mixed layer (ML) due to internal & intertial waves 285 /286 !------------------------------------------------------------------------287 &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionally:288 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb")289 332 / 290 333 !----------------------------------------------------------------------- -
trunk/NEMOGCM/CONFIG/GYRE_BFM/EXP00/namelist_top_cfg
r4152 r5836 23 23 !----------------------------------------------------------------------- 24 24 &namtrc_adv ! advection scheme for passive tracer 25 !----------------------------------------------------------------------- 26 ln_trcadv_tvd = .true. ! TVD scheme 27 ln_trcadv_muscl = .false. ! MUSCL scheme 25 !----------------------------------------------------------------------- 26 ln_trcadv_fct = .true. ! FCT scheme 27 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 28 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 29 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping 30 ! ! (number of sub-timestep = nn_fct_zts) 28 31 / 29 32 !----------------------------------------------------------------------- -
trunk/NEMOGCM/CONFIG/GYRE_BFM/cpp_GYRE_BFM.fcm
r4230 r5836 1 bld::tool::fppkeys key_dynspg_flt key_ ldfslp key_zdftke key_vectopt_loopkey_top key_my_trc key_mpp_mpi key_iomput1 bld::tool::fppkeys key_dynspg_flt key_zdftke key_top key_my_trc key_mpp_mpi key_iomput 2 2 inc $BFMDIR/src/nemo/bfm.fcm -
trunk/NEMOGCM/CONFIG/GYRE_PISCES/EXP00/namelist_cfg
r5102 r5836 104 104 / 105 105 !----------------------------------------------------------------------- 106 &namcla ! cross land advection107 !-----------------------------------------------------------------------108 /109 !-----------------------------------------------------------------------110 106 &nambfr ! bottom friction 111 107 !----------------------------------------------------------------------- … … 143 139 &namtra_adv ! advection scheme for tracer 144 140 !----------------------------------------------------------------------- 145 ln_traadv_msc_ups= .false. ! use upstream scheme within muscl 141 ln_traadv_fct = .true. ! FCT scheme 142 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 143 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 144 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping 145 ! ! (number of sub-timestep = nn_fct_zts) 146 146 / 147 147 !---------------------------------------------------------------------------------- 148 148 &namtra_ldf ! lateral diffusion scheme for tracers 149 149 !---------------------------------------------------------------------------------- 150 rn_aeiv_0 = 0. ! eddy induced velocity coefficient [m2/s] 151 rn_aht_0 = 1000. ! horizontal eddy diffusivity for tracers [m2/s] 150 ! ! Operator type: 151 ln_traldf_lap = .true. ! laplacian operator 152 ln_traldf_blp = .false. ! bilaplacian operator 153 ! ! Direction of action: 154 ln_traldf_lev = .false. ! iso-level 155 ln_traldf_hor = .false. ! horizontal (geopotential) 156 ln_traldf_iso = .true. ! iso-neutral (standard operator) 157 ln_traldf_triad = .false. ! iso-neutral (triad operator) 158 ! 159 ! ! iso-neutral options: 160 ln_traldf_msc = .true. ! Method of Stabilizing Correction (both operators) 161 rn_slpmax = 0.01 ! slope limit (both operators) 162 ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) 163 rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) 164 ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) 165 ! 166 ! ! Coefficients: 167 nn_aht_ijk_t = 0 ! space/time variation of eddy coef 168 ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file 169 ! ! = 0 constant 170 ! ! = 10 F(k) =ldf_c1d 171 ! ! = 20 F(i,j) =ldf_c2d 172 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 173 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d 174 ! ! = 31 F(i,j,k,t)=F(local velocity) 175 rn_aht_0 = 1000. ! lateral eddy diffusivity (lap. operator) [m2/s] 176 rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] 177 / 178 !---------------------------------------------------------------------------------- 179 &namtra_ldfeiv ! eddy induced velocity param. 180 !---------------------------------------------------------------------------------- 181 ln_ldfeiv =.false. ! use eddy induced velocity parameterization 152 182 / 153 183 !----------------------------------------------------------------------- … … 163 193 &namdyn_vor ! option of physics/algorithm (not control by CPP keys) 164 194 !----------------------------------------------------------------------- 165 ln_dynvor_ene = .true. ! energy conserving scheme 166 ln_dynvor_ens = .false. ! enstrophy conserving scheme 195 ln_dynvor_ene = .true. ! enstrophy conserving scheme 196 ln_dynvor_ens = .false. ! energy conserving scheme 197 ln_dynvor_mix = .false. ! mixed scheme 167 198 ln_dynvor_een = .false. ! energy & enstrophy scheme 199 nn_een_e3f = 1 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 168 200 / 169 201 !----------------------------------------------------------------------- … … 176 208 &namdyn_ldf ! lateral diffusion on momentum 177 209 !----------------------------------------------------------------------- 210 ! ! Type of the operator : 211 ! ! no diffusion: set ln_dynldf_lap=..._blp=F 212 ln_dynldf_lap = .true. ! laplacian operator 213 ln_dynldf_blp = .false. ! bilaplacian operator 214 ! ! Direction of action : 215 ln_dynldf_lev = .true. ! iso-level 216 ln_dynldf_hor = .false. ! horizontal (geopotential) 217 ln_dynldf_iso = .false. ! iso-neutral 218 ! ! Coefficient 219 nn_ahm_ijk_t = 0 ! space/time variation of eddy coef 220 ! ! =-30 read in eddy_viscosity_3D.nc file 221 ! ! =-20 read in eddy_viscosity_2D.nc file 222 ! ! = 0 constant 223 ! ! = 10 F(k)=c1d 224 ! ! = 20 F(i,j)=F(grid spacing)=c2d 225 ! ! = 30 F(i,j,k)=c2d*c1d 226 ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) 227 rn_ahm_0 = 100000. ! horizontal laplacian eddy viscosity [m2/s] 228 rn_ahm_b = 0. ! background eddy viscosity for ldf_iso [m2/s] 229 rn_bhm_0 = 0. ! horizontal bilaplacian eddy viscosity [m4/s] 230 ! 231 ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 232 / 178 233 rn_ahm_0_lap = 100000. ! horizontal laplacian eddy viscosity [m2/s] 179 234 / -
trunk/NEMOGCM/CONFIG/GYRE_PISCES/EXP00/namelist_top_cfg
r4340 r5836 24 24 &namtrc_adv ! advection scheme for passive tracer 25 25 !----------------------------------------------------------------------- 26 ln_trcadv_tvd = .true. ! TVD scheme 27 ln_trcadv_muscl = .false. ! MUSCL scheme 26 ln_trcadv_fct = .true. ! FCT scheme 27 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 28 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 29 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping 30 ! ! (number of sub-timestep = nn_fct_zts) 28 31 / 29 32 !----------------------------------------------------------------------- -
trunk/NEMOGCM/CONFIG/GYRE_PISCES/cpp_GYRE_PISCES.fcm
r4990 r5836 1 bld::tool::fppkeys key_dynspg_flt key_ ldfslp key_zdftke key_top key_pisces_reduced key_iomputkey_mpp_mpi1 bld::tool::fppkeys key_dynspg_flt key_zdftke key_top key_pisces_reduced key_mpp_mpi -
trunk/NEMOGCM/CONFIG/GYRE_XIOS/EXP00/namelist_cfg
r5407 r5836 150 150 / 151 151 !----------------------------------------------------------------------- 152 &namcla ! cross land advection153 !-----------------------------------------------------------------------154 /155 !-----------------------------------------------------------------------156 &namobc ! open boundaries parameters ("key_obc")157 !-----------------------------------------------------------------------158 /159 !-----------------------------------------------------------------------160 152 &namagrif ! AGRIF zoom ("key_agrif") 161 153 !----------------------------------------------------------------------- … … 200 192 &namtra_adv ! advection scheme for tracer 201 193 !----------------------------------------------------------------------- 202 ln_traadv_msc_ups= .false. ! use upstream scheme within muscl 194 ln_traadv_fct = .true. ! FCT scheme 195 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 196 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 197 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping 198 ! ! (number of sub-timestep = nn_fct_zts) 203 199 / 204 200 !----------------------------------------------------------------------- … … 209 205 &namtra_ldf ! lateral diffusion scheme for tracers 210 206 !---------------------------------------------------------------------------------- 211 rn_aeiv_0 = 0. ! eddy induced velocity coefficient [m2/s] 212 rn_aht_0 = 1000. ! horizontal eddy diffusivity for tracers [m2/s] 207 ! ! Operator type: 208 ln_traldf_lap = .true. ! laplacian operator 209 ln_traldf_blp = .false. ! bilaplacian operator 210 ! ! Direction of action: 211 ln_traldf_lev = .false. ! iso-level 212 ln_traldf_hor = .false. ! horizontal (geopotential) 213 ln_traldf_iso = .true. ! iso-neutral (standard operator) 214 ln_traldf_triad = .false. ! iso-neutral (triad operator) 215 ! 216 ! ! iso-neutral options: 217 ln_traldf_msc = .true. ! Method of Stabilizing Correction (both operators) 218 rn_slpmax = 0.01 ! slope limit (both operators) 219 ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) 220 rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) 221 ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) 222 ! 223 ! ! Coefficients: 224 nn_aht_ijk_t = 0 ! space/time variation of eddy coef 225 ! ! =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 226 ! ! = 0 constant 227 ! ! = 10 F(k) =ldf_c1d 228 ! ! = 20 F(i,j) =ldf_c2d 229 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 230 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d 231 ! ! = 31 F(i,j,k,t)=F(local velocity) 232 rn_aht_0 = 1000. ! lateral eddy diffusivity (lap. operator) [m2/s] 233 rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] 234 / 235 !---------------------------------------------------------------------------------- 236 &namtra_ldfeiv ! eddy induced velocity param. 237 !---------------------------------------------------------------------------------- 238 ln_ldfeiv =.false. ! use eddy induced velocity parameterization 213 239 / 214 240 !----------------------------------------------------------------------- -
trunk/NEMOGCM/CONFIG/GYRE_XIOS/cpp_GYRE_XIOS.fcm
r4373 r5836 1 bld::tool::fppkeys key_dynspg_flt key_ ldfslp key_zdftke key_iomput key_mpp_mpi1 bld::tool::fppkeys key_dynspg_flt key_zdftke key_iomput key_mpp_mpi -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/1_namelist_cfg
r5656 r5836 4 4 !----------------------------------------------------------------------- 5 5 &namrun ! parameters of the run 6 nn_it000=1 6 7 !----------------------------------------------------------------------- 7 8 cn_exp = "Agulhas" ! experience name 8 nn_itend = 10950 ! last time step9 nn_itend = 10950 9 10 nn_stock = 10950 ! frequency of creation of a restart file (modulo referenced to 1) 10 11 nn_write = 10950 ! frequency of write in the output file (modulo referenced to nn_it000) 11 ln_clobber = .true. ! clobber (overwrite) an existing file12 ln_clobber = .true. 12 13 / 13 14 !----------------------------------------------------------------------- … … 116 117 / 117 118 !----------------------------------------------------------------------- 118 &namcla ! cross land advection119 !-----------------------------------------------------------------------120 /121 !-----------------------------------------------------------------------122 119 &namagrif ! AGRIF zoom ("key_agrif") 123 120 !----------------------------------------------------------------------- … … 145 142 / 146 143 !----------------------------------------------------------------------- 147 &namtra_adv ! advection scheme for tracer 148 !----------------------------------------------------------------------- 149 / 150 !----------------------------------------------------------------------- 151 &namtra_ldf ! lateral diffusion scheme for tracers 152 !----------------------------------------------------------------------- 153 154 rn_aht_0 = 1000. ! horizontal eddy diffusivity for tracers [m2/s] 155 rn_aeiv_0 = 0. ! eddy induced velocity coefficient [m2/s] (require "key_traldf_eiv") 144 &namtra_adv ! advection scheme for tracer 145 !----------------------------------------------------------------------- 146 ln_traadv_fct = .true. ! FCT scheme 147 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 148 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 149 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping 150 ! ! (number of sub-timestep = nn_fct_zts) 151 / 152 !----------------------------------------------------------------------- 153 &namtra_ldf ! lateral diffusion scheme for tracers 154 !----------------------------------------------------------------------- 155 ! ! Operator type: 156 ln_traldf_lap = .true. ! laplacian operator 157 ln_traldf_blp = .false. ! bilaplacian operator 158 ! ! Direction of action: 159 ln_traldf_lev = .false. ! iso-level 160 ln_traldf_hor = .false. ! horizontal (geopotential) 161 ln_traldf_iso = .true. ! iso-neutral (standard operator) 162 ln_traldf_triad = .false. ! iso-neutral (triad operator) 163 ! 164 ! ! iso-neutral options: 165 ln_traldf_msc = .true. ! Method of Stabilizing Correction (both operators) 166 rn_slpmax = 0.01 ! slope limit (both operators) 167 ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) 168 rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) 169 ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) 170 ! 171 ! ! Coefficients: 172 nn_aht_ijk_t = 20 ! space/time variation of eddy coef 173 ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file 174 ! ! = 0 constant 175 ! ! = 10 F(k) =ldf_c1d 176 ! ! = 20 F(i,j) =ldf_c2d 177 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 178 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d 179 ! ! = 31 F(i,j,k,t)=F(local velocity) 180 rn_aht_0 = 1000. ! lateral eddy diffusivity (lap. operator) [m2/s] 181 rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] 182 / 183 !---------------------------------------------------------------------------------- 184 &namtra_ldfeiv ! eddy induced velocity param. 185 !---------------------------------------------------------------------------------- 186 ln_ldfeiv =.false. ! use eddy induced velocity parameterization 156 187 / 157 188 !----------------------------------------------------------------------- … … 162 193 &namdyn_vor ! option of physics/algorithm (not control by CPP keys) 163 194 !----------------------------------------------------------------------- 195 ln_dynvor_ene = .false. ! enstrophy conserving scheme 196 ln_dynvor_ens = .false. ! energy conserving scheme 197 ln_dynvor_mix = .false. ! mixed scheme 198 ln_dynvor_een = .true. ! energy & enstrophy scheme 199 nn_een_e3f = 1 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 200 ln_dynvor_msk = .true. ! vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) 164 201 / 165 202 !----------------------------------------------------------------------- … … 170 207 &namdyn_ldf ! lateral diffusion on momentum 171 208 !----------------------------------------------------------------------- 172 ! ! Type of the operator : 173 ln_dynldf_lap = .false. ! laplacian operator 174 ln_dynldf_bilap = .true. ! bilaplacian operator 175 rn_ahm_0_blp = -8.5e+11 ! horizontal bilaplacian eddy viscosity [m4/s] 209 ! ! Type of the operator : 210 ! ! no diffusion: set ln_dynldf_lap=..._blp=F 211 ln_dynldf_lap = .false. ! laplacian operator 212 ln_dynldf_blp = .true. ! bilaplacian operator 213 ! ! Direction of action : 214 ln_dynldf_lev = .true. ! iso-level 215 ln_dynldf_hor = .false. ! horizontal (geopotential) 216 ln_dynldf_iso = .false. ! iso-neutral 217 ! ! Coefficient 218 nn_ahm_ijk_t = 0 ! space/time variation of eddy coef 219 ! ! =-30 read in eddy_viscosity_3D.nc file 220 ! ! =-20 read in eddy_viscosity_2D.nc file 221 ! ! = 0 constant 222 ! ! = 10 F(k)=c1d 223 ! ! = 20 F(i,j)=F(grid spacing)=c2d 224 ! ! = 30 F(i,j,k)=c2d*c1d 225 ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) 226 rn_ahm_0 = 40000. ! horizontal laplacian eddy viscosity [m2/s] 227 rn_ahm_b = 0. ! background eddy viscosity for ldf_iso [m2/s] 228 rn_bhm_0 = 8.5e+11 ! horizontal bilaplacian eddy viscosity [m4/s] 176 229 / 177 230 !----------------------------------------------------------------------- -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist_cfg
r4990 r5836 93 93 / 94 94 !----------------------------------------------------------------------- 95 &namcla ! cross land advection96 !-----------------------------------------------------------------------97 /98 !-----------------------------------------------------------------------99 95 &nambfr ! bottom friction 100 96 !----------------------------------------------------------------------- … … 115 111 &namtra_adv ! advection scheme for tracer 116 112 !----------------------------------------------------------------------- 113 ln_traadv_fct = .true. ! FCT scheme 114 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 115 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 116 nn_fct_zts = 0 ! > 1 , 2nd order FCT scheme with vertical sub-timestepping 117 ! ! (number of sub-timestep = nn_fct_zts) 117 118 / 118 119 !----------------------------------------------------------------------- … … 120 121 !----------------------------------------------------------------------- 121 122 / 122 !-----------------------------------------------------------------------123 &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param)124 !-----------------------------------------------------------------------125 /126 123 !---------------------------------------------------------------------------------- 127 124 &namtra_ldf ! lateral diffusion scheme for tracers 128 125 !---------------------------------------------------------------------------------- 126 ! ! Operator type: 127 ln_traldf_lap = .true. ! laplacian operator 128 ln_traldf_blp = .false. ! bilaplacian operator 129 ! ! Direction of action: 130 ln_traldf_lev = .false. ! iso-level 131 ln_traldf_hor = .false. ! horizontal (geopotential) 132 ln_traldf_iso = .true. ! iso-neutral (standard operator) 133 ln_traldf_triad = .false. ! iso-neutral (triad operator) 134 ! 135 ! ! iso-neutral options: 136 ln_traldf_msc = .true. ! Method of Stabilizing Correction (both operators) 137 rn_slpmax = 0.01 ! slope limit (both operators) 138 ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) 139 rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) 140 ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) 141 ! 142 ! ! Coefficients: 143 nn_aht_ijk_t = 20 ! space/time variation of eddy coef 144 ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file 145 ! ! = 0 constant 146 ! ! = 10 F(k) =ldf_c1d 147 ! ! = 20 F(i,j) =ldf_c2d 148 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 149 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d 150 ! ! = 31 F(i,j,k,t)=F(local velocity) 151 rn_aht_0 = 2000. ! lateral eddy diffusivity (lap. operator) [m2/s] 152 rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] 153 / 154 !---------------------------------------------------------------------------------- 155 &namtra_ldfeiv ! eddy induced velocity param. 156 !---------------------------------------------------------------------------------- 157 ln_ldfeiv =.true. ! use eddy induced velocity parameterization 158 ln_ldfeiv_dia =.true. ! diagnose eiv stream function and velocities 159 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] 160 nn_aei_ijk_t = 21 ! space/time variation of the eiv coeficient 161 ! ! =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 162 ! ! = 0 constant 163 ! ! = 10 F(k) =ldf_c1d 164 ! ! = 20 F(i,j) =ldf_c2d 165 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 166 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d 129 167 / 130 168 !----------------------------------------------------------------------- … … 139 177 &namdyn_vor ! option of physics/algorithm (not control by CPP keys) 140 178 !----------------------------------------------------------------------- 179 ln_dynvor_ene = .false. ! enstrophy conserving scheme 180 ln_dynvor_ens = .false. ! energy conserving scheme 181 ln_dynvor_mix = .false. ! mixed scheme 182 ln_dynvor_een = .true. ! energy & enstrophy scheme 183 nn_een_e3f = 1 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 141 184 / 142 185 !----------------------------------------------------------------------- … … 147 190 &namdyn_ldf ! lateral diffusion on momentum 148 191 !----------------------------------------------------------------------- 192 ! ! Type of the operator : 193 ! ! no diffusion: set ln_dynldf_lap=..._blp=F 194 ln_dynldf_lap = .true. ! laplacian operator 195 ln_dynldf_blp = .false. ! bilaplacian operator 196 ! ! Direction of action : 197 ln_dynldf_lev = .true. ! iso-level 198 ln_dynldf_hor = .false. ! horizontal (geopotential) 199 ln_dynldf_iso = .false. ! iso-neutral 200 ! ! Coefficient 201 nn_ahm_ijk_t = -30 ! space/time variation of eddy coef 202 ! ! =-30 read in eddy_viscosity_3D.nc file 203 ! ! =-20 read in eddy_viscosity_2D.nc file 204 ! ! = 0 constant 205 ! ! = 10 F(k)=c1d 206 ! ! = 20 F(i,j)=F(grid spacing)=c2d 207 ! ! = 30 F(i,j,k)=c2d*c1d 208 ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) 209 rn_ahm_0 = 40000. ! horizontal laplacian eddy viscosity [m2/s] 210 rn_ahm_b = 0. ! background eddy viscosity for ldf_iso [m2/s] 211 rn_bhm_0 = 1.e+12 ! horizontal bilaplacian eddy viscosity [m4/s] 149 212 / 150 213 !----------------------------------------------------------------------- -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/cpp_ORCA2_LIM.fcm
r5385 r5836 1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_ diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi key_diaobs key_asminc1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_zdftke key_zdfddm key_zdftmx key_mpp_mpi key_iomput -
trunk/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/1_namelist_cfg
r5499 r5836 116 116 / 117 117 !----------------------------------------------------------------------- 118 &namcla ! cross land advection119 !-----------------------------------------------------------------------120 /121 !-----------------------------------------------------------------------122 118 &namagrif ! AGRIF zoom ("key_agrif") 123 119 !----------------------------------------------------------------------- … … 145 141 / 146 142 !----------------------------------------------------------------------- 147 &namtra_adv ! advection scheme for tracer 148 !----------------------------------------------------------------------- 149 / 150 !----------------------------------------------------------------------- 151 &namtra_ldf ! lateral diffusion scheme for tracers 152 !----------------------------------------------------------------------- 153 154 rn_aht_0 = 1000. ! horizontal eddy diffusivity for tracers [m2/s] 155 rn_aeiv_0 = 0. ! eddy induced velocity coefficient [m2/s] (require "key_traldf_eiv") 143 &namtra_adv ! advection scheme for tracer 144 !----------------------------------------------------------------------- 145 ln_traadv_fct = .true. ! FCT scheme 146 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 147 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 148 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping 149 ! ! (number of sub-timestep = nn_fct_zts) 150 / 151 !----------------------------------------------------------------------- 152 &namtra_ldf ! lateral diffusion scheme for tracers 153 !---------------------------------------------------------------------------------- 154 ! ! Operator type: 155 ln_traldf_lap = .true. ! laplacian operator 156 ln_traldf_blp = .false. ! bilaplacian operator 157 ! ! Direction of action: 158 ln_traldf_lev = .false. ! iso-level 159 ln_traldf_hor = .false. ! horizontal (geopotential) 160 ln_traldf_iso = .true. ! iso-neutral (Standard operator) 161 ln_traldf_triad = .false. ! iso-neutral (Triads operator) 162 ! 163 ! ! iso-neutral options: 164 ln_traldf_msc = .true. ! Method of Stabilizing Correction (both operators) 165 rn_slpmax = 0.01 ! slope limit (both operators) 166 ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) 167 rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) 168 ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) 169 ! 170 ! ! Coefficients: 171 nn_aht_ijk_t = 20 ! space/time variation of eddy coef 172 ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file 173 ! ! = 0 constant 174 ! ! = 10 F(k) =ldf_c1d 175 ! ! = 20 F(i,j) =ldf_c2d 176 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 177 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d 178 ! ! = 31 F(i,j,k,t)=F(local velocity) 179 rn_aht_0 = 1000. ! lateral eddy diffusivity (lap. operator) [m2/s] 180 rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] 181 / 182 !---------------------------------------------------------------------------------- 183 &namtra_ldfeiv ! eddy induced velocity param. 184 !---------------------------------------------------------------------------------- 185 ln_ldfeiv =.false. ! use eddy induced velocity parameterization 156 186 / 157 187 !----------------------------------------------------------------------- … … 162 192 &namdyn_vor ! option of physics/algorithm (not control by CPP keys) 163 193 !----------------------------------------------------------------------- 194 ln_dynvor_ene = .false. ! enstrophy conserving scheme 195 ln_dynvor_ens = .false. ! energy conserving scheme 196 ln_dynvor_mix = .false. ! mixed scheme 197 ln_dynvor_een = .true. ! energy & enstrophy scheme 198 nn_een_e3f = 0 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 199 ln_dynvor_msk = .false. ! vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) 164 200 / 165 201 !----------------------------------------------------------------------- … … 170 206 &namdyn_ldf ! lateral diffusion on momentum 171 207 !----------------------------------------------------------------------- 172 ! ! Type of the operator : 173 ln_dynldf_lap = .false. ! laplacian operator 174 ln_dynldf_bilap = .true. ! bilaplacian operator 175 rn_ahm_0_blp = -8.5e+11 ! horizontal bilaplacian eddy viscosity [m4/s] 208 ! ! Type of the operator : 209 ! ! no diffusion: set ln_dynldf_lap=..._blp=F 210 ln_dynldf_lap = .false. ! laplacian operator 211 ln_dynldf_blp = .true. ! bilaplacian operator 212 ! ! Direction of action : 213 ln_dynldf_lev = .true. ! iso-level 214 ln_dynldf_hor = .false. ! horizontal (geopotential) 215 ln_dynldf_iso = .false. ! iso-neutral 216 ! ! Coefficient 217 nn_ahm_ijk_t = 20 ! space/time variation of eddy coef 218 ! ! =-30 read in eddy_viscosity_3D.nc file 219 ! ! =-20 read in eddy_viscosity_2D.nc file 220 ! ! = 0 constant 221 ! ! = 10 F(k)=c1d 222 ! ! = 20 F(i,j)=F(grid spacing)=c2d 223 ! ! = 30 F(i,j,k)=c2d*c1d 224 ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) 225 rn_ahm_0 = 40000. ! horizontal laplacian eddy viscosity [m2/s] 226 rn_ahm_b = 0. ! background eddy viscosity for ldf_iso [m2/s] 227 rn_bhm_0 = 8.5e+11 ! horizontal bilaplacian eddy viscosity [m4/s] 176 228 / 177 229 !----------------------------------------------------------------------- -
trunk/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/namelist_cfg
r4995 r5836 93 93 / 94 94 !----------------------------------------------------------------------- 95 &namcla ! cross land advection96 !-----------------------------------------------------------------------97 /98 !-----------------------------------------------------------------------99 95 &nambfr ! bottom friction 100 96 !----------------------------------------------------------------------- … … 115 111 &namtra_adv ! advection scheme for tracer 116 112 !----------------------------------------------------------------------- 113 ln_traadv_fct = .true. ! FCT scheme 114 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 115 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 116 nn_fct_zts = 0 ! > 1 , 2nd order FCT scheme with vertical sub-timestepping 117 ! ! (number of sub-timestep = nn_fct_zts) 117 118 / 118 119 !----------------------------------------------------------------------- … … 123 124 &namtra_ldf ! lateral diffusion scheme for tracers 124 125 !---------------------------------------------------------------------------------- 125 / 126 !----------------------------------------------------------------------- 126 ! ! Operator type: 127 ln_traldf_lap = .true. ! laplacian operator 128 ln_traldf_blp = .false. ! bilaplacian operator 129 ! ! Direction of action: 130 ln_traldf_lev = .false. ! iso-level 131 ln_traldf_hor = .false. ! horizontal (geopotential) 132 ln_traldf_iso = .true. ! iso-neutral (Standard operator) 133 ln_traldf_triad = .false. ! iso-neutral (Triads operator) 134 ! 135 ! ! iso-neutral options: 136 ln_traldf_msc = .true. ! Method of Stabilizing Correction (both operators) 137 rn_slpmax = 0.01 ! slope limit (both operators) 138 ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) 139 rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) 140 ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) 141 ! 142 ! ! Coefficients: 143 nn_aht_ijk_t = 20 ! space/time variation of eddy coef 144 ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file 145 ! ! = 0 constant 146 ! ! = 10 F(k) =ldf_c1d 147 ! ! = 20 F(i,j) =ldf_c2d 148 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 149 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d 150 ! ! = 31 F(i,j,k,t)=F(local velocity) 151 rn_aht_0 = 2000. ! lateral eddy diffusivity (lap. operator) [m2/s] 152 rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] 153 / 154 !---------------------------------------------------------------------------------- 155 &namtra_ldfeiv ! eddy induced velocity param. 156 !---------------------------------------------------------------------------------- 157 ln_ldfeiv =.true. ! use eddy induced velocity parameterization 158 ln_ldfeiv_dia =.true. ! diagnose eiv stream function and velocities 159 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] 160 nn_aei_ijk_t = 21 ! space/time variation of the eiv coeficient 161 ! ! =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 162 ! ! = 0 constant 163 ! ! = 10 F(k) =ldf_c1d 164 ! ! = 20 F(i,j) =ldf_c2d 165 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 166 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d 167 /!----------------------------------------------------------------------- 127 168 &namtra_dmp ! tracer: T & S newtonian damping 128 169 !----------------------------------------------------------------------- … … 135 176 &namdyn_vor ! option of physics/algorithm (not control by CPP keys) 136 177 !----------------------------------------------------------------------- 178 ln_dynvor_ene = .false. ! enstrophy conserving scheme 179 ln_dynvor_ens = .false. ! energy conserving scheme 180 ln_dynvor_mix = .false. ! mixed scheme 181 ln_dynvor_een = .true. ! energy & enstrophy scheme 182 nn_een_e3f = 0 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 137 183 / 138 184 !----------------------------------------------------------------------- … … 140 186 !----------------------------------------------------------------------- 141 187 ln_hpg_zco = .false. ! z-coordinate - full steps 142 ln_hpg_zps = .false. ! z-coordinate - partial steps (interpolation)188 ln_hpg_zps = .false. ! z-coordinate - partial steps (interpolation) 143 189 ln_hpg_sco = .false. ! s-coordinate (standard jacobian formulation) 144 190 ln_hpg_djc = .false. ! s-coordinate (Density Jacobian with Cubic polynomial) 145 ln_hpg_prj = .true. 191 ln_hpg_prj = .true. ! s-coordinate (Pressure Jacobian scheme) 146 192 ln_dynhpg_imp = .false. ! time stepping: semi-implicit time scheme (T) 147 193 ! centered time scheme (F) … … 150 196 &namdyn_ldf ! lateral diffusion on momentum 151 197 !----------------------------------------------------------------------- 198 ! ! Type of the operator : 199 ! ! no diffusion: set ln_dynldf_lap=..._blp=F 200 ln_dynldf_lap = .true. ! laplacian operator 201 ln_dynldf_blp = .false. ! bilaplacian operator 202 ! ! Direction of action : 203 ln_dynldf_lev = .true. ! iso-level 204 ln_dynldf_hor = .false. ! horizontal (geopotential) 205 ln_dynldf_iso = .false. ! iso-neutral 206 ! ! Coefficient 207 nn_ahm_ijk_t = -30 ! space/time variation of eddy coef 208 ! ! =-30 read in eddy_viscosity_3D.nc file 209 ! ! =-20 read in eddy_viscosity_2D.nc file 210 ! ! = 0 constant 211 ! ! = 10 F(k)=c1d 212 ! ! = 20 F(i,j)=F(grid spacing)=c2d 213 ! ! = 30 F(i,j,k)=c2d*c1d 214 ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) 215 rn_ahm_0 = 40000. ! horizontal laplacian eddy viscosity [m2/s] 216 rn_ahm_b = 0. ! background eddy viscosity for ldf_iso [m2/s] 217 rn_bhm_0 = 1.e+12 ! horizontal bilaplacian eddy viscosity [m4/s] 218 ! 219 ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 152 220 / 153 221 !----------------------------------------------------------------------- … … 191 259 !----------------------------------------------------------------------- 192 260 / 261 !----------------------------------------------------------------------- 262 &namobs ! observation usage ('key_diaobs') 263 !----------------------------------------------------------------------- 264 / 265 !----------------------------------------------------------------------- 266 &nam_asminc ! assimilation increments ('key_asminc') 267 !----------------------------------------------------------------------- 268 / -
trunk/NEMOGCM/CONFIG/ORCA2_LIM3/cpp_ORCA2_LIM3.fcm
r5385 r5836 1 bld::tool::fppkeys key_trabbl key_lim3 key_vvl key_dynspg_ts key_ diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi key_diaobs key_asminc1 bld::tool::fppkeys key_trabbl key_lim3 key_vvl key_dynspg_ts key_zdftke key_zdfddm key_zdftmx key_mpp_mpi key_diaobs key_asminc key_iomput key_nosignedzero -
trunk/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/1_namelist_cfg
r5407 r5836 5 5 !! namsbc_cpl, namtra_qsr, namsbc_rnf, 6 6 !! namsbc_apr, namsbc_ssr, namsbc_alb) 7 !! 4 - lateral boundary (namlbc, nam cla, namobc, namagrif, nambdy, nambdy_tide)7 !! 4 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) 8 8 !! 5 - bottom boundary (nambfr, nambbc, nambbl) 9 9 !! 6 - Tracer (nameos, namtra_adv, namtra_ldf, namtra_dmp) … … 302 302 !!====================================================================== 303 303 !! namlbc lateral momentum boundary condition 304 !! namcla cross land advection305 !! namobc open boundaries parameters ("key_obc")306 304 !! namagrif agrif nested grid ( read by child model only ) ("key_agrif") 307 305 !! nambdy Unstructured open boundaries ("key_bdy") … … 314 312 rn_shlat = 2. ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat 315 313 ! free slip ! partial slip ! no slip ! strong slip 316 /317 !-----------------------------------------------------------------------318 &namcla ! cross land advection319 !-----------------------------------------------------------------------320 nn_cla = 0 ! advection between 2 ocean pts separates by land321 /322 !-----------------------------------------------------------------------323 &namobc ! open boundaries parameters ("key_obc")324 !-----------------------------------------------------------------------325 ln_obc_clim = .false. ! climatological obc data files (T) or not (F)326 ln_vol_cst = .true. ! impose the total volume conservation (T) or not (F)327 ln_obc_fla = .false. ! Flather open boundary condition328 nn_obcdta = 1 ! = 0 the obc data are equal to the initial state329 ! = 1 the obc data are read in 'obc.dta' files330 cn_obcdta = 'annual' ! set to annual if obc datafile hold 1 year of data331 ! set to monthly if obc datafile hold 1 month of data332 rn_dpein = 1. ! damping time scale for inflow at east open boundary333 rn_dpwin = 1. ! - - - west - -334 rn_dpnin = 1. ! - - - north - -335 rn_dpsin = 1. ! - - - south - -336 rn_dpeob = 3000. ! time relaxation (days) for the east open boundary337 rn_dpwob = 15. ! - - - west - -338 rn_dpnob = 3000. ! - - - north - -339 rn_dpsob = 15. ! - - - south - -340 rn_volemp = 1. ! = 0 the total volume change with the surface flux (E-P-R)341 ! = 1 the total volume remains constant342 314 / 343 315 !----------------------------------------------------------------------- … … 394 366 &nambdy_tide ! tidal forcing at open boundaries 395 367 !----------------------------------------------------------------------- 396 filtide= 'bdydta/amm12_bdytide_' ! file name root of tidal forcing files368 filtide = 'bdydta/amm12_bdytide_' ! file name root of tidal forcing files 397 369 tide_cpt(1) ='Q1' ! names of tidal components used 398 370 tide_cpt(2) ='O1' ! names of tidal components used … … 485 457 / 486 458 !----------------------------------------------------------------------- 487 &namtra_adv ! advection scheme for tracer 488 !----------------------------------------------------------------------- 489 ln_traadv_cen2 = .false. ! 2nd order centered scheme 490 ln_traadv_tvd = .true. ! TVD scheme 491 ln_traadv_muscl = .false. ! MUSCL scheme 492 ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 493 ln_traadv_ubs = .false. ! UBS scheme 494 ln_traadv_qck = .false. ! QUCIKEST scheme 459 &namtra_adv ! advection scheme for tracer 460 !----------------------------------------------------------------------- 461 ln_traadv_fct = .true. ! FCT scheme 462 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 463 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 464 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping 465 ! ! (number of sub-timestep = nn_fct_zts) 495 466 / 496 467 !----------------------------------------------------------------------- 497 468 &namtra_ldf ! lateral diffusion scheme for tracer 498 469 !----------------------------------------------------------------------- 499 ! ! Type of the operator : 500 ln_traldf_lap = .true. ! laplacian operator 501 ln_traldf_bilap = .false. ! bilaplacian operator 502 ! ! Direction of action : 503 ln_traldf_level = .false. ! iso-level 504 ln_traldf_hor = .false. ! horizontal (geopotential) (require "key_ldfslp" when ln_sco=T) 505 ln_traldf_iso = .true. ! iso-neutral (require "key_ldfslp") 506 ln_traldf_grif = .false. ! griffies skew flux formulation (require "key_ldfslp") ! UNDER TEST, DO NOT USE 507 ln_traldf_gdia = .false. ! griffies operator strfn diagnostics (require "key_ldfslp") ! UNDER TEST, DO NOT USE 508 ! ! Coefficient 509 rn_aht_0 = 1000. ! horizontal eddy diffusivity for tracers [m2/s] 510 rn_ahtb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 511 rn_aeiv_0 = 0. ! eddy induced velocity coefficient [m2/s] (require "key_traldf_eiv") 470 ! ! Operator type: 471 ! ! no diffusion: set ln_traldf_lap=..._blp=F 472 ln_traldf_lap = .true. ! laplacian operator 473 ln_traldf_blp = .false. ! bilaplacian operator 474 ! ! Direction of action: 475 ln_traldf_lev = .false. ! iso-level 476 ln_traldf_hor = .false. ! horizontal (geopotential) 477 ln_traldf_iso = .true. ! iso-neutral (standard operator) 478 ln_traldf_triad = .false. ! iso-neutral (triad operator) 479 ! 480 ! ! iso-neutral options: 481 ln_traldf_msc = .true. ! Method of Stabilizing Correction (both operators) 482 rn_slpmax = 0.01 ! slope limit (both operators) 483 ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) 484 rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) 485 ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) 486 ! 487 ! ! Coefficients: 488 nn_aht_ijk_t = 0 ! space/time variation of eddy coef 489 ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file 490 ! ! = 0 constant 491 ! ! = 10 F(k) =ldf_c1d 492 ! ! = 20 F(i,j) =ldf_c2d 493 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 494 ! ! = 30 F(i,j,k) =ldf_c2d * ldf_c1d 495 ! ! = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 496 rn_aht_0 = 1000. ! lateral eddy diffusivity (lap. operator) [m2/s] 497 rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] 498 / 499 !---------------------------------------------------------------------------------- 500 &namtra_ldfeiv ! eddy induced velocity param. 501 !---------------------------------------------------------------------------------- 502 ln_ldfeiv =.false. ! use eddy induced velocity parameterization 512 503 / 513 504 !----------------------------------------------------------------------- … … 546 537 &namdyn_vor ! option of physics/algorithm (not control by CPP keys) 547 538 !----------------------------------------------------------------------- 548 ln_dynvor_ene = .false. ! enstrophy conserving scheme 549 ln_dynvor_ens = .false. ! energy conserving scheme 550 ln_dynvor_mix = .false. ! mixed scheme 551 ln_dynvor_een = .true. ! energy & enstrophy scheme 539 ln_dynvor_ene = .false. ! enstrophy conserving scheme 540 ln_dynvor_ens = .false. ! energy conserving scheme 541 ln_dynvor_mix = .false. ! mixed scheme 542 ln_dynvor_een = .true. ! energy & enstrophy scheme 543 nn_een_e3f = 1 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 544 ln_dynvor_msk = .true. ! vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) 552 545 / 553 546 !----------------------------------------------------------------------- … … 572 565 &namdyn_ldf ! lateral diffusion on momentum 573 566 !----------------------------------------------------------------------- 574 ! ! Type of the operator : 575 ln_dynldf_lap = .false. ! laplacian operator 576 ln_dynldf_bilap = .true. ! bilaplacian operator 577 ! ! Direction of action : 578 ln_dynldf_level = .false. ! iso-level 579 ln_dynldf_hor = .true. ! horizontal (geopotential) (require "key_ldfslp" in s-coord.) 580 ln_dynldf_iso = .false. ! iso-neutral (require "key_ldfslp") 567 ! ! Type of the operator : 568 ! ! no diffusion: set ln_dynldf_lap=..._blp=F 569 ln_dynldf_lap = .false. ! laplacian operator 570 ln_dynldf_blp = .true. ! bilaplacian operator 571 ! ! Direction of action : 572 ln_dynldf_lev = .true. ! iso-level 573 ln_dynldf_hor = .false. ! horizontal (geopotential) 574 ln_dynldf_iso = .false. ! iso-neutral 581 575 ! ! Coefficient 582 rn_ahm_0_lap = 40000. ! horizontal laplacian eddy viscosity [m2/s] 583 rn_ahmb_0 = 0. ! background eddy viscosity for ldf_iso [m2/s] 584 rn_ahm_0_blp = -8.5e+11 ! horizontal bilaplacian eddy viscosity [m4/s] 576 nn_ahm_ijk_t = 0 ! space/time variation of eddy coef 577 ! ! =-30 read in eddy_viscosity_3D.nc file 578 ! ! =-20 read in eddy_viscosity_2D.nc file 579 ! ! = 0 constant 580 ! ! = 10 F(k)=c1d 581 ! ! = 20 F(i,j)=F(grid spacing)=c2d 582 ! ! = 30 F(i,j,k)=c2d*c1d 583 ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) 584 rn_ahm_0 = 40000. ! horizontal laplacian eddy viscosity [m2/s] 585 rn_ahm_b = 0. ! background eddy viscosity for ldf_iso [m2/s] 586 rn_bhm_0 = 8.5e+11 ! horizontal bilaplacian eddy viscosity [m4/s] 585 587 / 586 588 -
trunk/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/1_namelist_top_cfg
r4147 r5836 39 39 &namtrc_adv ! advection scheme for passive tracer 40 40 !----------------------------------------------------------------------- 41 ln_trcadv_cen2 = .false. ! 2nd order centered scheme 42 ln_trcadv_tvd = .true. ! TVD scheme 43 ln_trcadv_muscl = .false. ! MUSCL scheme 44 ln_trcadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 45 ln_trcadv_ubs = .false. ! UBS scheme 46 ln_trcadv_qck = .false. ! QUICKEST scheme 41 ln_trcadv_fct = .true. ! FCT scheme 42 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 43 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 44 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping 45 ! ! (number of sub-timestep = nn_fct_zts) 47 46 / 48 47 !----------------------------------------------------------------------- -
trunk/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/namelist_cfg
r5407 r5836 132 132 / 133 133 !----------------------------------------------------------------------- 134 &namcla ! cross land advection135 !-----------------------------------------------------------------------136 /137 !-----------------------------------------------------------------------138 &namobc ! open boundaries parameters ("key_obc")139 !-----------------------------------------------------------------------140 /141 !-----------------------------------------------------------------------142 134 &namagrif ! AGRIF zoom ("key_agrif") 143 135 !----------------------------------------------------------------------- … … 176 168 / 177 169 !----------------------------------------------------------------------- 178 &namtra_adv ! advection scheme for tracer 179 !----------------------------------------------------------------------- 170 &namtra_adv ! advection scheme for tracer 171 !----------------------------------------------------------------------- 172 ln_traadv_fct = .true. ! FCT scheme 173 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 174 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 175 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping 176 ! ! (number of sub-timestep = nn_fct_zts) 180 177 / 181 178 !---------------------------------------------------------------------------------- 182 179 &namtra_ldf ! lateral diffusion scheme for tracers 183 180 !---------------------------------------------------------------------------------- 181 ! ! Operator type: 182 ! ! no diffusion: set ln_traldf_lap=..._blp=F 183 ln_traldf_lap = .true. ! laplacian operator 184 ln_traldf_blp = .false. ! bilaplacian operator 185 ! ! Direction of action: 186 ln_traldf_lev = .false. ! iso-level 187 ln_traldf_hor = .false. ! horizontal (geopotential) 188 ln_traldf_iso = .true. ! iso-neutral (standard operator) 189 ln_traldf_triad = .false. ! iso-neutral (triad operator) 190 ! 191 ! ! iso-neutral options: 192 ln_traldf_msc = .true. ! Method of Stabilizing Correction (both operators) 193 rn_slpmax = 0.01 ! slope limit (both operators) 194 ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) 195 rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) 196 ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) 197 ! 198 ! ! Coefficients: 199 nn_aht_ijk_t = 20 ! space/time variation of eddy coef 200 ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file 201 ! ! = 0 constant 202 ! ! = 10 F(k) =ldf_c1d 203 ! ! = 20 F(i,j) =ldf_c2d 204 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 205 ! ! = 30 F(i,j,k) =ldf_c2d * ldf_c1d 206 ! ! = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 207 rn_aht_0 = 2000. ! lateral eddy diffusivity (lap. operator) [m2/s] 208 rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] 209 / 210 !---------------------------------------------------------------------------------- 211 &namtra_ldfeiv ! eddy induced velocity param. 212 !---------------------------------------------------------------------------------- 213 ln_ldfeiv =.true. ! use eddy induced velocity parameterization 214 ln_ldfeiv_dia =.false. ! diagnose eiv stream function and velocities 215 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] 216 nn_aei_ijk_t = 21 ! space/time variation of the eiv coeficient 217 ! ! =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 218 ! ! = 0 constant 219 ! ! = 10 F(k) =ldf_c1d 220 ! ! = 20 F(i,j) =ldf_c2d 221 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 222 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d 184 223 / 185 224 !----------------------------------------------------------------------- … … 188 227 / 189 228 !----------------------------------------------------------------------- 229 &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param) 230 !----------------------------------------------------------------------- 231 / 232 !----------------------------------------------------------------------- 190 233 &namdyn_adv ! formulation of the momentum advection 191 234 !----------------------------------------------------------------------- 192 235 / 193 236 !----------------------------------------------------------------------- 194 &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param)195 !-----------------------------------------------------------------------196 /197 !-----------------------------------------------------------------------198 237 &namdyn_vor ! option of physics/algorithm (not control by CPP keys) 199 238 !----------------------------------------------------------------------- 239 ln_dynvor_ene = .false. ! enstrophy conserving scheme 240 ln_dynvor_ens = .false. ! energy conserving scheme 241 ln_dynvor_mix = .false. ! mixed scheme 242 ln_dynvor_een = .true. ! energy & enstrophy scheme 243 nn_een_e3f = 1 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 200 244 / 201 245 !----------------------------------------------------------------------- … … 209 253 &namdyn_ldf ! lateral diffusion on momentum 210 254 !----------------------------------------------------------------------- 255 ! ! Type of the operator : 256 ! ! no diffusion: set ln_dynldf_lap=..._blp=F 257 ln_dynldf_lap = .true. ! laplacian operator 258 ln_dynldf_blp = .false. ! bilaplacian operator 259 ! ! Direction of action : 260 ln_dynldf_lev = .true. ! iso-level 261 ln_dynldf_hor = .false. ! horizontal (geopotential) 262 ln_dynldf_iso = .false. ! iso-neutral 263 ! ! Coefficient 264 nn_ahm_ijk_t = -30 ! space/time variation of eddy coef 265 ! ! =-30 read in eddy_viscosity_3D.nc file 266 ! ! =-20 read in eddy_viscosity_2D.nc file 267 ! ! = 0 constant 268 ! ! = 10 F(k)=c1d 269 ! ! = 20 F(i,j)=F(grid spacing)=c2d 270 ! ! = 30 F(i,j,k)=c2d*c1d 271 ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) 272 rn_ahm_0 = 40000. ! horizontal laplacian eddy viscosity [m2/s] 273 rn_ahm_b = 0. ! background eddy viscosity for ldf_iso [m2/s] 274 rn_bhm_0 = 1.e+12 ! horizontal bilaplacian eddy viscosity [m4/s] 211 275 / 212 276 !----------------------------------------------------------------------- … … 221 285 &namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion ("key_zdftke") 222 286 !----------------------------------------------------------------------- 223 /224 !------------------------------------------------------------------------225 &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionally:226 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb")227 287 / 228 288 !----------------------------------------------------------------------- -
trunk/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/namelist_top_cfg
r4340 r5836 24 24 &namtrc_adv ! advection scheme for passive tracer 25 25 !----------------------------------------------------------------------- 26 ln_trcadv_tvd = .true. ! TVD scheme 27 ln_trcadv_muscl = .false. ! MUSCL scheme 26 ln_trcadv_fct = .true. ! FCT scheme 27 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 28 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 29 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping 30 ! ! (number of sub-timestep = nn_fct_zts) 28 31 / 29 32 !----------------------------------------------------------------------- -
trunk/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/cpp_ORCA2_LIM_CFC_C14b.fcm
r4523 r5836 1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_ diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_top key_cfc key_c14b key_iomput key_mpp_mpi1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_zdftke key_zdfddm key_zdftmx key_top key_cfc key_c14b key_iomput key_mpp_mpi -
trunk/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_cfg
r4370 r5836 93 93 / 94 94 !----------------------------------------------------------------------- 95 &namcla ! cross land advection96 !-----------------------------------------------------------------------97 /98 !-----------------------------------------------------------------------99 95 &nambfr ! bottom friction 100 96 !----------------------------------------------------------------------- … … 115 111 &namtra_adv ! advection scheme for tracer 116 112 !----------------------------------------------------------------------- 113 ln_traadv_fct = .true. ! FCT scheme 114 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 115 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 116 nn_fct_zts = 0 ! > 1 , 2nd order FCT scheme with vertical sub-timestepping 117 ! ! (number of sub-timestep = nn_fct_zts) 117 118 / 118 119 !----------------------------------------------------------------------- … … 120 121 !----------------------------------------------------------------------- 121 122 / 122 !-----------------------------------------------------------------------123 &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param)124 !-----------------------------------------------------------------------125 /126 123 !---------------------------------------------------------------------------------- 127 124 &namtra_ldf ! lateral diffusion scheme for tracers 128 125 !---------------------------------------------------------------------------------- 126 ! ! Operator type: 127 ln_traldf_lap = .true. ! laplacian operator 128 ln_traldf_blp = .false. ! bilaplacian operator 129 ! ! Direction of action: 130 ln_traldf_lev = .false. ! iso-level 131 ln_traldf_hor = .false. ! horizontal (geopotential) 132 ln_traldf_iso = .true. ! iso-neutral (standard operator) 133 ln_traldf_triad = .false. ! iso-neutral (triad operator) 134 ! 135 ! ! iso-neutral options: 136 ln_traldf_msc = .true. ! Method of Stabilizing Correction (both operators) 137 rn_slpmax = 0.01 ! slope limit (both operators) 138 ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) 139 rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) 140 ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) 141 ! 142 ! ! Coefficients: 143 nn_aht_ijk_t = 20 ! space/time variation of eddy coef 144 ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file 145 ! ! = 0 constant 146 ! ! = 10 F(k) =ldf_c1d 147 ! ! = 20 F(i,j) =ldf_c2d 148 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 149 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d 150 ! ! = 31 F(i,j,k,t)=F(local velocity) 151 rn_aht_0 = 2000. ! lateral eddy diffusivity (lap. operator) [m2/s] 152 rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] 153 / 154 !---------------------------------------------------------------------------------- 155 &namtra_ldfeiv ! eddy induced velocity param. 156 !---------------------------------------------------------------------------------- 157 ln_ldfeiv =.true. ! use eddy induced velocity parameterization 158 ln_ldfeiv_dia =.true. ! diagnose eiv stream function and velocities 159 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] 160 nn_aei_ijk_t = 21 ! space/time variation of the eiv coeficient 161 ! ! =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 162 ! ! = 0 constant 163 ! ! = 10 F(k) =ldf_c1d 164 ! ! = 20 F(i,j) =ldf_c2d 165 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 166 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d 129 167 / 130 168 !----------------------------------------------------------------------- … … 139 177 &namdyn_vor ! option of physics/algorithm (not control by CPP keys) 140 178 !----------------------------------------------------------------------- 179 ln_dynvor_ene = .false. ! enstrophy conserving scheme 180 ln_dynvor_ens = .false. ! energy conserving scheme 181 ln_dynvor_mix = .false. ! mixed scheme 182 ln_dynvor_een = .true. ! energy & enstrophy scheme 183 nn_een_e3f = 0 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 141 184 / 142 185 !----------------------------------------------------------------------- … … 147 190 &namdyn_ldf ! lateral diffusion on momentum 148 191 !----------------------------------------------------------------------- 192 ! ! Type of the operator : 193 ! ! no diffusion: set ln_dynldf_lap=..._blp=F 194 ln_dynldf_lap = .true. ! laplacian operator 195 ln_dynldf_blp = .false. ! bilaplacian operator 196 ! ! Direction of action : 197 ln_dynldf_lev = .true. ! iso-level 198 ln_dynldf_hor = .false. ! horizontal (geopotential) 199 ln_dynldf_iso = .false. ! iso-neutral 200 ! ! Coefficient 201 nn_ahm_ijk_t = -30 ! space/time variation of eddy coef 202 ! ! =-30 read in eddy_viscosity_3D.nc file 203 ! ! =-20 read in eddy_viscosity_2D.nc file 204 ! ! = 0 constant 205 ! ! = 10 F(k)=c1d 206 ! ! = 20 F(i,j)=F(grid spacing)=c2d 207 ! ! = 30 F(i,j,k)=c2d*c1d 208 ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) 209 rn_ahm_0 = 40000. ! horizontal laplacian eddy viscosity [m2/s] 210 rn_ahm_b = 0. ! background eddy viscosity for ldf_iso [m2/s] 211 rn_bhm_0 = 1.e+12 ! horizontal bilaplacian eddy viscosity [m4/s] 149 212 / 150 213 !----------------------------------------------------------------------- -
trunk/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_top_cfg
r5385 r5836 62 62 &namtrc_adv ! advection scheme for passive tracer 63 63 !----------------------------------------------------------------------- 64 ln_trcadv_ tvd = .false. ! TVDscheme65 ln_trcadv_muscl = .true. ! MUSCL scheme64 ln_trcadv_mus = .true. ! MUSCL scheme 65 ln_mus_ups = .false. ! use upstream scheme near river mouths 66 66 / 67 67 !----------------------------------------------------------------------- -
trunk/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/cpp_ORCA2_LIM_PISCES.fcm
r4523 r5836 1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_ diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_top key_pisces key_iomput key_mpp_mpi1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_zdftke key_zdfddm key_zdftmx key_top key_pisces key_mpp_mpi key_iomput -
trunk/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist_cfg
r5385 r5836 78 78 !----------------------------------------------------------------------- 79 79 / 80 81 !-----------------------------------------------------------------------82 &namcla ! cross land advection83 !-----------------------------------------------------------------------84 /85 80 !----------------------------------------------------------------------- 86 81 &nambbl ! bottom boundary layer scheme … … 91 86 !----------------------------------------------------------------------- 92 87 / 93 !----------------------------------------------------------------------- 94 &namtra_ldf ! lateral diffusion scheme for tracer 95 !----------------------------------------------------------------------- 96 ln_triad_iso = .true. ! griffies operator calculates triads twice => pure lateral mixing in ML (require "key_ldfslp") 88 !---------------------------------------------------------------------------------- 89 &namtra_ldf ! lateral diffusion scheme for tracers 90 !---------------------------------------------------------------------------------- 91 ! ! Operator type: 92 ln_traldf_lap = .true. ! laplacian operator 93 ln_traldf_blp = .false. ! bilaplacian operator 94 ! ! Direction of action: 95 ln_traldf_lev = .false. ! iso-level 96 ln_traldf_hor = .false. ! horizontal (geopotential) 97 ln_traldf_iso = .true. ! iso-neutral 98 ln_traldf_triad = .false. ! iso-neutral using Griffies triads 99 ! 100 ! ! iso-neutral options: 101 ln_traldf_msc = .true. ! Method of Stabilizing Correction (both operators) 102 rn_slpmax = 0.01 ! slope limit (both operators) 103 ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) 104 rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) 105 ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) 106 ! 107 ! ! Coefficients: 108 nn_aht_ijk_t = 0 ! space/time variation of eddy coef 109 ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file 110 ! ! = 0 constant 111 ! ! = 10 F(k) =ldf_c1d 112 ! ! = 20 F(i,j) =ldf_c2d 113 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 114 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d 115 ! ! = 31 F(i,j,k,t)=F(local velocity) 116 rn_aht_0 = 2000. ! lateral eddy diffusivity (lap. operator) [m2/s] 117 rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] 118 / 119 !---------------------------------------------------------------------------------- 120 &namtra_ldfeiv ! eddy induced velocity param. 121 !---------------------------------------------------------------------------------- 122 ln_ldfeiv =.true. ! use eddy induced velocity parameterization 123 ln_ldfeiv_dia =.false. ! diagnose eiv stream function and velocities 124 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] 125 nn_aei_ijk_t = 0 ! space/time variation of the eiv coeficient 126 ! ! =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 127 ! ! = 0 constant 128 ! ! = 10 F(k) =ldf_c1d 129 ! ! = 20 F(i,j) =ldf_c2d 130 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 131 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d 97 132 / 98 133 !----------------------------------------------------------------------- … … 115 150 sn_ubl = 'dyna_grid_U' , 120 , 'sobblcox' , .true. , .true. , 'yearly' , '' , '' , '' 116 151 sn_vbl = 'dyna_grid_V' , 120 , 'sobblcoy' , .true. , .true. , 'yearly' , '' , '' , '' 117 sn_ahu = 'dyna_grid_U' , 120 , 'vozoahtu' , .true. , .true. , 'yearly' , '' , '' , ''118 sn_ahv = 'dyna_grid_V' , 120 , 'vomeahtv' , .true. , .true. , 'yearly' , '' , '' , ''119 sn_ahw = 'dyna_grid_W' , 120 , 'voveahtz' , .true. , .true. , 'yearly' , '' , '' , ''120 sn_eiu = 'dyna_grid_U' , 120 , 'vozoaeiu' , .true. , .true. , 'yearly' , '' , '' , ''121 sn_eiv = 'dyna_grid_V' , 120 , 'vomeaeiv' , .true. , .true. , 'yearly' , '' , '' , ''122 sn_eiw = 'dyna_grid_W' , 120 , 'soleaeiw' , .true. , .true. , 'yearly' , '' , '' , ''123 152 ! 124 153 cn_dir = './' ! root directory for the location of the dynamical files 125 ln_degrad = .false. ! flag for degradation - requires ("key_degrad")126 154 ln_dynwzv = .true. ! computation of vertical velocity instead of using the one read in file 127 155 ln_dynbbl = .true. ! bbl coef are in files, so read them - requires ("key_trabbl") -
trunk/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist_top_cfg
r5385 r5836 65 65 &namtrc_adv ! advection scheme for passive tracer 66 66 !----------------------------------------------------------------------- 67 ln_trcadv_ tvd = .false. ! TVDscheme68 ln_trcadv_muscl = .true. ! MUSCL scheme67 ln_trcadv_mus = .true. ! MUSCL scheme 68 ln_mus_ups = .false. ! use upstream scheme near river mouths 69 69 / 70 70 !----------------------------------------------------------------------- -
trunk/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/cpp_ORCA2_OFF_PISCES.fcm
r4523 r5836 1 bld::tool::fppkeys key_trabbl key_ ldfslp key_traldf_c2d key_traldf_eiv key_top key_offline key_pisces key_iomput key_mpp_mpi1 bld::tool::fppkeys key_trabbl key_top key_offline key_pisces key_iomput key_mpp_mpi -
trunk/NEMOGCM/CONFIG/ORCA2_SAS_LIM/EXP00/namelist_cfg
r4370 r5836 87 87 / 88 88 !----------------------------------------------------------------------- 89 &namcla ! cross land advection90 !-----------------------------------------------------------------------91 /92 !-----------------------------------------------------------------------93 89 &nameos ! ocean physical parameters 94 90 !----------------------------------------------------------------------- -
trunk/NEMOGCM/CONFIG/ORCA2_SAS_LIM/cpp_ORCA2_SAS_LIM.fcm
r4523 r5836 1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_ diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi -
trunk/NEMOGCM/CONFIG/SHARED/domain_def.xml
r5426 r5836 163 163 164 164 <domain id="grid_V" long_name="grid V"/> 165 165 166 <domain_group id="grid_W"> 166 167 <domain id="grid_W" long_name="grid W"/> … … 168 169 <domain id="EqW" zoom_ibegin="1" zoom_jbegin="0000" zoom_ni="0000" zoom_nj="1" /> 169 170 </domain_group> 171 172 <domain id="grid_F" long_name="grid F"/> 170 173 171 174 <domain_group id="scalarpoint"> -
trunk/NEMOGCM/CONFIG/SHARED/field_def.xml
r5517 r5836 25 25 26 26 <field id="toce" long_name="temperature" standard_name="sea_water_potential_temperature" unit="degC" grid_ref="grid_T_3D"/> 27 <field id="toce_e3t" long_name="temperature * e3t" unit="degC*m"grid_ref="grid_T_3D" > toce * e3t </field >27 <field id="toce_e3t" long_name="temperature (thickness weighted)" unit="degC" grid_ref="grid_T_3D" > toce * e3t </field > 28 28 <field id="soce" long_name="salinity" standard_name="sea_water_practical_salinity" unit="1e-3" grid_ref="grid_T_3D"/> 29 <field id="soce_e3t" long_name="salinity * e3t" unit="1e-3*m" grid_ref="grid_T_3D" > soce * e3t </field > 29 <field id="soce_e3t" long_name="salinity (thickness weighted)" unit="1e-3" grid_ref="grid_T_3D" > soce * e3t </field > 30 31 <!-- t-eddy viscosity coefficients (ldfdyn) --> 32 <field id="ahmt_2d" long_name=" surface t-eddy viscosity coefficient" unit="m2/s or m4/s" /> 33 <field id="ahmt_3d" long_name=" 3D t-eddy viscosity coefficient" unit="m2/s or m4/s" grid_ref="grid_T_3D"/> 30 34 31 35 <field id="sst" long_name="sea surface temperature" standard_name="sea_surface_temperature" unit="degC" /> … … 367 371 <field id="utau" long_name="Wind Stress along i-axis" standard_name="surface_downward_x_stress" unit="N/m2" /> 368 372 <field id="uoce" long_name="ocean current along i-axis" standard_name="sea_water_x_velocity" unit="m/s" grid_ref="grid_U_3D" /> 369 <field id="uoce_e3u" long_name="ocean current along i-axis * e3u" unit="m2/s"grid_ref="grid_U_3D" > uoce * e3u </field>373 <field id="uoce_e3u" long_name="ocean current along i-axis (thickness weighted)" unit="m/s" grid_ref="grid_U_3D" > uoce * e3u </field> 370 374 <field id="ssu" long_name="ocean surface current along i-axis" unit="m/s" /> 371 375 <field id="sbu" long_name="ocean bottom current along i-axis" unit="m/s" /> … … 375 379 <field id="uoces" long_name="ocean transport along i-axis times salinity (CRS)" unit="1e-3*m/s" grid_ref="grid_U_3D" /> 376 380 381 <!-- u-eddy coefficients (ldftra) --> 382 <field id="ahtu_2d" long_name=" surface u-eddy diffusivity coefficient" unit="m2/s or m4/s" /> 383 <field id="ahtu_3d" long_name=" 3D u-EIV coefficient" unit="m2/s or m4/s" grid_ref="grid_U_3D"/> 384 <field id="aeiu_2d" long_name=" surface u-EIV coefficient" unit="m2/s" /> 385 <field id="aeiu_3d" long_name=" 3D u-EIV coefficient" unit="m2/s" grid_ref="grid_U_3D"/> 386 377 387 <!-- variables available with MLE --> 378 388 <field id="psiu_mle" long_name="MLE streamfunction along i-axis" unit="m3/s" grid_ref="grid_U_3D" /> 379 389 380 <!-- uoce_eiv: available with key_traldf_eiv and key_diaeiv-->390 <!-- uoce_eiv: available EIV --> 381 391 <field id="uoce_eiv" long_name="EIV ocean current along i-axis" standard_name="bolus_sea_water_x_velocity" unit="m/s" grid_ref="grid_U_3D" /> 382 392 … … 402 412 <field id="vtau" long_name="Wind Stress along j-axis" standard_name="surface_downward_y_stress" unit="N/m2" /> 403 413 <field id="voce" long_name="ocean current along j-axis" standard_name="sea_water_y_velocity" unit="m/s" grid_ref="grid_V_3D" /> 404 <field id="voce_e3v" long_name="ocean current along j-axis * e3v" unit="m2/s"grid_ref="grid_V_3D" > voce * e3v </field>414 <field id="voce_e3v" long_name="ocean current along j-axis (thickness weighted)" unit="m/s" grid_ref="grid_V_3D" > voce * e3v </field> 405 415 <field id="ssv" long_name="ocean surface current along j-axis" unit="m/s" /> 406 416 <field id="sbv" long_name="ocean bottom current along j-axis" unit="m/s" /> … … 410 420 <field id="voces" long_name="ocean transport along j-axis times salinity (CRS)" unit="1e-3*m/s" grid_ref="grid_V_3D" /> 411 421 422 <!-- v-eddy coefficients (ldftra, ldfdyn) --> 423 <field id="ahtv_2d" long_name=" surface v-eddy diffusivity coefficient" unit="m2/s or (m4/s)^1/2" /> 424 <field id="ahtv_3d" long_name=" 3D v-eddy diffusivity coefficient" unit="m2/s or (m4/s)^1/2" grid_ref="grid_V_3D"/> 425 <field id="aeiv_2d" long_name=" surface v-EIV coefficient" unit="m2/s" /> 426 <field id="aeiv_3d" long_name=" 3D v-EIV coefficient" unit="m2/s" grid_ref="grid_V_3D" /> 427 412 428 <!-- variables available with MLE --> 413 429 <field id="psiv_mle" long_name="MLE streamfunction along j-axis" unit="m3/s" grid_ref="grid_V_3D" /> 414 430 415 <!-- voce_eiv: available with key_traldf_eiv and key_diaeiv-->431 <!-- voce_eiv: available with EIV --> 416 432 <field id="voce_eiv" long_name="EIV ocean current along j-axis" standard_name="bolus_sea_water_y_velocity" unit="m/s" grid_ref="grid_V_3D" /> 417 433 … … 438 454 <field id="wocetr_eff" long_name="effective ocean vertical transport" unit="m3/s" /> 439 455 440 <!-- woce_eiv: available with key_traldf_eiv and key_diaeiv-->456 <!-- woce_eiv: available with EIV --> 441 457 <field id="woce_eiv" long_name="EIV ocean vertical velocity" standard_name="bolus_upward_sea_water_velocity" unit="m/s" /> 442 458 … … 449 465 450 466 <!-- avt_evd and avm_evd: available with ln_zdfevd --> 451 <field id="avt_evd" long_name="convective enhancement tovertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_convection" unit="m2/s" />452 <field id="avm_evd" long_name="convective enhancement tovertical viscosity" standard_name="ocean_vertical_momentum_diffusivity_due_to_convection" unit="m2/s" />467 <field id="avt_evd" long_name="convective enhancement of vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_convection" unit="m2/s" /> 468 <field id="avm_evd" long_name="convective enhancement of vertical viscosity" standard_name="ocean_vertical_momentum_diffusivity_due_to_convection" unit="m2/s" /> 453 469 454 470 <!-- avt_tide: available with key_zdftmx --> … … 459 475 <field id="w_masstr2" long_name="square of vertical mass transport" standard_name="square_of_upward_ocean_mass_transport" unit="kg2/s2" /> 460 476 461 <!-- aht2d and aht2d_eiv : available with key_traldf_eiv and key_traldf_c2d-->477 <!-- aht2d and aht2d_eiv --> 462 478 <field id="aht2d" long_name="lateral eddy diffusivity" standard_name="ocean_tracer_xy_laplacian_diffusivity" unit="m2/s" grid_ref="grid_W_2D" /> 463 479 <field id="aht2d_eiv" long_name="EIV lateral eddy diffusivity" standard_name="ocean_tracer_bolus_laplacian_diffusivity" unit="m2/s" grid_ref="grid_W_2D" /> 464 480 </field_group> 465 481 482 <!-- F grid --> 483 <!-- f-eddy viscosity coefficients (ldfdyn) --> 484 <field id="ahmf_2d" long_name=" surface f-eddy viscosity coefficient" unit="m2/s or m4/s" /> 485 <field id="ahmf_3d" long_name=" 3D f-eddy viscosity coefficient" unit="m2/s or m4/s" grid_ref="grid_T_3D"/> 486 466 487 <!-- scalar variables available with key_diaar5 --> 467 488 -
trunk/NEMOGCM/CONFIG/SHARED/namelist_ref
r5656 r5836 1 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 !! namelist_ref 1 3 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 4 !! NEMO/OPA : 1 - run manager (namrun) … … 5 7 !! namsbc_cpl, namtra_qsr, namsbc_rnf, 6 8 !! namsbc_apr, namsbc_ssr, namsbc_alb) 7 !! 4 - lateral boundary (namlbc, nam cla, namobc, namagrif, nambdy, nambdy_tide)9 !! 4 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) 8 10 !! 5 - bottom boundary (nambfr, nambbc, nambbl) 9 !! 6 - Tracer (nameos, namtra_adv, namtra_ldf, namtra_ dmp)11 !! 6 - Tracer (nameos, namtra_adv, namtra_ldf, namtra_ldfeiv, namtra_dmp) 10 12 !! 7 - dynamics (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 11 !! 8 - Verical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_ kpp, namzdf_ddm, namzdf_tmx)13 !! 8 - Verical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_ddm, namzdf_tmx) 12 14 !! 9 - diagnostics (namnc4, namtrd, namspr, namflo, namhsb, namsto) 13 15 !! 10 - miscellaneous (namsol, nammpp, namctl) … … 200 202 !----------------------------------------------------------------------- 201 203 &namtsd ! data : Temperature & Salinity 202 !-----------------------------------------------------------------------203 204 !----------------------------------------------------------------------- 204 205 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! … … 545 546 !!====================================================================== 546 547 !! namlbc lateral momentum boundary condition 547 !! namcla cross land advection548 548 !! namobc open boundaries parameters ("key_obc") 549 549 !! namagrif agrif nested grid ( read by child model only ) ("key_agrif") … … 558 558 ! free slip ! partial slip ! no slip ! strong slip 559 559 ln_vorlat = .false. ! consistency of vorticity boundary condition with analytical eqs. 560 /561 !-----------------------------------------------------------------------562 &namcla ! cross land advection563 !-----------------------------------------------------------------------564 nn_cla = 0 ! advection between 2 ocean pts separates by land565 /566 !-----------------------------------------------------------------------567 &namobc ! open boundaries parameters ("key_obc")568 !-----------------------------------------------------------------------569 ln_obc_clim = .false. ! climatological obc data files (T) or not (F)570 ln_vol_cst = .true. ! impose the total volume conservation (T) or not (F)571 ln_obc_fla = .false. ! Flather open boundary condition572 nn_obcdta = 1 ! = 0 the obc data are equal to the initial state573 ! = 1 the obc data are read in 'obc.dta' files574 cn_obcdta = 'annual' ! set to annual if obc datafile hold 1 year of data575 ! set to monthly if obc datafile hold 1 month of data576 rn_dpein = 1. ! damping time scale for inflow at east open boundary577 rn_dpwin = 1. ! - - - west - -578 rn_dpnin = 1. ! - - - north - -579 rn_dpsin = 1. ! - - - south - -580 rn_dpeob = 3000. ! time relaxation (days) for the east open boundary581 rn_dpwob = 15. ! - - - west - -582 rn_dpnob = 3000. ! - - - north - -583 rn_dpsob = 15. ! - - - south - -584 rn_volemp = 1. ! = 0 the total volume change with the surface flux (E-P-R)585 ! = 1 the total volume remains constant586 560 / 587 561 !----------------------------------------------------------------------- … … 711 685 ! = 2 variable flux (read in geothermal_heating.nc in mW/m2) 712 686 rn_geoflx_cst = 86.4e-3 ! Constant value of geothermal heat flux [W/m2] 713 714 687 / 715 688 !----------------------------------------------------------------------- … … 725 698 !! Tracer (T & S ) namelists 726 699 !!====================================================================== 727 !! nameos equation of state728 !! namtra_adv advection scheme700 !! nameos equation of state 701 !! namtra_adv advection scheme 729 702 !! namtra_adv_mle mixed layer eddy param. (Fox-Kemper param.) 730 !! namtra_ldf lateral diffusion scheme 731 !! namtra_dmp T & S newtonian damping 703 !! namtra_ldf lateral diffusion scheme 704 !! namtra_ldfeiv eddy induced velocity param. 705 !! namtra_dmp T & S newtonian damping 732 706 !!====================================================================== 733 707 ! … … 740 714 ! = 1, S-EOS (simplified eos) 741 715 ln_useCT = .true. ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 742 !!716 ! 743 717 ! ! S-EOS coefficients : 744 !! rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS718 ! rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 745 719 rn_a0 = 1.6550e-1 ! thermal expension coefficient (nn_eos= 1) 746 720 rn_b0 = 7.6554e-1 ! saline expension coefficient (nn_eos= 1) … … 754 728 &namtra_adv ! advection scheme for tracer 755 729 !----------------------------------------------------------------------- 756 ln_traadv_cen2 = .false. ! 2nd order centered scheme 757 ln_traadv_tvd = .true. ! TVD scheme 758 ln_traadv_muscl = .false. ! MUSCL scheme 759 ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 760 ln_traadv_ubs = .false. ! UBS scheme 761 ln_traadv_qck = .false. ! QUICKEST scheme 762 ln_traadv_msc_ups= .false. ! use upstream scheme within muscl 763 ln_traadv_tvd_zts= .false. ! TVD scheme with sub-timestepping of vertical tracer advection 730 ln_traadv_cen = .false. ! 2nd order centered scheme 731 nn_cen_h = 4 ! =2/4, horizontal 2nd order CEN / 4th order CEN 732 nn_cen_v = 4 ! =2/4, vertical 2nd order CEN / 4th order COMPACT 733 ln_traadv_fct = .false. ! FCT scheme 734 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 735 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 736 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping 737 ! ! (number of sub-timestep = nn_fct_zts) 738 ln_traadv_mus = .false. ! MUSCL scheme 739 ln_mus_ups = .false. ! use upstream scheme near river mouths 740 ln_traadv_ubs = .false. ! UBS scheme 741 nn_ubs_v = 2 ! =2 , vertical 2nd order FCT 742 ln_traadv_qck = .false. ! QUICKEST scheme 764 743 / 765 744 !----------------------------------------------------------------------- 766 745 &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param) 767 746 !----------------------------------------------------------------------- 768 ln_mle = . true. ! (T) use the Mixed Layer Eddy (MLE) parameterisation747 ln_mle = .false. ! (T) use the Mixed Layer Eddy (MLE) parameterisation 769 748 rn_ce = 0.06 ! magnitude of the MLE (typical value: 0.06 to 0.08) 770 749 nn_mle = 1 ! MLE type: =0 standard Fox-Kemper ; =1 new formulation … … 780 759 !---------------------------------------------------------------------------------- 781 760 ! ! Operator type: 782 ln_traldf_lap = .true. ! laplacian operator 783 ln_traldf_bilap = .false. ! bilaplacian operator 761 ! ! no diffusion: set ln_traldf_lap=..._blp=F 762 ln_traldf_lap = .false. ! laplacian operator 763 ln_traldf_blp = .false. ! bilaplacian operator 784 764 ! ! Direction of action: 785 ln_traldf_level = .false. ! iso-level 786 ln_traldf_hor = .false. ! horizontal (geopotential) (needs "key_ldfslp" when ln_sco=T) 787 ln_traldf_iso = .true. ! iso-neutral (needs "key_ldfslp") 788 ! ! Griffies parameters (all need "key_ldfslp") 789 ln_traldf_grif = .false. ! use griffies triads 790 ln_traldf_gdia = .false. ! output griffies eddy velocities 791 ln_triad_iso = .false. ! pure lateral mixing in ML 792 ln_botmix_grif = .false. ! lateral mixing on bottom 793 ! ! Coefficients 794 ! Eddy-induced (GM) advection always used with Griffies; otherwise needs "key_traldf_eiv" 795 ! Value rn_aeiv_0 is ignored unless = 0 with Held-Larichev spatially varying aeiv 796 ! (key_traldf_c2d & key_traldf_eiv & key_orca_r2, _r1 or _r05) 797 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] 798 rn_aht_0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 799 rn_ahtb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 800 ! (normally=0; not used with Griffies) 801 rn_slpmax = 0.01 ! slope limit 802 rn_chsmag = 1. ! multiplicative factor in Smagorinsky diffusivity 803 rn_smsh = 1. ! Smagorinsky diffusivity: = 0 - use only sheer 804 rn_aht_m = 2000. ! upper limit or stability criteria for lateral eddy diffusivity (m2/s) 765 ln_traldf_lev = .false. ! iso-level 766 ln_traldf_hor = .false. ! horizontal (geopotential) 767 ln_traldf_iso = .false. ! iso-neutral (standard operator) 768 ln_traldf_triad = .false. ! iso-neutral (triad operator) 769 ! 770 ! ! iso-neutral options: 771 ln_traldf_msc = .false. ! Method of Stabilizing Correction (both operators) 772 rn_slpmax = 0.01 ! slope limit (both operators) 773 ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) 774 rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) 775 ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) 776 ! 777 ! ! Coefficients: 778 nn_aht_ijk_t = 0 ! space/time variation of eddy coef 779 ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file 780 ! ! = 0 constant 781 ! ! = 10 F(k) =ldf_c1d 782 ! ! = 20 F(i,j) =ldf_c2d 783 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 784 ! ! = 30 F(i,j,k) =ldf_c2d * ldf_c1d 785 ! ! = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 786 rn_aht_0 = 2000. ! lateral eddy diffusivity (lap. operator) [m2/s] 787 rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] 788 / 789 !---------------------------------------------------------------------------------- 790 &namtra_ldfeiv ! eddy induced velocity param. 791 !---------------------------------------------------------------------------------- 792 ln_ldfeiv =.false. ! use eddy induced velocity parameterization 793 ln_ldfeiv_dia =.false. ! diagnose eiv stream function and velocities 794 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] 795 nn_aei_ijk_t = 21 ! space/time variation of the eiv coeficient 796 ! ! =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 797 ! ! = 0 constant 798 ! ! = 10 F(k) =ldf_c1d 799 ! ! = 20 F(i,j) =ldf_c2d 800 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 801 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d 805 802 / 806 803 !----------------------------------------------------------------------- … … 853 850 ln_dynvor_ens = .false. ! energy conserving scheme 854 851 ln_dynvor_mix = .false. ! mixed scheme 855 ln_dynvor_een = .true. ! energy & enstrophy scheme 856 ln_dynvor_een_old = .false. ! energy & enstrophy scheme - original formulation 852 ln_dynvor_een = .false. ! energy & enstrophy scheme 853 nn_een_e3f = 1 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 854 ln_dynvor_msk = .false. ! vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) ! PLEASE DO NOT USE 857 855 / 858 856 !----------------------------------------------------------------------- … … 879 877 !----------------------------------------------------------------------- 880 878 ! ! Type of the operator : 881 ln_dynldf_lap = .true. ! laplacian operator 882 ln_dynldf_bilap = .false. ! bilaplacian operator 879 ! ! no diffusion: set ln_dynldf_lap=..._blp=F 880 ln_dynldf_lap = .false. ! laplacian operator 881 ln_dynldf_blp = .false. ! bilaplacian operator 883 882 ! ! Direction of action : 884 ln_dynldf_lev el = .false.! iso-level885 ln_dynldf_hor = .true. ! horizontal (geopotential) (require "key_ldfslp" in s-coord.)886 ln_dynldf_iso = .false. ! iso-neutral (require "key_ldfslp")883 ln_dynldf_lev = .false. ! iso-level 884 ln_dynldf_hor = .false. ! horizontal (geopotential) 885 ln_dynldf_iso = .false. ! iso-neutral 887 886 ! ! Coefficient 888 rn_ahm_0_lap = 40000. ! horizontal laplacian eddy viscosity [m2/s] 889 rn_ahmb_0 = 0. ! background eddy viscosity for ldf_iso [m2/s] 890 rn_ahm_0_blp = 0. ! horizontal bilaplacian eddy viscosity [m4/s] 891 rn_cmsmag_1 = 3. ! constant in laplacian Smagorinsky viscosity 892 rn_cmsmag_2 = 3 ! constant in bilaplacian Smagorinsky viscosity 893 rn_cmsh = 1. ! 1 or 0 , if 0 -use only shear for Smagorinsky viscosity 894 rn_ahm_m_blp = -1.e12 ! upper limit for bilap abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp) 895 rn_ahm_m_lap = 40000. ! upper limit for lap ahm < min(dx^2/16rdt, rn_ahm_m_lap) 887 nn_ahm_ijk_t = 0 ! space/time variation of eddy coef 888 ! ! =-30 read in eddy_viscosity_3D.nc file 889 ! ! =-20 read in eddy_viscosity_2D.nc file 890 ! ! = 0 constant 891 ! ! = 10 F(k)=c1d 892 ! ! = 20 F(i,j)=F(grid spacing)=c2d 893 ! ! = 30 F(i,j,k)=c2d*c1d 894 ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) 895 rn_ahm_0 = 40000. ! horizontal laplacian eddy viscosity [m2/s] 896 rn_ahm_b = 0. ! background eddy viscosity for ldf_iso [m2/s] 897 rn_bhm_0 = 1.e+12 ! horizontal bilaplacian eddy viscosity [m4/s] 898 ! 899 ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 896 900 / 897 901 … … 902 906 !! namzdf_ric richardson number dependent vertical mixing ("key_zdfric") 903 907 !! namzdf_tke TKE dependent vertical mixing ("key_zdftke") 904 !! namzdf_kpp KPP dependent vertical mixing ("key_zdfkpp")905 908 !! namzdf_ddm double diffusive mixing parameterization ("key_zdfddm") 906 909 !! namzdf_tmx tidal mixing parameterization ("key_zdftmx") … … 963 966 ! = 0 constant 10 m length scale 964 967 ! = 1 0.5m at the equator to 30m poleward of 40 degrees 965 /966 !------------------------------------------------------------------------967 &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionally:968 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb")969 ln_kpprimix = .true. ! shear instability mixing970 rn_difmiw = 1.0e-04 ! constant internal wave viscosity [m2/s]971 rn_difsiw = 0.1e-04 ! constant internal wave diffusivity [m2/s]972 rn_riinfty = 0.8 ! local Richardson Number limit for shear instability973 rn_difri = 0.0050 ! maximum shear mixing at Rig = 0 [m2/s]974 rn_bvsqcon = -0.01e-07 ! Brunt-Vaisala squared for maximum convection [1/s2]975 rn_difcon = 1. ! maximum mixing in interior convection [m2/s]976 nn_avb = 0 ! horizontal averaged (=1) or not (=0) on avt and amv977 nn_ave = 1 ! constant (=0) or profile (=1) background on avt978 968 / 979 969 !----------------------------------------------------------------------- -
trunk/NEMOGCM/CONFIG/SHARED/namelist_top_ref
r5416 r5836 41 41 &namtrc_adv ! advection scheme for passive tracer 42 42 !----------------------------------------------------------------------- 43 ln_trcadv_cen2 = .false. ! 2nd order centered scheme 44 ln_trcadv_tvd = .true. ! TVD scheme 45 ln_trcadv_muscl = .false. ! MUSCL scheme 46 ln_trcadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 47 ln_trcadv_ubs = .false. ! UBS scheme 48 ln_trcadv_qck = .false. ! QUICKEST scheme 49 ln_trcadv_msc_ups = .false. ! use upstream scheme within muscl 43 ln_trcadv_cen = .false. ! 2nd order centered scheme 44 nn_cen_h = 4 ! =2/4, horizontal 2nd order CEN / 4th order CEN 45 nn_cen_v = 4 ! =2/4, vertical 2nd order CEN / 4th order COMPACT 46 ln_trcadv_fct = .false. ! FCT scheme 47 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 48 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 49 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping 50 ! ! (number of sub-timestep = nn_fct_zts) 51 ln_trcadv_mus = .false. ! MUSCL scheme 52 ln_mus_ups = .false. ! use upstream scheme near river mouths 53 ln_trcadv_ubs = .false. ! UBS scheme 54 nn_ubs_v = 2 ! =2 , vertical 2nd order FCT 55 ln_trcadv_qck = .false. ! QUICKEST scheme 50 56 / 51 57 !----------------------------------------------------------------------- 52 58 &namtrc_ldf ! lateral diffusion scheme for passive tracer 53 59 !----------------------------------------------------------------------- 54 ! ! Type of the operator : 55 ln_trcldf_lap = .true. ! laplacian operator 56 ln_trcldf_bilap = .false. ! bilaplacian operator 57 ! Direction of action : 58 ln_trcldf_level = .false. ! iso-level 59 ln_trcldf_hor = .false. ! horizontal (geopotential) (require "key_ldfslp" when ln_sco=T) 60 ln_trcldf_iso = .true. ! iso-neutral (require "key_ldfslp") 61 ! ! Coefficient 62 rn_ahtrc_0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 63 rn_ahtrb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 60 ! ! Type of the operator: 61 ln_trcldf_lap = .true. ! laplacian operator 62 ln_trcldf_blp = .false. ! bilaplacian operator 63 ! ! Direction of action: 64 ln_trcldf_lev = .false. ! iso-level 65 ln_trcldf_hor = .false. ! horizontal (geopotential) 66 ln_trcldf_iso = .true. ! iso-neutral (standard operator) 67 ln_trcldf_triad = .false. ! iso-neutral (triad operator) 68 ! ! Coefficient 69 rn_ahtrc_0 = 2000. ! lateral eddy diffusivity (lap. operator) [m2/s] 70 rn_bhtrc_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] 64 71 / 65 72 !----------------------------------------------------------------------- -
trunk/NEMOGCM/NEMO/LIM_SRC_2/limdmp_2.F90
r4624 r5836 71 71 CALL fld_read( kt, nn_fsbc, sf_icedmp ) 72 72 ! 73 !CDIR COLLAPSE74 73 hicif(:,:) = MAX( 0._wp, & ! h >= 0 avoid spurious out of physical range 75 74 & hicif(:,:) - rdt_ice * resto_ice(:,:,1) * ( hicif(:,:) - sf_icedmp(jp_hicif)%fnow(:,:,1) ) ) 76 !CDIR COLLAPSE77 75 frld (:,:) = MAX( 0._wp, MIN( 1._wp, & ! 0<= frld<=1 values which blow the run up 78 76 & frld (:,:) - rdt_ice * resto_ice(:,:,1) * ( frld (:,:) - sf_icedmp(jp_frld )%fnow(:,:,1) ) ) ) -
trunk/NEMOGCM/NEMO/LIM_SRC_2/limrhg_2.F90
r5123 r5836 160 160 !------------------------------------------------------------------- 161 161 162 !CDIR NOVERRCHK163 162 DO jj = k_j1 , k_jpj-1 164 !CDIR NOVERRCHK165 163 DO ji = 1 , jpi 166 164 ! only the sinus changes its sign with the hemisphere … … 245 243 ! Computation of free drift field for free slip boundary conditions. 246 244 247 !CDIR NOVERRCHK248 245 DO jj = k_j1, k_jpj-1 249 !CDIR NOVERRCHK250 246 DO ji = 1, fs_jpim1 251 247 !- Rate of strain tensor. … … 401 397 iflag: DO jter = 1 , nbitdr ! Relaxation ! 402 398 ! ! ================ ! 403 !CDIR NOVERRCHK404 399 DO jj = k_j1+1, k_jpj-1 405 !CDIR NOVERRCHK406 400 DO ji = 2, fs_jpim1 ! NO vector opt. 407 401 ! -
trunk/NEMOGCM/NEMO/LIM_SRC_2/limsbc_2.F90
r5407 r5836 319 319 ! 320 320 IF( MOD( kt-1, nn_fsbc ) == 0 ) THEN !== Ice time-step only ==! (i.e. surface module time-step) 321 !CDIR NOVERRCHK 321 ! 322 322 DO jj = 1, jpj !* modulus of ice-ocean relative velocity at I-point 323 !CDIR NOVERRCHK324 323 DO ji = 1, jpi 325 324 zu_i = u_ice(ji,jj) - u_oce(ji,jj) ! ice-ocean relative velocity at I-point … … 328 327 END DO 329 328 END DO 330 !CDIR NOVERRCHK331 329 DO jj = 1, jpjm1 !* update the modulus of stress at ocean surface (T-point) 332 !CDIR NOVERRCHK333 330 DO ji = 1, jpim1 ! NO vector opt. 334 331 ! ! modulus of U_ice-U_oce at T-point … … 383 380 ! 384 381 IF( MOD( kt-1, nn_fsbc ) == 0 ) THEN !== Ice time-step only ==! (i.e. surface module time-step) 385 !CDIR NOVERRCHK 382 ! 386 383 DO jj = 2, jpjm1 !* modulus of the ice-ocean velocity at T-point 387 !CDIR NOVERRCHK388 384 DO ji = fs_2, fs_jpim1 389 385 zu_t = u_ice(ji,jj) + u_ice(ji-1,jj) - u_oce(ji,jj) - u_oce(ji-1,jj) ! 2*(U_ice-U_oce) at T-point -
trunk/NEMOGCM/NEMO/LIM_SRC_2/limthd_2.F90
r5407 r5836 196 196 !-------------------------------------------------------------------------- 197 197 198 !CDIR NOVERRCHK199 198 DO jj = 1, jpj 200 !CDIR NOVERRCHK201 199 DO ji = 1, jpi 202 200 zthsnice = hsnif(ji,jj) + hicif(ji,jj) -
trunk/NEMOGCM/NEMO/LIM_SRC_2/limthd_lac_2.F90
r3625 r5836 134 134 !--------------------------------------------------------------------- 135 135 136 !CDIR NOVERRCHK137 136 DO ji = kideb , kiut 138 137 iicefr = 1 - MAX( 0, INT( SIGN( 1.5 * zone , zfrl_old(ji) - 1.0 + epsi13 ) ) ) -
trunk/NEMOGCM/NEMO/LIM_SRC_3/limadv.F90
r5429 r5836 97 97 98 98 ! Initialize volumes of boxes (=area if adv_x first called, =psm otherwise) 99 psm (:,:) = MAX( pcrh * e1 2t(:,:) + ( 1.0 - pcrh ) * psm(:,:) , epsi20 )99 psm (:,:) = MAX( pcrh * e1e2t(:,:) + ( 1.0 - pcrh ) * psm(:,:) , epsi20 ) 100 100 101 101 ! Calculate fluxes and moments between boxes i<-->i+1 … … 282 282 283 283 ! Initialize volumes of boxes (=area if adv_x first called, =psm otherwise) 284 psm(:,:) = MAX( pcrh * e1 2t(:,:) + ( 1.0 - pcrh ) * psm(:,:) , epsi20 )284 psm(:,:) = MAX( pcrh * e1e2t(:,:) + ( 1.0 - pcrh ) * psm(:,:) , epsi20 ) 285 285 286 286 ! Calculate fluxes and moments between boxes j<-->j+1 -
trunk/NEMOGCM/NEMO/LIM_SRC_3/limcons.F90
r5183 r5836 185 185 zfs_b = glob_sum( ( sfx_bri(:,:) + sfx_bog(:,:) + sfx_bom(:,:) + sfx_sum(:,:) + sfx_sni(:,:) + & 186 186 & sfx_opw(:,:) + sfx_res(:,:) + sfx_dyn(:,:) & 187 & ) * e1 2t(:,:) * tmask(:,:,1) * zconv )187 & ) * e1e2t(:,:) * tmask(:,:,1) * zconv ) 188 188 189 189 ! water flux 190 190 zfw_b = glob_sum( -( wfx_bog(:,:) + wfx_bom(:,:) + wfx_sum(:,:) + wfx_sni(:,:) + wfx_opw(:,:) + & 191 191 & wfx_res(:,:) + wfx_dyn(:,:) + wfx_snw(:,:) + wfx_sub(:,:) + wfx_spr(:,:) & 192 & ) * e1 2t(:,:) * tmask(:,:,1) * zconv )192 & ) * e1e2t(:,:) * tmask(:,:,1) * zconv ) 193 193 194 194 ! heat flux 195 195 zft_b = glob_sum( ( hfx_sum(:,:) + hfx_bom(:,:) + hfx_bog(:,:) + hfx_dif(:,:) + hfx_opw(:,:) + hfx_snw(:,:) & 196 196 & - hfx_thd(:,:) - hfx_dyn(:,:) - hfx_res(:,:) - hfx_sub(:,:) - hfx_spr(:,:) & 197 & ) * e1 2t(:,:) * tmask(:,:,1) * zconv )198 199 zvi_b = glob_sum( SUM( v_i * rhoic + v_s * rhosn, dim=3 ) * e1 2t * tmask(:,:,1) * zconv )200 201 zsmv_b = glob_sum( SUM( smv_i * rhoic , dim=3 ) * e1 2t * tmask(:,:,1) * zconv )197 & ) * e1e2t(:,:) * tmask(:,:,1) * zconv ) 198 199 zvi_b = glob_sum( SUM( v_i * rhoic + v_s * rhosn, dim=3 ) * e1e2t * tmask(:,:,1) * zconv ) 200 201 zsmv_b = glob_sum( SUM( smv_i * rhoic , dim=3 ) * e1e2t * tmask(:,:,1) * zconv ) 202 202 203 203 zei_b = glob_sum( ( SUM( SUM( e_i(:,:,1:nlay_i,:), dim=4 ), dim=3 ) + & 204 204 & SUM( SUM( e_s(:,:,1:nlay_s,:), dim=4 ), dim=3 ) & 205 ) * e1 2t * tmask(:,:,1) * zconv )205 ) * e1e2t * tmask(:,:,1) * zconv ) 206 206 207 207 ELSEIF( icount == 1 ) THEN … … 210 210 zfs = glob_sum( ( sfx_bri(:,:) + sfx_bog(:,:) + sfx_bom(:,:) + sfx_sum(:,:) + sfx_sni(:,:) + & 211 211 & sfx_opw(:,:) + sfx_res(:,:) + sfx_dyn(:,:) & 212 & ) * e1 2t(:,:) * tmask(:,:,1) * zconv ) - zfs_b212 & ) * e1e2t(:,:) * tmask(:,:,1) * zconv ) - zfs_b 213 213 214 214 ! water flux 215 215 zfw = glob_sum( -( wfx_bog(:,:) + wfx_bom(:,:) + wfx_sum(:,:) + wfx_sni(:,:) + wfx_opw(:,:) + & 216 216 & wfx_res(:,:) + wfx_dyn(:,:) + wfx_snw(:,:) + wfx_sub(:,:) + wfx_spr(:,:) & 217 & ) * e1 2t(:,:) * tmask(:,:,1) * zconv ) - zfw_b217 & ) * e1e2t(:,:) * tmask(:,:,1) * zconv ) - zfw_b 218 218 219 219 ! heat flux 220 220 zft = glob_sum( ( hfx_sum(:,:) + hfx_bom(:,:) + hfx_bog(:,:) + hfx_dif(:,:) + hfx_opw(:,:) + hfx_snw(:,:) & 221 221 & - hfx_thd(:,:) - hfx_dyn(:,:) - hfx_res(:,:) - hfx_sub(:,:) - hfx_spr(:,:) & 222 & ) * e1 2t(:,:) * tmask(:,:,1) * zconv ) - zft_b222 & ) * e1e2t(:,:) * tmask(:,:,1) * zconv ) - zft_b 223 223 224 224 ! outputs 225 225 zvi = ( ( glob_sum( SUM( v_i * rhoic + v_s * rhosn, dim=3 ) & 226 & * e1 2t * tmask(:,:,1) * zconv ) - zvi_b ) * r1_rdtice - zfw ) * rday226 & * e1e2t * tmask(:,:,1) * zconv ) - zvi_b ) * r1_rdtice - zfw ) * rday 227 227 228 228 zsmv = ( ( glob_sum( SUM( smv_i * rhoic , dim=3 ) & 229 & * e1 2t * tmask(:,:,1) * zconv ) - zsmv_b ) * r1_rdtice + zfs ) * rday229 & * e1e2t * tmask(:,:,1) * zconv ) - zsmv_b ) * r1_rdtice + zfs ) * rday 230 230 231 231 zei = glob_sum( ( SUM( SUM( e_i(:,:,1:nlay_i,:), dim=4 ), dim=3 ) + & 232 232 & SUM( SUM( e_s(:,:,1:nlay_s,:), dim=4 ), dim=3 ) & 233 & ) * e1 2t * tmask(:,:,1) * zconv ) * r1_rdtice - zei_b * r1_rdtice + zft233 & ) * e1e2t * tmask(:,:,1) * zconv ) * r1_rdtice - zei_b * r1_rdtice + zft 234 234 235 235 ! zvtrp and zetrp must be close to 0 if the advection scheme is conservative 236 zvtrp = glob_sum( ( diag_trp_vi * rhoic + diag_trp_vs * rhosn ) * e1 2t * tmask(:,:,1) * zconv ) * rday237 zetrp = glob_sum( ( diag_trp_ei + diag_trp_es ) * e1 2t * tmask(:,:,1) * zconv )236 zvtrp = glob_sum( ( diag_trp_vi * rhoic + diag_trp_vs * rhosn ) * e1e2t * tmask(:,:,1) * zconv ) * rday 237 zetrp = glob_sum( ( diag_trp_ei + diag_trp_es ) * e1e2t * tmask(:,:,1) * zconv ) 238 238 239 239 zvmin = glob_min( v_i ) … … 242 242 243 243 ! set threshold values and calculate the ice area (+epsi10 to set a threshold > 0 when there is no ice) 244 zarea = glob_sum( SUM( a_i + epsi10, dim=3 ) * e1 2t * zconv ) ! in 1.e9 m2244 zarea = glob_sum( SUM( a_i + epsi10, dim=3 ) * e1e2t * zconv ) ! in 1.e9 m2 245 245 zv_sill = zarea * 2.5e-5 246 246 zs_sill = zarea * 25.e-5 … … 286 286 #if ! defined key_bdy 287 287 ! heat flux 288 zhfx = glob_sum( ( hfx_in - hfx_out - diag_heat - diag_trp_ei - diag_trp_es - hfx_sub ) * e1 2t * tmask(:,:,1) * zconv )288 zhfx = glob_sum( ( hfx_in - hfx_out - diag_heat - diag_trp_ei - diag_trp_es - hfx_sub ) * e1e2t * tmask(:,:,1) * zconv ) 289 289 ! salt flux 290 zsfx = glob_sum( ( sfx + diag_smvi ) * e1 2t * tmask(:,:,1) * zconv ) * rday290 zsfx = glob_sum( ( sfx + diag_smvi ) * e1e2t * tmask(:,:,1) * zconv ) * rday 291 291 ! water flux 292 zvfx = glob_sum( ( wfx_ice + wfx_snw + wfx_spr + wfx_sub + diag_vice + diag_vsnw ) * e1 2t * tmask(:,:,1) * zconv ) * rday292 zvfx = glob_sum( ( wfx_ice + wfx_snw + wfx_spr + wfx_sub + diag_vice + diag_vsnw ) * e1e2t * tmask(:,:,1) * zconv ) * rday 293 293 294 294 ! set threshold values and calculate the ice area (+epsi10 to set a threshold > 0 when there is no ice) 295 zarea = glob_sum( SUM( a_i + epsi10, dim=3 ) * e1 2t * zconv ) ! in 1.e9 m2295 zarea = glob_sum( SUM( a_i + epsi10, dim=3 ) * e1e2t * zconv ) ! in 1.e9 m2 296 296 zv_sill = zarea * 2.5e-5 297 297 zs_sill = zarea * 25.e-5 -
trunk/NEMOGCM/NEMO/LIM_SRC_3/limctl.F90
r5167 r5836 306 306 WRITE(numout,*) ' - Cell values ' 307 307 WRITE(numout,*) ' ~~~~~~~~~~~ ' 308 WRITE(numout,*) ' cell area : ', e1 2t(ji,jj)308 WRITE(numout,*) ' cell area : ', e1e2t(ji,jj) 309 309 WRITE(numout,*) ' at_i : ', at_i(ji,jj) 310 310 WRITE(numout,*) ' vt_i : ', vt_i(ji,jj) … … 350 350 WRITE(numout,*) ' - Cell values ' 351 351 WRITE(numout,*) ' ~~~~~~~~~~~ ' 352 WRITE(numout,*) ' cell area : ', e1 2t(ji,jj)352 WRITE(numout,*) ' cell area : ', e1e2t(ji,jj) 353 353 WRITE(numout,*) ' at_i : ', at_i(ji,jj) 354 354 WRITE(numout,*) ' vt_i : ', vt_i(ji,jj) -
trunk/NEMOGCM/NEMO/LIM_SRC_3/limdiahsb.F90
r5215 r5836 71 71 72 72 ! 1/area 73 z1_area = 1._wp / MAX( glob_sum( e1 2t(:,:) * tmask(:,:,1) ), epsi06 )74 75 rswitch = MAX( 0._wp , SIGN( 1._wp , glob_sum( e1 2t(:,:) * tmask(:,:,1) ) - epsi06 ) )73 z1_area = 1._wp / MAX( glob_sum( e1e2t(:,:) * tmask(:,:,1) ), epsi06 ) 74 75 rswitch = MAX( 0._wp , SIGN( 1._wp , glob_sum( e1e2t(:,:) * tmask(:,:,1) ) - epsi06 ) ) 76 76 ! ----------------------- ! 77 77 ! 1 - Content variations ! 78 78 ! ----------------------- ! 79 zbg_ivo = glob_sum( vt_i(:,:) * e1 2t(:,:) * tmask(:,:,1) ) ! volume ice80 zbg_svo = glob_sum( vt_s(:,:) * e1 2t(:,:) * tmask(:,:,1) ) ! volume snow81 zbg_are = glob_sum( at_i(:,:) * e1 2t(:,:) * tmask(:,:,1) ) ! area82 zbg_sal = glob_sum( SUM( smv_i(:,:,:), dim=3 ) * e1 2t(:,:) * tmask(:,:,1) ) ! mean salt content83 zbg_tem = glob_sum( ( tm_i(:,:) - rt0 ) * vt_i(:,:) * e1 2t(:,:) * tmask(:,:,1) ) ! mean temp content84 85 !zbg_ihc = glob_sum( et_i(:,:) * e1 2t(:,:) * tmask(:,:,1) ) / MAX( zbg_ivo,epsi06 ) ! ice heat content86 !zbg_shc = glob_sum( et_s(:,:) * e1 2t(:,:) * tmask(:,:,1) ) / MAX( zbg_svo,epsi06 ) ! snow heat content79 zbg_ivo = glob_sum( vt_i(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! volume ice 80 zbg_svo = glob_sum( vt_s(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! volume snow 81 zbg_are = glob_sum( at_i(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! area 82 zbg_sal = glob_sum( SUM( smv_i(:,:,:), dim=3 ) * e1e2t(:,:) * tmask(:,:,1) ) ! mean salt content 83 zbg_tem = glob_sum( ( tm_i(:,:) - rt0 ) * vt_i(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! mean temp content 84 85 !zbg_ihc = glob_sum( et_i(:,:) * e1e2t(:,:) * tmask(:,:,1) ) / MAX( zbg_ivo,epsi06 ) ! ice heat content 86 !zbg_shc = glob_sum( et_s(:,:) * e1e2t(:,:) * tmask(:,:,1) ) / MAX( zbg_svo,epsi06 ) ! snow heat content 87 87 88 88 ! Volume 89 89 ztmp = rswitch * z1_area * r1_rau0 * rday 90 zbg_vfx = ztmp * glob_sum( emp(:,:) * e1 2t(:,:) * tmask(:,:,1) )91 zbg_vfx_bog = ztmp * glob_sum( wfx_bog(:,:) * e1 2t(:,:) * tmask(:,:,1) )92 zbg_vfx_opw = ztmp * glob_sum( wfx_opw(:,:) * e1 2t(:,:) * tmask(:,:,1) )93 zbg_vfx_sni = ztmp * glob_sum( wfx_sni(:,:) * e1 2t(:,:) * tmask(:,:,1) )94 zbg_vfx_dyn = ztmp * glob_sum( wfx_dyn(:,:) * e1 2t(:,:) * tmask(:,:,1) )95 zbg_vfx_bom = ztmp * glob_sum( wfx_bom(:,:) * e1 2t(:,:) * tmask(:,:,1) )96 zbg_vfx_sum = ztmp * glob_sum( wfx_sum(:,:) * e1 2t(:,:) * tmask(:,:,1) )97 zbg_vfx_res = ztmp * glob_sum( wfx_res(:,:) * e1 2t(:,:) * tmask(:,:,1) )98 zbg_vfx_spr = ztmp * glob_sum( wfx_spr(:,:) * e1 2t(:,:) * tmask(:,:,1) )99 zbg_vfx_snw = ztmp * glob_sum( wfx_snw(:,:) * e1 2t(:,:) * tmask(:,:,1) )100 zbg_vfx_sub = ztmp * glob_sum( wfx_sub(:,:) * e1 2t(:,:) * tmask(:,:,1) )90 zbg_vfx = ztmp * glob_sum( emp(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 91 zbg_vfx_bog = ztmp * glob_sum( wfx_bog(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 92 zbg_vfx_opw = ztmp * glob_sum( wfx_opw(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 93 zbg_vfx_sni = ztmp * glob_sum( wfx_sni(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 94 zbg_vfx_dyn = ztmp * glob_sum( wfx_dyn(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 95 zbg_vfx_bom = ztmp * glob_sum( wfx_bom(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 96 zbg_vfx_sum = ztmp * glob_sum( wfx_sum(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 97 zbg_vfx_res = ztmp * glob_sum( wfx_res(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 98 zbg_vfx_spr = ztmp * glob_sum( wfx_spr(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 99 zbg_vfx_snw = ztmp * glob_sum( wfx_snw(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 100 zbg_vfx_sub = ztmp * glob_sum( wfx_sub(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 101 101 102 102 ! Salt 103 zbg_sfx = ztmp * glob_sum( sfx(:,:) * e1 2t(:,:) * tmask(:,:,1) )104 zbg_sfx_bri = ztmp * glob_sum( sfx_bri(:,:) * e1 2t(:,:) * tmask(:,:,1) )105 zbg_sfx_res = ztmp * glob_sum( sfx_res(:,:) * e1 2t(:,:) * tmask(:,:,1) )106 zbg_sfx_dyn = ztmp * glob_sum( sfx_dyn(:,:) * e1 2t(:,:) * tmask(:,:,1) )107 108 zbg_sfx_bog = ztmp * glob_sum( sfx_bog(:,:) * e1 2t(:,:) * tmask(:,:,1) )109 zbg_sfx_opw = ztmp * glob_sum( sfx_opw(:,:) * e1 2t(:,:) * tmask(:,:,1) )110 zbg_sfx_sni = ztmp * glob_sum( sfx_sni(:,:) * e1 2t(:,:) * tmask(:,:,1) )111 zbg_sfx_bom = ztmp * glob_sum( sfx_bom(:,:) * e1 2t(:,:) * tmask(:,:,1) )112 zbg_sfx_sum = ztmp * glob_sum( sfx_sum(:,:) * e1 2t(:,:) * tmask(:,:,1) )103 zbg_sfx = ztmp * glob_sum( sfx(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 104 zbg_sfx_bri = ztmp * glob_sum( sfx_bri(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 105 zbg_sfx_res = ztmp * glob_sum( sfx_res(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 106 zbg_sfx_dyn = ztmp * glob_sum( sfx_dyn(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 107 108 zbg_sfx_bog = ztmp * glob_sum( sfx_bog(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 109 zbg_sfx_opw = ztmp * glob_sum( sfx_opw(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 110 zbg_sfx_sni = ztmp * glob_sum( sfx_sni(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 111 zbg_sfx_bom = ztmp * glob_sum( sfx_bom(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 112 zbg_sfx_sum = ztmp * glob_sum( sfx_sum(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 113 113 114 114 ! Heat budget 115 zbg_ihc = glob_sum( et_i(:,:) * e1 2t(:,:) * 1.e-20 )! ice heat content [1.e20 J]116 zbg_shc = glob_sum( et_s(:,:) * e1 2t(:,:) * 1.e-20 )! snow heat content [1.e20 J]117 zbg_hfx_dhc = glob_sum( diag_heat(:,:) * e1 2t(:,:) * tmask(:,:,1) ) ! [in W]118 zbg_hfx_spr = glob_sum( hfx_spr(:,:) * e1 2t(:,:) * tmask(:,:,1) ) ! [in W]119 120 zbg_hfx_thd = glob_sum( hfx_thd(:,:) * e1 2t(:,:) * tmask(:,:,1) ) ! [in W]121 zbg_hfx_dyn = glob_sum( hfx_dyn(:,:) * e1 2t(:,:) * tmask(:,:,1) ) ! [in W]122 zbg_hfx_res = glob_sum( hfx_res(:,:) * e1 2t(:,:) * tmask(:,:,1) ) ! [in W]123 zbg_hfx_sub = glob_sum( hfx_sub(:,:) * e1 2t(:,:) * tmask(:,:,1) ) ! [in W]124 zbg_hfx_snw = glob_sum( hfx_snw(:,:) * e1 2t(:,:) * tmask(:,:,1) ) ! [in W]125 zbg_hfx_sum = glob_sum( hfx_sum(:,:) * e1 2t(:,:) * tmask(:,:,1) ) ! [in W]126 zbg_hfx_bom = glob_sum( hfx_bom(:,:) * e1 2t(:,:) * tmask(:,:,1) ) ! [in W]127 zbg_hfx_bog = glob_sum( hfx_bog(:,:) * e1 2t(:,:) * tmask(:,:,1) ) ! [in W]128 zbg_hfx_dif = glob_sum( hfx_dif(:,:) * e1 2t(:,:) * tmask(:,:,1) ) ! [in W]129 zbg_hfx_opw = glob_sum( hfx_opw(:,:) * e1 2t(:,:) * tmask(:,:,1) ) ! [in W]130 zbg_hfx_out = glob_sum( hfx_out(:,:) * e1 2t(:,:) * tmask(:,:,1) ) ! [in W]131 zbg_hfx_in = glob_sum( hfx_in(:,:) * e1 2t(:,:) * tmask(:,:,1) ) ! [in W]115 zbg_ihc = glob_sum( et_i(:,:) * e1e2t(:,:) ) * 1.e-20 ! ice heat content [1.e20 J] 116 zbg_shc = glob_sum( et_s(:,:) * e1e2t(:,:) ) * 1.e-20 ! snow heat content [1.e20 J] 117 zbg_hfx_dhc = glob_sum( diag_heat(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 118 zbg_hfx_spr = glob_sum( hfx_spr(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 119 120 zbg_hfx_thd = glob_sum( hfx_thd(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 121 zbg_hfx_dyn = glob_sum( hfx_dyn(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 122 zbg_hfx_res = glob_sum( hfx_res(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 123 zbg_hfx_sub = glob_sum( hfx_sub(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 124 zbg_hfx_snw = glob_sum( hfx_snw(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 125 zbg_hfx_sum = glob_sum( hfx_sum(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 126 zbg_hfx_bom = glob_sum( hfx_bom(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 127 zbg_hfx_bog = glob_sum( hfx_bog(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 128 zbg_hfx_dif = glob_sum( hfx_dif(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 129 zbg_hfx_opw = glob_sum( hfx_opw(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 130 zbg_hfx_out = glob_sum( hfx_out(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 131 zbg_hfx_in = glob_sum( hfx_in(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 132 132 133 133 ! --------------------------------------------- ! 134 134 ! 2 - Trends due to forcing and ice growth/melt ! 135 135 ! --------------------------------------------- ! 136 z_frc_vol = r1_rau0 * glob_sum( - emp(:,:) * e1 2t(:,:) * tmask(:,:,1) ) ! volume fluxes137 z_frc_sal = r1_rau0 * glob_sum( sfx(:,:) * e1 2t(:,:) * tmask(:,:,1) ) ! salt fluxes136 z_frc_vol = r1_rau0 * glob_sum( - emp(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! volume fluxes 137 z_frc_sal = r1_rau0 * glob_sum( sfx(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! salt fluxes 138 138 z_bg_grme = glob_sum( - ( wfx_bog(:,:) + wfx_opw(:,:) + wfx_sni(:,:) + wfx_dyn(:,:) + & 139 139 & wfx_bom(:,:) + wfx_sum(:,:) + wfx_res(:,:) + wfx_snw(:,:) + & 140 & wfx_sub(:,:) ) * e1 2t(:,:) * tmask(:,:,1) ) ! volume fluxes140 & wfx_sub(:,:) ) * e1e2t(:,:) * tmask(:,:,1) ) ! volume fluxes 141 141 ! 142 142 frc_vol = frc_vol + z_frc_vol * rdt_ice -
trunk/NEMOGCM/NEMO/LIM_SRC_3/limdyn.F90
r5123 r5836 191 191 CALL prt_ctl(tab2d_1=delta_i , clinfo1=' lim_dyn : delta_i :') 192 192 CALL prt_ctl(tab2d_1=strength , clinfo1=' lim_dyn : strength :') 193 CALL prt_ctl(tab2d_1=e1 2t, clinfo1=' lim_dyn : cell area :')193 CALL prt_ctl(tab2d_1=e1e2t , clinfo1=' lim_dyn : cell area :') 194 194 CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_dyn : at_i :') 195 195 CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_dyn : vt_i :') -
trunk/NEMOGCM/NEMO/LIM_SRC_3/limhdf.F90
r5429 r5836 76 76 DO jj = 2, jpjm1 77 77 DO ji = fs_2 , fs_jpim1 ! vector opt. 78 efact(ji,jj) = ( e2u(ji,jj) + e2u(ji-1,jj) + e1v(ji,jj) + e1v(ji,jj-1) ) * r1_e1 2t(ji,jj)78 efact(ji,jj) = ( e2u(ji,jj) + e2u(ji-1,jj) + e1v(ji,jj) + e1v(ji,jj-1) ) * r1_e1e2t(ji,jj) 79 79 END DO 80 80 END DO … … 107 107 DO jj= 2, jpjm1 ! diffusive trend : divergence of the fluxes 108 108 DO ji = fs_2 , fs_jpim1 ! vector opt. 109 zdiv(ji,jj) = ( zflu(ji,jj) - zflu(ji-1,jj) + zflv(ji,jj) - zflv(ji,jj-1) ) * r1_e1 2t(ji,jj)109 zdiv(ji,jj) = ( zflu(ji,jj) - zflu(ji-1,jj) + zflv(ji,jj) - zflv(ji,jj-1) ) * r1_e1e2t(ji,jj) 110 110 END DO 111 111 END DO … … 149 149 DO jj= 2, jpjm1 ! diffusive trend : divergence of the fluxes 150 150 DO ji = fs_2 , fs_jpim1 ! vector opt. 151 zdiv(ji,jj) = ( zflu(ji,jj) - zflu(ji-1,jj) + zflv(ji,jj) - zflv(ji,jj-1) ) * r1_e1 2t(ji,jj)151 zdiv(ji,jj) = ( zflu(ji,jj) - zflu(ji-1,jj) + zflv(ji,jj) - zflv(ji,jj-1) ) * r1_e1e2t(ji,jj) 152 152 ptab(ji,jj) = ztab0(ji,jj) + 0.5 * ( zdiv(ji,jj) + zdiv0(ji,jj) ) 153 153 END DO -
trunk/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90
r5202 r5836 377 377 CALL prt_ctl_info(' - Cell values : ') 378 378 CALL prt_ctl_info(' ~~~~~~~~~~~~~ ') 379 CALL prt_ctl(tab2d_1=e1 2t, clinfo1=' lim_itd_me : cell area :')379 CALL prt_ctl(tab2d_1=e1e2t, clinfo1=' lim_itd_me : cell area :') 380 380 CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_itd_me : at_i :') 381 381 CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_itd_me : vt_i :') -
trunk/NEMOGCM/NEMO/LIM_SRC_3/limrhg.F90
r5429 r5836 355 355 divu_i(ji,jj) = ( e2u(ji,jj) * u_ice(ji,jj) - e2u(ji-1,jj) * u_ice(ji-1,jj) & 356 356 & + e1v(ji,jj) * v_ice(ji,jj) - e1v(ji,jj-1) * v_ice(ji,jj-1) & 357 & ) * r1_e1 2t(ji,jj)357 & ) * r1_e1e2t(ji,jj) 358 358 359 359 zdt(ji,jj) = ( ( u_ice(ji,jj) * r1_e2u(ji,jj) - u_ice(ji-1,jj) * r1_e2u(ji-1,jj) ) * e2t(ji,jj) * e2t(ji,jj) & 360 360 & - ( v_ice(ji,jj) * r1_e1v(ji,jj) - v_ice(ji,jj-1) * r1_e1v(ji,jj-1) ) * e1t(ji,jj) * e1t(ji,jj) & 361 & ) * r1_e1 2t(ji,jj)361 & ) * r1_e1e2t(ji,jj) 362 362 363 363 ! 364 364 zds(ji,jj) = ( ( u_ice(ji,jj+1) * r1_e1u(ji,jj+1) - u_ice(ji,jj) * r1_e1u(ji,jj) ) * e1f(ji,jj) * e1f(ji,jj) & 365 365 & + ( v_ice(ji+1,jj) * r1_e2v(ji+1,jj) - v_ice(ji,jj) * r1_e2v(ji,jj) ) * e2f(ji,jj) * e2f(ji,jj) & 366 & ) * r1_e1 2f(ji,jj) * ( 2._wp - fmask(ji,jj,1) ) &366 & ) * r1_e1e2f(ji,jj) * ( 2._wp - fmask(ji,jj,1) ) & 367 367 & * zmask(ji,jj) * zmask(ji,jj+1) * zmask(ji+1,jj) * zmask(ji+1,jj+1) 368 368 … … 386 386 zdst = ( e2u(ji,jj) * v_ice1(ji,jj) - e2u(ji-1,jj ) * v_ice1(ji-1,jj ) & 387 387 & + e1v(ji,jj) * u_ice2(ji,jj) - e1v(ji ,jj-1) * u_ice2(ji ,jj-1) & 388 & ) * r1_e1 2t(ji,jj)388 & ) * r1_e1e2t(ji,jj) 389 389 390 390 delta = SQRT( divu_i(ji,jj)**2 + ( zdt(ji,jj)**2 + zdst**2 ) * usecc2 ) … … 394 394 zddc = ( ( v_ice1(ji,jj+1) * r1_e1u(ji,jj+1) - v_ice1(ji,jj) * r1_e1u(ji,jj) ) * e1f(ji,jj) * e1f(ji,jj) & 395 395 & + ( u_ice2(ji+1,jj) * r1_e2v(ji+1,jj) - u_ice2(ji,jj) * r1_e2v(ji,jj) ) * e2f(ji,jj) * e2f(ji,jj) & 396 & ) * r1_e1 2f(ji,jj)396 & ) * r1_e1e2f(ji,jj) 397 397 398 398 zdtc = (- ( v_ice1(ji,jj+1) * r1_e1u(ji,jj+1) - v_ice1(ji,jj) * r1_e1u(ji,jj) ) * e1f(ji,jj) * e1f(ji,jj) & 399 399 & + ( u_ice2(ji+1,jj) * r1_e2v(ji+1,jj) - u_ice2(ji,jj) * r1_e2v(ji,jj) ) * e2f(ji,jj) * e2f(ji,jj) & 400 & ) * r1_e1 2f(ji,jj)400 & ) * r1_e1e2f(ji,jj) 401 401 402 402 zddc = SQRT( zddc**2 + ( zdtc**2 + zds(ji,jj)**2 ) * usecc2 ) + rn_creepl … … 423 423 & + ( zs2(ji+1,jj) * e2t(ji+1,jj)**2 - zs2(ji,jj) * e2t(ji,jj)**2 ) * r1_e2u(ji,jj) & 424 424 & + 2.0 * ( zs12(ji,jj) * e1f(ji,jj)**2 - zs12(ji,jj-1) * e1f(ji,jj-1)**2 ) * r1_e1u(ji,jj) & 425 & ) * r1_e1 2u(ji,jj)425 & ) * r1_e1e2u(ji,jj) 426 426 ! contribution of zs1, zs2 and zs12 to zf2 427 427 zf2(ji,jj) = 0.5 * ( ( zs1(ji,jj+1) - zs1(ji,jj) ) * e1v(ji,jj) & 428 428 & - ( zs2(ji,jj+1) * e1t(ji,jj+1)**2 - zs2(ji,jj) * e1t(ji,jj)**2 ) * r1_e1v(ji,jj) & 429 429 & + 2.0 * ( zs12(ji,jj) * e2f(ji,jj)**2 - zs12(ji-1,jj) * e2f(ji-1,jj)**2 ) * r1_e2v(ji,jj) & 430 & ) * r1_e1 2v(ji,jj)430 & ) * r1_e1e2v(ji,jj) 431 431 END DO 432 432 END DO … … 607 607 divu_i(ji,jj) = ( e2u(ji,jj) * u_ice(ji,jj) - e2u(ji-1,jj ) * u_ice(ji-1,jj ) & 608 608 & + e1v(ji,jj) * v_ice(ji,jj) - e1v(ji ,jj-1) * v_ice(ji ,jj-1) & 609 & ) * r1_e1 2t(ji,jj)609 & ) * r1_e1e2t(ji,jj) 610 610 611 611 zdt(ji,jj) = ( ( u_ice(ji,jj) * r1_e2u(ji,jj) - u_ice(ji-1,jj) * r1_e2u(ji-1,jj) ) * e2t(ji,jj) * e2t(ji,jj) & 612 612 & -( v_ice(ji,jj) * r1_e1v(ji,jj) - v_ice(ji,jj-1) * r1_e1v(ji,jj-1) ) * e1t(ji,jj) * e1t(ji,jj) & 613 & ) * r1_e1 2t(ji,jj)613 & ) * r1_e1e2t(ji,jj) 614 614 ! 615 615 ! SB modif because ocean has no slip boundary condition 616 616 zds(ji,jj) = ( ( u_ice(ji,jj+1) * r1_e1u(ji,jj+1) - u_ice(ji,jj) * r1_e1u(ji,jj) ) * e1f(ji,jj) * e1f(ji,jj) & 617 617 & +( v_ice(ji+1,jj) * r1_e2v(ji+1,jj) - v_ice(ji,jj) * r1_e2v(ji,jj) ) * e2f(ji,jj) * e2f(ji,jj) & 618 & ) * r1_e1 2f(ji,jj) * ( 2.0 - fmask(ji,jj,1) ) &618 & ) * r1_e1e2f(ji,jj) * ( 2.0 - fmask(ji,jj,1) ) & 619 619 & * zmask(ji,jj) * zmask(ji,jj+1) * zmask(ji+1,jj) * zmask(ji+1,jj+1) 620 620 621 621 zdst = ( e2u(ji,jj) * v_ice1(ji,jj) - e2u(ji-1,jj ) * v_ice1(ji-1,jj ) & 622 & + e1v(ji,jj) * u_ice2(ji,jj) - e1v(ji ,jj-1) * u_ice2(ji ,jj-1) ) * r1_e1 2t(ji,jj)622 & + e1v(ji,jj) * u_ice2(ji,jj) - e1v(ji ,jj-1) * u_ice2(ji ,jj-1) ) * r1_e1e2t(ji,jj) 623 623 624 624 delta = SQRT( divu_i(ji,jj)**2 + ( zdt(ji,jj)**2 + zdst**2 ) * usecc2 ) … … 637 637 DO ji = fs_2, fs_jpim1 638 638 zdst = ( e2u(ji,jj) * v_ice1(ji,jj) - e2u( ji-1, jj ) * v_ice1(ji-1,jj) & 639 & + e1v(ji,jj) * u_ice2(ji,jj) - e1v( ji , jj-1 ) * u_ice2(ji,jj-1) ) * r1_e1 2t(ji,jj)639 & + e1v(ji,jj) * u_ice2(ji,jj) - e1v( ji , jj-1 ) * u_ice2(ji,jj-1) ) * r1_e1e2t(ji,jj) 640 640 shear_i(ji,jj) = SQRT( zdt(ji,jj) * zdt(ji,jj) + zdst * zdst ) 641 641 END DO -
trunk/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90
r5407 r5836 325 325 CALL prt_ctl_info(' - Cell values : ') 326 326 CALL prt_ctl_info(' ~~~~~~~~~~~~~ ') 327 CALL prt_ctl(tab2d_1=e1 2t, clinfo1=' lim_thd : cell area :')327 CALL prt_ctl(tab2d_1=e1e2t, clinfo1=' lim_thd : cell area :') 328 328 CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_thd : at_i :') 329 329 CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_thd : vt_i :') … … 382 382 CALL prt_ctl_info(' - Cell values : ') 383 383 CALL prt_ctl_info(' ~~~~~~~~~~~~~ ') 384 CALL prt_ctl(tab2d_1=e1 2t, clinfo1=' lim_itd_th : cell area :')384 CALL prt_ctl(tab2d_1=e1e2t, clinfo1=' lim_itd_th : cell area :') 385 385 CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_itd_th : at_i :') 386 386 CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_itd_th : vt_i :') -
trunk/NEMOGCM/NEMO/LIM_SRC_3/limtrp.F90
r5202 r5836 95 95 ENDIF 96 96 97 zsm(:,:) = e1 2t(:,:)97 zsm(:,:) = e1e2t(:,:) 98 98 99 99 ! !-------------------------------------! … … 162 162 ! transported fields 163 163 !------------------------- 164 z0opw(:,:,1) = ato_i(:,:) * e1 2t(:,:) ! Open water area165 DO jl = 1, jpl 166 z0snw (:,:,jl) = v_s (:,:, jl) * e12t(:,:)! Snow volume167 z0ice(:,:,jl) = v_i (:,:, jl) * e12t(:,:)! Ice volume168 z0ai (:,:,jl) = a_i (:,:, jl) * e12t(:,:)! Ice area169 z0smi (:,:,jl) = smv_i(:,:, jl) * e12t(:,:)! Salt content170 z0oi (:,:,jl) = oa_i (:,:, jl) * e12t(:,:)! Age content171 z0es (:,:,jl) = e_s (:,:,1,jl) * e1 2t(:,:) ! Snow heat content164 z0opw(:,:,1) = ato_i(:,:) * e1e2t(:,:) ! Open water area 165 DO jl = 1, jpl 166 z0snw (:,:,jl) = v_s (:,:, jl) * e1e2t(:,:) ! Snow volume 167 z0ice(:,:,jl) = v_i (:,:, jl) * e1e2t(:,:) ! Ice volume 168 z0ai (:,:,jl) = a_i (:,:, jl) * e1e2t(:,:) ! Ice area 169 z0smi (:,:,jl) = smv_i(:,:, jl) * e1e2t(:,:) ! Salt content 170 z0oi (:,:,jl) = oa_i (:,:, jl) * e1e2t(:,:) ! Age content 171 z0es (:,:,jl) = e_s (:,:,1,jl) * e1e2t(:,:) ! Snow heat content 172 172 DO jk = 1, nlay_i 173 z0ei (:,:,jk,jl) = e_i (:,:,jk,jl) * e1 2t(:,:) ! Ice heat content173 z0ei (:,:,jk,jl) = e_i (:,:,jk,jl) * e1e2t(:,:) ! Ice heat content 174 174 END DO 175 175 END DO … … 263 263 ! Recover the properties from their contents 264 264 !------------------------------------------- 265 ato_i(:,:) = z0opw(:,:,1) * r1_e1 2t(:,:)266 DO jl = 1, jpl 267 v_i (:,:, jl) = z0ice(:,:,jl) * r1_e12t(:,:)268 v_s (:,:, jl) = z0snw(:,:,jl) * r1_e12t(:,:)269 smv_i(:,:, jl) = z0smi(:,:,jl) * r1_e12t(:,:)270 oa_i (:,:, jl) = z0oi (:,:,jl) * r1_e12t(:,:)271 a_i (:,:, jl) = z0ai (:,:,jl) * r1_e12t(:,:)272 e_s (:,:,1,jl) = z0es (:,:,jl) * r1_e1 2t(:,:)265 ato_i(:,:) = z0opw(:,:,1) * r1_e1e2t(:,:) 266 DO jl = 1, jpl 267 v_i (:,:, jl) = z0ice(:,:,jl) * r1_e1e2t(:,:) 268 v_s (:,:, jl) = z0snw(:,:,jl) * r1_e1e2t(:,:) 269 smv_i(:,:, jl) = z0smi(:,:,jl) * r1_e1e2t(:,:) 270 oa_i (:,:, jl) = z0oi (:,:,jl) * r1_e1e2t(:,:) 271 a_i (:,:, jl) = z0ai (:,:,jl) * r1_e1e2t(:,:) 272 e_s (:,:,1,jl) = z0es (:,:,jl) * r1_e1e2t(:,:) 273 273 DO jk = 1, nlay_i 274 e_i(:,:,jk,jl) = z0ei(:,:,jk,jl) * r1_e1 2t(:,:)274 e_i(:,:,jk,jl) = z0ei(:,:,jk,jl) * r1_e1e2t(:,:) 275 275 END DO 276 276 END DO -
trunk/NEMOGCM/NEMO/LIM_SRC_3/limupdate1.F90
r5215 r5836 146 146 CALL prt_ctl_info(' - Cell values : ') 147 147 CALL prt_ctl_info(' ~~~~~~~~~~~~~ ') 148 CALL prt_ctl(tab2d_1=e1 2t, clinfo1=' lim_update1 : cell area :')148 CALL prt_ctl(tab2d_1=e1e2t , clinfo1=' lim_update1 : cell area :') 149 149 CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_update1 : at_i :') 150 150 CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_update1 : vt_i :') -
trunk/NEMOGCM/NEMO/LIM_SRC_3/limupdate2.F90
r5410 r5836 191 191 CALL prt_ctl_info(' - Cell values : ') 192 192 CALL prt_ctl_info(' ~~~~~~~~~~~~~ ') 193 CALL prt_ctl(tab2d_1=e1 2t, clinfo1=' lim_update2 : cell area :')193 CALL prt_ctl(tab2d_1=e1e2t , clinfo1=' lim_update2 : cell area :') 194 194 CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_update2 : at_i :') 195 195 CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_update2 : vt_i :') -
trunk/NEMOGCM/NEMO/NST_SRC/agrif_opa_sponge.F90
r5656 r5836 210 210 DO jj = j1,j2-1 211 211 DO ji = i1,i2-1 212 zabe1 = fsaht_spu(ji,jj) * umask(ji,jj,jk) * re2u_e1u(ji,jj) * fse3u_n(ji,jj,jk)213 zabe2 = fsaht_spv(ji,jj) * vmask(ji,jj,jk) * re1v_e2v(ji,jj) * fse3v_n(ji,jj,jk)212 zabe1 = fsaht_spu(ji,jj) * umask(ji,jj,jk) * e2_e1u(ji,jj) * fse3u_n(ji,jj,jk) 213 zabe2 = fsaht_spv(ji,jj) * vmask(ji,jj,jk) * e1_e2v(ji,jj) * fse3v_n(ji,jj,jk) 214 214 ztu(ji,jj,jk) = zabe1 * ( tsbdiff(ji+1,jj ,jk,jn) - tsbdiff(ji,jj,jk,jn) ) 215 215 ztv(ji,jj,jk) = zabe2 * ( tsbdiff(ji ,jj+1,jk,jn) - tsbdiff(ji,jj,jk,jn) ) … … 239 239 240 240 IF (.NOT. tabspongedone_tsn(ji,jj)) THEN 241 zbtr = r1_e1 2t(ji,jj) / fse3t_n(ji,jj,jk)241 zbtr = r1_e1e2t(ji,jj) / fse3t_n(ji,jj,jk) 242 242 ! horizontal diffusive trends 243 243 ztsa = zbtr * ( ztu(ji,jj,jk) - ztu(ji-1,jj,jk) + ztv(ji,jj,jk) - ztv(ji ,jj-1,jk) ) … … 290 290 DO jj = j1,j2 291 291 DO ji = i1+1,i2 ! vector opt. 292 zbtr = r1_e1 2t(ji,jj) / fse3t_n(ji,jj,jk) * fsahm_spt(ji,jj)292 zbtr = r1_e1e2t(ji,jj) / fse3t_n(ji,jj,jk) * fsahm_spt(ji,jj) 293 293 hdivdiff(ji,jj,jk) = ( e2u(ji ,jj)*fse3u_n(ji ,jj,jk) * ubdiff(ji ,jj,jk) & 294 294 & -e2u(ji-1,jj)*fse3u_n(ji-1,jj,jk) * ubdiff(ji-1,jj,jk) ) * zbtr … … 298 298 DO jj = j1,j2-1 299 299 DO ji = i1,i2 ! vector opt. 300 zbtr = r1_e1 2f(ji,jj) * fse3f_n(ji,jj,jk) * fsahm_spf(ji,jj)300 zbtr = r1_e1e2f(ji,jj) * fse3f_n(ji,jj,jk) * fsahm_spf(ji,jj) 301 301 rotdiff(ji,jj,jk) = (-e1u(ji,jj+1) * ubdiff(ji,jj+1,jk) & 302 302 +e1u(ji,jj ) * ubdiff(ji,jj ,jk) & … … 396 396 DO jj = j1+1,j2 397 397 DO ji = i1,i2 ! vector opt. 398 zbtr = r1_e1 2t(ji,jj) / fse3t_n(ji,jj,jk) * fsahm_spt(ji,jj)398 zbtr = r1_e1e2t(ji,jj) / fse3t_n(ji,jj,jk) * fsahm_spt(ji,jj) 399 399 hdivdiff(ji,jj,jk) = ( e1v(ji,jj ) * fse3v(ji,jj ,jk) * vbdiff(ji,jj ,jk) & 400 400 & -e1v(ji,jj-1) * fse3v(ji,jj-1,jk) * vbdiff(ji,jj-1,jk) ) * zbtr … … 403 403 DO jj = j1,j2 404 404 DO ji = i1,i2-1 ! vector opt. 405 zbtr = r1_e1 2f(ji,jj) * fse3f_n(ji,jj,jk) * fsahm_spf(ji,jj)405 zbtr = r1_e1e2f(ji,jj) * fse3f_n(ji,jj,jk) * fsahm_spf(ji,jj) 406 406 rotdiff(ji,jj,jk) = ( e2v(ji+1,jj) * vbdiff(ji+1,jj,jk) & 407 407 & -e2v(ji ,jj) * vbdiff(ji ,jj,jk) & -
trunk/NEMOGCM/NEMO/NST_SRC/agrif_top_sponge.F90
r5656 r5836 74 74 DO jj = j1,j2-1 75 75 DO ji = i1,i2-1 76 zabe1 = fsaht_spu(ji,jj) * umask(ji,jj,jk) * re2u_e1u(ji,jj) * fse3u_n(ji,jj,jk)77 zabe2 = fsaht_spv(ji,jj) * vmask(ji,jj,jk) * re1v_e2v(ji,jj) * fse3v_n(ji,jj,jk)76 zabe1 = fsaht_spu(ji,jj) * umask(ji,jj,jk) * e2_e1u(ji,jj) * fse3u_n(ji,jj,jk) 77 zabe2 = fsaht_spv(ji,jj) * vmask(ji,jj,jk) * e1_e2v(ji,jj) * fse3v_n(ji,jj,jk) 78 78 ztu(ji,jj) = zabe1 * ( trbdiff(ji+1,jj ,jk,jn) - trbdiff(ji,jj,jk,jn) ) 79 79 ztv(ji,jj) = zabe2 * ( trbdiff(ji ,jj+1,jk,jn) - trbdiff(ji,jj,jk,jn) ) … … 85 85 86 86 IF (.NOT. tabspongedone_trn(ji,jj)) THEN 87 zbtr = r1_e1 2t(ji,jj) / fse3t(ji,jj,jk)87 zbtr = r1_e1e2t(ji,jj) / fse3t(ji,jj,jk) 88 88 ! horizontal diffusive trends 89 89 ztra = zbtr * ( ztu(ji,jj) - ztu(ji-1,jj ) + ztv(ji,jj) - ztv(ji ,jj-1) ) -
trunk/NEMOGCM/NEMO/OFF_SRC/domrea.F90
r5504 r5836 4 4 !! Ocean initialization : domain initialization 5 5 !!============================================================================== 6 !! History : OPA ! 1990-10 (C. Levy - G. Madec) Original code 7 !! ! 1992-01 (M. Imbard) insert time step initialization 8 !! ! 1996-06 (G. Madec) generalized vertical coordinate 9 !! ! 1997-02 (G. Madec) creation of domwri.F 10 !! ! 2001-05 (E.Durand - G. Madec) insert closed sea 11 !! NEMO 1.0 ! 2002-08 (G. Madec) F90: Free form and module 12 !!---------------------------------------------------------------------- 6 13 7 14 !!---------------------------------------------------------------------- … … 10 17 !! dom_ctl : control print for the ocean domain 11 18 !!---------------------------------------------------------------------- 12 !! * Modules used13 19 USE oce ! 20 USE trc_oce ! shared ocean/biogeochemical variables 14 21 USE dom_oce ! ocean space and time domain 15 22 USE phycst ! physical constants 23 USE domstp ! domain: set the time-step 24 ! 16 25 USE in_out_manager ! I/O manager 17 26 USE lib_mpp ! distributed memory computing library 18 19 USE domstp ! domain: set the time-step20 21 27 USE lbclnk ! lateral boundary condition - MPP exchanges 22 USE trc_oce ! shared ocean/biogeochemical variables23 28 USE wrk_nemo 24 29 … … 26 31 PRIVATE 27 32 28 !! * Routine accessibility 29 PUBLIC dom_rea ! called by opa.F90 33 PUBLIC dom_rea ! called by nemogcm.F90 30 34 31 35 !! * Substitutions … … 33 37 # include "vectopt_loop_substitute.h90" 34 38 !!---------------------------------------------------------------------- 35 !! NEMO/OFF 3. 3 , NEMO Consortium (2010)39 !! NEMO/OFF 3.7 , NEMO Consortium (2015) 36 40 !! $Id$ 37 41 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 38 42 !!---------------------------------------------------------------------- 39 40 43 CONTAINS 41 44 … … 51 54 !! - dom_stp: defined the model time step 52 55 !! - dom_rea: read the meshmask file if nmsh=1 53 !! 54 !! History : 55 !! ! 90-10 (C. Levy - G. Madec) Original code 56 !! ! 91-11 (G. Madec) 57 !! ! 92-01 (M. Imbard) insert time step initialization 58 !! ! 96-06 (G. Madec) generalized vertical coordinate 59 !! ! 97-02 (G. Madec) creation of domwri.F 60 !! ! 01-05 (E.Durand - G. Madec) insert closed sea 61 !! 8.5 ! 02-08 (G. Madec) F90: Free form and module 62 !!---------------------------------------------------------------------- 63 !! * Local declarations 64 INTEGER :: jk ! dummy loop argument 65 INTEGER :: iconf = 0 ! temporary integers 66 !!---------------------------------------------------------------------- 67 56 !!---------------------------------------------------------------------- 57 INTEGER :: jk ! dummy loop index 58 INTEGER :: iconf = 0 ! local integers 59 !!---------------------------------------------------------------------- 60 ! 68 61 IF(lwp) THEN 69 62 WRITE(numout,*) … … 71 64 WRITE(numout,*) '~~~~~~~~' 72 65 ENDIF 73 74 CALL dom_nam ! read namelist ( namrun, namdom , namcla)66 ! 67 CALL dom_nam ! read namelist ( namrun, namdom ) 75 68 CALL dom_zgr ! Vertical mesh and bathymetry option 76 69 CALL dom_grd ! Create a domain file 77 78 !79 ! - ML - Used in dom_vvl_sf_nxt and lateral diffusion routines80 ! but could be usefull in many other routines81 e12t (:,:) = e1t(:,:) * e2t(:,:)82 e1e2t (:,:) = e1t(:,:) * e2t(:,:)83 e12u (:,:) = e1u(:,:) * e2u(:,:)84 e12v (:,:) = e1v(:,:) * e2v(:,:)85 e1 2f (:,:) = e1f(:,:) * e2f(:,:)86 r1_e12t (:,:) = 1._wp / e12t(:,:)87 r1_e12u (:,:) = 1._wp / e12u(:,:)88 r1_e12v (:,:) = 1._wp / e12v(:,:)89 r1_e12f (:,:) = 1._wp / e12f(:,:)90 re2u_e1u(:,:) = e2u(:,:) / e1u(:,:)91 re1v_e2v(:,:) = e1v(:,:) / e2v(:,:)92 ! 93 hu(:,:) = 0._wp 70 ! 71 ! ! associated horizontal metrics 72 ! 73 r1_e1t(:,:) = 1._wp / e1t(:,:) ; r1_e2t (:,:) = 1._wp / e2t(:,:) 74 r1_e1u(:,:) = 1._wp / e1u(:,:) ; r1_e2u (:,:) = 1._wp / e2u(:,:) 75 r1_e1v(:,:) = 1._wp / e1v(:,:) ; r1_e2v (:,:) = 1._wp / e2v(:,:) 76 r1_e1f(:,:) = 1._wp / e1f(:,:) ; r1_e2f (:,:) = 1._wp / e2f(:,:) 77 ! 78 e1e2t (:,:) = e1t(:,:) * e2t(:,:) ; r1_e1e2t(:,:) = 1._wp / e1e2t(:,:) 79 e1e2u (:,:) = e1u(:,:) * e2u(:,:) ; r1_e1e2u(:,:) = 1._wp / e1e2u(:,:) 80 e1e2v (:,:) = e1v(:,:) * e2v(:,:) ; r1_e1e2v(:,:) = 1._wp / e1e2v(:,:) 81 e1e2f (:,:) = e1f(:,:) * e2f(:,:) ; r1_e1e2f(:,:) = 1._wp / e1e2f(:,:) 82 ! 83 e2_e1u(:,:) = e2u(:,:) / e1u(:,:) 84 e1_e2v(:,:) = e1v(:,:) / e2v(:,:) 85 ! 86 hu(:,:) = 0._wp ! Ocean depth at U- and V-points 94 87 hv(:,:) = 0._wp 95 88 DO jk = 1, jpk … … 100 93 hur(:,:) = 1._wp / ( hu(:,:) + 1._wp - umask(:,:,1) ) * umask(:,:,1) 101 94 hvr(:,:) = 1._wp / ( hv(:,:) + 1._wp - vmask(:,:,1) ) * vmask(:,:,1) 102 95 ! 103 96 CALL dom_stp ! Time step 104 97 CALL dom_msk ! Masks 105 98 CALL dom_ctl ! Domain control 106 99 ! 107 100 END SUBROUTINE dom_rea 101 108 102 109 103 SUBROUTINE dom_nam … … 115 109 !! ** input : - namrun namelist 116 110 !! - namdom namelist 117 !! - namcla namelist118 111 !!---------------------------------------------------------------------- 119 112 USE ioipsl 120 INTEGER :: ios ! Local integer output status for namelist read 113 INTEGER :: ios ! Local integer output status for namelist read 114 ! 121 115 NAMELIST/namrun/ cn_ocerst_indir, cn_ocerst_outdir, nn_stocklist, ln_rst_list, & 122 116 & nn_no , cn_exp , cn_ocerst_in, cn_ocerst_out, ln_rstart , nn_rstctl, & … … 130 124 & ppsur, ppa0, ppa1, ppkth, ppacr, ppdzmin, pphmax, ldbletanh, & 131 125 & ppa2, ppkth2, ppacr2 132 NAMELIST/namcla/ nn_cla133 126 #if defined key_netcdf4 134 127 NAMELIST/namnc4/ nn_nchunks_i, nn_nchunks_j, nn_nchunks_k, ln_nc4zip … … 178 171 nstocklist = nn_stocklist 179 172 nwrite = nn_write 180 181 173 ! 182 174 ! ! control of output frequency 183 175 IF ( nstock == 0 .OR. nstock > nitend ) THEN … … 275 267 rdth = rn_rdth 276 268 277 REWIND( numnam_ref ) ! Namelist namcla in reference namelist : Cross land advection278 READ ( numnam_ref, namcla, IOSTAT = ios, ERR = 905)279 905 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcla in reference namelist', lwp )280 281 REWIND( numnam_cfg ) ! Namelist namcla in configuration namelist : Cross land advection282 READ ( numnam_cfg, namcla, IOSTAT = ios, ERR = 906 )283 906 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcla in configuration namelist', lwp )284 IF(lwm) WRITE( numond, namcla )285 286 IF(lwp) THEN287 WRITE(numout,*)288 WRITE(numout,*) ' Namelist namcla'289 WRITE(numout,*) ' cross land advection nn_cla = ', nn_cla290 ENDIF291 292 269 #if defined key_netcdf4 293 270 ! ! NetCDF 4 case ("key_netcdf4" defined) … … 321 298 END SUBROUTINE dom_nam 322 299 300 323 301 SUBROUTINE dom_zgr 324 302 !!---------------------------------------------------------------------- … … 374 352 END SUBROUTINE dom_zgr 375 353 354 376 355 SUBROUTINE dom_ctl 377 356 !!---------------------------------------------------------------------- … … 382 361 !! ** Method : compute and print extrema of masked scale factors 383 362 !! 384 !! History : 385 !! 8.5 ! 02-08 (G. Madec) Original code 386 !!---------------------------------------------------------------------- 387 !! * Local declarations 363 !!---------------------------------------------------------------------- 388 364 INTEGER :: iimi1, ijmi1, iimi2, ijmi2, iima1, ijma1, iima2, ijma2 389 365 INTEGER, DIMENSION(2) :: iloc ! … … 421 397 ijma2 = iloc(2) + njmpp - 1 422 398 ENDIF 423 399 ! 424 400 IF(lwp) THEN 425 401 WRITE(numout,"(14x,'e1t maxi: ',1f10.2,' at i = ',i5,' j= ',i5)") ze1max, iima1, ijma1 … … 428 404 WRITE(numout,"(14x,'e2t mini: ',1f10.2,' at i = ',i5,' j= ',i5)") ze2min, iimi2, ijmi2 429 405 ENDIF 430 406 ! 431 407 END SUBROUTINE dom_ctl 408 432 409 433 410 SUBROUTINE dom_grd … … 538 515 CALL iom_get( inum2, jpdom_data, 'facvolt', facvol ) 539 516 #endif 540 541 517 ! ! horizontal mesh (inum3) 542 518 CALL iom_get( inum3, jpdom_data, 'glamt', glamt ) … … 756 732 !! (min value = 1 over land) 757 733 !!---------------------------------------------------------------------- 758 !759 734 INTEGER :: ji, jj ! dummy loop indices 760 735 REAL(wp), POINTER, DIMENSION(:,:) :: zmbk … … 785 760 END SUBROUTINE zgr_bot_level 786 761 762 787 763 SUBROUTINE dom_msk 788 764 !!--------------------------------------------------------------------- … … 799 775 !! tpol : ??? 800 776 !!---------------------------------------------------------------------- 801 ! 802 INTEGER :: ji, jj, jk ! dummy loop indices 803 INTEGER :: iif, iil, ijf, ijl ! local integers 777 INTEGER :: ji, jj, jk ! dummy loop indices 778 INTEGER :: iif, iil, ijf, ijl ! local integers 804 779 INTEGER, POINTER, DIMENSION(:,:) :: imsk 805 !806 780 !!--------------------------------------------------------------------- 807 781 … … 853 827 ! 3. Ocean/land mask at wu-, wv- and w points 854 828 !---------------------------------------------- 855 wmask (:,:,1) = tmask(:,:,1) ! ????????856 wumask(:,:,1) = umask(:,:,1) ! ????????857 wvmask(:,:,1) = vmask(:,:,1) ! ????????858 DO jk =2,jpk859 wmask (:,:,jk) =tmask(:,:,jk) * tmask(:,:,jk-1)860 wumask(:,:,jk) =umask(:,:,jk) * umask(:,:,jk-1)861 wvmask(:,:,jk) =vmask(:,:,jk) * vmask(:,:,jk-1)829 wmask (:,:,1) = tmask(:,:,1) ! surface value 830 wumask(:,:,1) = umask(:,:,1) 831 wvmask(:,:,1) = vmask(:,:,1) 832 DO jk = 2, jpk ! deeper value 833 wmask (:,:,jk) = tmask(:,:,jk) * tmask(:,:,jk-1) 834 wumask(:,:,jk) = umask(:,:,jk) * umask(:,:,jk-1) 835 wvmask(:,:,jk) = vmask(:,:,jk) * vmask(:,:,jk-1) 862 836 END DO 863 837 ! -
trunk/NEMOGCM/NEMO/OFF_SRC/dtadyn.F90
r5768 r5836 26 26 USE trc_oce ! share ocean/biogeo variables 27 27 USE phycst ! physical constants 28 USE ldftra ! lateral diffusivity coefficients 28 29 USE trabbl ! active tracer: bottom boundary layer 29 30 USE ldfslp ! lateral diffusion: iso-neutral slopes 30 USE ldfeiv ! eddy induced velocity coef.31 USE ldftra_oce ! ocean tracer lateral physics32 31 USE zdfmxl ! vertical physics: mixed layer depth 33 32 USE eosbn2 ! equation of state - Brunt Vaisala frequency … … 40 39 USE fldread ! read input fields 41 40 USE timing ! Timing 41 USE wrk_nemo 42 42 43 43 IMPLICIT NONE … … 50 50 LOGICAL :: ln_dynwzv !: vertical velocity read in a file (T) or computed from u/v (F) 51 51 LOGICAL :: ln_dynbbl !: bbl coef read in a file (T) or computed (F) 52 LOGICAL :: ln_degrad !: degradation option enabled or not53 52 LOGICAL :: ln_dynrnf !: read runoff data in file (T) or set to zero (F) 54 53 55 INTEGER , PARAMETER :: jpfld = 21! maximum number of fields to read54 INTEGER , PARAMETER :: jpfld = 15 ! maximum number of fields to read 56 55 INTEGER , SAVE :: jf_tem ! index of temperature 57 56 INTEGER , SAVE :: jf_sal ! index of salinity … … 68 67 INTEGER , SAVE :: jf_ubl ! index of u-bbl coef 69 68 INTEGER , SAVE :: jf_vbl ! index of v-bbl coef 70 INTEGER , SAVE :: jf_ahu ! index of u-diffusivity coef71 INTEGER , SAVE :: jf_ahv ! index of v-diffusivity coef72 INTEGER , SAVE :: jf_ahw ! index of w-diffusivity coef73 INTEGER , SAVE :: jf_eiu ! index of u-eiv74 INTEGER , SAVE :: jf_eiv ! index of v-eiv75 INTEGER , SAVE :: jf_eiw ! index of w-eiv76 69 INTEGER , SAVE :: jf_fmf ! index of downward salt flux 77 70 … … 112 105 !! - interpolates data if needed 113 106 !!---------------------------------------------------------------------- 114 ! 115 USE oce, ONLY: zts => tsa 107 USE oce, ONLY: zts => tsa 116 108 USE oce, ONLY: zuslp => ua , zvslp => va 117 USE oce, ONLY: zwslpi => rotb , zwslpj => rotn118 USE oce, ONLY: zu => ub , zv => vb, zw => hdivb109 USE oce, ONLY: zwslpi => ua_sv , zwslpj => va_sv 110 USE oce, ONLY: zu => ub , zv => vb, zw => rke 119 111 ! 120 112 INTEGER, INTENT(in) :: kt ! ocean time-step index 113 ! 114 ! REAL(wp), DIMENSION(jpi,jpj,jpk,jpts) :: zts 115 ! REAL(wp), DIMENSION(jpi,jpj,jpk ) :: zuslp, zvslp, zwslpi, zwslpj 116 ! REAL(wp), DIMENSION(jpi,jpj,jpk ) :: zu, zv, zw 117 ! 121 118 ! 122 119 INTEGER :: ji, jj ! dummy loop indices … … 138 135 CALL fld_read( kt, 1, sf_dyn ) !== read data at kt time step ==! 139 136 ! 140 IF( l k_ldfslp .AND. .NOT.lk_c1d .AND. sf_dyn(jf_tem)%ln_tint ) THEN ! Computes slopes (here avt is used as workspace)137 IF( l_ldfslp .AND. .NOT.lk_c1d .AND. sf_dyn(jf_tem)%ln_tint ) THEN ! Computes slopes (here avt is used as workspace) 141 138 zts(:,:,:,jp_tem) = sf_dyn(jf_tem)%fdta(:,:,:,1) * tmask(:,:,:) ! temperature 142 139 zts(:,:,:,jp_sal) = sf_dyn(jf_sal)%fdta(:,:,:,1) * tmask(:,:,:) ! salinity … … 162 159 ENDIF 163 160 ! 164 IF( l k_ldfslp .AND. .NOT.lk_c1d ) THEN ! Computes slopes (here avt is used as workspace)161 IF( l_ldfslp .AND. .NOT.lk_c1d ) THEN ! Computes slopes (here avt is used as workspace) 165 162 iswap_tem = 0 166 163 IF( kt /= nit000 .AND. ( sf_dyn(jf_tem)%nrec_a(2) - nrecprev_tem ) /= 0 ) iswap_tem = 1 … … 267 264 rnf (:,:) = sf_dyn(jf_rnf)%fnow(:,:,1) * tmask(:,:,1) ! river runoffs 268 265 266 ! ! update eddy diffusivity coeff. and/or eiv coeff. at kt 267 IF( l_ldftra_time .OR. l_ldfeiv_time ) CALL ldf_tra( kt ) 269 268 ! ! bbl diffusive coef 270 269 #if defined key_trabbl && ! defined key_c1d … … 276 275 CALL bbl( kt, nit000, 'TRC') 277 276 END IF 278 #endif279 #if ( ! defined key_degrad && defined key_traldf_c2d && defined key_traldf_eiv ) && ! defined key_c1d280 aeiw(:,:) = sf_dyn(jf_eiw)%fnow(:,:,1) * tmask(:,:,1) ! w-eiv281 ! ! Computes the horizontal values from the vertical value282 DO jj = 2, jpjm1283 DO ji = fs_2, fs_jpim1 ! vector opt.284 aeiu(ji,jj) = .5 * ( aeiw(ji,jj) + aeiw(ji+1,jj ) ) ! Average the diffusive coefficient at u- v- points285 aeiv(ji,jj) = .5 * ( aeiw(ji,jj) + aeiw(ji ,jj+1) ) ! at u- v- points286 END DO287 END DO288 CALL lbc_lnk( aeiu, 'U', 1. ) ; CALL lbc_lnk( aeiv, 'V', 1. ) ! lateral boundary condition289 #endif290 291 #if defined key_degrad && ! defined key_c1d292 ! ! degrad option : diffusive and eiv coef are 3D293 ahtu(:,:,:) = sf_dyn(jf_ahu)%fnow(:,:,:) * umask(:,:,:)294 ahtv(:,:,:) = sf_dyn(jf_ahv)%fnow(:,:,:) * vmask(:,:,:)295 ahtw(:,:,:) = sf_dyn(jf_ahw)%fnow(:,:,:) * tmask(:,:,:)296 # if defined key_traldf_eiv297 aeiu(:,:,:) = sf_dyn(jf_eiu)%fnow(:,:,:) * umask(:,:,:)298 aeiv(:,:,:) = sf_dyn(jf_eiv)%fnow(:,:,:) * vmask(:,:,:)299 aeiw(:,:,:) = sf_dyn(jf_eiw)%fnow(:,:,:) * tmask(:,:,:)300 # endif301 277 #endif 302 278 ! … … 339 315 TYPE(FLD_N), DIMENSION(jpfld) :: slf_d ! array of namelist informations on the fields to read 340 316 TYPE(FLD_N) :: sn_tem, sn_sal, sn_mld, sn_emp, sn_ice, sn_qsr, sn_wnd, sn_rnf ! informations about the fields to be read 341 TYPE(FLD_N) :: sn_uwd, sn_vwd, sn_wwd, sn_avt, sn_ubl, sn_vbl ! " " 342 TYPE(FLD_N) :: sn_ahu, sn_ahv, sn_ahw, sn_eiu, sn_eiv, sn_eiw, sn_fmf ! " " 343 !!---------------------------------------------------------------------- 344 ! 345 NAMELIST/namdta_dyn/cn_dir, ln_dynwzv, ln_dynbbl, ln_degrad, ln_dynrnf, & 317 TYPE(FLD_N) :: sn_uwd, sn_vwd, sn_wwd, sn_avt, sn_ubl, sn_vbl, sn_fmf ! " " 318 NAMELIST/namdta_dyn/cn_dir, ln_dynwzv, ln_dynbbl, ln_dynrnf, & 346 319 & sn_tem, sn_sal, sn_mld, sn_emp, sn_ice, sn_qsr, sn_wnd, sn_rnf, & 347 & sn_uwd, sn_vwd, sn_wwd, sn_avt, sn_ubl, sn_vbl, &348 & sn_ahu, sn_ahv, sn_ahw, sn_eiu, sn_eiv, sn_eiw, sn_fmf320 & sn_uwd, sn_vwd, sn_wwd, sn_avt, sn_ubl, sn_vbl, sn_fmf 321 !!---------------------------------------------------------------------- 349 322 ! 350 323 REWIND( numnam_ref ) ! Namelist namdta_dyn in reference namelist : Offline: init. of dynamical data … … 365 338 WRITE(numout,*) ' vertical velocity read from file (T) or computed (F) ln_dynwzv = ', ln_dynwzv 366 339 WRITE(numout,*) ' bbl coef read from file (T) or computed (F) ln_dynbbl = ', ln_dynbbl 367 WRITE(numout,*) ' degradation option enabled (T) or not (F) ln_degrad = ', ln_degrad368 340 WRITE(numout,*) ' river runoff option enabled (T) or not (F) ln_dynrnf = ', ln_dynrnf 369 341 WRITE(numout,*) 370 342 ENDIF 371 343 ! 372 IF( ln_degrad .AND. .NOT.lk_degrad ) THEN373 CALL ctl_warn( 'dta_dyn_init: degradation option requires key_degrad activated ; force ln_degrad to false' )374 ln_degrad = .FALSE.375 ENDIF376 344 IF( ln_dynbbl .AND. ( .NOT.lk_trabbl .OR. lk_c1d ) ) THEN 377 345 CALL ctl_warn( 'dta_dyn_init: bbl option requires key_trabbl activated ; force ln_dynbbl to false' ) … … 395 363 ENDIF 396 364 397 ! 398 IF( .NOT.ln_degrad ) THEN ! no degrad option 399 IF( lk_traldf_eiv .AND. ln_dynbbl ) THEN ! eiv & bbl 400 jf_ubl = jfld + 1 ; jf_vbl = jfld + 2 ; jf_eiw = jfld + 3 ; jfld = jf_eiw 401 slf_d(jf_ubl) = sn_ubl ; slf_d(jf_vbl) = sn_vbl ; slf_d(jf_eiw) = sn_eiw 402 ENDIF 403 IF( .NOT.lk_traldf_eiv .AND. ln_dynbbl ) THEN ! no eiv & bbl 365 IF( ln_dynbbl ) THEN ! eiv & bbl 404 366 jf_ubl = jfld + 1 ; jf_vbl = jfld + 2 ; jfld = jf_vbl 405 367 slf_d(jf_ubl) = sn_ubl ; slf_d(jf_vbl) = sn_vbl 406 ENDIF 407 IF( lk_traldf_eiv .AND. .NOT.ln_dynbbl ) THEN ! eiv & no bbl 408 jf_eiw = jfld + 1 ; jfld = jf_eiw ; slf_d(jf_eiw) = sn_eiw 409 ENDIF 410 ELSE 411 jf_ahu = jfld + 1 ; jf_ahv = jfld + 2 ; jf_ahw = jfld + 3 ; jfld = jf_ahw 412 slf_d(jf_ahu) = sn_ahu ; slf_d(jf_ahv) = sn_ahv ; slf_d(jf_ahw) = sn_ahw 413 IF( lk_traldf_eiv .AND. ln_dynbbl ) THEN ! eiv & bbl 414 jf_ubl = jfld + 1 ; jf_vbl = jfld + 2 ; 415 slf_d(jf_ubl) = sn_ubl ; slf_d(jf_vbl) = sn_vbl 416 jf_eiu = jfld + 3 ; jf_eiv = jfld + 4 ; jf_eiw = jfld + 5 ; jfld = jf_eiw 417 slf_d(jf_eiu) = sn_eiu ; slf_d(jf_eiv) = sn_eiv ; slf_d(jf_eiw) = sn_eiw 418 ENDIF 419 IF( .NOT.lk_traldf_eiv .AND. ln_dynbbl ) THEN ! no eiv & bbl 420 jf_ubl = jfld + 1 ; jf_vbl = jfld + 2 ; jfld = jf_vbl 421 slf_d(jf_ubl) = sn_ubl ; slf_d(jf_vbl) = sn_vbl 422 ENDIF 423 IF( lk_traldf_eiv .AND. .NOT.ln_dynbbl ) THEN ! eiv & no bbl 424 jf_eiu = jfld + 1 ; jf_eiv = jfld + 2 ; jf_eiw = jfld + 3 ; jfld = jf_eiw 425 slf_d(jf_eiu) = sn_eiu ; slf_d(jf_eiv) = sn_eiv ; slf_d(jf_eiw) = sn_eiw 426 ENDIF 427 ENDIF 428 368 ENDIF 369 370 429 371 ALLOCATE( sf_dyn(jfld), STAT=ierr ) ! set sf structure 430 372 IF( ierr > 0 ) THEN … … 452 394 END DO 453 395 ! 454 IF( l k_ldfslp .AND. .NOT.lk_c1d ) THEN ! slopes396 IF( l_ldfslp .AND. .NOT.lk_c1d ) THEN ! slopes 455 397 IF( sf_dyn(jf_tem)%ln_tint ) THEN ! time interpolation 456 398 ALLOCATE( uslpdta (jpi,jpj,jpk,2), vslpdta (jpi,jpj,jpk,2), & … … 511 453 zv = pv(ji ,jj ,jk) * vmask(ji ,jj ,jk) * e1v(ji ,jj ) * fse3v(ji ,jj ,jk) 512 454 zv1 = pv(ji ,jj-1,jk) * vmask(ji ,jj-1,jk) * e1v(ji ,jj-1) * fse3v(ji ,jj-1,jk) 513 zet = 1. / ( e1 t(ji,jj) *e2t(ji,jj) * fse3t(ji,jj,jk) )455 zet = 1. / ( e1e2t(ji,jj) * fse3t(ji,jj,jk) ) 514 456 zhdiv(ji,jj,jk) = ( zu - zu1 + zv - zv1 ) * zet 515 457 END DO 516 458 END DO 517 459 END DO 460 ! ! update the horizontal divergence with the runoff inflow 461 IF( ln_dynrnf ) zhdiv(:,:,1) = zhdiv(:,:,1) - rnf(:,:) * r1_rau0 / fse3t(:,:,1) 462 ! 518 463 CALL lbc_lnk( zhdiv, 'T', 1. ) ! Lateral boundary conditions on zhdiv 519 !520 464 ! computation of vertical velocity from the bottom 521 465 pw(:,:,jpk) = 0._wp … … 540 484 REAL(wp), DIMENSION(jpi,jpj,jpk) , INTENT(out) :: pwslpj ! meridional diapycnal slopes 541 485 !!--------------------------------------------------------------------- 542 #if defined key_ldfslp && ! defined key_c1d 543 CALL eos ( pts, rhd, rhop, gdept_0(:,:,:) )544 CALL eos_rab( pts, rab_n ) ! now local thermal/haline expension ratio at T-points545 CALL bn2 ( pts, rab_n, rn2 ) ! now Brunt-Vaisala546 547 ! Partial steps: before Horizontal DErivative548 IF( ln_zps .AND. .NOT. ln_isfcav) &549 & CALL zps_hde ( kt, jpts, pts, gtsu, gtsv, & ! Partial steps: before horizontal gradient550 & rhd, gru , grv ) ! of t, s, rd at the last ocean level551 IF( ln_zps .AND. ln_isfcav) &552 & CALL zps_hde_isf( kt, jpts, pts, gtsu, gtsv, & ! Partial steps for top cell (ISF)553 & rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv , &554 & gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi ) ! of t, s, rd at the first ocean level555 556 rn2b(:,:,:) = rn2(:,:,:) ! need for zdfmxl557 CALL zdf_mxl( kt ) ! mixed layer depth558 CALL ldf_slp( kt, rhd, rn2 ) ! slopes559 puslp (:,:,:) = uslp (:,:,:)560 pvslp (:,:,:) = vslp (:,:,:)561 pwslpi(:,:,:) = wslpi(:,:,:)562 pwslpj(:,:,:) = wslpj(:,:,:)563 #else 564 puslp (:,:,:) = 0. ! to avoid warning when compiling565 pvslp (:,:,:) = 0.566 pwslpi(:,:,:) = 0.567 pwslpj(:,:,:) = 0.568 #endif 486 IF( l_ldfslp .AND. .NOT.lk_c1d ) THEN ! Computes slopes (here avt is used as workspace) 487 CALL eos ( pts, rhd, rhop, gdept_0(:,:,:) ) 488 CALL eos_rab( pts, rab_n ) ! now local thermal/haline expension ratio at T-points 489 CALL bn2 ( pts, rab_n, rn2 ) ! now Brunt-Vaisala 490 491 ! Partial steps: before Horizontal DErivative 492 IF( ln_zps .AND. .NOT. ln_isfcav) & 493 & CALL zps_hde ( kt, jpts, pts, gtsu, gtsv, & ! Partial steps: before horizontal gradient 494 & rhd, gru , grv ) ! of t, s, rd at the last ocean level 495 IF( ln_zps .AND. ln_isfcav) & 496 & CALL zps_hde_isf( kt, jpts, pts, gtsu, gtsv, & ! Partial steps for top cell (ISF) 497 & rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv , & 498 & gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi ) ! of t, s, rd at the first ocean level 499 500 rn2b(:,:,:) = rn2(:,:,:) ! need for zdfmxl 501 CALL zdf_mxl( kt ) ! mixed layer depth 502 CALL ldf_slp( kt, rhd, rn2 ) ! slopes 503 puslp (:,:,:) = uslp (:,:,:) 504 pvslp (:,:,:) = vslp (:,:,:) 505 pwslpi(:,:,:) = wslpi(:,:,:) 506 pwslpj(:,:,:) = wslpj(:,:,:) 507 ELSE 508 puslp (:,:,:) = 0. ! to avoid warning when compiling 509 pvslp (:,:,:) = 0. 510 pwslpi(:,:,:) = 0. 511 pwslpj(:,:,:) = 0. 512 ENDIF 569 513 ! 570 514 END SUBROUTINE dta_dyn_slp -
trunk/NEMOGCM/NEMO/OFF_SRC/nemogcm.F90
r5504 r5836 26 26 USE traqsr ! solar radiation penetration (tra_qsr_init routine) 27 27 USE trabbl ! bottom boundary layer (tra_bbl_init routine) 28 USE traldf ! lateral physics (tra_ldf_init routine) 28 29 USE zdfini ! vertical physics: initialization 29 30 USE sbcmod ! surface boundary condition (sbc_init routine) … … 283 284 CALL sbc_init ! Forcings : surface module 284 285 285 #if ! defined key_degrad286 286 CALL ldf_tra_init ! Lateral ocean tracer physics 287 #endif 288 IF( lk_ldfslp ) CALL ldf_slp_init ! slope of lateral mixing 287 CALL ldf_eiv_init ! Eddy induced velocity param 288 CALL tra_ldf_init ! lateral mixing 289 IF( l_ldfslp ) CALL ldf_slp_init ! slope of lateral mixing 289 290 290 291 CALL tra_qsr_init ! penetrative solar radiation qsr … … 444 445 USE dom_oce, ONLY: dom_oce_alloc 445 446 USE zdf_oce, ONLY: zdf_oce_alloc 446 USE ldftra_oce, ONLY: ldftra_oce_alloc447 447 USE trc_oce, ONLY: trc_oce_alloc 448 448 ! … … 453 453 ierr = ierr + dia_wri_alloc () 454 454 ierr = ierr + dom_oce_alloc () ! ocean domain 455 ierr = ierr + ldftra_oce_alloc() ! ocean lateral physics : tracers456 455 ierr = ierr + zdf_oce_alloc () ! ocean vertical physics 457 456 ! -
trunk/NEMOGCM/NEMO/OPA_SRC/ASM/asmbkg.F90
r5215 r5836 18 18 19 19 !!---------------------------------------------------------------------- 20 !! 'key_asminc' : Switch on the assimilation increment interface21 !!----------------------------------------------------------------------22 20 !! asm_bkg_wri : Write out the background state 23 21 !! asm_trj_wri : Write out the model state trajectory (used with 4D-Var) … … 27 25 USE zdf_oce ! Vertical mixing variables 28 26 USE zdfddm ! Double diffusion mixing parameterization 29 USE ldftra _oce ! Lateral tracer mixing coefficient defined in memory30 USE ldfslp ! Slopes of neutral surfaces27 USE ldftra ! Lateral diffusion: eddy diffusivity coefficients 28 USE ldfslp ! Lateral diffusion: slopes of neutral surfaces 31 29 USE tradmp ! Tracer damping 32 30 #if defined key_zdftke … … 41 39 USE asmpar ! Parameters for the assmilation interface 42 40 USE zdfmxl ! mixed layer depth 43 #if defined key_traldf_c2d44 USE ldfeiv ! eddy induced velocity coef. (ldf_eiv routine)45 #endif46 41 #if defined key_lim2 47 42 USE ice_2 … … 155 150 CALL iom_rstput( kt, nitdin_r, inum, 'sshn' , sshn ) 156 151 #if defined key_lim2 || defined key_lim3 157 IF( ( nn_ice == 2 ) .OR. ( nn_ice == 3 )) THEN158 IF(ALLOCATED(frld)) THEN159 CALL iom_rstput( kt, nitdin_r, inum, 'iceconc', 1. 0- frld(:,:) )152 IF( nn_ice == 2 .OR. nn_ice == 3 ) THEN 153 IF( ALLOCATED(frld) ) THEN 154 CALL iom_rstput( kt, nitdin_r, inum, 'iceconc', 1._wp - frld(:,:) ) 160 155 ELSE 161 CALL ctl_warn('Ice concentration not written to background as ice variable frld not allocated on this timestep')162 ENDIF156 CALL ctl_warn('Ice concentration not written to background as ice variable frld not allocated on this timestep') 157 ENDIF 163 158 ENDIF 164 159 #endif -
trunk/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90
r5541 r5836 14 14 15 15 !!---------------------------------------------------------------------- 16 !! 'key_asminc' : Switch on the assimilation increment interface17 !!----------------------------------------------------------------------18 16 !! asm_inc_init : Initialize the increment arrays and IAU weights 19 17 !! calc_date : Compute the calendar date YYYYMMDD on a given step … … 28 26 USE domvvl ! domain: variable volume level 29 27 USE oce ! Dynamics and active tracers defined in memory 30 USE ldfdyn _oce ! ocean dynamics: lateral physics28 USE ldfdyn ! lateral diffusion: eddy viscosity coefficients 31 29 USE eosbn2 ! Equation of state - in situ and potential density 32 30 USE zpshde ! Partial step : Horizontal Derivative … … 56 54 LOGICAL, PUBLIC, PARAMETER :: lk_asminc = .FALSE. !: No assimilation increments 57 55 #endif 58 LOGICAL, PUBLIC :: ln_bkgwri = .FALSE.!: No output of the background state fields59 LOGICAL, PUBLIC :: ln_asmiau = .FALSE.!: No applying forcing with an assimilation increment60 LOGICAL, PUBLIC :: ln_asmdin = .FALSE.!: No direct initialization61 LOGICAL, PUBLIC :: ln_trainc = .FALSE.!: No tracer (T and S) assimilation increments62 LOGICAL, PUBLIC :: ln_dyninc = .FALSE.!: No dynamics (u and v) assimilation increments63 LOGICAL, PUBLIC :: ln_sshinc = .FALSE.!: No sea surface height assimilation increment64 LOGICAL, PUBLIC :: ln_seaiceinc 65 LOGICAL, PUBLIC :: ln_salfix = .FALSE.!: Apply minimum salinity check56 LOGICAL, PUBLIC :: ln_bkgwri !: No output of the background state fields 57 LOGICAL, PUBLIC :: ln_asmiau !: No applying forcing with an assimilation increment 58 LOGICAL, PUBLIC :: ln_asmdin !: No direct initialization 59 LOGICAL, PUBLIC :: ln_trainc !: No tracer (T and S) assimilation increments 60 LOGICAL, PUBLIC :: ln_dyninc !: No dynamics (u and v) assimilation increments 61 LOGICAL, PUBLIC :: ln_sshinc !: No sea surface height assimilation increment 62 LOGICAL, PUBLIC :: ln_seaiceinc !: No sea ice concentration increment 63 LOGICAL, PUBLIC :: ln_salfix !: Apply minimum salinity check 66 64 LOGICAL, PUBLIC :: ln_temnofreeze = .FALSE. !: Don't allow the temperature to drop below freezing 67 INTEGER, PUBLIC :: nn_divdmp 65 INTEGER, PUBLIC :: nn_divdmp !: Apply divergence damping filter nn_divdmp times 68 66 69 67 REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: t_bkg , s_bkg !: Background temperature and salinity … … 90 88 !! * Substitutions 91 89 # include "domzgr_substitute.h90" 92 # include "ldfdyn_substitute.h90"93 90 # include "vectopt_loop_substitute.h90" 94 91 !!---------------------------------------------------------------------- … … 139 136 ! Read Namelist nam_asminc : assimilation increment interface 140 137 !----------------------------------------------------------------------- 141 ln_seaiceinc = .FALSE.138 ln_seaiceinc = .FALSE. 142 139 ln_temnofreeze = .FALSE. 143 140 … … 428 425 429 426 IF ( ln_dyninc .AND. nn_divdmp > 0 ) THEN 430 431 CALL wrk_alloc( jpi,jpj,hdiv)432 433 DO 434 427 ! 428 CALL wrk_alloc( jpi,jpj, hdiv ) 429 ! 430 DO jt = 1, nn_divdmp 431 ! 435 432 DO jk = 1, jpkm1 436 437 433 hdiv(:,:) = 0._wp 438 439 434 DO jj = 2, jpjm1 440 435 DO ji = fs_2, fs_jpim1 ! vector opt. … … 444 439 + e1v(ji ,jj ) * fse3v(ji ,jj ,jk) * v_bkginc(ji ,jj ,jk) & 445 440 - e1v(ji ,jj-1) * fse3v(ji ,jj-1,jk) * v_bkginc(ji ,jj-1,jk) ) & 446 / ( e1 t(ji,jj) *e2t(ji,jj) * fse3t(ji,jj,jk) )441 / ( e1e2t(ji,jj) * fse3t(ji,jj,jk) ) 447 442 END DO 448 443 END DO 449 450 444 CALL lbc_lnk( hdiv, 'T', 1. ) ! lateral boundary cond. (no sign change) 451 445 ! 452 446 DO jj = 2, jpjm1 453 447 DO ji = fs_2, fs_jpim1 ! vector opt. 454 u_bkginc(ji,jj,jk) = u_bkginc(ji,jj,jk) + 0.2_wp * ( e1 t(ji+1,jj)*e2t(ji+1,jj) * hdiv(ji+1,jj) &455 - e1t(ji ,jj)*e2t(ji ,jj) * hdiv(ji ,jj) ) &456 /e1u(ji,jj) * umask(ji,jj,jk)457 v_bkginc(ji,jj,jk) = v_bkginc(ji,jj,jk) + 0.2_wp * ( e1 t(ji,jj+1)*e2t(ji,jj+1) * hdiv(ji,jj+1) &458 - e1t(ji,jj )*e2t(ji,jj ) * hdiv(ji,jj ) ) &459 /e2v(ji,jj) * vmask(ji,jj,jk)448 u_bkginc(ji,jj,jk) = u_bkginc(ji,jj,jk) + 0.2_wp * ( e1e2t(ji+1,jj) * hdiv(ji+1,jj) & 449 & - e1e2t(ji ,jj) * hdiv(ji ,jj) ) & 450 & * r1_e1u(ji,jj) * umask(ji,jj,jk) 451 v_bkginc(ji,jj,jk) = v_bkginc(ji,jj,jk) + 0.2_wp * ( e1e2t(ji,jj+1) * hdiv(ji,jj+1) & 452 & - e1e2t(ji,jj ) * hdiv(ji,jj ) ) & 453 & * r1_e2v(ji,jj) * vmask(ji,jj,jk) 460 454 END DO 461 455 END DO 462 463 456 END DO 464 457 ! 465 458 END DO 466 467 CALL wrk_dealloc( jpi,jpj,hdiv)468 459 ! 460 CALL wrk_dealloc( jpi,jpj, hdiv ) 461 ! 469 462 ENDIF 470 471 472 463 473 464 !----------------------------------------------------------------------- … … 476 467 477 468 IF ( ln_asmdin ) THEN 478 469 ! 479 470 ALLOCATE( t_bkg(jpi,jpj,jpk) ) 480 471 ALLOCATE( s_bkg(jpi,jpj,jpk) ) … … 482 473 ALLOCATE( v_bkg(jpi,jpj,jpk) ) 483 474 ALLOCATE( ssh_bkg(jpi,jpj) ) 484 485 t_bkg(:,:,:) = 0. 0486 s_bkg(:,:,:) = 0. 0487 u_bkg(:,:,:) = 0. 0488 v_bkg(:,:,:) = 0. 0489 ssh_bkg(:,:) = 0. 0490 475 ! 476 t_bkg(:,:,:) = 0._wp 477 s_bkg(:,:,:) = 0._wp 478 u_bkg(:,:,:) = 0._wp 479 v_bkg(:,:,:) = 0._wp 480 ssh_bkg(:,:) = 0._wp 481 ! 491 482 !-------------------------------------------------------------------- 492 483 ! Read from file the background state at analysis time 493 484 !-------------------------------------------------------------------- 494 485 ! 495 486 CALL iom_open( c_asmdin, inum ) 496 487 ! 497 488 CALL iom_get( inum, 'rdastp', zdate_bkg ) 498 489 ! 499 490 IF(lwp) THEN 500 491 WRITE(numout,*) 501 WRITE(numout,*) 'asm_inc_init : Assimilation background state valid at : ', & 502 & NINT( zdate_bkg ) 492 WRITE(numout,*) 'asm_inc_init : Assimilation background state valid at : ', NINT( zdate_bkg ) 503 493 WRITE(numout,*) '~~~~~~~~~~~~' 504 494 ENDIF 505 495 ! 506 496 IF ( NINT( zdate_bkg ) /= iitdin_date ) & 507 497 & CALL ctl_warn( ' Validity time of assimilation background state does', & 508 498 & ' not agree with Direct Initialization time' ) 509 499 ! 510 500 IF ( ln_trainc ) THEN 511 501 CALL iom_get( inum, jpdom_autoglo, 'tn', t_bkg ) … … 514 504 s_bkg(:,:,:) = s_bkg(:,:,:) * tmask(:,:,:) 515 505 ENDIF 516 506 ! 517 507 IF ( ln_dyninc ) THEN 518 508 CALL iom_get( inum, jpdom_autoglo, 'un', u_bkg ) … … 521 511 v_bkg(:,:,:) = v_bkg(:,:,:) * vmask(:,:,:) 522 512 ENDIF 523 513 ! 524 514 IF ( ln_sshinc ) THEN 525 515 CALL iom_get( inum, jpdom_autoglo, 'sshn', ssh_bkg ) 526 516 ssh_bkg(:,:) = ssh_bkg(:,:) * tmask(:,:,1) 527 517 ENDIF 528 518 ! 529 519 CALL iom_close( inum ) 530 520 ! 531 521 ENDIF 532 522 ! … … 574 564 ! If kt = kit000 - 1 then set the date to the restart date 575 565 IF ( kt == kit000 - 1 ) THEN 576 577 566 kdate = ndastp 578 567 RETURN 579 580 568 ENDIF 581 569 … … 646 634 !! ** Action : 647 635 !!---------------------------------------------------------------------- 648 INTEGER, INTENT(IN) :: kt! Current time step649 ! 650 INTEGER :: ji,jj,jk651 INTEGER :: it636 INTEGER, INTENT(IN) :: kt ! Current time step 637 ! 638 INTEGER :: ji, jj, jk 639 INTEGER :: it 652 640 REAL(wp) :: zincwgt ! IAU weight for current time step 653 641 REAL (wp), DIMENSION(jpi,jpj,jpk) :: fzptnz ! 3d freezing point values 654 642 !!---------------------------------------------------------------------- 655 643 ! 656 644 ! freezing point calculation taken from oc_fz_pt (but calculated for all depths) 657 645 ! used to prevent the applied increments taking the temperature below the local freezing point 658 659 646 DO jk = 1, jpkm1 660 647 CALL eos_fzp( tsn(:,:,jk,jp_sal), fzptnz(:,:,jk), fsdept(:,:,jk) ) 661 648 END DO 662 663 IF ( ln_asmiau ) THEN 664 665 !-------------------------------------------------------------------- 666 ! Incremental Analysis Updating 667 !-------------------------------------------------------------------- 668 649 ! 650 ! !-------------------------------------- 651 IF ( ln_asmiau ) THEN ! Incremental Analysis Updating 652 ! !-------------------------------------- 653 ! 669 654 IF ( ( kt >= nitiaustr_r ).AND.( kt <= nitiaufin_r ) ) THEN 670 655 ! 671 656 it = kt - nit000 + 1 672 657 zincwgt = wgtiau(it) / rdt ! IAU weight for the current time step 673 658 ! 674 659 IF(lwp) THEN 675 660 WRITE(numout,*) … … 677 662 WRITE(numout,*) '~~~~~~~~~~~~' 678 663 ENDIF 679 664 ! 680 665 ! Update the tracer tendencies 681 666 DO jk = 1, jpkm1 … … 700 685 ENDIF 701 686 END DO 702 703 ENDIF 704 687 ! 688 ENDIF 689 ! 705 690 IF ( kt == nitiaufin_r + 1 ) THEN ! For bias crcn to work 706 691 DEALLOCATE( t_bkginc ) 707 692 DEALLOCATE( s_bkginc ) 708 693 ENDIF 709 710 711 ELSEIF ( ln_asmdin ) THEN 712 713 !-------------------------------------------------------------------- 714 ! Direct Initialization 715 !-------------------------------------------------------------------- 716 694 ! !-------------------------------------- 695 ELSEIF ( ln_asmdin ) THEN ! Direct Initialization 696 ! !-------------------------------------- 697 ! 717 698 IF ( kt == nitdin_r ) THEN 718 699 ! 719 700 neuler = 0 ! Force Euler forward step 720 701 ! 721 702 ! Initialize the now fields with the background + increment 722 703 IF (ln_temnofreeze) THEN … … 745 726 !!gm 746 727 747 748 IF( ln_zps .AND. .NOT. lk_c1d .AND. .NOT. ln_isfcav) & 749 & CALL zps_hde ( kt, jpts, tsb, gtsu, gtsv, & ! Partial steps: before horizontal gradient 750 & rhd, gru , grv ) ! of t, s, rd at the last ocean level 751 IF( ln_zps .AND. .NOT. lk_c1d .AND. ln_isfcav) & 752 & CALL zps_hde_isf( nit000, jpts, tsb, gtsu, gtsv, & ! Partial steps for top cell (ISF) 753 & rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv , & 754 & gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi ) ! of t, s, rd at the last ocean level 755 756 #if defined key_zdfkpp 757 CALL eos( tsn, rhd, fsdept_n(:,:,:) ) ! Compute rhd 758 !!gm fabien CALL eos( tsn, rhd ) ! Compute rhd 759 #endif 760 728 IF( ln_zps .AND. .NOT. lk_c1d ) THEN ! Partial steps: before horizontal gradient 729 IF(ln_isfcav) THEN ! ocean cavities: top and bottom cells (ISF) 730 CALL zps_hde_isf( nit000, jpts, tsb, gtsu, gtsv, gtui, gtvi, & 731 & rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv , & 732 & grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi ) 733 ELSE ! no ocean cavities: bottom cells 734 CALL zps_hde ( kt, jpts, tsb, gtsu, gtsv, & ! 735 & rhd, gru , grv ) ! of t, s, rd at the last ocean level 736 ENDIF 737 ENDIF 738 ! 761 739 DEALLOCATE( t_bkginc ) 762 740 DEALLOCATE( s_bkginc ) … … 767 745 ENDIF 768 746 ! Perhaps the following call should be in step 769 IF 747 IF ( ln_seaiceinc ) CALL seaice_asm_inc ( kt ) ! apply sea ice concentration increment 770 748 ! 771 749 END SUBROUTINE tra_asm_inc … … 788 766 REAL(wp) :: zincwgt ! IAU weight for current time step 789 767 !!---------------------------------------------------------------------- 790 791 IF ( ln_asmiau ) THEN 792 793 !-------------------------------------------------------------------- 794 ! Incremental Analysis Updating 795 !-------------------------------------------------------------------- 796 768 ! 769 ! !-------------------------------------------- 770 IF ( ln_asmiau ) THEN ! Incremental Analysis Updating 771 ! !-------------------------------------------- 772 ! 797 773 IF ( ( kt >= nitiaustr_r ).AND.( kt <= nitiaufin_r ) ) THEN 798 774 ! 799 775 it = kt - nit000 + 1 800 776 zincwgt = wgtiau(it) / rdt ! IAU weight for the current time step 801 777 ! 802 778 IF(lwp) THEN 803 779 WRITE(numout,*) 804 WRITE(numout,*) 'dyn_asm_inc : Dynamics IAU at time step = ', & 805 & kt,' with IAU weight = ', wgtiau(it) 780 WRITE(numout,*) 'dyn_asm_inc : Dynamics IAU at time step = ', kt,' with IAU weight = ', wgtiau(it) 806 781 WRITE(numout,*) '~~~~~~~~~~~~' 807 782 ENDIF 808 783 ! 809 784 ! Update the dynamic tendencies 810 785 DO jk = 1, jpkm1 … … 812 787 va(:,:,jk) = va(:,:,jk) + v_bkginc(:,:,jk) * zincwgt 813 788 END DO 814 789 ! 815 790 IF ( kt == nitiaufin_r ) THEN 816 791 DEALLOCATE( u_bkginc ) 817 792 DEALLOCATE( v_bkginc ) 818 793 ENDIF 819 820 ENDIF 821 822 ELSEIF ( ln_asmdin ) THEN 823 824 !-------------------------------------------------------------------- 825 ! Direct Initialization 826 !-------------------------------------------------------------------- 827 794 ! 795 ENDIF 796 ! !----------------------------------------- 797 ELSEIF ( ln_asmdin ) THEN ! Direct Initialization 798 ! !----------------------------------------- 799 ! 828 800 IF ( kt == nitdin_r ) THEN 829 801 ! 830 802 neuler = 0 ! Force Euler forward step 831 803 ! 832 804 ! Initialize the now fields with the background + increment 833 805 un(:,:,:) = u_bkg(:,:,:) + u_bkginc(:,:,:) 834 806 vn(:,:,:) = v_bkg(:,:,:) + v_bkginc(:,:,:) 835 807 ! 836 808 ub(:,:,:) = un(:,:,:) ! Update before fields 837 809 vb(:,:,:) = vn(:,:,:) 838 810 ! 839 811 DEALLOCATE( u_bkg ) 840 812 DEALLOCATE( v_bkg ) … … 864 836 REAL(wp) :: zincwgt ! IAU weight for current time step 865 837 !!---------------------------------------------------------------------- 866 867 IF ( ln_asmiau ) THEN 868 869 !-------------------------------------------------------------------- 870 ! Incremental Analysis Updating 871 !-------------------------------------------------------------------- 872 838 ! 839 ! !----------------------------------------- 840 IF ( ln_asmiau ) THEN ! Incremental Analysis Updating 841 ! !----------------------------------------- 842 ! 873 843 IF ( ( kt >= nitiaustr_r ).AND.( kt <= nitiaufin_r ) ) THEN 874 844 ! 875 845 it = kt - nit000 + 1 876 846 zincwgt = wgtiau(it) / rdt ! IAU weight for the current time step 877 847 ! 878 848 IF(lwp) THEN 879 849 WRITE(numout,*) … … 882 852 WRITE(numout,*) '~~~~~~~~~~~~' 883 853 ENDIF 884 854 ! 885 855 ! Save the tendency associated with the IAU weighted SSH increment 886 856 ! (applied in dynspg.*) … … 891 861 DEALLOCATE( ssh_bkginc ) 892 862 ENDIF 893 894 ENDIF 895 896 ELSEIF ( ln_asmdin ) THEN 897 898 !-------------------------------------------------------------------- 899 ! Direct Initialization 900 !-------------------------------------------------------------------- 901 863 ! 864 ENDIF 865 ! !----------------------------------------- 866 ELSEIF ( ln_asmdin ) THEN ! Direct Initialization 867 ! !----------------------------------------- 868 ! 902 869 IF ( kt == nitdin_r ) THEN 903 904 neuler = 0 ! Force Euler forward step 905 906 ! Initialize the now fields the background + increment 907 sshn(:,:) = ssh_bkg(:,:) + ssh_bkginc(:,:) 908 909 ! Update before fields 910 sshb(:,:) = sshn(:,:) 911 870 ! 871 neuler = 0 ! Force Euler forward step 872 ! 873 sshn(:,:) = ssh_bkg(:,:) + ssh_bkginc(:,:) ! Initialize the now fields the background + increment 874 ! 875 sshb(:,:) = sshn(:,:) ! Update before fields 876 ! 912 877 IF( lk_vvl ) THEN 913 878 DO jk = 1, jpk … … 915 880 END DO 916 881 ENDIF 917 882 ! 918 883 DEALLOCATE( ssh_bkg ) 919 884 DEALLOCATE( ssh_bkginc ) 920 885 ! 921 886 ENDIF 922 887 ! … … 937 902 !! 938 903 !!---------------------------------------------------------------------- 939 IMPLICIT NONE 940 ! 941 INTEGER, INTENT(in) :: kt ! Current time step 904 INTEGER, INTENT(in) :: kt ! Current time step 942 905 INTEGER, INTENT(in), OPTIONAL :: kindic ! flag for disabling the deallocation 943 906 !