Subversion Repositories lagranto.ecmwf

Rev

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

Rev 29 Rev 31
Line 598... Line 598...
598 
      varname = dimname(3)
 598 
      varname = dimname(3)
599 
      ierr = NF90_INQ_VARID(fid,varname,varid)
 599 
      ierr = NF90_INQ_VARID(fid,varname,varid)
600 
      IF(ierr /= nf90_NoErr) PRINT *,NF90_STRERROR(ierr)
 600 
      IF(ierr /= nf90_NoErr) PRINT *,NF90_STRERROR(ierr)
601 
      ierr = nf90_get_var(fid,varid,lev)
 601 
      ierr = nf90_get_var(fid,varid,lev)
602 
      IF(ierr /= nf90_NoErr) PRINT *,NF90_STRERROR(ierr)
 602 
      IF(ierr /= nf90_NoErr) PRINT *,NF90_STRERROR(ierr)
- 
 
 603 
      varmin(1) = lon(1)
- 
 
 604 
      varmax(1) = lon( vardim(1) )
- 
 
 605 
      varmin(2) = lat(1)
- 
 
 606 
      varmax(2) = lat( vardim(2) )
603 
 
 607 
 
604 
c     Get ak and bk
 608 
c     Get ak and bk
605 
      varname='hyam'
 609 
      varname='hyam'
606 
      ierr = NF90_INQ_VARID(fid,varname,varid)
 610 
      ierr = NF90_INQ_VARID(fid,varname,varid)
607 
      IF(ierr /= nf90_NoErr) PRINT *,NF90_STRERROR(ierr)
 611 
      IF(ierr /= nf90_NoErr) PRINT *,NF90_STRERROR(ierr)
Line 675... Line 679...
675 
      if (abs(delta+dx).lt.eps) then
 679 
      if (abs(delta+dx).lt.eps) then
676 
          closear = 1
 680 
          closear = 1
677 
      else
 681 
      else
678 
          closear = 0
 682 
          closear = 0
679 
      endif
 683 
      endif
- 
 
 684 
      print*,'2: closear=',closear,varmax(1),varmin(1)
680 
 
 685 
 
681 
c     Save output array - close array and swap on demand
 686 
c     Save output array - close array and swap on demand
682 
      do j=1,vardim(2)
 687 
      do j=1,vardim(2)
683 
        do k=1,vardim(3)
 688 
        do k=1,vardim(3)
684 
          do i=1,vardim(1)
 689 
          do i=1,vardim(1)