Changeset 10473 for branches/UKMO/AMM15_v3_6_STABLE_package_collate_coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90

Timestamp:

2019-01-08T18:02:36+01:00 (5 years ago)

Author:

jcastill

Message:

Merged branch UKMO/r6232_INGV1_WAVE-coupling@7620

File:

• branches/UKMO/AMM15_v3_6_STABLE_package_collate_coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90 (modified) (5 diffs)

branches/UKMO/AMM15_v3_6_STABLE_package_collate_coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90

@@ -89 +89 @@
          &             ln_blk_mfs, ln_apr_dyn, nn_ice, nn_ice_embd, ln_dm2dc   , ln_rnf   ,   &
          &             ln_ssr    , nn_isf    , nn_fwb, ln_cdgw    , ln_wave    , ln_sdw   ,   &
-         &             nn_lsm    , nn_limflx , nn_components, ln_cpl
+         &             ln_tauoc  , ln_stcor  , nn_lsm, nn_limflx , nn_components, ln_cpl
       INTEGER  ::   ios
       INTEGER  ::   ierr, ierr0, ierr1, ierr2, ierr3, jpm
@@ -132 +132 @@
          WRITE(numout,*) '              ocean-atmosphere coupled formulation       ln_cpl      = ', ln_cpl
          WRITE(numout,*) '              forced-coupled mixed formulation           ln_mixcpl   = ', ln_mixcpl
+         WRITE(numout,*) '              wave physics                               ln_wave     = ', ln_wave  
+         WRITE(numout,*) '                 Stokes drift corr. to vert. velocity    ln_sdw      = ', ln_sdw  
+         WRITE(numout,*) '                 wave modified ocean stress              ln_tauoc    = ', ln_tauoc  
+         WRITE(numout,*) '                 Stokes coriolis term                    ln_stcor    = ', ln_stcor  
+         WRITE(numout,*) '                 neutral drag coefficient (CORE, MFS)    ln_cdgw     = ', ln_cdgw
          WRITE(numout,*) '              OASIS coupling (with atm or sas)           lk_oasis    = ', lk_oasis
          WRITE(numout,*) '              components of your executable              nn_components = ', nn_components
@@ -216 +221 @@
       
       IF ( ln_wave ) THEN
-      !Activated wave module but neither drag nor stokes drift activated
-         IF ( .NOT.(ln_cdgw .OR. ln_sdw) )   THEN
-            CALL ctl_warn( 'Ask for wave coupling but nor drag coefficient (ln_cdgw=F) neither stokes drift activated (ln_sdw=F)' )
-      !drag coefficient read from wave model definable only with mfs bulk formulae and core 
+         !Activated wave module but neither drag nor stokes drift activated  
+         IF ( .NOT.(ln_cdgw .OR. ln_sdw .OR. ln_tauoc .OR. ln_stcor ) )   THEN   
+             CALL ctl_warn( 'Ask for wave coupling but ln_cdgw=F, ln_sdw=F, ln_tauoc=F, ln_stcor=F')  
+         !drag coefficient read from wave model definable only with mfs bulk formulae and core
          ELSEIF (ln_cdgw .AND. .NOT.(ln_blk_mfs .OR. ln_blk_core) )       THEN       
              CALL ctl_stop( 'drag coefficient read from wave model definable only with mfs bulk formulae and core')
+         ELSEIF (ln_stcor .AND. .NOT. ln_sdw) THEN   
+             CALL ctl_stop( 'Stokes-Coriolis term calculated only if activated Stokes Drift ln_sdw=T')
          ENDIF
       ELSE
-      IF ( ln_cdgw .OR. ln_sdw  )                                         & 
-         &   CALL ctl_stop('Not Activated Wave Module (ln_wave=F) but     &
-         & asked coupling with drag coefficient (ln_cdgw =T) or Stokes drift (ln_sdw=T) ')
+         IF ( ln_cdgw .OR. ln_sdw .OR. ln_tauoc .OR. ln_stcor ) &   
+            &   CALL ctl_stop( 'Not Activated Wave Module (ln_wave=F) but asked coupling ',    &  
+            &                  'with drag coefficient (ln_cdgw =T) '  ,                        &  
+            &                  'or Stokes Drift (ln_sdw=T) ' ,                                 &  
+            &                  'or ocean stress modification due to waves (ln_tauoc=T) ',      &  
+            &                  'or Stokes-Coriolis term (ln_stcori=T)'  )
       ENDIF 
       !                          ! Choice of the Surface Boudary Condition (set nsbc)
@@ -307 +317 @@
 
       IF( nn_ice == 4      )   CALL cice_sbc_init( nsbc )      ! CICE initialisation
+      !  
+      IF( ln_wave          )   CALL sbc_wave_init              ! surface wave initialisation  
+      !
       
    END SUBROUTINE sbc_init
@@ -386 +399 @@
       IF( ln_mixcpl )        CALL sbc_cpl_rcv ( kt, nn_fsbc, nn_ice )   ! forced-coupled mixed formulation after forcing
 
+      IF ( ln_wave .AND. ln_tauoc) THEN                 ! Wave stress subctracted  
+            utau(:,:) = utau(:,:)*tauoc_wave(:,:)  
+            vtau(:,:) = vtau(:,:)*tauoc_wave(:,:)  
+            taum(:,:) = taum(:,:)*tauoc_wave(:,:)  
+      !  
+            SELECT CASE( nsbc )  
+            CASE(  0,1,2,3,5,-1 )  ;  
+                IF(lwp .AND. kt == nit000 ) WRITE(numout,*) 'WARNING: You are subtracting the wave stress to the ocean. &  
+                        & If not requested select ln_tauoc=.false'  
+            END SELECT  
+      !  
+      END IF
 
       !                                            !==  Misc. Options  ==!