New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 11227 – NEMO

Changeset 11227


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

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

Location:
NEMO/branches/UKMO/NEMO_4.0_GO8_package
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0_GO8_package/cfgs/SHARED/field_def_nemo-oce.xml

    r11110 r11227  
    529529   <field id="ws"           long_name="product_of_upward_sea_water_velocity_and_salinity"              unit="PSU m/s"       grid_ref="grid_W_3D"   /> 
    530530   <field id="wrhop"        long_name="product_of_upward_sea_water_velocity_and_potential_density"     unit="(kg/m3).(m/s)" grid_ref="grid_W_3D"   /> 
     531   <field id="uv"           long_name="product_of_sea_water_x_velocity_and_sea_water_y_velocity"       unit="m2/s2   "      grid_ref="grid_T_3D"   /> 
     532   <field id="uw"           long_name="product_of_upward_sea_water_velocity_and_sea_water_x_velocity"  unit="m2/s2   "      grid_ref="grid_W_3D"   /> 
     533   <field id="vw"           long_name="product_of_upward_sea_water_velocity_and_sea_water_y_velocity"  unit="m2/s2"         grid_ref="grid_W_3D"   /> 
    531534      </field_group> 
    532535 
  • 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.