Changeset 32 for trunk/NEMO/OPA_SRC/DYN
- Timestamp:
- 2004-02-17T10:20:15+01:00 (20 years ago)
- Location:
- trunk/NEMO/OPA_SRC/DYN
- Files:
-
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DYN/dynhpg.F90
r3 r32 31 31 !! 'key_autotasking' : j-k-i loop (j-slab) 32 32 !!---------------------------------------------------------------------- 33 LOGICAL, PUBLIC :: l_dyn_hpg_tsk = .TRUE. ! ???34 LOGICAL, PUBLIC :: l_dyn_hpg = .FALSE. ! ???33 LOGICAL, PUBLIC, PARAMETER :: lk_dynhpg_tsk = .TRUE. !: autotasked hpg flag 34 LOGICAL, PUBLIC, PARAMETER :: lk_dynhpg = .FALSE. !: vector hpg flag 35 35 #else 36 36 !!---------------------------------------------------------------------- 37 37 !! default case : k-j-i loop (vector opt.) 38 38 !!---------------------------------------------------------------------- 39 LOGICAL, PUBLIC :: l_dyn_hpg_tsk = .FALSE. ! ???40 LOGICAL, PUBLIC :: l_dyn_hpg = .TRUE. ! ???39 LOGICAL, PUBLIC, PARAMETER :: lk_dynhpg_tsk = .FALSE. !: autotasked hpg flag 40 LOGICAL, PUBLIC, PARAMETER :: lk_dynhpg = .TRUE. !: vector hpg flag 41 41 #endif 42 42 … … 67 67 !! level. s-coordinates ('key_s_coord'): a corrective term is added 68 68 !! to the horizontal pressure gradient : 69 !! zhpi = g ..... + 1/e1u mi(rhd) di[ gdep3w ]70 !! zhpj = g ..... + 1/e2v mj(rhd) dj[ gdep3w ]69 !! zhpi = grav ..... + 1/e1u mi(rhd) di[ grav dep3w ] 70 !! zhpj = grav ..... + 1/e2v mj(rhd) dj[ grav dep3w ] 71 71 !! add it to the general momentum trend (ua,va). 72 72 !! ua = ua - 1/e1u * zhpi … … 104 104 ! 0. Local constant initialization 105 105 ! -------------------------------- 106 zcoef0 = - g* 0.5106 zcoef0 = - grav * 0.5 107 107 zuap = 0.e0 108 108 zvap = 0.e0 … … 187 187 !! jk is computed by taking the vertical integral of the in-situ 188 188 !! density gradient along the model level from the suface to that 189 !! level: zhpi = g .....190 !! zhpj = g .....189 !! level: zhpi = grav ..... 190 !! zhpj = grav ..... 191 191 !! add it to the general momentum trend (ua,va). 192 192 !! ua = ua - 1/e1u * zhpi … … 223 223 ! 0. Local constant initialization 224 224 ! -------------------------------- 225 zcoef0 = - g* 0.5225 zcoef0 = - grav * 0.5 226 226 zuap = 0.e0 227 227 zvap = 0.e0 … … 342 342 !! jk is computed by taking the vertical integral of the in-situ 343 343 !! density gradient along the model level from the suface to that 344 !! level: zhpi = g .....345 !! zhpj = g .....344 !! level: zhpi = grav ..... 345 !! zhpj = grav ..... 346 346 !! add it to the general momentum trend (ua,va). 347 347 !! ua = ua - 1/e1u * zhpi … … 379 379 ! 0. Local constant initialization 380 380 ! -------------------------------- 381 zcoef0 = - g* 0.5381 zcoef0 = - grav * 0.5 382 382 zuap = 0.e0 383 383 zvap = 0.e0 -
trunk/NEMO/OPA_SRC/DYN/dynhpg_atsk.F90
r3 r32 56 56 !! level. s-coordinate case ('key_s_coord'): a corrective term is 57 57 !! added to the horizontal pressure gradient : 58 !! zhpi = g ..... + 1/e1u mi(rhd) di[ gdep3w ]59 !! zhpj = g ..... + 1/e2v mj(rhd) dj[ gdep3w ]58 !! zhpi = grav ..... + 1/e1u mi(rhd) di[ grav dep3w ] 59 !! zhpj = grav ..... + 1/e2v mj(rhd) dj[ grav dep3w ] 60 60 !! add it to the general momentum trend (ua,va). 61 61 !! ua = ua - 1/e1u * zhpi … … 92 92 ! 0. Local constant initialization 93 93 ! -------------------------------- 94 zcoef0 = - g* 0.594 zcoef0 = - grav * 0.5 95 95 96 96 ! ! =============== … … 174 174 !! jk is computed by taking the vertical integral of the in-situ 175 175 !! density gradient along the model level from the suface to that 176 !! level: zhpi = g .....177 !! zhpj = g .....176 !! level: zhpi = grav ..... 177 !! zhpj = grav ..... 178 178 !! add it to the general momentum trend (ua,va). 179 179 !! ua = ua - 1/e1u * zhpi … … 210 210 ! 0. Local constant initialization 211 211 ! -------------------------------- 212 zcoef0 = - g* 0.5212 zcoef0 = - grav * 0.5 213 213 zuap = 0.e0 214 214 zvap = 0.e0 … … 313 313 !! jk is computed by taking the vertical integral of the in-situ 314 314 !! density gradient along the model level from the suface to that 315 !! level: zhpi = g .....316 !! zhpj = g .....315 !! level: zhpi = grav ..... 316 !! zhpj = grav ..... 317 317 !! add it to the general momentum trend (ua,va). 318 318 !! ua = ua - 1/e1u * zhpi … … 351 351 ! 0. Local constant initialization 352 352 ! -------------------------------- 353 zcoef0 = - g* 0.5353 zcoef0 = - grav * 0.5 354 354 355 355 ! ! =============== -
trunk/NEMO/OPA_SRC/DYN/dynldf_bilapg.F90
r3 r32 465 465 CONTAINS 466 466 SUBROUTINE dyn_ldf_bilapg( kt ) ! Dummy routine 467 WRITE(*,*) kt467 WRITE(*,*) 'dyn_ldf_bilapg: You should not have seen this print! error?', kt 468 468 END SUBROUTINE dyn_ldf_bilapg 469 469 #endif -
trunk/NEMO/OPA_SRC/DYN/dynldf_iso.F90
r3 r32 251 251 CONTAINS 252 252 SUBROUTINE dyn_ldf_iso( kt ) ! Empty routine 253 WRITE(*,*) kt253 WRITE(*,*) 'dyn_ldf_iso: You should not have seen this print! error?', kt 254 254 END SUBROUTINE dyn_ldf_iso 255 255 #endif -
trunk/NEMO/OPA_SRC/DYN/dynnxt.F90
r3 r32 35 35 !! After velocity is compute using a leap-frog scheme environment: 36 36 !! (ua,va) = (ub,vb) + 2 rdt (ua,va) 37 !! Note that if "key_dynspg_fsc" defined, the time stepping38 !! has already been performed in dynspg.F routine37 !! Note that if lk_dynspg_fsc=T, the time stepping has already been 38 !! performed in dynspg module 39 39 !! Time filter applied on now horizontal velocity to avoid the 40 40 !! divergence of two consecutive time-steps and swap of dynamics -
trunk/NEMO/OPA_SRC/DYN/dynvor.F90
r3 r32 27 27 28 28 !! * Shared module variables 29 LOGICAL, PUBLIC :: ln_dynvor_ene = .FALSE. ! energy conserving scheme30 LOGICAL, PUBLIC :: ln_dynvor_ens = .TRUE. ! enstrophy conserving scheme31 LOGICAL, PUBLIC :: ln_dynvor_mix = .FALSE. ! mixed scheme29 LOGICAL, PUBLIC :: ln_dynvor_ene = .FALSE. !: energy conserving scheme 30 LOGICAL, PUBLIC :: ln_dynvor_ens = .TRUE. !: enstrophy conserving scheme 31 LOGICAL, PUBLIC :: ln_dynvor_mix = .FALSE. !: mixed scheme 32 32 33 33 !! * Substitutions -
trunk/NEMO/OPA_SRC/DYN/dynzdf_iso.F90
r3 r32 387 387 CONTAINS 388 388 SUBROUTINE dyn_zdf_iso( kt ) ! Dummy routine 389 WRITE(*,*) kt389 WRITE(*,*) 'dyn_zdf_iso: You should not have seen this print! error?', kt 390 390 END SUBROUTINE dyn_zdf_iso 391 391 #endif
Note: See TracChangeset
for help on using the changeset viewer.