Changeset 2220 for LMDZ5/branches/testing/libf/phylmd/thermcell_plume.F90
- Timestamp:
- Mar 3, 2015, 2:41:13 PM (10 years ago)
- Location:
- LMDZ5/branches/testing
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/branches/testing
- Property svn:mergeinfo changed
/LMDZ5/trunk merged: 2188-2195,2197-2202,2204-2210,2213-2216
- Property svn:mergeinfo changed
-
LMDZ5/branches/testing/libf/phylmd/thermcell_plume.F90
r2187 r2220 359 359 zbuoyjam(ig,l)=fact_shell*RG*(ztva_est(ig,l)-ztv_est(ig,l))/ztv_est(ig,l) & 360 360 & +(1.-fact_shell)*zbuoy(ig,l) 361 print*,'on est pass?? par l??1',l,lt,ztv1,ztv2,ztv(ig,lt),ztv_est(ig,l),ztva_est(ig,l),ztv(ig,l), &362 & zinv,zlmelup,zbuoy(ig,l),zbuoyjam(ig,l)363 361 elseif (zlmelup.ge.zinv) then 364 362 ztv_est2=atv2*0.5*(zlmelup+zinv)+btv2 … … 370 368 & ztv_est1)/ztv_est1)+(1.-fact_shell)*zbuoy(ig,l) 371 369 372 print*,'on est pass?? par l??2',l,lt,ztv_est1,ztv_est2,ztv(ig,lt),ztv_est(ig,l),ztva_est(ig,l),ztv(ig,l), &373 & zinv,zlmelup,zbuoy(ig,l),zbuoyjam(ig,l)374 370 else 375 371 ztv_est(ig,l)=atv1*zlmel+btv1 376 372 zbuoyjam(ig,l)=fact_shell*RG*(ztva_est(ig,l)-ztv_est(ig,l))/ztv_est(ig,l) & 377 373 & +(1.-fact_shell)*zbuoy(ig,l) 378 print*,'on est pass?? par l??3',l,lt,ztv1,ztv2,ztv(ig,lt),ztv_est(ig,l),ztva_est(ig,l),ztv(ig,l), &379 & zinv,zlmelup,zbuoy(ig,l),zbuoyjam(ig,l)380 374 endif 381 375 … … 392 386 endif 393 387 394 print*,'on est pass?? par l??4',l,lt,ztv1,ztv2,ztv(ig,lt),ztv(ig,l),ztva_est(ig,l), &395 & zlmelup,zbuoy(ig,l),zbuoyjam(ig,l)396 388 ! zbuoyjam(ig,l)=fact_shell*RG*(((zlmelup-zltdwn)/zdz)*(ztva_est(ig,l)- & 397 389 ! & ztv1)/ztv1+((zltdwn-zlmeldwn)/zdz)*(ztva_est(ig,l)- & … … 417 409 & ztv(ig,lt))/ztv(ig,lt)+((zdz2-lmel)/zdz3)*(ztva_est(ig,l)- & 418 410 & ztv(ig,lt-1))/ztv(ig,lt-1))+0.*zbuoy(ig,l) 419 print*,'on est pass?? par l??',l,lt,ztv(ig,lt),ztva_est(ig,l),ztv(ig,l), &420 & zbuoy(ig,l),zbuoyjam(ig,l)421 411 endif ! if (iflag_thermals_ed.lt.8) then 422 412
Note: See TracChangeset
for help on using the changeset viewer.