- Timestamp:
- 2015-02-11T16:15:11+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limthd_lac.F90
r5070 r5078 362 362 DO ji = 1, nbpac 363 363 364 zEi = - ze_newice(ji) / rhoic! specific enthalpy of forming ice [J/kg]365 366 zEw = rcp * ( t_bo_1d(ji) - rt0 ) 364 zEi = - ze_newice(ji) * r1_rhoic ! specific enthalpy of forming ice [J/kg] 365 366 zEw = rcp * ( t_bo_1d(ji) - rt0 ) ! specific enthalpy of seawater at t_bo_1d [J/kg] 367 367 ! clem: we suppose we are already at the freezing point (condition qlead<0 is satisfyied) 368 368 … … 371 371 zfmdt = - qlead_1d(ji) / zdE ! Fm.dt [kg/m2] (<0) 372 372 ! clem: we use qlead instead of zqld (limthd) because we suppose we are at the freezing point 373 zv_newice(ji) = - zfmdt /rhoic373 zv_newice(ji) = - zfmdt * r1_rhoic 374 374 375 375 zQm = zfmdt * zEw ! heat to the ocean >0 associated with mass flux … … 458 458 DO jk = 1, nlay_i 459 459 DO ji = 1, nbpac 460 h_i_old (ji,jk) = zv_i_1d(ji,jl) / REAL( nlay_i )460 h_i_old (ji,jk) = zv_i_1d(ji,jl) * r1_nlay_i 461 461 qh_i_old(ji,jk) = ze_i_1d(ji,jk,jl) * h_i_old(ji,jk) 462 462 END DO … … 525 525 DO ji = 1, jpi 526 526 ! heat content in J/m2 527 e_i(ji,jj,jk,jl) = e_i(ji,jj,jk,jl) * v_i(ji,jj,jl) / REAL( nlay_i ,wp )527 e_i(ji,jj,jk,jl) = e_i(ji,jj,jk,jl) * v_i(ji,jj,jl) * r1_nlay_i 528 528 END DO 529 529 END DO
Note: See TracChangeset
for help on using the changeset viewer.