Changeset 6673


Ignore:
Timestamp:
2016-06-07T17:07:47+02:00 (4 years ago)
Author:
timgraham
Message:

Bugfix for heat transport from vT

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/v3_6_extra_CMIP6_diagnostics/NEMOGCM/NEMO/OPA_SRC/DIA/diaprod.F90

    r6672 r6673  
    112112            DO jj = 2, jpjm1 
    113113               DO ji = fs_2, fs_jpim1   ! vector opt. 
    114                   z3d(ji,jj,jk) = vn(ji,jj,jk) * fse3v(ji,jj,jk) * e1v(ji,jj) 
     114                  z3d(ji,jj,jk) = z3d(ji,jj,jk) * fse3v(ji,jj,jk) * e1v(ji,jj) 
    115115               END DO 
    116116            END DO 
     
    148148      ENDIF 
    149149 
    150       IF( iom_use("vs") ) THEN 
     150      IF( iom_use("vs") .OR. iom_use("sopst_vs") ) THEN 
    151151         z3d(:,:,:) = 0.e0  
    152152         DO jk = 1, jpkm1 
     
    158158         END DO 
    159159         CALL iom_put( "vs", z3d )                  ! product of salinity and meridional velocity at V points 
     160         DO jk = 1, jpkm1 
     161            DO jj = 2, jpjm1 
     162               DO ji = fs_2, fs_jpim1   ! vector opt. 
     163                  z3d(ji,jj,jk) = z3d(ji,jj,jk) * fse3v(ji,jj,jk) * e1v(ji,jj) 
     164               END DO 
     165            END DO 
     166         END DO 
     167         CALL dia_ptr_ohst_components( jp_sal, 'vts', z3d) 
    160168      ENDIF 
    161169 
Note: See TracChangeset for help on using the changeset viewer.