- Timestamp:
- Aug 2, 2024, 2:12:03 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_dry.F90
r5117 r5158 38 38 39 39 !initialisations 40 doig=1,ngrid41 dol=1,nlay+140 DO ig=1,ngrid 41 DO l=1,nlay+1 42 42 zw2(ig,l)=0. 43 43 wa_moy(ig,l)=0. 44 44 enddo 45 45 enddo 46 doig=1,ngrid47 dol=1,nlay46 DO ig=1,ngrid 47 DO l=1,nlay 48 48 ztva(ig,l)=ztv(ig,l) 49 49 enddo 50 50 enddo 51 doig=1,ngrid51 DO ig=1,ngrid 52 52 wmax(ig)=0. 53 53 wmaxa(ig)=0. … … 58 58 ! Calcul des F^*, integrale verticale de E^* 59 59 f_star(:,1)=0. 60 dol=1,nlay60 DO l=1,nlay 61 61 f_star(:,l+1)=f_star(:,l)+alim_star(:,l) 62 62 enddo … … 71 71 ! lmax 72 72 73 dol=1,nlay-274 doig=1,ngrid73 DO l=1,nlay-2 74 DO ig=1,ngrid 75 75 IF (l==lmin(ig).AND.lalim(ig)>1) THEN 76 76 !------------------------------------------------------------------------ … … 136 136 137 137 ! Determination de zw2 max 138 doig=1,ngrid138 DO ig=1,ngrid 139 139 wmax(ig)=0. 140 140 enddo 141 141 142 dol=1,nlay143 doig=1,ngrid142 DO l=1,nlay 143 DO ig=1,ngrid 144 144 IF (l<=lmax(ig)) THEN 145 145 zw2(ig,l)=sqrt(zw2(ig,l)) … … 152 152 153 153 ! Longueur caracteristique correspondant a la hauteur des thermiques. 154 doig=1,ngrid154 DO ig=1,ngrid 155 155 zmax(ig)=0. 156 156 zlevinter(ig)=zlev(ig,1) 157 157 enddo 158 doig=1,ngrid158 DO ig=1,ngrid 159 159 ! calcul de zlevinter 160 160 zlevinter(ig)=zlev(ig,lmax(ig)) + &
Note: See TracChangeset
for help on using the changeset viewer.