Changeset 6620 for branches/UKMO/dev_r5518_sshinc_with_VVL
- Timestamp:
- 2016-05-25T16:09:21+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_sshinc_with_VVL/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90
r6605 r6620 33 33 USE timing ! Timing 34 34 USE eosbn2 35 #if defined key_asminc 36 USE asminc ! Assimilation increment 37 #endif 35 38 36 39 IMPLICIT NONE … … 287 290 END DO 288 291 ENDIF 292 293 #if defined key_asminc 294 ! WARNING: THIS MAY WELL NOT BE REQUIRED - WE DON'T WANT TO CHANGE T&S BUT THIS MAY COMPENSATE ANOTHER TERM... 295 ! Rate of change in e3t for each level is ssh_iau*e3t_0/ht_0 296 ! Contribution to tsa should be rate of change in level / per m of ocean? (hence the division by fse3t_n) 297 IF( ln_sshinc ) THEN ! input of heat and salt due to assimilation 298 DO jj = 2, jpj 299 DO ji = fs_2, fs_jpim1 300 zdep = ssh_iau(ji,jj) / ( ht_0(ji,jj) + 1.0 - ssmask(ji, jj) ) 301 DO jk = 1, jpkm1 302 tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) & 303 & * ( 1.0 + zdep * ( e3t_0(ji,jj,jk) / fse3t_n(ji,jj,jk) ) ) 304 tsa(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal) & 305 & * ( 1.0 + zdep * ( e3t_0(ji,jj,jk) / fse3t_n(ji,jj,jk) ) ) 306 END DO 307 END DO 308 END DO 309 ENDIF 310 #endif 289 311 290 312 IF( l_trdtra ) THEN ! send trends for further diagnostics
Note: See TracChangeset
for help on using the changeset viewer.