Changeset 1705 for trunk/NEMO/OPA_SRC/SBC/sbcblk_core.F90

Timestamp:

2009-11-03T17:37:00+01:00 (14 years ago)

Author:

smasson

Message:

impact of HF winds in TKE, see ticket:585

File:

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

trunk/NEMO/OPA_SRC/SBC/sbcblk_core.F90

@@ -43 +43 @@
    PUBLIC   blk_ice_core       ! routine called in sbc_ice_lim module
       
-   INTEGER , PARAMETER ::   jpfld   = 8           ! maximum number of files to read 
+   INTEGER , PARAMETER ::   jpfld   = 9           ! maximum number of files to read 
    INTEGER , PARAMETER ::   jp_wndi = 1           ! index of 10m wind velocity (i-component) (m/s)    at T-point
    INTEGER , PARAMETER ::   jp_wndj = 2           ! index of 10m wind velocity (j-component) (m/s)    at T-point
@@ -52 +52 @@
    INTEGER , PARAMETER ::   jp_prec = 7           ! index of total precipitation (rain+snow) (Kg/m2/s)
    INTEGER , PARAMETER ::   jp_snow = 8           ! index of snow (solid prcipitation)       (kg/m2/s)
+   INTEGER , PARAMETER ::   jp_tdif = 9           ! index of tau diff associated to HF tau   (N/m2)   at T-point
    TYPE(FLD), ALLOCATABLE, DIMENSION(:) ::   sf   ! structure of input fields (file informations, fields read)
          
@@ -63 +64 @@
 
    !                                !!* Namelist namsbc_core : CORE bulk parameters
-   LOGICAL  ::   ln_2m = .FALSE.     ! logical flag for height of air temp. and hum
-   REAL(wp) ::   rn_pfac = 1.        ! multiplication factor for precipitation
+   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_pfac   = 1.          ! multiplication factor for precipitation
 
    !! * Substitutions
@@ -93 +95 @@
       !!      the total precipitation (rain+snow) (Kg/m2/s)
       !!      the snow (solid prcipitation)       (kg/m2/s)
+      !!   OPTIONAL parameter (see ln_taudif namelist flag):
+      !!      the tau diff associated to HF tau   (N/m2)   at T-point 
       !!              (2) CALL blk_oce_core
       !!
@@ -110 +114 @@
       INTEGER  ::   ierror   ! return error code
       INTEGER  ::   ifpr     ! dummy loop indice
+      INTEGER  ::   jfld     ! dummy loop arguments
       !!
       CHARACTER(len=100) ::  cn_dir   !   Root directory for location of core files
@@ -115 +120 @@
       TYPE(FLD_N) ::   sn_wndi, sn_wndj, sn_humi, sn_qsr       ! informations about the fields to be read
       TYPE(FLD_N) ::   sn_qlw , sn_tair, sn_prec, sn_snow      !   "                                 "
-      NAMELIST/namsbc_core/ cn_dir, ln_2m, rn_pfac, sn_wndi, sn_wndj, sn_humi, sn_qsr,   &
-         &                                                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 ,           &
+         &                  sn_qlw , sn_tair, sn_prec  , sn_snow, sn_tdif
       !!---------------------------------------------------------------------
 
@@ -136 +143 @@
          sn_prec = FLD_N( 'precip'  ,    -1.    ,  'precip'  ,  .true.    , .false. ,   'yearly'  , ''       , ''         )
          sn_snow = FLD_N( 'snow'    ,    -1.    ,  'snow'    ,  .true.    , .false. ,   'yearly'  , ''       , ''         )
+         sn_tdif = FLD_N( 'taudif'  ,    24.    ,  'taudif'  ,  .true.    , .false. ,   'yearly'  , ''       , ''         )
          !
          REWIND( numnam )                    ! ... read in namlist namsbc_core
@@ -145 +153 @@
          slf_i(jp_tair) = sn_tair   ;   slf_i(jp_humi) = sn_humi
          slf_i(jp_prec) = sn_prec   ;   slf_i(jp_snow) = sn_snow
+         slf_i(jp_tdif) = sn_tdif
+         !
+         ! do we use HF tau information?
+         lhftau = ln_taudif
+         jfld = jpfld - COUNT( (/.NOT. lhftau/) )
          !
          ! set sf structure
-         ALLOCATE( sf(jpfld), STAT=ierror )
+         ALLOCATE( sf(jfld), STAT=ierror )
          IF( ierror > 0 ) THEN
             CALL ctl_stop( 'sbc_blk_core: unable to allocate sf structure' )   ;   RETURN
          ENDIF
-         DO ifpr= 1, jpfld
+         DO ifpr= 1, jfld
             ALLOCATE( sf(ifpr)%fnow(jpi,jpj) )
             ALLOCATE( sf(ifpr)%fdta(jpi,jpj,2) )
@@ -291 +304 @@
          END DO
       END DO
+
+      ! ... add the HF tau contribution to the wind stress module?
+      IF( lhftau ) THEN 
+!CDIR COLLAPSE
+         taum(:,:) = taum(:,:) + sf(jp_tdif)%fnow(:,:)
+      ENDIF
+      CALL iom_put( "taum_oce", taum )   ! output wind stress module
+
       ! ... utau, vtau at U- and V_points, resp.
       !     Note the use of 0.5*(2-umask) in order to unmask the stress along coastlines