Ignore:
Timestamp:
2019-07-09T12:05:00+02:00 (15 months ago)
Author:
davestorkey
Message:

branches/UKMO/NEMO_4.0_GO8_package : Add uv, uw, vw diagnostics to diaprod.F90.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0_GO8_package/src/OCE/DIA/diaprod.F90

    r11110 r11227  
    163163      ENDIF 
    164164 
     165      IF( iom_use("uv") ) THEN 
     166         z3d(:,:,:) = 0.e0  
     167         DO jk = 1, jpkm1 
     168            DO jj = 2, jpjm1 
     169               DO ji = fs_2, fs_jpim1   ! vector opt. 
     170                  z3d(ji,jj,jk) = 0.25 * ( un(ji-1,jj,jk) + un(ji,jj,jk) ) * ( vn(ji,jj-1,jk) + vn(ji,jj,jk) )  
     171               END DO 
     172            END DO 
     173         END DO 
     174         CALL iom_put( "uv", z3d )                  ! product of zonal velocity and meridional velocity at T points 
     175      ENDIF 
     176 
     177      IF( iom_use("uw") ) THEN 
     178         z3d(:,:,:) = 0.e0  
     179         DO jj = 2, jpjm1 
     180            DO ji = fs_2, fs_jpim1   ! vector opt. 
     181               z3d(ji,jj,1) = 0.5 * ( wn(ji,jj,1) + wn(ji+1,jj,1) ) * un(ji,jj,1)  
     182            END DO 
     183         END DO 
     184         DO jk = 2, jpkm1 
     185            DO jj = 2, jpjm1 
     186               DO ji = fs_2, fs_jpim1   ! vector opt. 
     187                  z3d(ji,jj,jk) = 0.25 * ( wn(ji,jj,jk) + wn(ji+1,jj,jk) ) * ( un(ji,jj,jk-1) + un(ji,jj,jk) )  
     188               END DO 
     189            END DO 
     190         END DO 
     191         CALL iom_put( "uw", z3d )                  ! product of zonal velocity and vertical velocity at UW points 
     192      ENDIF 
     193 
     194      IF( iom_use("vw") ) THEN 
     195         z3d(:,:,:) = 0.e0  
     196         DO jj = 2, jpjm1 
     197            DO ji = fs_2, fs_jpim1   ! vector opt. 
     198               z3d(ji,jj,1) = 0.5 * ( wn(ji,jj,1) + wn(ji,jj+1,1) ) * vn(ji,jj,1)  
     199            END DO 
     200         END DO 
     201         DO jk = 2, jpkm1 
     202            DO jj = 2, jpjm1 
     203               DO ji = fs_2, fs_jpim1   ! vector opt. 
     204                  z3d(ji,jj,jk) = 0.25 * ( wn(ji,jj,jk) + wn(ji,jj+1,jk) ) * ( vn(ji,jj,jk-1) + vn(ji,jj,jk) )  
     205               END DO 
     206            END DO 
     207         END DO 
     208         CALL iom_put( "vw", z3d )                  ! product of meriodional velocity and vertical velocity at VW points 
     209      ENDIF 
     210 
    165211      IF( iom_use("urhop") ) THEN 
    166212         z3d(:,:,:) = 0.e0  
Note: See TracChangeset for help on using the changeset viewer.