228,6 → 228,7 |
if ( ( umode.ne.'m' ).and. |
> ( umode.ne.'m,agl' ).and. |
> ( umode.ne.'hPa' ).and. |
> ( umode.ne.'hPa,agl' ).and. |
> ( umode.ne.'K' ).and. |
> ( umode.ne.'PVU' ).and. |
> ( umode.ne.'INDEX' ) ) |
566,7 → 567,7 |
|
c ------ Pressure -------------------------------------------------- |
|
if ( umode.eq.'hPa' ) then |
if ( (umode.eq.'hPa').or.(umode.eq.'hPa,agl') ) then |
|
c Read pressure from first data file (pfile0) on U-grid |
call input_open (fid,pfile0) |
616,6 → 617,16 |
enddo |
enddo |
|
c Convert the pressure: Pa -> hPa |
do i=1,nx |
do j=1,ny |
do k=1,nz |
pr(i,j,k) = 0.01 * pr(i,j,k) |
enddo |
prs(i,j) = 0.01 * prs(i,j) |
enddo |
enddo |
|
endif |
|
c ------ Potential temperature ------------------------------------- |
1173,6 → 1184,21 |
start_hgt(i) = tmp1 |
enddo |
|
c Vertical mode <hPa,agl> |
elseif ( umode.eq.'hPa,agl' ) then |
|
do i=1,start_n |
call get_index3(rid,rjd,rkd,start_lon(i),start_lat(i), |
> start_lev(i),4,pr,prs,nx,ny,nz,xmin,ymin,dx,dy) |
tmp1 = int_index3 (prs,nx,ny,1,rid,rjd,1.,mdv) |
start_lev(i) = tmp1 + start_lev(i) |
call get_index3(rid,rjd,rkd,start_lon(i),start_lat(i), |
> start_lev(i),4,pr,prs,nx,ny,nz,xmin,ymin,dx,dy) |
tmp1 = int_index3 (z,nx,ny,nz,rid,rjd,rkd,mdv) |
start_hgt(i) = tmp1 |
enddo |
|
|
c Vertical mode <K> |
elseif ( umode.eq.'K' ) then |
|