- Timestamp:
- 2015-12-07T16:11:45+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_MetOffice_merge_2015/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90
r6006 r6012 119 119 INTEGER :: ikt, ikb 120 120 REAL(wp) :: zfact, z1_e3t, zdep 121 REAL(wp) :: z alpha, zhk121 REAL(wp) :: zt_frz, zpress 122 122 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztrdt, ztrds 123 123 !!---------------------------------------------------------------------- … … 218 218 !---------------------------------------- 219 219 ! 220 IF( nn_isf > 0) THEN220 IF( ln_isf ) THEN 221 221 zfact = 0.5_wp 222 222 DO jj = 2, jpj … … 227 227 228 228 ! level fully include in the ice shelf boundary layer 229 ! if isfdiv, we have to remove heat flux due to inflow at 0oC (as in rnf when you add rnf at sst)230 229 ! sign - because fwf sign of evapo (rnf sign of precip) 231 230 DO jk = ikt, ikb - 1 232 231 ! compute trend 233 tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) & 234 & + zfact * (risf_tsc_b(ji,jj,jp_tem) + risf_tsc(ji,jj,jp_tem)) * r1_hisf_tbl(ji,jj) 235 tsa(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal) & 236 & + zfact * (risf_tsc_b(ji,jj,jp_sal) + risf_tsc(ji,jj,jp_sal)) * r1_hisf_tbl(ji,jj) 232 tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) & 233 & + zfact * ( risf_tsc_b(ji,jj,jp_tem) + risf_tsc(ji,jj,jp_tem) ) & 234 & * r1_hisf_tbl(ji,jj) 237 235 END DO 238 236 239 237 ! level partially include in ice shelf boundary layer 240 238 ! compute trend 241 tsa(ji,jj,ikb,jp_tem) = tsa(ji,jj,ikb,jp_tem) &242 & + zfact * ( risf_tsc_b(ji,jj,jp_tem) + risf_tsc(ji,jj,jp_tem)) * r1_hisf_tbl(ji,jj) * ralpha(ji,jj)243 tsa(ji,jj,ikb,jp_sal) = tsa(ji,jj,ikb,jp_sal) &244 & + zfact * (risf_tsc_b(ji,jj,jp_sal) + risf_tsc(ji,jj,jp_sal)) * r1_hisf_tbl(ji,jj) * ralpha(ji,jj) 239 tsa(ji,jj,ikb,jp_tem) = tsa(ji,jj,ikb,jp_tem) & 240 & + zfact * ( risf_tsc_b(ji,jj,jp_tem) + risf_tsc(ji,jj,jp_tem) ) & 241 & * r1_hisf_tbl(ji,jj) * ralpha(ji,jj) 242 245 243 END DO 246 244 END DO
Note: See TracChangeset
for help on using the changeset viewer.