Changeset 9966


Ignore:
Timestamp:
2018-07-18T16:26:36+02:00 (2 years ago)
Author:
lovato
Message:

trunk : correct peak wave number computation using wave model frequency data (Stokes drift case nn_sdrift=2)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/SBC/sbcwave.F90

    r9821 r9966  
    131131         END DO 
    132132      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 
    133138         DO jj = 1, jpjm1 
    134139            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) ) / grav 
    136                zk_v(ji,jj) = 0.5_wp * ( wfreq(ji,jj)*wfreq(ji,jj) + wfreq(ji,jj+1)*wfreq(ji,jj+1) ) / grav 
     140               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) ) 
    137142               ! 
    138143               zu0_sd(ji,jj) = 0.5_wp * ( ut0sd(ji,jj) + ut0sd(ji+1,jj) ) 
Note: See TracChangeset for help on using the changeset viewer.