Changeset 7549 for branches/UKMO/dev_r5518_GO6_package/NEMOGCM
- Timestamp:
- 2017-01-11T17:15:08+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90
r7061 r7549 136 136 CALL trd_tra( kt, 'TRA', jp_sal, jptra_zdfp, ztrds ) 137 137 ENDIF 138 ! total trend (note slightly inaccurate because tsb is time-filtered and tsa isn't)138 ! total trend for the non-time-filtered variables. 139 139 ztrdt(:,:,:) = 0._wp 140 140 ztrds(:,:,:) = 0._wp 141 141 DO jk = 1, jpkm1 142 ztrdt(:,:,jk) = ( tsa(:,:,jk,jp_tem) - tsb(:,:,jk,jp_tem) ) / r2dtra(jk) 143 ztrds(:,:,jk) = ( tsa(:,:,jk,jp_sal) - tsb(:,:,jk,jp_sal) ) / r2dtra(jk) 142 zfact = 1.0 / rdttra(jk) 143 ztrdt(:,:,jk) = ( tsa(:,:,jk,jp_tem) - tsn(:,:,jk,jp_tem) ) * zfact 144 ztrds(:,:,jk) = ( tsa(:,:,jk,jp_sal) - tsn(:,:,jk,jp_sal) ) * zfact 144 145 END DO 145 146 CALL trd_tra( kt, 'TRA', jp_tem, jptra_tot, ztrdt )
Note: See TracChangeset
for help on using the changeset viewer.