var/mri4dvar/nc_vplocal.ncl

   1 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
   2 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
   3 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
   4
   5 begin
   6
   7
   8    filename = "vp_0020"
   9
  10    varnames = (/"psi","chi_u","t_u","rh","ps_u"/)
  11    nvar = dimsizes(varnames)
  12
  13    setfileoption("bin","ReadByteOrder","BigEndian")
  14    dims = fbinrecread(filename,0,15,"integer")
  15    ni=dims(12)
  16    nj=dims(13)
  17    nk=dims(14)
  18    print("ni, nj, nk = "+ni+", "+nj+", "+nk)
  19
  20    v = fbinrecread(filename,1,(/5,nk,nj,ni/),"double")
  21
  22    system("/bin/rm -f "+filename+".nc")   ; remove any pre-existing file
  23    ncdf = addfile(filename+".nc" ,"c")  ; open output netCDF file
  24
  25     ;===================================================================
  26     ; make time an UNLIMITED dimension; recommended  for most applications
  27     ;===================================================================
  28     ;   filedimdef(ncdf,"time",-1,True)
  29
  30    ;v1 = new((/nk, nj,ni/),double)
  31    ;v2 = new((/nk, nj,ni/),double)
  32    ;v3 = new((/nk, nj,ni/),double)
  33    v4 = new((/nk, nj,ni/),double)
  34    v5 = new((/nk, nj,ni/),double)
  35
  36    ;do k = 0, nk-1
  37    ;do j = 0, nj-1
  38    ;do i = 0, ni-1
  39    ;    v1(k,j,i) = v(i,j,k,1)
  40    ;    v2(k,j,i) = v(i,j,k,2)
  41    ;    v3(k,j,i) = v(i,j,k,3)
  42    ;    v4(k,j,i) = v(i,j,k,4)
  43    ;    ;v5(k,j,i) = v(i,j,k,5)
  44    ;end do
  45    ;end do
  46    ;end do
  47
  48    ;v1(:,:,:) = v(1,:,:,:)
  49    ;v2(:,:,:) = v(2,:,:,:)
  50    ;v3(:,:,:) = v(3,:,:,:)
  51    v4(:,:,:) = v(4,:,:,:)
  52    ;v5(:,:,:) = v(5,:,:,:)
  53
  54    ;ncdf->v1 = v1
  55    ;ncdf->v2 = v2
  56    ;ncdf->v3 = v3
  57    ncdf->v4 = v4
  58    ;ncdf->v5 = v5
  59
  60    ;************************************************
  61    ; end of reading be.dat data
  62    ;************************************************
  63
  64 end