Please someone help me... I waiting for answer for my last message also... looks like all are busy in userlist
Hi I got working on some data in which I want to add/remove zaxis for some purpose How this is possible
I tried to add zaxis from define axis and define grid command I could not succeed
Please help me., for first data set I want to add zaxis and for second dataset I want to remove z axis
Here is scenario
yes? use Quickscat_Wind_2006_2013.nc
yes? use Ascat_Wind_2007_13_1.nc
yes? sh d
currently SET data sets:
1> ./Quickscat_Wind_2006_2013.nc
name
title I J K L M N
WIND_SPEED
wind speed module 1:180 1:120 ... 1:454 ... ...
WIND_STRESS
wind stress magnitude
1:180 1:120 ... 1:454 ... ...
WIND_STRESS_CURL
wind stress curl 1:180 1:120 ... 1:454 ... ...
U zonal wind speed component 1:180 1:120 ... 1:454 ... ...
V meridional wind speed component 1:180
1:120 ... 1:454 ... ...
TAUX zonal wind stress component 1:180 1:120 ... 1:454 ... ...
TAUY meridional wind stress componen 1:180 1:120 ... 1:454 ... ...
2> ./Ascat_Wind_2007_13_1.nc (default)
name title
I J K L M N
WIND_SPEED
wind speed 1:480 1:240 1:1 1:2252 ... ...
WIND_STRESS
wind stress 1:480 1:240 1:1
1:2252 ... ...
WIND_STRESS_CURL
TAUY[X=@DDC]-TAUX[Y=@DDC] 1:480 1:240 1:1 1:2252 ... ...
U eastward wind speed 1:480 1:240 1:1 1:2252 ... ...
V northward wind speed
1:480 1:240 1:1 1:2252 ... ...
TAUX eastward wind stress 1:480 1:240 1:1 1:2252 ... ...
TAUY northward wind stress 1:480 1:240 1:1 1:2252 ... ...
yes? sh g u[d=1]
GRID
GBE1
name axis # pts start end
LONGITUDE421_600 LONGITUDE 180mr 30.25E 119.75E
LATITUDE101_220 LATITUDE 120 r 29.75S 29.75N
normal Z
TIME TIME 454 i 01-JAN-2006 12:00 31-MAR-2007 12:00
normal E
normal F
yes? sh g u[d=2]
GRID
GQF1
name axis # pts start end
LONGITUDE LONGITUDE 480mr 20.125E 139.88E
LATITUDE LATITUDE 240 r 29.875S 29.875N
DEPTH DEPTH (m) 1 r- 10
10
TIME1 TIME 2252 i 01-APR-2007 12:00 30-JUN-2013 12:00
normal E
normal F
yes?
ncdump -h Ascat_Wind_2007_13_1.nc
netcdf Ascat_Wind_2007_13_1 {
dimensions:
LONGITUDE = 480 ;
LATITUDE = 240 ;
DEPTH = 1 ;
TIME = UNLIMITED ; // (2252 currently)
bnds = 2 ;
variables:
double LONGITUDE(LONGITUDE) ;
LONGITUDE:axis = "X" ;
LONGITUDE:point_spacing = "even"
;
LONGITUDE:units = "degrees_east" ;
LONGITUDE:long_name = "longitude" ;
LONGITUDE:modulo = "yes" ;
LONGITUDE:standard_name = "longitude" ;
double LATITUDE(LATITUDE) ;
LATITUDE:axis = "Y" ;
LATITUDE:long_name = "latitude" ;
LATITUDE:modulo = "false" ;
LATITUDE:point_spacing = "even" ;
LATITUDE:units = "degrees_north" ;
LATITUDE:standard_name = "latitude" ;
double DEPTH(DEPTH) ;
DEPTH:axis = "Z" ;
DEPTH:point_spacing = "even"
;
DEPTH:units = "m" ;
DEPTH:long_name = "depth" ;
DEPTH:modulo = "false" ;
DEPTH:positive = "down" ;
DEPTH:standard_name = "depth" ;
double TIME(TIME) ;
TIME:axis = "T" ;
TIME:time_origin = "1900-01-01 00:00:00" ;
TIME:units = "hours since 1900-01-01 00:00:00" ;
TIME:bounds = "TIME_bnds" ;
TIME:long_name = "time" ;
TIME:modulo = "false" ;
TIME:calendar = "GREGORIAN" ;
TIME:standard_name = "time"
;
double TIME_bnds(TIME, bnds) ;
float WIND_SPEED(TIME, DEPTH, LATITUDE, LONGITUDE) ;
WIND_SPEED:_FillValue = -327.68f ;
WIND_SPEED:long_name = "wind speed" ;
WIND_SPEED:units = "m/s" ;
float WIND_STRESS(TIME, DEPTH, LATITUDE, LONGITUDE) ;
WIND_STRESS:_FillValue = -3.2768f ;
WIND_STRESS:long_name = "wind stress" ;
WIND_STRESS:units = "Pa" ;
double WIND_STRESS_CURL(TIME, DEPTH, LATITUDE, LONGITUDE) ;
WIND_STRESS_CURL:missing_value = -1.e+34 ;
WIND_STRESS_CURL:_FillValue = -1.e+34 ;
WIND_STRESS_CURL:long_name = "TAUY[X=@DDC]-TAUX[Y=@DDC]" ;
float U(TIME, DEPTH, LATITUDE, LONGITUDE) ;
U:missing_value = -1.e+34f ;
U:_FillValue = -1.e+34f ;
U:long_name = "eastward wind speed" ;
U:units = "m/s" ;
float V(TIME, DEPTH, LATITUDE, LONGITUDE) ;
V:missing_value = -1.e+34f ;
V:_FillValue = -1.e+34f ;
V:long_name = "northward wind speed" ;
V:units = "m/s" ;
float TAUX(TIME, DEPTH, LATITUDE, LONGITUDE) ;
TAUX:missing_value = -1.e+34f ;
TAUX:_FillValue = -1.e+34f
;
TAUX:long_name = "eastward wind stress" ;
TAUX:units = "Pa" ;
float TAUY(TIME, DEPTH, LATITUDE, LONGITUDE) ;
TAUY:missing_value = -1.e+34f ;
TAUY:_FillValue = -1.e+34f ;
TAUY:long_name = "northward wind stress" ;
TAUY:units = "Pa" ;
// global attributes:
:Conventions = "CF-1.0" ;
}
$ ncdump -h Quickscat_Wind_2006_2013.nc
netcdf Quickscat_Wind_2006_2013 {
dimensions:
LONGITUDE421_600 = 180 ;
LATITUDE101_220 = 120
;
TIME = UNLIMITED ; // (454 currently)
bnds = 2 ;
variables:
double LONGITUDE421_600(LONGITUDE421_600) ;
LONGITUDE421_600:axis = "X" ;
LONGITUDE421_600:point_spacing = "even" ;
LONGITUDE421_600:units = "degrees_east" ;
LONGITUDE421_600:long_name = "longitude" ;
LONGITUDE421_600:modulo = "yes" ;
LONGITUDE421_600:standard_name = "longitude" ;
double LATITUDE101_220(LATITUDE101_220) ;
LATITUDE101_220:axis = "Y" ;
LATITUDE101_220:long_name = "latitude" ;
LATITUDE101_220:modulo = "false"
;
LATITUDE101_220:point_spacing = "even" ;
LATITUDE101_220:units = "degrees_north" ;
LATITUDE101_220:standard_name = "latitude" ;
double TIME(TIME) ;
TIME:time_origin = "1900-01-01 00:00:00" ;
TIME:axis = "T" ;
TIME:units = "hours since 1900-01-01 00:00:00" ;
TIME:bounds = "TIME_bnds" ;
TIME:long_name = "time" ;
TIME:modulo = "false" ;
TIME:calendar = "GREGORIAN" ;
TIME:standard_name = "time" ;
double TIME_bnds(TIME, bnds) ;
float
WIND_SPEED(TIME, LATITUDE101_220, LONGITUDE421_600) ;
WIND_SPEED:_FillValue = 327.67f ;
WIND_SPEED:long_name = "wind speed module" ;
WIND_SPEED:units = "m/s" ;
float WIND_STRESS(TIME, LATITUDE101_220, LONGITUDE421_600) ;
WIND_STRESS:_FillValue = 32.767f ;
WIND_STRESS:long_name = "wind stress magnitude" ;
WIND_STRESS:units = "Pa" ;
float WIND_STRESS_CURL(TIME, LATITUDE101_220, LONGITUDE421_600) ;
WIND_STRESS_CURL:_FillValue = 3.2767e-05f ;
WIND_STRESS_CURL:long_name = "wind stress curl" ;
WIND_STRESS_CURL:units = "Pa/m" ;
float
U(TIME, LATITUDE101_220, LONGITUDE421_600) ;
U:missing_value = -1.e+34f ;
U:_FillValue = -1.e+34f ;
U:long_name = "zonal wind speed component" ;
U:units = "m/s" ;
float V(TIME, LATITUDE101_220, LONGITUDE421_600) ;
V:missing_value = -1.e+34f ;
V:_FillValue = -1.e+34f ;
V:long_name = "meridional wind speed component" ;
V:units = "m/s" ;
float TAUX(TIME, LATITUDE101_220, LONGITUDE421_600) ;
TAUX:missing_value = -1.e+34f ;
TAUX:_FillValue = -1.e+34f ;
TAUX:long_name =
"zonal wind stress component" ;
TAUX:units = "Pa" ;
float TAUY(TIME, LATITUDE101_220, LONGITUDE421_600) ;
TAUY:missing_value = -1.e+34f ;
TAUY:_FillValue = -1.e+34f ;
TAUY:long_name = "meridional wind stress component" ;
TAUY:units = "Pa" ;
// global attributes:
:Conventions = "CF-1.0" ;
}
- Dona Clara