Changeset 83 for trunk

Timestamp:

12/17/13 23:02:01 (10 years ago)

Author:

smasson

Message:

introduce rn_usecrt (1 by default)

File:

• trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90 (modified) (3 diffs)

trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90

@@ -74 +74 @@
    LOGICAL  ::   ln_2m     = .FALSE.   ! logical flag for height of air temp. and hum
    LOGICAL  ::   ln_taudif = .FALSE.   ! logical flag to use the "mean of stress module - module of mean stress" data
+   REAL(wp) ::   rn_usecrt = 1.        ! weighted use surface currents to compute relative wind speed
    REAL(wp) ::   rn_pfac   = 1.        ! multiplication factor for precipitation
 
@@ -128 +129 @@
       TYPE(FLD_N) ::   sn_qlw , sn_tair, sn_prec, sn_snow      !   "                                 "
       TYPE(FLD_N) ::   sn_tdif                                 !   "                                 "
-      NAMELIST/namsbc_core/ cn_dir , ln_2m  , ln_taudif, rn_pfac,           &
-         &                  sn_wndi, sn_wndj, sn_humi  , sn_qsr ,           &
+      NAMELIST/namsbc_core/ cn_dir , ln_2m  , ln_taudif, rn_pfac, rn_usecrt, &
+         &                  sn_wndi, sn_wndj, sn_humi  , sn_qsr ,            &
          &                  sn_qlw , sn_tair, sn_prec  , sn_snow, sn_tdif
       !!---------------------------------------------------------------------
@@ -263 +264 @@
 
       ! ... components ( U10m - U_oce ) at T-point (unmasked)
-      zwnd_i(:,:) = 0.e0  
-      zwnd_j(:,:) = 0.e0
+      IF( rn_usecrt /= 0. ) THEN
 #if defined key_vectopt_loop
 !CDIR COLLAPSE
 #endif
-      DO jj = 2, jpjm1
-         DO ji = fs_2, fs_jpim1   ! vect. opt.
-            zwnd_i(ji,jj) = (  sf(jp_wndi)%fnow(ji,jj,1) - 0.5 * ( pu(ji-1,jj  ) + pu(ji,jj) )  )
-            zwnd_j(ji,jj) = (  sf(jp_wndj)%fnow(ji,jj,1) - 0.5 * ( pv(ji  ,jj-1) + pv(ji,jj) )  )
+         DO jj = 2, jpjm1
+            DO ji = fs_2, fs_jpim1   ! vect. opt.
+               zwnd_i(ji,jj) = (  sf(jp_wndi)%fnow(ji,jj,1) - 0.5 * rn_usecrt * ( pu(ji-1,jj  ) + pu(ji,jj) )  )
+               zwnd_j(ji,jj) = (  sf(jp_wndj)%fnow(ji,jj,1) - 0.5 * rn_usecrt * ( pv(ji  ,jj-1) + pv(ji,jj) )  )
+            END DO
          END DO
-      END DO
-      CALL lbc_lnk( zwnd_i(:,:) , 'T', -1. )
-      CALL lbc_lnk( zwnd_j(:,:) , 'T', -1. )
-      ! ... scalar wind ( = | U10m - U_oce | ) at T-point (masked)
+         CALL lbc_lnk( zwnd_i(:,:) , 'T', -1. )
+         CALL lbc_lnk( zwnd_j(:,:) , 'T', -1. )
+      ELSE
+         zwnd_i(:,:) = sf(jp_wndi)%fnow(:,:,1)
+         zwnd_j(:,:) = sf(jp_wndj)%fnow(:,:,1)
+      END IF
+      ! ... scalar wind module at T-point (masked)
 !CDIR NOVERRCHK
 !CDIR COLLAPSE