Changeset 8553

Timestamp:

2017-09-21T13:09:26+02:00 (7 years ago)

Author:

jcastill

Message:

Cap the value of tauoc as recieved via coupling or via forcing

Location:

branches/UKMO/r6232_HZG_WAVE-coupling/NEMOGCM/NEMO/OPA_SRC/SBC

Files:

• sbccpl.F90 (modified) (1 diff)
• sbcwave.F90 (modified) (1 diff)

branches/UKMO/r6232_HZG_WAVE-coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

@@ -1200 +1200 @@
       !                                                      ! Stress adsorbed by waves  ! 
       !                                                      ! ========================= !  
-      IF( srcv(jpr_tauoc)%laction .AND. ln_tauoc ) tauoc_wave(:,:) = frcv(jpr_tauoc)%z3(:,:,1) 
+      IF( srcv(jpr_tauoc)%laction .AND. ln_tauoc ) THEN
+         tauoc_wave(:,:) = frcv(jpr_tauoc)%z3(:,:,1)
+         ! cap the value of tauoc
+         WHERE(tauoc_wave <   0.0 ) tauoc_wave = 1.0
+         WHERE(tauoc_wave > 100.0 ) tauoc_wave = 1.0
+      ENDIF
       
       !                                                      ! ========================= !  

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

@@ -261 +261 @@
          CALL fld_read( kt, nn_fsbc, sf_tauoc )          ! read wave norm stress from external forcing
          tauoc_wave(:,:) = sf_tauoc(1)%fnow(:,:,1)
-         WHERE( tauoc_wave < -100.0 ) tauoc_wave = 1.0
-         WHERE( tauoc_wave >  100.0 ) tauoc_wave = 1.0
+         WHERE( tauoc_wave <   0.0 ) tauoc_wave = 1.0
+         WHERE( tauoc_wave > 100.0 ) tauoc_wave = 1.0
       ENDIF