Changes in src/readwind_gfs.f90 [db91eb7:5d74ed9] in flexpart.git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/readwind_gfs.f90
rdb91eb7 r5d74ed9 83 83 84 84 ! NCEP 85 integer :: numpt,numpu,numpv,numpw,numprh 85 integer :: numpt,numpu,numpv,numpw,numprh,numpclwch 86 86 real :: help, temp, ew 87 87 real :: elev … … 134 134 numpw=0 135 135 numprh=0 136 numpclwch=0 136 137 ifield=0 137 138 10 ifield=ifield+1 … … 557 558 endif 558 559 ! SEC & IP 12/2018 read GFS clouds 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((isec1(6).eq.153).and.(isec1(7).eq.100)) then !! CLWCR Cloud liquid water content [kg/kg] 561 if((i.eq.0).and.(j.eq.0)) then 562 do ii=1,nuvz 563 if ((isec1(8)*100.0).eq.akz(ii)) numpclwch=ii 564 end do 565 endif 566 help=zsec4(nxfield*(ny-j-1)+i+1) 567 if(i.le.i180) then 568 clwch(i179+i,j,numpclwch,n)=help 569 else 570 clwch(i-i181,j,numpclwch,n)=help 571 endif 561 572 readclouds=.true. 562 573 sumclouds=.true. 574 ! readclouds=.false. 575 ! sumclouds=.false. 563 576 endif 564 577
Note: See TracChangeset
for help on using the changeset viewer.