Problem with saving data without long-lag into NetCDF file

From: <Arnold.Sullivan_at_nyahnyahspammersnyahnyah>
Date: Tue Dec 22 2009 - 22:28:21 MST

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