- Timestamp:
- 2017-12-01T09:43:23+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_METO_2017/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90
r7788 r8868 27 27 USE trd_oce ! trends: ocean variables 28 28 USE trdtra ! trends manager: tracers 29 #if defined key_asminc 30 USE asminc ! Assimilation increment 31 #endif 29 32 ! 30 33 USE in_out_manager ! I/O manager … … 72 75 INTEGER, INTENT(in) :: kt ! ocean time-step index 73 76 ! 74 INTEGER :: ji, jj, jk, jn ! dummy loop indices75 INTEGER :: ikt, ikb ! local integers76 REAL(wp) :: zfact, z1_e3t, zdep ! local scalar77 INTEGER :: ji, jj, jk, jn ! dummy loop indices 78 INTEGER :: ikt, ikb ! local integers 79 REAL(wp) :: zfact, z1_e3t, zdep, ztim ! local scalar 77 80 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztrdt, ztrds 78 81 !!---------------------------------------------------------------------- … … 208 211 IF( iom_use('rnf_x_sss') ) CALL iom_put( "rnf_x_sss", rnf*tsn(:,:,1,jp_sal) ) ! runoff term on sss 209 212 213 #if defined key_asminc 214 ! 215 !---------------------------------------- 216 ! Assmilation effects 217 !---------------------------------------- 218 ! 219 IF( ln_sshinc ) THEN ! input of heat and salt due to assimilation 220 ! 221 IF( ln_linssh ) THEN 222 DO jj = 2, jpj 223 DO ji = fs_2, fs_jpim1 224 ztim = ssh_iau(ji,jj) / e3t_n(ji,jj,1) 225 tsa(ji,jj,1,jp_tem) = tsa(ji,jj,1,jp_tem) + tsn(ji,jj,1,jp_tem) * ztim 226 tsa(ji,jj,1,jp_sal) = tsa(ji,jj,1,jp_sal) + tsn(ji,jj,1,jp_sal) * ztim 227 END DO 228 END DO 229 ELSE 230 DO jj = 2, jpj 231 DO ji = fs_2, fs_jpim1 232 ztim = ssh_iau(ji,jj) / ( ht_n(ji,jj) + 1. - ssmask(ji, jj) ) 233 tsa(ji,jj,:,jp_tem) = tsa(ji,jj,:,jp_tem) + tsn(ji,jj,:,jp_tem) * ztim 234 tsa(ji,jj,:,jp_sal) = tsa(ji,jj,:,jp_sal) + tsn(ji,jj,:,jp_sal) * ztim 235 END DO 236 END DO 237 ENDIF 238 ! 239 ENDIF 240 ! 241 #endif 242 210 243 ! 211 244 !----------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.