[Thread Prev][Thread Next][Index]
[ferret_users] CRU ascii file to netcdf
- To: ferret_users@xxxxxxxx
- Subject: [ferret_users] CRU ascii file to netcdf
- From: Peter Szabo <szabpet83@xxxxxxxxx>
- Date: Thu, 13 Mar 2008 09:26:40 +0100
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:received:received:message-id:date:from:to:subject:mime-version:content-type; bh=3lv/vCs2Mlig3fv1c1kgeO13mgP5M5r7QeKCWXARW/k=; b=PSeKFFPUQzUgk8EgWzS3RyMN6zu79nOtsjoIOhVIsGEgb88WjREhel9k2Txt6snkFy7/a9eODN8PuNNc/zqIr76EleK86zTru8hCc1A8zNG3idLRLAvS0pNo6c8/+CZpY8W/5TgXI0dWE6/UzOjFm+829Aukfo34IaCM7fn44v0=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=message-id:date:from:to:subject:mime-version:content-type; b=qbOcB9OX+y45OtfD/CU7n2dFrnkBvNSE8GoIoPQ3x3DJskQWUr/gWap6vFTStlbX3lPQ0o4GxNVvP9Q/7KN5Sw0KVI6zzfBgaWIEUV0kmyzMNKe9i+LFQ/DBPcT6fxkj3axss6O7a/3ADNMn3OsD6HKKqsJl+mBJnpKniDKVV90=
- Sender: owner-ferret_users@xxxxxxxx
Dear Ferret Users,
Its me, again...
I have data set of CRU 10" in ASCII format like this:
4 1
-999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999
-999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999
-999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999
-999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999
-999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999
-999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999
-999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999
-999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999
-999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999
-999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999
4 2
-999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999
-999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999
-999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999
-999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999
-999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999
-999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999
-999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999
-999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999
-999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999
-999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 -999
The 1st, 12th, 23th... lines are the coordinate numbers of the data (the longitude goes from 4:258, the latitude goes from 1:228)
The "real" first row means a year 1961... The last one 1970. The columns mean months of the year. So 10*12 matrix i have after the grid defining numbers.
I would like to make a netcdf of it. Here is how i was trying to solve this problem.
--------------------------
define axis/t=1:120:1/np=120 time
define grid/t=time grid
repeat/range=1:58140:1/name=m (file/ez/var="lon,lat"/skip=`2+11*(m-1)`:`11+11*(m-1)` "myfile" ;\
file/var=pre/skip=`1+11*(m-1)`/columns=12/grid=grid "myfile" ;\
save/file="output.nc"/append/clobber/ilimits=1:258/jlimits=1:228/i=`lon`/j=`lat` pre)
----------------------------
Has anyone faced this problem?
Thank You,
Peter Szabo, HMS
[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Privacy Policy | Disclaimer | Accessibility Statement