Subversion Repositories lagranto.um

Rev

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

Rev 7 Rev 11
Line 119... Line 119...
119 
      real         tmp(200)
 119 
      real         tmp(200)
120 
      character*80 varname
 120 
      character*80 varname
121 
      real         rtime
 121 
      real         rtime
122 
      integer      varid
 122 
      integer      varid
123 
      integer      cdfid
 123 
      integer      cdfid
- 
 
 124 
      real         rmax,rmin
124 
 
 125 
 
125 
c     Set file identifier
 126 
c     Set file identifier
126 
      if (fid.lt.0) then
 127 
      if (fid.lt.0) then
127 
        cdfid = -fid
 128 
        cdfid = -fid
128 
      else
129 
      else
Line 279... Line 280...
279 
      varname='PS'
 280 
      varname='PS'
280 
      ierr = NF90_INQ_VARID(cdfid,varname,varid)
 281 
      ierr = NF90_INQ_VARID(cdfid,varname,varid)
281 
      IF(ierr /= nf90_NoErr) PRINT *,NF90_STRERROR(ierr)
 282 
      IF(ierr /= nf90_NoErr) PRINT *,NF90_STRERROR(ierr)
282 
      ierr = nf90_get_var(cdfid,varid,ps)
 283 
      ierr = nf90_get_var(cdfid,varid,ps)
283 
      IF(ierr /= nf90_NoErr) PRINT *,NF90_STRERROR(ierr)
 284 
      IF(ierr /= nf90_NoErr) PRINT *,NF90_STRERROR(ierr)
- 
 
 285 
 
- 
 
 286 
c     Check that PS is in hPa - quick and dirty
- 
 
 287 
      rmax = -1.e9
- 
 
 288 
      rmin = +1.e9
- 
 
 289 
      do i=1,nx
- 
 
 290 
       do j=1,ny
- 
 
 291 
        if ( (ps(i,j).gt.rmax).and.abs(ps(i,j)-mdv).gt.eps)  then
- 
 
 292 
          rmax = ps(i,j)
- 
 
 293 
        endif
- 
 
 294 
        if ( (ps(i,j).lt.rmin).and.abs(ps(i,j)-mdv).gt.eps)  then
- 
 
 295 
          rmin = ps(i,j)
- 
 
 296 
        endif
- 
 
 297 
       enddo
- 
 
 298 
      enddo
- 
 
 299 
      if ( (rmin.lt.200.).or.(rmax.gt.1500.) ) then
- 
 
 300 
        print*,' ERROR: PS must be in hPa... Stop'
- 
 
 301 
        stop
- 
 
 302 
      endif
284
303
285 
c     Calculate layer and level pressures
 304 
c     Calculate layer and level pressures
286 
      do i=1,nx
 305 
      do i=1,nx
287 
         do j=1,ny
 306 
         do j=1,ny
288 
               do k=1,nz
 307 
               do k=1,nz