Changeset 7346


Ignore:
Timestamp:
2016-11-28T11:26:34+01:00 (4 years ago)
Author:
emanuelaclementi
Message:

#1643 removal of surface wn modification due to Stokes Drift in sshsshwzv.F90 and minor modifications in development branch 2015/dev_r5936_INGV1_WAVE

Location:
branches/2015/dev_r5936_INGV1_WAVE/NEMOGCM/NEMO/OPA_SRC
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5936_INGV1_WAVE/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90

    r7340 r7346  
    99   !!             -   !  2010-09  (D.Storkey and E.O'Dea) bug fixes for BDY module 
    1010   !!            3.3  !  2011-10  (M. Leclair) split former ssh_wzv routine and remove all vvl related work 
    11    !!            3.6  !  2014-10  (E. Clementi, P. Oddo) add wave contribution to surface vertical velocity  
    1211   !!---------------------------------------------------------------------- 
    1312 
     
    3938   USE wrk_nemo        ! Memory Allocation 
    4039   USE timing          ! Timing 
    41    USE sbcwave         ! Stokes velocities 
    4240 
    4341   IMPLICIT NONE 
     
    162160      REAL(wp), POINTER, DIMENSION(:,:  ) ::  z2d 
    163161      REAL(wp), POINTER, DIMENSION(:,:,:) ::  z3d, zhdiv 
    164       REAL(wp) ::  zdsshu, zdsshv 
    165162      !!---------------------------------------------------------------------- 
    166163      ! 
     
    210207      ENDIF 
    211208 
    212       IF( ln_wave .AND. ln_sdw ) THEN 
    213          ! Compute d(ssh)/dx  and d(ssh)/dy   
    214          ! Compute the surface vertical velocity accounting for the Stokes Drift 
    215          DO jj = 1 , jpjm1 
    216             DO ji = 1 , fs_jpim1 
    217               zdsshu = ( sshn(ji+1,jj) - sshn(ji,jj) ) / e1u(ji,jj) 
    218               zdsshv = ( sshn(ji,jj+1) - sshn(ji,jj) ) / e2v(ji,jj) 
    219               wn(ji,jj,1) = wn(ji,jj,1) + ( usd3dt(ji,jj,1) * zdsshu   & 
    220                  &                        + vsd3dt(ji,jj,1) * zdsshv   & 
    221                  &        - wsd3d (ji,jj,1) ) * tmask(ji,jj,1) 
    222             END DO 
    223          END DO 
    224       ENDIF 
    225       ! 
    226209#if defined key_bdy 
    227210      IF( lk_bdy ) THEN 
  • branches/2015/dev_r5936_INGV1_WAVE/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90

    r7343 r7346  
    377377            SELECT CASE( nsbc ) 
    378378            CASE(  0,1,2,3,5,-1 )  ; 
    379                 IF(lwp .AND. kt == nit000 ) WRITE(numout,*) 'WARNING: You are subtracting the wave stress to the ocean. & 
     379                IF(lwp) WRITE(numout,*) 'WARNING: You are subtracting the wave stress to the ocean. & 
    380380                        & If not requested select ln_tauoc=.false' 
    381381            END SELECT 
  • branches/2015/dev_r5936_INGV1_WAVE/NEMOGCM/NEMO/OPA_SRC/SBC/sbcwave.F90

    r7343 r7346  
    77   !!         :  3.4  !   2012-10  (M. Adani)  Stokes Drift  
    88   !!            3.6  !   2014-09  (E. Clementi,P. Oddo) New Stokes Drift Computation 
    9    !!            3.6  !   2014-09  (Clementi E, Oddo P)New Stokes Drift Computation 
    109   !!---------------------------------------------------------------------- 
    1110 
     
    5453   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:)       :: tauoc_wave 
    5554   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:)       :: tsd2d 
    56    REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:)       :: usd2d, vsd2d 
    5755   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:)       :: zusd2dt, zvsd2dt 
    5856   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:,:)     :: usd3d, vsd3d, wsd3d  
     
    281279               CALL fld_fill( sf_sd, slf_i, cn_dir, 'sbc_wave', 'Wave module ', 'namsbc_wave' ) 
    282280            ENDIF 
    283             ALLOCATE( usd2d(jpi,jpj),vsd2d(jpi,jpj) ) 
    284281            ALLOCATE( usd3d(jpi,jpj,jpk),vsd3d(jpi,jpj,jpk),wsd3d(jpi,jpj,jpk) ) 
    285282            ALLOCATE( usd3dt(jpi,jpj,jpk),vsd3dt(jpi,jpj,jpk) ) 
  • branches/2015/dev_r5936_INGV1_WAVE/NEMOGCM/NEMO/OPA_SRC/step.F90

    r7343 r7346  
    220220                         CALL dyn_adv       ( kstp )  ! advection (vector or flux form) 
    221221                         CALL dyn_vor       ( kstp )  ! vorticity term including Coriolis 
    222       IF( ln_wave .AND. ln_sdw .AND. ln_stcor    )  & 
     222      IF( ln_wave .AND. ln_sdw .AND. ln_stcor)           & 
    223223               &         CALL dyn_stcor     ( kstp )  ! Stokes-Coriolis forcing 
    224224                         CALL dyn_ldf       ( kstp )  ! lateral mixing 
    225       IF( ln_wave .AND. ln_sdw .AND. ln_stcor)           &  
    226                &         CALL dyn_stcor     ( kstp )  ! Stokes-Coriolis forcing 
    227225                         CALL dyn_hpg       ( kstp )  ! horizontal gradient of Hydrostatic pressure 
    228226                         CALL dyn_spg       ( kstp )  ! surface pressure gradient 
Note: See TracChangeset for help on using the changeset viewer.