Changeset 2189
- Timestamp:
- Jan 31, 2015, 10:31:49 AM (10 years ago)
- Location:
- LMDZ5/trunk/libf/phylmd
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/ozonecm_m.F90
r1955 r2189 91 91 ozonecm = max(ozonecm, 1e-12) 92 92 93 print*,'ozonecm Version2'93 ! print*,'ozonecm Version2' 94 94 95 95 END function ozonecm -
LMDZ5/trunk/libf/phylmd/pbl_surface_mod.F90
r2188 r2189 1700 1700 1701 1701 CASE(is_ter) 1702 print*,"DEBUGTS",yts(knon/2),ylwdown(knon/2)1702 ! print*,"DEBUGTS",yts(knon/2),ylwdown(knon/2) 1703 1703 CALL surf_land(itap, dtime, date0, jour, knon, ni,& 1704 1704 rlon, rlat, & -
LMDZ5/trunk/libf/phylmd/thermcell_plume.F90
r2159 r2189 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.