- Timestamp:
- 2015-12-04T17:56:07+01:00 (8 years ago)
- Location:
- branches/2015/dev_MetOffice_merge_2015/NEMOGCM/NEMO/OPA_SRC/TRA
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_MetOffice_merge_2015/NEMOGCM/NEMO/OPA_SRC/TRA/tradmp.F90
r5836 r6006 100 100 ! 101 101 CALL wrk_alloc( jpi,jpj,jpk,jpts, zts_dta ) 102 !103 102 IF( l_trdtra ) THEN !* Save ta and sa trends 104 103 CALL wrk_alloc( jpi,jpj,jpk,jpts, ztrdts ) -
branches/2015/dev_MetOffice_merge_2015/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90
r5643 r6006 22 22 USE sbcrnf ! River runoff 23 23 USE sbcisf ! Ice shelf 24 USE iscplini ! Ice sheet coupling 24 25 USE traqsr ! solar radiation penetration 25 26 USE trd_oce ! trends: ocean variables … … 117 118 INTEGER :: ji, jj, jk, jn ! dummy loop indices 118 119 INTEGER :: ikt, ikb 119 INTEGER :: nk_isf120 120 REAL(wp) :: zfact, z1_e3t, zdep 121 121 REAL(wp) :: zalpha, zhk … … 219 219 ! 220 220 IF( nn_isf > 0 ) THEN 221 zfact = 0.5 e0221 zfact = 0.5_wp 222 222 DO jj = 2, jpj 223 223 DO ji = fs_2, fs_jpim1 … … 230 230 ! sign - because fwf sign of evapo (rnf sign of precip) 231 231 DO jk = ikt, ikb - 1 232 ! compute tfreez for the temperature correction (we add water at freezing temperature)233 232 ! compute trend 234 233 tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) & … … 239 238 240 239 ! level partially include in ice shelf boundary layer 241 ! compute tfreez for the temperature correction (we add water at freezing temperature)242 240 ! compute trend 243 241 tsa(ji,jj,ikb,jp_tem) = tsa(ji,jj,ikb,jp_tem) & … … 279 277 ENDIF 280 278 281 IF( l_trdtra ) THEN ! send trends for further diagnostics 279 !---------------------------------------- 280 ! Ice Sheet coupling imbalance correction to have conservation 281 !---------------------------------------- 282 ! 283 IF( ln_iscpl .AND. ln_hsb) THEN ! input of heat and salt due to river runoff 284 DO jk = 1,jpk 285 DO jj = 2, jpj 286 DO ji = fs_2, fs_jpim1 287 zdep = 1._wp / fse3t_n(ji,jj,jk) 288 tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) - htsc_iscpl(ji,jj,jk,jp_tem) & 289 & * zdep 290 tsa(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal) - htsc_iscpl(ji,jj,jk,jp_sal) & 291 & * zdep 292 END DO 293 END DO 294 END DO 295 ENDIF 296 297 IF( l_trdtra ) THEN ! save the horizontal diffusive trends for further diagnostics 282 298 ztrdt(:,:,:) = tsa(:,:,:,jp_tem) - ztrdt(:,:,:) 283 299 ztrds(:,:,:) = tsa(:,:,:,jp_sal) - ztrds(:,:,:)
Note: See TracChangeset
for help on using the changeset viewer.