Subversion Repositories lagranto.ecmwf

Rev

Rev 29 | Rev 36 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 29 Rev 31
Line 416... Line 416...
416 
         if ( leveltype.eq.'hybrid_sigma_pressure' ) then
 416 
         if ( leveltype.eq.'hybrid_sigma_pressure' ) then
417 
           do k=1,vardim(3)
 417 
           do k=1,vardim(3)
418 
              if ( vertical_swap.eq.1 ) then
 418 
              if ( vertical_swap.eq.1 ) then
419 
                 ak(k) = aktmp(vardim(3)-k+1)
 419 
                 ak(k) = aktmp(vardim(3)-k+1)
420 
                 bk(k) = bktmp(vardim(3)-k+1)
 420 
                 bk(k) = bktmp(vardim(3)-k+1)
- 
 
 421 
              else
- 
 
 422 
                 ak(k) = aktmp(k)
- 
 
 423 
                 bk(k) = bktmp(k)
421 
              endif
 424 
              endif
422 
           enddo
 425 
           enddo
423 
         elseif (leveltype.eq.'air_pressure' ) then
 426 
         elseif (leveltype.eq.'air_pressure' ) then
424 
           do k=1,vardim(3)
 427 
           do k=1,vardim(3)
425 
              if ( vertical_swap.eq.1 ) then
 428 
              if ( vertical_swap.eq.1 ) then
Line 598... Line 601...
598 
      varname = dimname(3)
 601 
      varname = dimname(3)
599 
      ierr = NF90_INQ_VARID(fid,varname,varid)
 602 
      ierr = NF90_INQ_VARID(fid,varname,varid)
600 
      IF(ierr /= nf90_NoErr) PRINT *,NF90_STRERROR(ierr)
 603 
      IF(ierr /= nf90_NoErr) PRINT *,NF90_STRERROR(ierr)
601 
      ierr = nf90_get_var(fid,varid,lev)
 604 
      ierr = nf90_get_var(fid,varid,lev)
602 
      IF(ierr /= nf90_NoErr) PRINT *,NF90_STRERROR(ierr)
 605 
      IF(ierr /= nf90_NoErr) PRINT *,NF90_STRERROR(ierr)
- 
 
 606 
      varmin(1) = lon(1)
- 
 
 607 
      varmax(1) = lon( vardim(1) )
- 
 
 608 
      varmin(2) = lat(1)
- 
 
 609 
      varmax(2) = lat( vardim(2) )
603 
 
 610 
 
604 
c     Get ak and bk
 611 
c     Get ak and bk
605 
      varname='hyam'
 612 
      varname='hyam'
606 
      ierr = NF90_INQ_VARID(fid,varname,varid)
 613 
      ierr = NF90_INQ_VARID(fid,varname,varid)
607 
      IF(ierr /= nf90_NoErr) PRINT *,NF90_STRERROR(ierr)
 614 
      IF(ierr /= nf90_NoErr) PRINT *,NF90_STRERROR(ierr)