Changeset 14153
- Timestamp:
- 2020-12-10T18:32:07+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/AMM15_v3_6_STABLE_package_collate_sshiau/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90
r9178 r14153 283 283 ENDIF 284 284 285 #if defined key_asminc286 ! WARNING: THIS MAY WELL NOT BE REQUIRED - WE DON'T WANT TO CHANGE T&S BUT THIS MAY COMPENSATE ANOTHER TERM...287 ! Rate of change in e3t for each level is ssh_iau*e3t_0/ht_0288 ! Contribution to tsa should be rate of change in level / per m of ocean? (hence the division by fse3t_n)289 IF( ln_sshinc ) THEN ! input of heat and salt due to assimilation290 DO jj = 2, jpj291 DO ji = fs_2, fs_jpim1292 zdep = ssh_iau(ji,jj) / ( ht_0(ji,jj) + 1.0 - ssmask(ji, jj) )293 DO jk = 1, jpkm1294 tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) &295 & + tsn(ji,jj,jk,jp_tem) * zdep * ( e3t_0(ji,jj,jk) / fse3t_n(ji,jj,jk) )296 tsa(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal) &297 & + tsn(ji,jj,jk,jp_sal) * zdep * ( e3t_0(ji,jj,jk) / fse3t_n(ji,jj,jk) )298 END DO299 END DO300 END DO301 ENDIF302 #endif303 304 285 IF( l_trdtra ) THEN ! send trends for further diagnostics 305 286 ztrdt(:,:,:) = tsa(:,:,:,jp_tem) - ztrdt(:,:,:)
Note: See TracChangeset
for help on using the changeset viewer.