Changeset 7717 for branches/UKMO/r6232_INGV1_WAVE-coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcwave.F90

Timestamp:

2017-02-22T17:49:40+01:00 (7 years ago)

Author:

jcastill

Message:

Add surface drag coefficient direct forcing and coupling

File:

• branches/UKMO/r6232_INGV1_WAVE-coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcwave.F90 (modified) (4 diffs)

branches/UKMO/r6232_INGV1_WAVE-coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcwave.F90

@@ -63 +63 @@
    REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:)   ::   ut0sd, vt0sd        !: surface Stokes drift velocities at t-point
    REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:,:) ::   usd  , vsd  , wsd   !: Stokes drift velocities at u-, v- & w-points, resp.
+
+   INTEGER, PUBLIC            ::   nn_drag       ! type of formula to calculate wind stress from wind components
+   INTEGER, PUBLIC, PARAMETER ::   jp_ukmo = 0        ! UKMO SHELF formulation
+   INTEGER, PUBLIC, PARAMETER ::   jp_std  = 1        ! standard formulation: 
 
    !! * Substitutions
@@ -210 +214 @@
          CALL fld_read( kt, nn_fsbc, sf_cd )             ! read from external forcing
          cdn_wave(:,:) = sf_cd(1)%fnow(:,:,1)
+         ! make sure that the interpolation does not write silly fill values in case it fails in a grid point
+         where( cdn_wave > 1.0 ) cdn_wave = 1.5e-3 
       ENDIF
 
@@ -265 +271 @@
                              &   sn_hsw, sn_wmp, sn_wnum, sn_tauoc      ! informations about the fields to be read
       !
-      NAMELIST/namsbc_wave/  sn_cdg, cn_dir, sn_usd, sn_vsd, sn_hsw, sn_wmp, sn_wnum, sn_tauoc
+      NAMELIST/namsbc_wave/  sn_cdg, cn_dir, sn_usd, sn_vsd, sn_hsw, sn_wmp, sn_wnum, sn_tauoc, nn_drag
       !!---------------------------------------------------------------------
       !
@@ -278 +284 @@
       !
       IF( ln_cdgw ) THEN
+         ! check the value of nn_drag
+         IF( .NOT. ( nn_drag == jp_ukmo .OR. nn_drag == jp_std )) THEN
+            CALL ctl_stop( 'sbc_wave_init: nn_drag does not specify a valid drag formulation' )
+            RETURN  
+         ENDIF
+           
          IF( .NOT. cpl_wdrag ) THEN
             ALLOCATE( sf_cd(1), STAT=ierror )           !* allocate and fill sf_wave with sn_cdg