- Timestamp:
- 12/17/13 23:02:01 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90
r41 r83 74 74 LOGICAL :: ln_2m = .FALSE. ! logical flag for height of air temp. and hum 75 75 LOGICAL :: ln_taudif = .FALSE. ! logical flag to use the "mean of stress module - module of mean stress" data 76 REAL(wp) :: rn_usecrt = 1. ! weighted use surface currents to compute relative wind speed 76 77 REAL(wp) :: rn_pfac = 1. ! multiplication factor for precipitation 77 78 … … 128 129 TYPE(FLD_N) :: sn_qlw , sn_tair, sn_prec, sn_snow ! " " 129 130 TYPE(FLD_N) :: sn_tdif ! " " 130 NAMELIST/namsbc_core/ cn_dir , ln_2m , ln_taudif, rn_pfac, 131 & sn_wndi, sn_wndj, sn_humi , sn_qsr , &131 NAMELIST/namsbc_core/ cn_dir , ln_2m , ln_taudif, rn_pfac, rn_usecrt, & 132 & sn_wndi, sn_wndj, sn_humi , sn_qsr , & 132 133 & sn_qlw , sn_tair, sn_prec , sn_snow, sn_tdif 133 134 !!--------------------------------------------------------------------- … … 263 264 264 265 ! ... components ( U10m - U_oce ) at T-point (unmasked) 265 zwnd_i(:,:) = 0.e0 266 zwnd_j(:,:) = 0.e0 266 IF( rn_usecrt /= 0. ) THEN 267 267 #if defined key_vectopt_loop 268 268 !CDIR COLLAPSE 269 269 #endif 270 DO jj = 2, jpjm1 271 DO ji = fs_2, fs_jpim1 ! vect. opt. 272 zwnd_i(ji,jj) = ( sf(jp_wndi)%fnow(ji,jj,1) - 0.5 * ( pu(ji-1,jj ) + pu(ji,jj) ) ) 273 zwnd_j(ji,jj) = ( sf(jp_wndj)%fnow(ji,jj,1) - 0.5 * ( pv(ji ,jj-1) + pv(ji,jj) ) ) 270 DO jj = 2, jpjm1 271 DO ji = fs_2, fs_jpim1 ! vect. opt. 272 zwnd_i(ji,jj) = ( sf(jp_wndi)%fnow(ji,jj,1) - 0.5 * rn_usecrt * ( pu(ji-1,jj ) + pu(ji,jj) ) ) 273 zwnd_j(ji,jj) = ( sf(jp_wndj)%fnow(ji,jj,1) - 0.5 * rn_usecrt * ( pv(ji ,jj-1) + pv(ji,jj) ) ) 274 END DO 274 275 END DO 275 END DO 276 CALL lbc_lnk( zwnd_i(:,:) , 'T', -1. ) 277 CALL lbc_lnk( zwnd_j(:,:) , 'T', -1. ) 278 ! ... scalar wind ( = | U10m - U_oce | ) at T-point (masked) 276 CALL lbc_lnk( zwnd_i(:,:) , 'T', -1. ) 277 CALL lbc_lnk( zwnd_j(:,:) , 'T', -1. ) 278 ELSE 279 zwnd_i(:,:) = sf(jp_wndi)%fnow(:,:,1) 280 zwnd_j(:,:) = sf(jp_wndj)%fnow(:,:,1) 281 END IF 282 ! ... scalar wind module at T-point (masked) 279 283 !CDIR NOVERRCHK 280 284 !CDIR COLLAPSE
Note: See TracChangeset
for help on using the changeset viewer.