Changeset 7807 for branches/UKMO/r6232_HZG_WAVE/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90

Timestamp:

2017-03-17T10:44:05+01:00 (7 years ago)

Author:

jcastill

Message:

Changes as in HZG wave forcing branch, but adapted to r6232

File:

• branches/UKMO/r6232_HZG_WAVE/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90 (modified) (6 diffs)

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

@@ -85 +85 @@
       !!----------------------------------------------------------------------
       INTEGER ::   icpt   ! local integer
-      !!
+      !! V.Alari 26.06.2016 Added logicals for waves, deleted ln_wave ln_sdw, which came
+      !with UKMO
       NAMELIST/namsbc/ nn_fsbc   , ln_ana    , ln_flx, ln_blk_clio, ln_blk_core, ln_mixcpl,   &
          &             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_ssr    , nn_isf    , nn_fwb,                                        &
+         &             nn_lsm    , nn_limflx , nn_components, ln_cpl,ln_stcor, ln_wavetke,    &
+         &             ln_tauoc, ln_cdgw
       INTEGER  ::   ios
       INTEGER  ::   ierr, ierr0, ierr1, ierr2, ierr3, jpm
@@ -120 +122 @@
           nn_ice      =   0
       ENDIF
-
+     ! V.Alari added wave logicalc (26.06.2016)
       IF(lwp) THEN               ! Control print
          WRITE(numout,*) '        Namelist namsbc (partly overwritten with CPP key setting)'
@@ -146 +148 @@
          WRITE(numout,*) '              closed sea (=0/1) (set in namdom)          nn_closea   = ', nn_closea
          WRITE(numout,*) '              n. of iterations if land-sea-mask applied  nn_lsm      = ', nn_lsm
-      ENDIF
+         WRITE(numout,*) '         Stokes-Coriolis forcing from wavemodel ln_stcor       = ', ln_stcor
+         WRITE(numout,*) '              TKE wave breaking BC from wave model ln_wavetke        = ', ln_wavetke
+         WRITE(numout,*) '              Wave-modified stress from wave model ln_tauoc          = ', ln_tauoc
+         WRITE(numout,*) '              Drag coefficient from wave model ln_cdgw              = ', ln_cdgw
+ 
+
+     ENDIF
 
       ! LIM3 Multi-category heat flux formulation
@@ -179 +187 @@
 
       !                          ! Checks:
-      IF( nn_isf .EQ. 0 ) THEN                      ! variable initialisation if no ice shelf 
+      IF( nn_isf .EQ. 0 ) THEN                      ! no specific treatment in vicinity of ice shelf 
          IF( sbc_isf_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'sbc_init : unable to allocate sbc_isf arrays' )
          fwfisf  (:,:)   = 0.0_wp ; fwfisf_b  (:,:)   = 0.0_wp
@@ -214 +222 @@
       IF( ln_dm2dc .AND. .NOT.( ln_flx .OR. ln_blk_core ) .AND. nn_components /= jp_iam_opa )   &
          &   CALL ctl_stop( 'diurnal cycle into qsr field from daily values requires a flux or core-bulk formulation' )
-      
-      IF ( ln_wave ) THEN
+   
+     ! V.Alari comment out UKMO wave stuff   
+    !  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)' )
+     !    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 
-         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')
-         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) ')
-      ENDIF 
+       !  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')
+       !  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) ')
+     ! ENDIF 
+      ! End V.Alari
       !                          ! Choice of the Surface Boudary Condition (set nsbc)
       ll_purecpl = ln_cpl .AND. .NOT. ln_mixcpl
@@ -352 +362 @@
       IF( nn_components /= jp_iam_sas )   CALL sbc_ssm( kt )   ! ocean sea surface variables (sst_m, sss_m, ssu_m, ssv_m)
       !                                                        ! averaged over nf_sbc time-step
-
-      IF (ln_wave) CALL sbc_wave( kt )
+            ! START V.Alari 26.06.2016
+     IF (ln_cdgw) CALL sbc_wave( kt )
+      ! Read wave parameters if Stokes-Coriolis forcing OR wave parameters for
+      ! TKE are needed
+      ! Also read it if drag coefficient to ensure that we use the neutral 10m
+      ! from WAM.
+      IF (ln_stcor .OR. ln_wavetke .OR. ln_tauoc) THEN
+            CALL sbc_wavepar( kt )
+      ENDIF
+
+       ! END
+
+!      IF (ln_wave) CALL sbc_wave( kt )
                                                    !==  sbc formulation  ==!