wake_popdyn_1 Subroutine

public subroutine wake_popdyn_1(klon, klev, dtime, cstar, tau_wk_inv, wgen, wdens, awdens, sigmaw, wdensmin, dtimesub, gfl, rad_wk, f_shear, drdt_pos, d_awdens, d_wdens, d_sigmaw, iflag_wk_act, wk_adv, cin, wape, drdt, d_dens_gen, d_dens_death, d_dens_col, d_dens_bnd, d_sig_gen, d_sig_death, d_sig_col, d_sig_spread, d_sig_bnd, d_wdens_targ, d_sigmaw_targ)

USE lmdz_wake_ini , ONLY : iflag_wk_pop_dyn, wdensmin tau_wk(i) = max(rad_wk(i)/(3.cstar(i))((cstar(i)/cstart)1.5 - 1), 100.) (1 - 2sigmaw(i)(1.-f_shear(i))) d_wdens(i) = ( wgen(i)(1.+2.(sigmaw(i)-sigmad)) & - wdens(i)tau_wk_inv_min & - 2.gfl(i)wdens(i)Cstar(i) )*dtimesub

      d_wdens(i) = ( wgen(i) - (wdens(i)-awdens(i))*tau_wk_inv_min -  &
                     2.*wdens(i)*gfl(i)*drdt_pos )*dtimesub

jyg+mlt d_dens_bnd(i) = d_dens_bnd(i) + d_wdens_targ - d_wdens(i) d_wdens(i) = max(d_wdens(i), wdensmin-wdens(i)) jyg d_sigmaw(i) = ( (1.-2f_shear(i)sigmaw(i))(gfl(i)Cstar(i)+wgen(i)sigmad/wdens(i)) & + 2.f_shear(i)wgen(i)sigmaw(i)2/wdens(i) & - sigmaw(i)tau_wk_inv_min )dtimesub

jyg+mlt d_sig_bnd(i) = d_sig_bnd(i) + d_sigmaw_targ - d_sigmaw(i) d_sig_bnd_provis(i) = d_sigmaw_targ - d_sigmaw(i) d_sigmaw(i) = max(d_sigmaw(i), sigmad-sigmaw(i)) jyg

Arguments