Changeset 7717 for branches/UKMO/r6232_INGV1_WAVE-coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcflx.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/sbcflx.F90 (modified) (5 diffs)

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

@@ -22 +22 @@
    USE lib_mpp         ! distribued memory computing library
    USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
+   USE sbcwave         ! wave physics
 
    IMPLICIT NONE
@@ -35 +36 @@
    INTEGER , PARAMETER ::   jp_emp  = 5   ! index of evaporation-precipation file
    TYPE(FLD), ALLOCATABLE, DIMENSION(:) ::   sf    ! structure of input fields (file informations, fields read)
+   LOGICAL , PUBLIC    ::   ln_shelf_flx = .FALSE. ! UKMO SHELF specific flux flag
+   INTEGER             ::   jpfld_local   ! maximum number of files to read (locally modified depending on ln_shelf_flx) 
 
    !! * Substitutions
@@ -81 +84 @@
       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
       !!
       CHARACTER(len=100) ::  cn_dir                               ! Root directory for location of flx files
+      NAMELIST/namsbc_flx/ ln_shelf_flx                           ! Put here to allow merging with another UKMO branch
       TYPE(FLD_N), DIMENSION(jpfld) ::   slf_i                    ! array of namelist information structures
       TYPE(FLD_N) ::   sn_utau, sn_vtau, sn_qtot, sn_qsr, sn_emp  ! informations about the fields to be read
@@ -140 +145 @@
             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 .AND. 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) = sf(jp_utau)%fnow(ji,jj,1)
+                     vtau(ji,jj) = sf(jp_vtau)%fnow(ji,jj,1)
+                  END DO
+               END DO
+            ELSE
+               DO jj = 1, jpj
+                  DO ji = 1, jpi
+                     totwind = sqrt((sf(jp_utau)%fnow(ji,jj,1))**2.0 + (sf(jp_vtau)%fnow(ji,jj,1))**2.0)
+                     utau(ji,jj) = zrhoa * cdn_wave(ji,jj) * sf(jp_utau)%fnow(ji,jj,1) * totwind
+                     vtau(ji,jj) = zrhoa * cdn_wave(ji,jj) * sf(jp_vtau)%fnow(ji,jj,1) * totwind
+                  END DO
+               END DO
+            ENDIF
+         ENDIF
          !                                                        ! add to qns the heat due to e-p
          qns(:,:) = qns(:,:) - emp(:,:) * sst_m(:,:) * rcp        ! mass flux is at SST
@@ -153 +180 @@
                zmod = 0.5 * SQRT( ztx * ztx + zty * zty )
                taum(ji,jj) = zmod
+               IF( ln_shelf_flx ) THEN
+                  ztx = sf(jp_utau)%fnow(ji-1,jj  ) + sf(jp_utau)%fnow(ji,jj)
+                  zty = sf(jp_vtau)%fnow(ji  ,jj-1) + sf(jp_vtau)%fnow(ji,jj)
+                  wndm(ji,jj) = 0.5 * SQRT( ztx * ztx + zty * zty )
+               ELSE
                wndm(ji,jj) = SQRT( zmod * zcoef )
+               ENDIF
             END DO
          END DO