Changeset 938 for LMDZ4/trunk/libf/phylmd/thermcell.F
- Timestamp:
- Mar 31, 2008, 6:14:45 PM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk/libf/phylmd/thermcell.F
r878 r938 991 991 #ifdef troisD 992 992 c if (sorties) then 993 print*,'Debut des wrgradsfi'993 c print*,'Debut des wrgradsfi' 994 994 995 995 c print*,'16 OK convect8' … … 1068 1068 c print*,'18 OK convect8' 1069 1069 c endif 1070 print*,'Fin des wrgradsfi'1070 c print*,'Fin des wrgradsfi' 1071 1071 #endif 1072 1072 … … 1167 1167 c detr_star(ig,l)=0. 1168 1168 endif 1169 1169 c print*,'ok detr_star' 1170 1170 cprise en compte du detrainement dans le calcul du flux 1171 1171 f_star(ig,l+1)=f_star(ig,l)+alim_star(ig,l)+entr_star(ig,l) … … 1188 1188 s -zw2(ig,l))/(zw2(ig,l+1)-zw2(ig,l)) 1189 1189 zw2(ig,l+1)=0. 1190 1190 c print*,'linter=',linter(ig) 1191 1191 else 1192 1192 wa_moy(ig,l+1)=sqrt(zw2(ig,l+1)) … … 1199 1199 enddo 1200 1200 enddo 1201 1201 c print*,'fin calcul zw2' 1202 1202 c 1203 1203 c Calcul de la couche correspondant a la hauteur du thermique … … 1230 1230 if (l.le.lmax(ig)) then 1231 1231 if (zw2(ig,l).lt.0.)then 1232 1232 c print*,'pb2 zw2<0' 1233 1233 endif 1234 1234 zw2(ig,l)=sqrt(zw2(ig,l))
Note: See TracChangeset
for help on using the changeset viewer.