- Timestamp:
- 2017-06-27T15:26:53+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package_inc_asm/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90
r7947 r8225 236 236 DO jk = ikt, ikb - 1 237 237 ! compute tfreez for the temperature correction (we add water at freezing temperature) 238 ! zpress = grav*rau0*fsdept(ji,jj,jk)*1.e-04239 zt_frz = -1.9 !eos_fzp( tsn(ji,jj,jk,jp_sal), zpress )240 238 ! compute trend 241 239 tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) & 242 & + zfact * (risf_tsc_b(ji,jj,jp_tem) + risf_tsc(ji,jj,jp_tem) & 243 & - rdivisf * (fwfisf(ji,jj) + fwfisf_b(ji,jj)) * zt_frz * r1_rau0) & 244 & * r1_hisf_tbl(ji,jj) 240 & + zfact * (risf_tsc_b(ji,jj,jp_tem) + risf_tsc(ji,jj,jp_tem)) * r1_hisf_tbl(ji,jj) 245 241 tsa(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal) & 246 242 & + zfact * (risf_tsc_b(ji,jj,jp_sal) + risf_tsc(ji,jj,jp_sal)) * r1_hisf_tbl(ji,jj) … … 249 245 ! level partially include in ice shelf boundary layer 250 246 ! compute tfreez for the temperature correction (we add water at freezing temperature) 251 ! zpress = grav*rau0*fsdept(ji,jj,ikb)*1.e-04252 zt_frz = -1.9 !eos_fzp( tsn(ji,jj,ikb,jp_sal), zpress )253 247 ! compute trend 254 248 tsa(ji,jj,ikb,jp_tem) = tsa(ji,jj,ikb,jp_tem) & 255 & + zfact * (risf_tsc_b(ji,jj,jp_tem) + risf_tsc(ji,jj,jp_tem) & 256 & - rdivisf * (fwfisf(ji,jj) + fwfisf_b(ji,jj)) * zt_frz * r1_rau0) & 257 & * r1_hisf_tbl(ji,jj) * ralpha(ji,jj) 249 & + zfact * (risf_tsc_b(ji,jj,jp_tem) + risf_tsc(ji,jj,jp_tem)) * r1_hisf_tbl(ji,jj) * ralpha(ji,jj) 258 250 tsa(ji,jj,ikb,jp_sal) = tsa(ji,jj,ikb,jp_sal) & 259 251 & + zfact * (risf_tsc_b(ji,jj,jp_sal) + risf_tsc(ji,jj,jp_sal)) * r1_hisf_tbl(ji,jj) * ralpha(ji,jj)
Note: See TracChangeset
for help on using the changeset viewer.