Changeset 6592 for branches/UKMO/dev_r5518_Surge_Modelling/NEMOGCM/NEMO

Timestamp:

2016-05-20T16:06:01+02:00 (8 years ago)

Author:

rfurner

Message:

changes for outputting and namelist changes

Location:

branches/UKMO/dev_r5518_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC/SBC

Files:

• sbc_oce.F90 (modified) (1 diff)
• sbcsurge.F90 (modified) (2 diffs)

branches/UKMO/dev_r5518_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_oce.F90

@@ -102 +102 @@
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   utau   , utau_b   !: sea surface i-stress (ocean referential)     [N/m2]
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   vtau   , vtau_b   !: sea surface j-stress (ocean referential)     [N/m2]
+   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   uwnd              !: wind speed i-component                       [m2]
+   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   vwnd              !: wind speed j-component                       [m2]
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   taum              !: module of sea surface stress (at T-point)    [N/m2] 
    !! wndm is used onmpute surface gases exchanges in ice-free ocean or leads

branches/UKMO/dev_r5518_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcsurge.F90

@@ -172 +172 @@
       DO jj = 2, jpjm1
          DO ji = fs_2, fs_jpim1   ! vect. opt.
-            zwnd_i(ji,jj) = (  sf(jp_wndi)%fnow(ji,jj,1) - rn_vfac * 0.5 * ( pu(ji-1,jj  ) + pu(ji,jj) )  )
-            zwnd_j(ji,jj) = (  sf(jp_wndj)%fnow(ji,jj,1) - rn_vfac * 0.5 * ( pv(ji  ,jj-1) + pv(ji,jj) )  )
+            uwnd(ji,jj) = (  sf(jp_wndi)%fnow(ji,jj,1) - rn_vfac * 0.5 * ( pu(ji-1,jj  ) + pu(ji,jj) )  )
+            vwnd(ji,jj) = (  sf(jp_wndj)%fnow(ji,jj,1) - rn_vfac * 0.5 * ( pv(ji  ,jj-1) + pv(ji,jj) )  )
          END DO
       END DO
+      zwnd_i(:,:) = uwnd(:,:)
+      zwnd_j(:,:) = vwnd(:,:)
       CALL lbc_lnk( zwnd_i(:,:) , 'T', -1. )
       CALL lbc_lnk( zwnd_j(:,:) , 'T', -1. )
@@ -217 +219 @@
 
       CALL iom_put( "taum_oce", taum )   ! output wind stress module
+      CALL iom_put( "uwnd", uwnd )   ! output wind stress module
+      CALL iom_put( "vwnd", vwnd )   ! output wind stress module
 
       ! ... utau, vtau at U- and V_points, resp.