- Timestamp:
- 2018-11-07T18:25:49+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r9866_HPC_03_globcom/src/OCE/SBC/sbcwave.F90
r9821 r10288 74 74 # include "vectopt_loop_substitute.h90" 75 75 !!---------------------------------------------------------------------- 76 !! NEMO/OCE 4.0 , NEMO Consortium (2018) 76 !! NEMO/OCE 4.0 , NEMO Consortium (2018) 77 77 !! $Id$ 78 !! Software governed by the CeCILL licen ce (./LICENSE)78 !! Software governed by the CeCILL license (see ./LICENSE) 79 79 !!---------------------------------------------------------------------- 80 80 CONTAINS … … 131 131 END DO 132 132 ELSE IF( ll_st_peakfr ) THEN ! peak wave number calculated from the peak frequency received by the wave model 133 DO jj = 1, jpj 134 DO ji = 1, jpi 135 zk_t(ji,jj) = ( 2.0_wp * rpi * wfreq(ji,jj) ) * ( 2.0_wp * rpi * wfreq(ji,jj) ) / grav 136 END DO 137 END DO 133 138 DO jj = 1, jpjm1 134 139 DO ji = 1, jpim1 135 zk_u(ji,jj) = 0.5_wp * ( wfreq(ji,jj)*wfreq(ji,jj) + wfreq(ji+1,jj)*wfreq(ji+1,jj) ) / grav136 zk_v(ji,jj) = 0.5_wp * ( wfreq(ji,jj)*wfreq(ji,jj) + wfreq(ji,jj+1)*wfreq(ji,jj+1) ) / grav140 zk_u(ji,jj) = 0.5_wp * ( zk_t(ji,jj) + zk_t(ji+1,jj) ) 141 zk_v(ji,jj) = 0.5_wp * ( zk_t(ji,jj) + zk_t(ji,jj+1) ) 137 142 ! 138 143 zu0_sd(ji,jj) = 0.5_wp * ( ut0sd(ji,jj) + ut0sd(ji+1,jj) )
Note: See TracChangeset
for help on using the changeset viewer.