Changeset 10478 for branches/UKMO/AMM15_v3_6_STABLE_package_collate_coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcflx.F90

Timestamp:

2019-01-09T12:11:55+01:00 (5 years ago)

Author:

jcastill

Message:

Merged branch UKMO/r6232_HZG_WAVE-coupling@9868

File:

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

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

@@ -20 +20 @@
    USE iom             ! IOM library
    USE in_out_manager  ! I/O manager
+   USE sbcwave         ! wave physics
    USE lib_mpp         ! distribued memory computing library
    USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
@@ -87 +88 @@
       REAL(wp) ::   zrhoa  = 1.22         ! Air density kg/m3
       REAL(wp) ::   zcdrag = 1.5e-3       ! drag coefficient
+      REAL(wp) ::   totwind               ! UKMO SHELF: Module of wind speed
       REAL(wp) ::   ztx, zty, zmod, zcoef ! temporary variables
       REAL     ::   cs                    ! UKMO SHELF: Friction co-efficient at surface
@@ -115 +117 @@
          IF(lwm) WRITE ( numond, namsbc_flx ) 
          !
+         IF(lwp) THEN                        ! Namelist print 
+            WRITE(numout,*) 
+            WRITE(numout,*) 'sbc_flx : Flux forcing' 
+            WRITE(numout,*) '~~~~~~~~~~~' 
+            WRITE(numout,*) '       Namelist namsbc_flx : shelf seas configuration (force with winds instead of momentum)' 
+            WRITE(numout,*) '          shelf seas configuration    ln_shelf_flx    = ', ln_shelf_flx 
+            WRITE(numout,*) '          relative wind speed         ln_rel_wind     = ', ln_rel_wind 
+            WRITE(numout,*) '          wind multiplication factor  rn_wfac         = ', rn_wfac 
+         ENDIF
          !                                         ! check: do we plan to use ln_dm2dc with non-daily forcing?
          IF( ln_dm2dc .AND. sn_qsr%nfreqh /= 24 )   &
@@ -210 +221 @@
             END DO
          END DO
+         !                                                        ! add modification due to drag coefficient read from wave forcing 
+         !                                                        ! this code is inefficient but put here to allow merging with another UKMO branch 
+         IF( ln_shelf_flx ) THEN 
+            ! calculate first the wind module, as it will be used later 
+            DO jj = 2, jpjm1 
+               DO ji = fs_2, fs_jpim1    ! vect. opt. 
+                  ztx = zwnd_i(ji-1,jj  ) + zwnd_i(ji,jj) 
+                  zty = zwnd_j(ji  ,jj-1) + zwnd_j(ji,jj) 
+                  wndm(ji,jj) = 0.5 * SQRT( ztx * ztx + zty * zty ) 
+               END DO 
+            END DO 
+            CALL lbc_lnk( wndm(:,:), 'T', 1. ) 
+        
+            IF( ln_cdgw .AND. nn_drag == jp_std ) THEN 
+               IF( cpl_wdrag ) THEN  
+                  ! reset utau and vtau to the wind components: the momentum will 
+                  ! be calculated from the coupled value of the drag coefficient 
+                  DO jj = 1, jpj 
+                     DO ji = 1, jpi 
+                        utau(ji,jj) = zwnd_i(ji,jj) 
+                        vtau(ji,jj) = zwnd_j(ji,jj) 
+                     END DO 
+                  END DO 
+               ELSE 
+                  DO jj = 1, jpjm1 
+                     DO ji = 1, jpim1 
+                        utau(ji,jj) = zrhoa * 0.5 * ( cdn_wave(ji,jj) + cdn_wave(ji+1,jj) ) * zwnd_i(ji,jj) * & 
+                                                                        0.5 * ( wndm(ji,jj) + wndm(ji+1,jj) ) 
+                        vtau(ji,jj) = zrhoa * 0.5 * ( cdn_wave(ji,jj) + cdn_wave(ji,jj+1) ) * zwnd_j(ji,jj) * & 
+                                                                        0.5 * ( wndm(ji,jj) + wndm(ji,jj+1) ) 
+                     END DO 
+                  END DO 
+                  CALL lbc_lnk_multi( utau(:,:), 'U', -1., vtau(:,:), 'V', -1. ) 
+               ENDIF 
+            ELSE IF( nn_drag == jp_const ) THEN 
+               DO jj = 1, jpjm1 
+                  DO ji = 1, jpim1 
+                     utau(ji,jj) = zrhoa * zcdrag * zwnd_i(ji,jj) * 0.5 * ( wndm(ji,jj) + wndm(ji+1,jj) ) 
+                     vtau(ji,jj) = zrhoa * zcdrag * zwnd_j(ji,jj) * 0.5 * ( wndm(ji,jj) + wndm(ji,jj+1) ) 
+                  END DO 
+               END DO 
+               CALL lbc_lnk_multi( utau(:,:), 'U', -1., vtau(:,:), 'V', -1. ) 
+            ENDIF 
+         ENDIF
          !                                                        ! add to qns the heat due to e-p
          qns(:,:) = qns(:,:) - emp(:,:) * sst_m(:,:) * rcp        ! mass flux is at SST
@@ -230 +285 @@
                zmod = 0.5 * SQRT( ztx * ztx + zty * zty )
                taum(ji,jj) = zmod
+               IF( .NOT. ln_shelf_flx ) THEN
                wndm(ji,jj) = SQRT( zmod * zcoef )
+               ENDIF
             END DO
          END DO