- Timestamp:
- 2018-10-29T15:20:26+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_AMM15_package/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90
r10249 r10251 25 25 USE trd_oce ! trends: ocean variables 26 26 USE trdtra ! trends manager: tracers 27 USE tradwl ! solar radiation penetration (downwell method)28 27 ! 29 28 USE in_out_manager ! I/O manager … … 34 33 USE timing ! Timing 35 34 USE eosbn2 36 #if defined key_asminc37 USE asminc ! Assimilation increment38 #endif39 35 40 36 IMPLICIT NONE … … 124 120 REAL(wp) :: zfact, z1_e3t, zdep 125 121 REAL(wp) :: zalpha, zhk 122 REAL(wp) :: zt_frz, zpress 126 123 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztrdt, ztrds 127 124 !!---------------------------------------------------------------------- … … 142 139 143 140 !!gm IF( .NOT.ln_traqsr ) qsr(:,:) = 0.e0 ! no solar radiation penetration 144 IF( .NOT.ln_traqsr .and. .NOT.ln_tradwl) THEN ! no solar radiation penetration141 IF( .NOT.ln_traqsr ) THEN ! no solar radiation penetration 145 142 qns(:,:) = qns(:,:) + qsr(:,:) ! total heat flux in qns 146 143 qsr(:,:) = 0.e0 ! qsr set to zero … … 235 232 DO jk = ikt, ikb - 1 236 233 ! compute tfreez for the temperature correction (we add water at freezing temperature) 234 ! zpress = grav*rau0*fsdept(ji,jj,jk)*1.e-04 235 zt_frz = -1.9 !eos_fzp( tsn(ji,jj,jk,jp_sal), zpress ) 237 236 ! compute trend 238 237 tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) & 239 & + zfact * (risf_tsc_b(ji,jj,jp_tem) + risf_tsc(ji,jj,jp_tem)) * r1_hisf_tbl(ji,jj) 238 & + zfact * (risf_tsc_b(ji,jj,jp_tem) + risf_tsc(ji,jj,jp_tem) & 239 & - rdivisf * (fwfisf(ji,jj) + fwfisf_b(ji,jj)) * zt_frz * r1_rau0) & 240 & * r1_hisf_tbl(ji,jj) 240 241 tsa(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal) & 241 242 & + zfact * (risf_tsc_b(ji,jj,jp_sal) + risf_tsc(ji,jj,jp_sal)) * r1_hisf_tbl(ji,jj) … … 244 245 ! level partially include in ice shelf boundary layer 245 246 ! compute tfreez for the temperature correction (we add water at freezing temperature) 247 ! zpress = grav*rau0*fsdept(ji,jj,ikb)*1.e-04 248 zt_frz = -1.9 !eos_fzp( tsn(ji,jj,ikb,jp_sal), zpress ) 246 249 ! compute trend 247 250 tsa(ji,jj,ikb,jp_tem) = tsa(ji,jj,ikb,jp_tem) & 248 & + zfact * (risf_tsc_b(ji,jj,jp_tem) + risf_tsc(ji,jj,jp_tem)) * r1_hisf_tbl(ji,jj) * ralpha(ji,jj) 251 & + zfact * (risf_tsc_b(ji,jj,jp_tem) + risf_tsc(ji,jj,jp_tem) & 252 & - rdivisf * (fwfisf(ji,jj) + fwfisf_b(ji,jj)) * zt_frz * r1_rau0) & 253 & * r1_hisf_tbl(ji,jj) * ralpha(ji,jj) 249 254 tsa(ji,jj,ikb,jp_sal) = tsa(ji,jj,ikb,jp_sal) & 250 255 & + zfact * (risf_tsc_b(ji,jj,jp_sal) + risf_tsc(ji,jj,jp_sal)) * r1_hisf_tbl(ji,jj) * ralpha(ji,jj) … … 282 287 END DO 283 288 ENDIF 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 289 304 290 IF( l_trdtra ) THEN ! send trends for further diagnostics
Note: See TracChangeset
for help on using the changeset viewer.