| [1403] | 1 | ! |
|---|
| 2 | ! $Header$ |
|---|
| 3 | ! |
|---|
| [972] | 4 | SUBROUTINE thermcell_init(ngrid,nlay,ztv,zlay,zlev, & |
|---|
| [878] | 5 | & lalim,lmin,alim_star,alim_star_tot,lev_out) |
|---|
| 6 | |
|---|
| 7 | !---------------------------------------------------------------------- |
|---|
| 8 | !thermcell_init: calcul du profil d alimentation du thermique |
|---|
| 9 | !---------------------------------------------------------------------- |
|---|
| [2408] | 10 | USE print_control_mod, ONLY: lunout |
|---|
| [878] | 11 | IMPLICIT NONE |
|---|
| [1026] | 12 | #include "thermcell.h" |
|---|
| [938] | 13 | |
|---|
| [878] | 14 | INTEGER l,ig |
|---|
| 15 | !arguments d entree |
|---|
| 16 | INTEGER ngrid,nlay |
|---|
| 17 | REAL ztv(ngrid,nlay) |
|---|
| [972] | 18 | REAL zlay(ngrid,nlay) |
|---|
| 19 | REAL zlev(ngrid,nlay+1) |
|---|
| [878] | 20 | !arguments de sortie |
|---|
| 21 | INTEGER lalim(ngrid) |
|---|
| 22 | INTEGER lmin(ngrid) |
|---|
| 23 | REAL alim_star(ngrid,nlay) |
|---|
| 24 | REAL alim_star_tot(ngrid) |
|---|
| 25 | integer lev_out ! niveau pour les print |
|---|
| 26 | |
|---|
| [972] | 27 | REAL zzalim(ngrid) |
|---|
| [878] | 28 | !CR: ponderation entrainement des couches instables |
|---|
| 29 | !def des alim_star tels que alim=f*alim_star |
|---|
| 30 | |
|---|
| 31 | |
|---|
| [1403] | 32 | write(lunout,*)'THERM INIT V20C ' |
|---|
| [1026] | 33 | |
|---|
| [1403] | 34 | alim_star_tot(:)=0. |
|---|
| 35 | alim_star(:,:)=0. |
|---|
| 36 | lmin(:)=1 |
|---|
| 37 | lalim(:)=1 |
|---|
| [878] | 38 | |
|---|
| [1403] | 39 | do l=1,nlay-1 |
|---|
| [878] | 40 | do ig=1,ngrid |
|---|
| [1403] | 41 | if (ztv(ig,l)> ztv(ig,l+1) .and. ztv(ig,1)>=ztv(ig,l) ) then |
|---|
| 42 | alim_star(ig,l)=MAX((ztv(ig,l)-ztv(ig,l+1)),0.) & |
|---|
| [878] | 43 | & *sqrt(zlev(ig,l+1)) |
|---|
| [1403] | 44 | lalim(:)=l+1 |
|---|
| 45 | alim_star_tot(ig)=alim_star_tot(ig)+alim_star(ig,l) |
|---|
| [878] | 46 | endif |
|---|
| 47 | enddo |
|---|
| 48 | enddo |
|---|
| 49 | do l=1,nlay |
|---|
| 50 | do ig=1,ngrid |
|---|
| [1403] | 51 | if (alim_star_tot(ig) > 1.e-10 ) then |
|---|
| [878] | 52 | alim_star(ig,l)=alim_star(ig,l)/alim_star_tot(ig) |
|---|
| 53 | endif |
|---|
| 54 | enddo |
|---|
| 55 | enddo |
|---|
| [1403] | 56 | alim_star_tot(:)=1. |
|---|
| [1026] | 57 | |
|---|
| [1403] | 58 | return |
|---|
| [878] | 59 | end |
|---|
