- Timestamp:
- 2020-12-10T18:29:37+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90
r10302 r14152 291 291 IF( iom_use('rnf_x_sss') ) CALL iom_put( "rnf_x_sss", rnf*tsn(:,:,1,jp_sal) ) ! runoff term on sss 292 292 293 #if defined key_asminc294 ! 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_0296 ! 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 assimilation298 DO jj = 2, jpj299 DO ji = fs_2, fs_jpim1300 zdep = ssh_iau(ji,jj) / ( ht_0(ji,jj) + 1.0 - ssmask(ji, jj) )301 DO jk = 1, jpkm1302 tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) &303 & + tsn(ji,jj,jk,jp_tem) * 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 & + tsn(ji,jj,jk,jp_sal) * zdep * ( e3t_0(ji,jj,jk) / fse3t_n(ji,jj,jk) )306 END DO307 END DO308 END DO309 ENDIF310 #endif293 !#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 ! & + tsn(ji,jj,jk,jp_tem) * 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 ! & + tsn(ji,jj,jk,jp_sal) * zdep * ( e3t_0(ji,jj,jk) / fse3t_n(ji,jj,jk) ) 306 ! END DO 307 ! END DO 308 ! END DO 309 ! ENDIF 310 !#endif 311 311 312 312 IF( l_trdtra ) THEN ! send trends for further diagnostics
Note: See TracChangeset
for help on using the changeset viewer.