- Timestamp:
- 2020-02-24T14:00:21+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/LDF/ldftra.F90
r12377 r12443 820 820 ! 821 821 IF( iom_use('weiv_masstr') ) THEN ! vertical mass transport & its square value 822 zw2d(:,:) = r au0 * e1e2t(:,:)822 zw2d(:,:) = rho0 * e1e2t(:,:) 823 823 DO jk = 1, jpk 824 824 zw3d(:,:,jk) = zw3d(:,:,jk) * zw2d(:,:) … … 830 830 zw3d(:,:,:) = 0.e0 831 831 DO jk = 1, jpkm1 832 zw3d(:,:,jk) = r au0 * ( psi_uw(:,:,jk+1) - psi_uw(:,:,jk) )832 zw3d(:,:,jk) = rho0 * ( psi_uw(:,:,jk+1) - psi_uw(:,:,jk) ) 833 833 END DO 834 834 CALL iom_put( "ueiv_masstr", zw3d ) ! mass transport in i-direction 835 835 ENDIF 836 836 ! 837 zztmp = 0.5_wp * r au0 * rcp837 zztmp = 0.5_wp * rho0 * rcp 838 838 IF( iom_use('ueiv_heattr') .OR. iom_use('ueiv_heattr3d') ) THEN 839 839 zw2d(:,:) = 0._wp … … 853 853 zw3d(:,:,:) = 0.e0 854 854 DO jk = 1, jpkm1 855 zw3d(:,:,jk) = r au0 * ( psi_vw(:,:,jk+1) - psi_vw(:,:,jk) )855 zw3d(:,:,jk) = rho0 * ( psi_vw(:,:,jk+1) - psi_vw(:,:,jk) ) 856 856 END DO 857 857 CALL iom_put( "veiv_masstr", zw3d ) ! mass transport in i-direction
Note: See TracChangeset
for help on using the changeset viewer.