Changeset 8918


Ignore:
Timestamp:
2017-12-06T13:27:23+01:00 (3 years ago)
Author:
jcastill
Message:

Fix the calculation of relative winds, as winds are received on the U,V grids and not on the T grid

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/AMM15_v3_6_STABLE_package/NEMOGCM/NEMO/OPA_SRC/SBC/sbcflx.F90

    r8416 r8918  
    157157 
    158158            IF( ln_rel_wind ) THEN 
    159                DO jj = 2, jpjm1 
    160                   DO ji = fs_2, fs_jpim1   ! vect. opt. 
    161                      zwnd_i(ji,jj) = ( sf(jp_utau)%fnow(ji,jj,1) - rn_wfac * 0.5 * ( ssu_m(ji-1,jj  ) + ssu_m(ji,jj) )) 
    162                      zwnd_j(ji,jj) = ( sf(jp_vtau)%fnow(ji,jj,1) - rn_wfac * 0.5 * ( ssv_m(ji  ,jj-1) + ssv_m(ji,jj) )) 
     159               DO jj = 1, jpj 
     160                  DO ji = 1, jpi 
     161                     zwnd_i(ji,jj) = sf(jp_utau)%fnow(ji,jj,1) - rn_wfac * ssu_m(ji,jj) 
     162                     zwnd_j(ji,jj) = sf(jp_vtau)%fnow(ji,jj,1) - rn_wfac * ssv_m(ji,jj) 
    163163                  END DO 
    164164               END DO 
    165                CALL lbc_lnk( zwnd_i(:,:) , 'T', -1. ) 
    166                CALL lbc_lnk( zwnd_j(:,:) , 'T', -1. ) 
    167165            ELSE 
    168166               zwnd_i(:,:) = sf(jp_utau)%fnow(:,:,1) 
Note: See TracChangeset for help on using the changeset viewer.