Re: addfile() for create new nc file

On 10 August 2011 13:59, gibies george <gibies@tropmet.res.in> wrote:

> Thank you Dennis and Mary for your support and suggestions.
>
> My problem is solved by doing a slight variation in my script according to
> your suggestions. Some how I had to retain the string information and so I
> made them to be attributes to the data variable.
>

But again there is a problem in writing the attributes. Here is the output
of ncdump

netcdf ond_ssn_multimodel_data {
dimensions:
time = 25 ;
lat = 180 ;
lon = 360 ;
season = 1 ;
component = 6 ;
ens = 1 ;
mdl = 9 ;
variables:
float multi_model_data(time, lat, lon, season, component, ens, mdl) ;
multi_model_data:unit =
"m/s\000\000\000\000m/s\000\000\000\000m/s\000\000\000\000mm/day\000hPa\000\000\000\000degC"
;
multi_model_data:units =
"m/s\000\000\000\000m/s\000\000\000\000m/s\000\000\000\000mm/day\000hPa\000\000\000\000degC"
;
multi_model_data:region = "global" ;
multi_model_data:ssnend = 12 ;
multi_model_data:ssnbegin = 10 ;
multi_model_data:dataset = "ond_ssn_1981_to_2005" ;
multi_model_data:datavar = "wind850_rain_mslp_tsfc" ;
multi_model_data:anal_start_year = 1981 ;
multi_model_data:anal_end_year = 2005 ;
multi_model_data:latsouth = -90 ;
multi_model_data:latnorth = 90 ;
multi_model_data:lonwest = 0 ;
multi_model_data:loneast = 360 ;
multi_model_data:mpCenterLonF = 180 ;
multi_model_data:novctr = 1 ;
multi_model_data:nosclr = 3 ;
multi_model_data:long_name = "wind850_rain_mslp_tsfc" ;
multi_model_data:ens = "ond_ssn" ;
multi_model_data:mdl =
"ond_ssn\000\000\000\000\000\000\000\000\000\000\000\000\000\000\000\000cfsv2_jul_init\000\000\000\000\000\000\000\000\000cfsv2_aug_init\000\000\000\000\000\000\000\000\000ECMWF_aug_init\000\000\000\000\000\000\000\000\000IFM_Geomer_aug_init\000\000\000\000Meteo_France_aug_init\000\000UK_MetOffice_aug_init\000\000CMCC_INGV_aug_init\000\000\000\000\000DePreSys_UKMO_aug_init"
;
multi_model_data:component =
"wind850_u\000\000\000\000\000wind850_v\000\000\000\000\000wind850_speed\000rain\000\000\000\000\000\000\000\000\000\000mslp\000\000\000\000\000\000\000\000\000\000tsfc"
;
multi_model_data:season = "ond" ;
multi_model_data:_FillValue = -9.96921e+36f ;
int time(time) ;
time:_FillValue = -999 ;
float lat(lat) ;
float lon(lon) ;
int season(season) ;
season:ssnbegin = 10 ;
season:ssnend = 12 ;
int component(component) ;
int ens(ens) ;
int mdl(mdl) ;

Can anybody suggest a solution.

> But I still feel it is better to have more direct way to do all these
> things. Please inform me if there is some simpler way.
>
> Thanking you once again
>
>
> *Here is my modified script.*
>
> ________________________________________________________________________________________________
>
> multidata =
> new((/dimsizes(time),dimsizes(lat),dimsizes(lon),dimsizes(season),dimsizes(component),\
> nofens,nofmdls/),float)
> multidata!0 = "time"
> multidata&time = time
> multidata!1 = "lat"
> multidata&lat = lat
> multidata!2 = "lon"
> multidata&lon = lon
> multidata!3 ="season"
> multidata&season=ind(.not.ismissing(season))
> multidata@season = season
> multidata!4 ="component"
> multidata&component=ind(.not.ismissing(component))
> multidata@component = component
> multidata!5 ="ens"
> multidata&ens = ispan(1,nofens,1)
> multidata@ens = ens
> multidata!6 ="mdl"
> multidata&mdl = ispan(1,nofmdls,1)
> multidata@mdl = mdl
> -------------------
> -----------
> ----------------------
> ----------
> multidata@dataset = mdl@analdset + "_" + predef@anal_start_year + "_to_" +
> predef@anal_end_year
> multidata@datavar = datavar
> multidata&ens = ind(.not.ismissing(ens))
> multidata&mdl = ind(.not.ismissing(mdl))
> print("---------------Multimodel Ensemble-------------------")
> print("_____________________________________________________")
> printVarSummary(multidata)
> print("_____________________________________________________")
>
> syout= systemfunc("rm " + mdl@analdset+"_multimodel_data.nc")
> outfile = addfile(mdl@analdset+"_multimodel_data.nc","c")
> outfile->multi_model_data = multidata
> print("Data is written. Press any key to continue.")
> response = systemfunc("read var; echo $var")
>
> ___________________________________________________________________________________________
>
> *
> *
> *Here is the output.*
>
> Variable: multidata
> Type: float
> Total Size: 12960000 bytes
> 3240000 values
> Number of Dimensions: 7
> Dimensions and sizes: [time | 25] x [lat | 180] x [lon | 360] x [season |
> 1] x [component | 1] x [ens | 2] x [mdl | 1]
> Coordinates:
> time: [1981..2005]
> lat: [-90..90]
> lon: [ 0..360]
> season: [0..0]
> component: [0..0]
> ens: [0..1]
> mdl: [0..0]
> Number Of Attributes: 22
> unit : mm/day
> units : mm/day
> long_name : _rain
> nosclr : 1
> novctr : 0
> mpCenterLonF : 180
> loneast : 360
> lonwest : 0
> latnorth : 90
> latsouth : -90
> anal_end_year : 2005
> anal_start_year : 1981
> datavar : _rain
> dataset : cfsv2_jjas_1981_to_2005
> ssnbegin : 6
> ssnend : 9
> region : global
> mdl : cfsv2
> ens : ( feb_initial , may_initial )
> component : rain
> season : jjas
> _FillValue : -9.96921e+36
> (0) _____________________________________________________
> (0) Data is written. Press any key to continue.
>
>
> ____________________________________________________________________________________________
>
>
>
>
> On 6 June 2011 21:53, Mary Haley <haley@ucar.edu> wrote:
>
>>
>> Hi,
>>
>> NetCDF doesn't allow you to have strings for coordinate arrays. You must
>> convert these to 2D character arrays.
>>
>> Please see the "Coordinate Systems" section of this NetCDF document:
>>
>> http://www.unidata.ucar.edu/software/netcdf/docs/BestPractices.html
>>
>> I don't know how you created "season" in your code below, but you can try
>> converting it to a character array, and use this as your coordinate
>> variable:
>>
>> cseason = tocharacter(season)
>> multidata&season = cseason
>>
>> --Mary
>>
>> On Jun 6, 2011, at 1:11 AM, gibies george wrote:
>>
>> Dear NCL users,
>>
>> Somebody please suggest a way to eliminate the following warning messages.
>>
>>
>> I was trying to write my dataset to an nc file.
>>
>> Thanking you in davance.
>>
>> *_______________________________________
>> ________________Script__________________
>> _______________________________________*
>> -------------------------
>> -----------------------
>> ----------------------------
>> ---------------------
>> --------------------------
>>
>> multidata@dataset = mdl@analdset + "_" + predef@anal_start_year +
>> "_to_" + predef@anal_end_year
>> multidata@datavar = datavar
>> multidata&ens = ens
>> multidata&mdl = mdl
>> print("---------------Multimodel Ensemble-------------------")
>> print("_____________________________________________________")
>> printVarSummary(multidata)
>> print("_____________________________________________________")
>>
>> syout= systemfunc("rm " + mdl@analdset+"_multimodel_data.nc")
>> outfile = addfile(mdl@analdset+"_multimodel_data.nc","c")
>> outfile->multi_model_data = multidata
>>
>> return(multidata)
>> end ;;; End of function definition get_multimodel_data() ;;;
>>
>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>> *_______________________________________*
>>
>> *_________________________________
>> ___________ Output _______________
>> _________________________________*
>>
>>
>> (0) ---------------Multimodel Ensemble-------------------
>> (0) _____________________________________________________
>>
>>
>> Variable: multidata
>> Type: float
>> Total Size: 272160000 bytes
>> 68040000 values
>> Number of Dimensions: 7
>> Dimensions and sizes: [time | 25] x [lat | 280] x [lon | 360] x [season
>> | 1] x [component | 3] x [ens | 1] x [mdl | 9]
>> Coordinates:
>> time: [1981..2005]
>> lat: [-35..35]
>> lon: [30..120]
>> season: [ond..ond]
>> component: [rain..tsfc]
>> ens: [..]
>> mdl: [cfsv2_jun_init..DePreSys_UKMO_aug_init]
>> Number Of Attributes: 18
>> unit : ( mm/day, hPa, degC )
>> units : ( mm/day, hPa, degC )
>> region : indotropic
>> ssnend : 12
>> ssnbegin : 10
>> dataset : ond_skill_1981_to_2005
>> datavar : _rain_mslp_tsfc
>> anal_start_year : 1981
>> anal_end_year : 2005
>> latsouth : -35
>> latnorth : 35
>> lonwest : 30
>> loneast : 120
>> mpCenterLonF : 75
>> novctr : 0
>> nosclr : 3
>> long_name : _rain_mslp_tsfc
>> _FillValue : -9.96921e+36
>> (0) _____________________________________________________
>> fatal:Attempting to write variable (season) of type (string) which is not
>> representable in the format of file (ond_skill_multimodel_data.nc)
>> fatal:Trying using a type conversion function
>> warning:FileWriteVarVar: Could not write coordinate variable (season) to
>> file (ond_skill_multimodel_data), continuing anyway
>> fatal:Attempting to write variable (component) of type (string) which is
>> not representable in the format of file (ond_skill_multimodel_data.nc)
>> fatal:Trying using a type conversion function
>> warning:FileWriteVarVar: Could not write coordinate variable (component)
>> to file (ond_skill_multimodel_data), continuing anyway
>> fatal:Attempting to write variable (ens) of type (string) which is not
>> representable in the format of file (ond_skill_multimodel_data.nc)
>> fatal:Trying using a type conversion function
>> warning:FileWriteVarVar: Could not write coordinate variable (ens) to file
>> (ond_skill_multimodel_data), continuing anyway
>> fatal:Attempting to write variable (mdl) of type (string) which is not
>> representable in the format of file (ond_skill_multimodel_data.nc)
>> fatal:Trying using a type conversion function
>> warning:FileWriteVarVar: Could not write coordinate variable (mdl) to file
>> (ond_skill_multimodel_data), continuing anyway
>> warning:Execute: Error occurred at or near line 3188 in file
>> /home/gibies/ncl/mylib/lib_gibies_read_data_module.ncl
>>
>>
>>
>> --
>> *Gibies George, CSIR-RF,*
>> *Monsoon Seasonal Prediction Group,
>> Climate and Global Modelling Division,
>> Indian Institute of Tropical Meteorology,
>> Dr. Homi Bhabha Road,
>> NCL (P. O.), Pashan,
>> Pune 411008, India.*
>>
>> *http://sites.google.com/site/gibiesge/*
>>
>> Please Think about the environment. Save paper; Save Trees; and don't
>> print this e-mail unless it is necessary.
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>>
>
>
> --
> *Gibies George, CSIR-RF,*
> *Monsoon Seasonal Prediction Group,
> Climate and Global Modelling Division,
> Indian Institute of Tropical Meteorology,
> Dr. Homi Bhabha Road,
> NCL (P. O.), Pashan,
> Pune 411008, India.*
>
> *http://sites.google.com/site/gibiesge/*
>
> Please Think about the environment. Save paper; Save Trees; and don't print
> this e-mail unless it is necessary.
>

-- 
*Gibies George,   CSIR-RF,*
*Monsoon Seasonal Prediction Group,
Climate and Global Modelling Division,
Indian Institute of Tropical Meteorology,
Dr. Homi Bhabha Road,
NCL (P. O.), Pashan,
Pune 411008, India.*
*http://sites.google.com/site/gibiesge/*
Please Think about the environment. Save paper; Save Trees; and don't print
this e-mail unless it is necessary.