Subversion Repositories lagranto.icon

Rev

Rev 3 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 3 Rev 10
Line 226... Line 226...
226 
c      Read units of vertical axis
 226 
c      Read units of vertical axis
227 
       read(10,*) umode
 227 
       read(10,*) umode
228 
       if ( ( umode.ne.'m'       ).and.
 228 
       if ( ( umode.ne.'m'       ).and.
229 
     >      ( umode.ne.'m,agl'   ).and.
 229 
     >      ( umode.ne.'m,agl'   ).and.
230 
     >      ( umode.ne.'hPa'     ).and.
 230 
     >      ( umode.ne.'hPa'     ).and.
- 
 
 231 
     >      ( umode.ne.'hPa,agl' ).and.
231 
     >      ( umode.ne.'K'       ).and.
 232 
     >      ( umode.ne.'K'       ).and.
232 
     >      ( umode.ne.'PVU'     ).and.
 233 
     >      ( umode.ne.'PVU'     ).and.
233 
     >      ( umode.ne.'INDEX'   )  )
 234 
     >      ( umode.ne.'INDEX'   )  )
234 
     > then
235 
     > then
235 
          print*,' ERROR: unit not supported ',trim(umode)
 236 
          print*,' ERROR: unit not supported ',trim(umode)
Line 564... Line 565...
564 
 
 565 
 
565 
      endif
 566 
      endif
566 
 
 567 
 
567 
c     ------ Pressure --------------------------------------------------
 568 
c     ------ Pressure --------------------------------------------------
568 
 
 569 
 
569 
      if ( umode.eq.'hPa' ) then
 570 
      if ( (umode.eq.'hPa').or.(umode.eq.'hPa,agl') ) then
570 
 
 571 
 
571 
c       Read pressure from first data file (pfile0) on U-grid
 572 
c       Read pressure from first data file (pfile0) on U-grid
572 
        call input_open (fid,pfile0)
 573 
        call input_open (fid,pfile0)
573 
        varname='P'                                      ! P
574 
        varname='P'                                      ! P
574 
        call input_wind
 575 
        call input_wind
Line 614... Line 615...
614 
           do j=1,ny
 615 
           do j=1,ny
615 
              prs(i,j) = pr(i,j,1)
 616 
              prs(i,j) = pr(i,j,1)
616 
           enddo
 617 
           enddo
617 
        enddo
 618 
        enddo
618 
 
 619 
 
- 
 
 620 
c       Convert the pressure: Pa -> hPa
- 
 
 621 
        do i=1,nx
- 
 
 622 
        	do j=1,ny
- 
 
 623 
        	   do k=1,nz
- 
 
 624 
        	      pr(i,j,k) = 0.01 * pr(i,j,k)
- 
 
 625 
        	   enddo
- 
 
 626 
        	   prs(i,j) = 0.01 * prs(i,j)
- 
 
 627 
        	enddo
- 
 
 628 
        enddo
- 
 
 629 
 
619 
      endif
 630 
      endif
620 
 
 631 
 
621 
c     ------ Potential temperature -------------------------------------
 632 
c     ------ Potential temperature -------------------------------------
622 
 
 633 
 
623 
       if ( (umode.eq.'K').or.(umode.eq.'PVU') ) then
 634 
       if ( (umode.eq.'K').or.(umode.eq.'PVU') ) then
Line 1171... Line 1182...
1171 
     >                start_lev(i),4,pr,prs,nx,ny,nz,xmin,ymin,dx,dy)
 1182 
     >                start_lev(i),4,pr,prs,nx,ny,nz,xmin,ymin,dx,dy)
1172 
            tmp1 = int_index3 (z,nx,ny,nz,rid,rjd,rkd,mdv)
 1183 
            tmp1 = int_index3 (z,nx,ny,nz,rid,rjd,rkd,mdv)
1173 
            start_hgt(i) = tmp1
 1184 
            start_hgt(i) = tmp1
1174 
         enddo
 1185 
         enddo
1175 
 
 1186 
 
- 
 
 1187 
c     Vertical mode <hPa,agl>
- 
 
 1188 
      elseif ( umode.eq.'hPa,agl' ) then
- 
 
 1189 
 
- 
 
 1190 
         do i=1,start_n
- 
 
 1191 
            call get_index3(rid,rjd,rkd,start_lon(i),start_lat(i),
- 
 
 1192 
     >                start_lev(i),4,pr,prs,nx,ny,nz,xmin,ymin,dx,dy)
- 
 
 1193 
            tmp1 = int_index3 (prs,nx,ny,1,rid,rjd,1.,mdv)
- 
 
 1194 
            start_lev(i) = tmp1 + start_lev(i)
- 
 
 1195 
            call get_index3(rid,rjd,rkd,start_lon(i),start_lat(i),
- 
 
 1196 
     >                start_lev(i),4,pr,prs,nx,ny,nz,xmin,ymin,dx,dy)
- 
 
 1197 
            tmp1 = int_index3 (z,nx,ny,nz,rid,rjd,rkd,mdv)
- 
 
 1198 
            start_hgt(i) = tmp1
- 
 
 1199 
         enddo
- 
 
 1200 
 
- 
 
 1201 
 
1176 
c     Vertical mode <K>
 1202 
c     Vertical mode <K>
1177 
      elseif ( umode.eq.'K' ) then
 1203 
      elseif ( umode.eq.'K' ) then
1178 
 
 1204 
 
1179 
         do i=1,start_n
 1205 
         do i=1,start_n
1180 
            call get_index3(rid,rjd,rkd,start_lon(i),start_lat(i),
 1206 
            call get_index3(rid,rjd,rkd,start_lon(i),start_lat(i),