Changeset 437c545 in flexpart.git
- Timestamp:
- Mar 20, 2019, 3:32:30 PM (5 years ago)
- Branches:
- master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug
- Children:
- 7873bf7
- Parents:
- 77783e3
- Location:
- src
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
src/readwind_gfs.f90
rdb91eb7 r437c545 558 558 ! SEC & IP 12/2018 read GFS clouds 559 559 if(isec1(6).eq.153) then !! CLWCR Cloud liquid water content [kg/kg] 560 clwch(i,j,nlev_ec-k+2,n)=zsec4(nxfield*(ny-j-1)+i+1) 560 if((i.eq.0).and.(j.eq.0)) then 561 do ii=1,nuvz 562 if ((isec1(8)*100.0).eq.akz(ii)) numpu=ii 563 end do 564 endif 565 help=zsec4(nxfield*(ny-j-1)+i+1) 566 if(i.le.i180) then 567 clwch(i179+i,j,numpu,n)=help 568 else 569 clwch(i-i181,j,numpu,n)=help 570 endif 561 571 readclouds=.true. 562 572 sumclouds=.true. 573 ! readclouds=.false. 574 ! sumclouds=.false. 563 575 endif 564 576 -
src/verttransform_ecmwf.f90
r79e0349 r437c545 654 654 if ((lsp.gt.0.01).or.(convp.gt.0.01)) then ! cloud and precipitation 655 655 656 do kz=nz, 1,-1 !go Bottom up!656 do kz=nz,2,-1 !go Bottom up! 657 657 if (clw(ix,jy,kz,n).gt. 0) then ! is in cloud 658 658 cloudsh(ix,jy,n)=cloudsh(ix,jy,n)+height(kz)-height(kz-1) -
src/verttransform_gfs.f90
rdb91eb7 r437c545 548 548 if ((lsp.gt.0.01).or.(convp.gt.0.01)) then ! cloud and precipitation 549 549 550 do kz=nz, 1,-1 !go Bottom up!550 do kz=nz,2,-1 !go Bottom up! 551 551 if (clw(ix,jy,kz,n).gt. 0) then ! is in cloud 552 552 cloudsh(ix,jy,n)=cloudsh(ix,jy,n)+height(kz)-height(kz-1) 553 553 clouds(ix,jy,kz,n)=1 ! is a cloud 554 554 if (lsp.ge.convp) then 555 clouds(ix,jy,kz,n)=3 ! lsp in-cloud555 clouds(ix,jy,kz,n)=3 ! lsp in-cloud 556 556 else 557 557 clouds(ix,jy,kz,n)=2 ! convp in-cloud
Note: See TracChangeset
for help on using the changeset viewer.