Dear NCL users,
I had probelm with saving data without long-lag into NetCDF file, and hope someone could help me.
**** code ************************************************
 diri = "/workst/all/var/"  ; directory where files reside
 fili = systemfunc ("ls "+diri+"*_var.nc")
 fin  = addfiles (fili, "r")   ; note the "s" of addfile
 ListSetType (fin, "cat")               ; concatenate or "merge"
 Nmodel    = 24
 Nseason   = 4
 SQ_VAR    = new((/Nmodel,Nseason/),float)
 SQ_VAR(:,0)  = fin[:]->SQ_VAR_DJF
 SQ_VAR(:,1)  = fin[:]->SQ_VAR_MAM
 SQ_VAR(:,2)  = fin[:]->SQ_VAR_JJA
 SQ_VAR(:,3)  = fin[:]->SQ_VAR_SON
 filedimdef(fout,"time",-1,True)
 fout@title = "enso var from model files"
~~~~~
 fout->SQ_VAR = (/SQ_VAR/)              ; write SQ_VAR to a file
*** original data looks like *****************************
Variable: SQ_VAR_SON (file variable)
Type: float
Total Size: 96 bytes
            24 values
Number of Dimensions: 1
Dimensions and sizes:   [time | 24]
Coordinates:
            time [   1..   0]
Number Of Attributes: 4
  missing_value :       -1e+34
  _FillValue :  -1e+34
**** the output data looks like *****************************
show from FERRET
    1> ./merge_out.nc  (default)
 name     title                             I         J         K         L
 SQ_VAR                              1:24       1:4      ...       ...
Which is not what I want. I want data save as
 name     title                             I         J         K         L
 SQ_VAR                                ...       ...      1:4       1:24
How can I do it?
Cheers,
Arnold
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Dec 22 22:28:29 2009
This archive was generated by hypermail 2.1.8 : Tue Dec 29 2009 - 10:29:16 MST