| 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 |
| |