Changeset 32 for trunk/NEMO/OPA_SRC/LDF/ldfslp.F90
- Timestamp:
- 2004-02-17T10:20:15+01:00 (20 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/LDF/ldfslp.F90
r3 r32 29 29 30 30 !! * Share module variables 31 LOGICAL , PUBLIC, PARAMETER :: lk_ldfslp = .TRUE. ! slopes flag32 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) :: & 33 uslp, wslpi, & ! i_slope at U- and W-points34 vslp, wslpj ! j-slope at V- and W-points31 LOGICAL , PUBLIC, PARAMETER :: lk_ldfslp = .TRUE. !: slopes flag 32 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) :: & !: 33 uslp, wslpi, & !: i_slope at U- and W-points 34 vslp, wslpj !: j-slope at V- and W-points 35 35 36 36 !! * Module variables … … 95 95 96 96 !! * Local declarations 97 INTEGER :: ji, jj, jk ! dummy loop indices 97 INTEGER :: ji, jj, jk ! dummy loop indices 98 INTEGER :: ii0, ii1, ij0, ij1 ! temporary integer 98 99 #if defined key_partial_steps 99 100 INTEGER :: iku, ikv ! temporary integers … … 120 121 121 122 zeps = 1.e-20 122 zmg = -1.0 / g 123 zm05g = -0.5 / g 123 zmg = -1.0 / grav 124 zm05g = -0.5 / grav 124 125 125 126 zww(:,:,:) = 0.e0 … … 177 178 178 179 ! Local vertical density gradient evaluated from N^2 179 ! zwy = d/dz(prd)= - ( prd ) / g * mk(pn2) -- at t point180 ! zwy = d/dz(prd)= - ( prd ) / grav * mk(pn2) -- at t point 180 181 181 182 DO jj = 1, jpj … … 323 324 324 325 ! Local vertical density gradient evaluated from N^2 325 ! zwy = d/dz(prd)= - mk ( prd ) / g * pn2 -- at w point326 ! zwy = d/dz(prd)= - mk ( prd ) / grav * pn2 -- at w point 326 327 DO jj = 1, jpj 327 328 DO ji = 1, jpi … … 450 451 ! 451 452 ! ! Gibraltar Strait 452 uslp ( mi0(69):mi1(71) , mj0(50):mj1(53) , jk ) = 0.0e0 453 vslp ( mi0(68):mi1(71) , mj0(51):mj1(53) , jk ) = 0.0e0 454 wslpi( mi0(69):mi1(71) , mj0(51):mj1(53) , jk ) = 0.0e0 455 wslpj( mi0(69):mi1(71) , mj0(51):mj1(53) , jk ) = 0.0e0 453 ij0 = 50 ; ij1 = 53 454 ii0 = 69 ; ii1 = 71 ; uslp ( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , jk ) = 0.e0 455 ij0 = 51 ; ij1 = 53 456 ii0 = 68 ; ii1 = 71 ; vslp ( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , jk ) = 0.e0 457 ii0 = 69 ; ii1 = 71 ; wslpi( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , jk ) = 0.e0 458 ii0 = 69 ; ii1 = 71 ; wslpj( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , jk ) = 0.e0 456 459 457 460 ! ! Mediterrannean Sea 458 uslp ( mi0(71):mi1(90) , mj0(49):mj1(56) , jk ) = 0.0e0 459 vslp ( mi0(70):mi1(90) , mj0(50):mj1(56) , jk ) = 0.0e0 460 wslpi( mi0(71):mi1(90) , mj0(50):mj1(56) , jk ) = 0.0e0 461 wslpj( mi0(71):mi1(90) , mj0(50):mj1(56) , jk ) = 0.0e0 461 ij0 = 49 ; ij1 = 56 462 ii0 = 71 ; ii1 = 90 ; uslp ( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , jk ) = 0.e0 463 ij0 = 50 ; ij1 = 56 464 ii0 = 70 ; ii1 = 90 ; vslp ( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , jk ) = 0.e0 465 ii0 = 71 ; ii1 = 90 ; wslpi( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , jk ) = 0.e0 466 ii0 = 71 ; ii1 = 90 ; wslpj( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , jk ) = 0.e0 462 467 ENDIF 463 468 ! ! =============== … … 527 532 528 533 zeps = 1.e-20 529 zmg = -1.0 / g 530 zm05g = -0.5 / g 534 zmg = -1.0 / grav 535 zm05g = -0.5 / grav 531 536 532 537 … … 567 572 568 573 ! Local vertical density gradient evaluated from N^2 569 ! zwy = d/dz(prd)= - ( prd ) / g * mk(pn2) -- at t point574 ! zwy = d/dz(prd)= - ( prd ) / grav * mk(pn2) -- at t point 570 575 571 576 !----------------------------------------------------------------------- … … 619 624 620 625 ! Local vertical density gradient evaluated from N^2 621 ! zwy = d/dz(prd)= - ( prd ) / g * mk(pn2) -- at t point626 ! zwy = d/dz(prd)= - ( prd ) / grav * mk(pn2) -- at t point 622 627 zwy ( :, jpj) = 0.0e0 623 628 zwy ( jpi, :) = 0.0e0 … … 673 678 674 679 ! Local vertical density gradient evaluated from N^2 675 ! zwy = d/dz(prd)= - mk ( prd ) / g * pn2 -- at w point680 ! zwy = d/dz(prd)= - mk ( prd ) / grav * pn2 -- at w point 676 681 # if defined key_vectopt_loop && ! defined key_autotasking 677 682 jj = 1 … … 809 814 !! Dummy module : NO Rotation of lateral mixing tensor 810 815 !!------------------------------------------------------------------------ 811 LOGICAL, PUBLIC, PARAMETER :: lk_ldfslp = .FALSE. ! slopes flag816 LOGICAL, PUBLIC, PARAMETER :: lk_ldfslp = .FALSE. !: slopes flag 812 817 CONTAINS 813 818 SUBROUTINE ldf_slp( kt, prd, pn2 ) ! Dummy routine 814 819 INTEGER, INTENT(in) :: kt 815 820 REAL,DIMENSION(:,:,:), INTENT(in) :: prd, pn2 816 WRITE(*,*) kt, prd, pn2821 WRITE(*,*) 'ldf_slp: You should not have seen this print! error?', kt, prd(1,1,1), pn2(1,1,1) 817 822 END SUBROUTINE ldf_slp 818 823 #endif
Note: See TracChangeset
for help on using the changeset viewer.