Changeset 1705 for trunk/NEMO/OPA_SRC/SBC

Timestamp:

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

Author:

smasson

Message:

impact of HF winds in TKE, see ticket:585

Location:

trunk/NEMO/OPA_SRC/SBC

Files:

• sbc_oce.F90 (modified) (2 diffs)
• sbcblk_core.F90 (modified) (9 diffs)
• sbccpl.F90 (modified) (2 diffs)
• sbcmod.F90 (modified) (1 diff)

trunk/NEMO/OPA_SRC/SBC/sbc_oce.F90

@@ -29 +29 @@
    !                                             !:  = 1 global mean of e-p-r set to zero at each nn_fsbc time step
    !                                             !:  = 2 annual global mean of e-p-r set to zero
-  INTEGER , PUBLIC ::   nn_ico_cpl  = 0          !: ice-ocean coupling indicator
+   INTEGER , PUBLIC ::   nn_ico_cpl  = 0          !: ice-ocean coupling indicator
    !                                             !:  = 0   LIM-3 old case
    !                                             !:  = 1   stresses computed using now ocean velocity
@@ -37 +37 @@
    !!              Ocean Surface Boundary Condition fields
    !!----------------------------------------------------------------------
+   LOGICAL , PUBLIC ::   lhftau = .FALSE.              !: HF tau contribution: mean of stress module - module of the mean stress
    REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   utau      !: sea surface i-stress (ocean referential)     [N/m2]
    REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   vtau      !: sea surface j-stress (ocean referential)     [N/m2]
    REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   taum      !: module of sea surface stress (at T-point)    [N/m2] 
-   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   wndm      !: wind speed module at T-point (=|U10m-Uoce|)  [m/s] Used only in PISCES
+   !! wndm is used only in PISCES to compute gases exchanges at the surface of the free ocean or in the leads in sea-ice parts
+   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   wndm      !: wind speed module at T-point (=|U10m-Uoce|)  [m/s] 
    REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   qsr       !: sea heat flux:     solar                     [W/m2]
    REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   qns       !: sea heat flux: non solar                     [W/m2]

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

trunk/NEMO/OPA_SRC/SBC/sbccpl.F90

@@ -419 +419 @@
       !                                                      ! ------------------------- !
       srcv(jpr_taum)%clname = 'O_TauMod'   ;   IF( TRIM(cn_rcv_taumod) == 'coupled' )   srcv(jpr_taum)%laction = .TRUE.
+      lhftau = srcv(jpr_taum)%laction
 
 #if defined key_cpl_carbon_cycle
@@ -693 +694 @@
          taum(:,:) = frcv(:,:,jpr_taum)
          wndm(:,:) = frcv(:,:,jpr_w10m)
+         CALL iom_put( "taum_oce", taum )   ! output wind stress module
          !  
       ENDIF

trunk/NEMO/OPA_SRC/SBC/sbcmod.F90

@@ -253 +253 @@
       CALL iom_put( "utau", utau )   ! i-wind stress   (stress can be updated at 
       CALL iom_put( "vtau", vtau )   ! j-wind stress    each time step in sea-ice)
-      CALL iom_put( "wspd", wndm )   ! wind speed module 
+      CALL iom_put( "taum", taum )   ! wind stress module 
+      CALL iom_put( "wspd", wndm )   ! wind speed  module 
       !
       IF(ln_ctl) THEN         ! print mean trends (used for debugging)