Changeset 616 for trunk/LMDZ.MARS/libf
- Timestamp:
- Apr 5, 2012, 3:43:29 PM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90
r615 r616 754 754 do ig=1,ngridmx 755 755 if (zw2(ig,l+1)>0. .and. zw2(ig,l+1).lt.1.e-10) then 756 IF (lwrite) THEN 756 757 print*,'On tombe sur le cas particulier de thermcell_plume' 758 ENDIF 757 759 zw2(ig,l+1)=0. 758 760 linter(ig)=l+1 … … 1029 1031 ncorecfm2=ncorecfm2+1 1030 1032 else 1033 IF (lwrite) THEN 1031 1034 print*,'fm(l+1)<0 : ig, l+1,lmax :',ig,l+1,lmax(ig),fm(ig,l+1) 1035 ENDIF 1032 1036 ncorecfm1=ncorecfm1+1 1033 1037 fm(ig,l+1)=fm(ig,l) … … 1229 1233 & fm(ig,k)*ztv(ig,k)+fm(ig,k)*ztva(ig,k)-fm(ig,k+1)*ztva(ig,k+1)) 1230 1234 if (ztv(ig,k) + ptimestep*zdthladj(ig,k) .le. 0.) then 1235 IF (lwrite) THEN 1231 1236 print*,'Teta<0 in thermcell_dTeta up: qenv .. dq : ', ztv(ig,k),ptimestep*zdthladj(ig,k) 1237 ENDIF 1232 1238 if(ztv(ig,k) .gt. 0.) then 1233 1239 zdthladj(ig,k)=0. … … 1315 1321 & fm_down(ig,k)*ztv(ig,k)+fm_down(ig,k)*ztvd(ig,k)-fm_down(ig,k+1)*ztvd(ig,k+1)) 1316 1322 if (ztv(ig,k) + ptimestep*zdthladj_down(ig,k) .le. 0.) then 1323 IF (lwrite) THEN 1317 1324 print*,'q<0 in thermcell_dTeta down: qenv .. dq : ', ztv(ig,k),ptimestep*zdthladj_down(ig,k) 1325 ENDIF 1318 1326 if(ztv(ig,k) .gt. 0.) then 1319 1327 zdthladj(ig,k)=0.
Note: See TracChangeset
for help on using the changeset viewer.