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