Changeset 356
- Timestamp:
- Nov 7, 2011, 3:32:08 PM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90
r337 r356 63 63 INTEGER ig,k,l,ll,iq 64 64 INTEGER lmax(ngridmx),lmin(ngridmx),lalim(ngridmx) 65 INTEGER lmix(ngridmx)66 INTEGER lmix_bis(ngridmx)67 65 REAL linter(ngridmx) 68 REAL zmix(ngridmx)69 66 REAL zmax(ngridmx) 70 67 REAL ztva(ngridmx,nlayermx),zw_est(ngridmx,nlayermx+1),ztva_est(ngridmx,nlayermx) … … 124 121 REAL denom(ngridmx) 125 122 REAL zlevinter(ngridmx) 126 127 ! =========================================128 129 ! ============= DRY VARIABLES =============130 131 REAL zw2_dry(ngridmx,nlayermx+1)132 REAL f_star_dry(ngridmx,nlayermx+1)133 REAL ztva_dry(ngridmx,nlayermx+1)134 REAL wmaxa_dry(ngridmx)135 REAL wa_moy_dry(ngridmx,nlayermx+1)136 REAL linter_dry(ngridmx),zlevinter_dry(ngridmx)137 INTEGER lmix_dry(ngridmx),lmax_dry(ngridmx)138 123 139 124 ! ========================================= … … 325 310 326 311 ! Initialisation des variables entieres 327 lmix(:)=1328 lmix_bis(:)=2329 312 wmaxa(:)=0. 330 313 lalim(:)=1 … … 627 610 if (wa_moy(ig,l+1).gt.wmaxa(ig)) then 628 611 ! lmix est le niveau de la couche ou w (wa_moy) est maximum 629 !on rajoute le calcul de lmix_bis630 lmix(ig)=l+1631 612 wmaxa(ig)=wa_moy(ig,l+1) 632 613 endif … … 668 649 669 650 !------------------------------------------------------------------------------- 670 ! Calcul des caracteristiques du thermique:zmax, zmix,wmax651 ! Calcul des caracteristiques du thermique:zmax,wmax 671 652 !------------------------------------------------------------------------------- 672 653 … … 738 719 endif 739 720 enddo 740 741 ! def de zmix continu (profil parabolique des vitesses)742 do ig=1,ngridmx743 if (lmix(ig).gt.1) then744 ! test745 if (((zw2(ig,lmix(ig)-1)-zw2(ig,lmix(ig))) &746 & *((zlev(ig,lmix(ig)))-(zlev(ig,lmix(ig)+1))) &747 & -(zw2(ig,lmix(ig))-zw2(ig,lmix(ig)+1)) &748 & *((zlev(ig,lmix(ig)-1))-(zlev(ig,lmix(ig))))).gt.1e-10) &749 & then750 !751 zmix(ig)=((zw2(ig,lmix(ig)-1)-zw2(ig,lmix(ig))) &752 & *((zlev(ig,lmix(ig))*zlev(ig,lmix(ig)))-(zlev(ig,lmix(ig)+1)*zlev(ig,lmix(ig)+1))) &753 & -(zw2(ig,lmix(ig))-zw2(ig,lmix(ig)+1)) &754 & *((zlev(ig,lmix(ig)-1)*zlev(ig,lmix(ig)-1))-(zlev(ig,lmix(ig))*zlev(ig,lmix(ig))))) &755 & /(2.*((zw2(ig,lmix(ig)-1)-zw2(ig,lmix(ig))) &756 & *((zlev(ig,lmix(ig)))-(zlev(ig,lmix(ig)+1))) &757 & -(zw2(ig,lmix(ig))-zw2(ig,lmix(ig)+1)) &758 & *((zlev(ig,lmix(ig)-1))-(zlev(ig,lmix(ig))))))759 else760 zmix(ig)=zlev(ig,lmix(ig))761 print*,'pb zmix'762 endif763 else764 zmix(ig)=0.765 endif766 !test767 if ((zmax(ig)-zmix(ig)).le.0.) then768 zmix(ig)=0.9*zmax(ig)769 endif770 enddo771 !772 ! calcul du nouveau lmix correspondant773 do ig=1,ngridmx774 do l=1,nlayermx775 if (zmix(ig).ge.zlev(ig,l).and. &776 & zmix(ig).lt.zlev(ig,l+1)) then777 lmix(ig)=l778 endif779 enddo780 enddo781 782 721 783 722 ! Attention, w2 est transforme en sa racine carree dans cette routine
Note: See TracChangeset
for help on using the changeset viewer.