grads-1.9b4 & gds-1.3

Thomas LOUBRIEU Thomas.Loubrieu at IFREMER.FR
Thu Jun 23 05:22:14 EDT 2005


Hi Jennifer,

Excuse me I forgot to thank you for your investigation in my
configuration files. It works now.
I have now to move on the ferret users list because of problems related
to the _FillValue management.

Thanks again,

Thomas

Jennifer Adams wrote:

> On Jun 10, 2005, at 12:47 PM, Thomas LOUBRIEU wrote:
>
>> Hi Jennifer,
>>
>> What you have foreseen is right. The first case fails while the second
>> one succeed.
>> I've tried other cases and it seems that :
>> 1) using sdfopen on one of my netcdf file : it always works
>> 2) using open on my ctl file :
>>    a) it works if the Y subset is below or equal 421 (for example 'set
>> y 40 412') and for any other subset on the other dimensions (x, t, z).
>>    b) it fails as soon as the Y subset contains a dimension above 421.
>>
>> The 421 value corresponds to the longitude (which is X ?) dimension
>> size
>> of my dataset.
>>
>> Does this help you to find out what happens ?
>
> Yes, it does.  Note in the output from ncdump:
>
>>         short temperature(time, depth, latitude, longitude) ;
>>         short pct_variance(time, depth, latitude, longitude) ;
>
> This tells you the "shape" of the variable, and when you list the
> varying dimensions in your descriptor file they need to be in the
> same order that they are given here. You need to change the variable
> declaration in your descriptor file from "t,z,x,y" to "t,z,y,x" .
> That should fix this problem.
>
> Jennifer
>
>
>
>>
>> Thank you again,
>>
>> Thomas
>>
>> PS. the ncdump -c output is in the message.
>>
>>
>>
>> Jennifer Adams wrote:
>>
>>
>>> OK, let's simplify. If I have read everything correctly, this should
>>> give you the error:
>>>
>>> open /home/tloubrie/dods/dodsG/OA_V302/coriolis-atlantic-temperature-
>>> weekly-oa_v3_rt.ctl
>>> set x 152 267
>>> set y 260 426
>>> set z 7 7
>>> set t 242 242
>>> set gxout fwrite
>>> set fwex
>>> d temperature
>>>
>>> While this should work:
>>>
>>> open /home/tloubrie/dods/dodsG/OA_V302/coriolis-atlantic-temperature-
>>> weekly-oa_v3_rt.ctl
>>> set x 101 301
>>> set y 260 401
>>> set z 7 7
>>> set t 242 242
>>> set gxout fwrite
>>> set fwex
>>> d temperature
>>>
>>> If that is the case, what is the output from ncdump -c for the data
>>> file that goes with time step 242?
>>> Jennifer
>>>
>>>
>>> On Jun 9, 2005, at 8:21 AM, Thomas LOUBRIEU wrote:
>>>
>>>
>>>> Hi Jennifer,
>>>>
>>>> Thanks for your quick reply.
>>>> Doing what you asked for, I've seen that there is a netcdf index
>>>> query
>>>> failure, but I do not understand where.
>>>> The elements  are :
>>>> 1) For the failing request :
>>>> GrADS session :
>>>>
>>>> ga-> /export/home/logiciels/gds/scripts/subset.gs /export/home/
>>>> logiciels/jakarta-tomcat-4.1.30/temp/gds/grads/dods/subset/coriolis-
>>>> atlantic-temperature-weekly-oa_v3_rt.temperature59646.subset /home/
>>>> tloubrie/dods/dodsG/OA_V302/coriolis-atlantic-temperature-weekly-
>>>> oa_v3_rt.ctl temperature 152 267 260 426 7 7 242 242
>>>>
>>>> output  = /export/home/logiciels/jakarta-tomcat-4.1.30/temp/gds/
>>>> grads/dods/subset/coriolis-atlantic-temperature-weekly-
>>>> oa_v3_rt.temperature59646.subset
>>>>
>>>> dataset = /home/tloubrie/dods/dodsG/OA_V302/coriolis-atlantic-
>>>> temperature-weekly-oa_v3_rt.ctl
>>>>
>>>> varname = temperature
>>>>
>>>> 152 -> 267 Varying
>>>>
>>>> 260 -> 426 Varying
>>>>
>>>> 7 -> 7
>>>>
>>>> 242 -> 242
>>>>
>>>> 2 Varying dimensions
>>>>
>>>>  Xsize = 421  Ysize = 597  Zsize = 59  Tsize = 255
>>>>
>>>> varname=temperature
>>>>
>>>> Output from q dims:
>>>>
>>>> Default file number is: 1
>>>>
>>>> X is varying   Lon = -100 to 40   X = 1 to 421
>>>>
>>>> Y is varying   Lat = -69.9662 to 69.9662   Y = 1 to 597
>>>>
>>>> Z is fixed     Lev = 5  Z = 1
>>>>
>>>> T is fixed     Time = 00Z05JAN2000  T = 1
>>>>
>>>>
>>>>
>>>> Output from q fwrite:
>>>>
>>>> FWrite file is closed
>>>>
>>>> FWrite file name is: /export/home/logiciels/jakarta-tomcat-4.1.30/
>>>> temp/gds/grads/dods/subset/coriolis-atlantic-temperature-weekly-
>>>> oa_v3_rt.temperature59646.subset
>>>>
>>>> FWrite byte order is big_endian ; machine byte order is
>>>> little_endian
>>>>
>>>>
>>>>
>>>> NetCDF Error (nc_get_vara_float): Index exceeds dimension bound
>>>>
>>>> Data Request Error:  Error for variable 'temperature'
>>>>
>>>>  Error ocurred at column 1
>>>>
>>>> DISPLAY error:  Invalid expression
>>>>
>>>>  Expression = temperature
>>>>
>>>> NetCDF Error (nc_get_vara_float): Index exceeds dimension bound
>>>>
>>>> Data Request Error:  Error for variable 'temperature'
>>>>
>>>>  Error ocurred at column 1
>>>>
>>>> DISPLAY error:  Invalid expression
>>>>
>>>>  Expression = temperature
>>>>
>>>>
>>>>
>>>> No hardcopy metafile open
>>>>
>>>> GX package terminated
>>>>
>>>>
>>>> no output file is created.
>>>>
>>>> 2) For a succesfull resquest :
>>>> GraDS session :
>>>>
>>>> ga-> /export/home/logiciels/gds/scripts/subset.gs /export/home/
>>>> logiciels/jakarta-tomcat-4.1.30/temp/gds/grads/dods/subset/coriolis-
>>>> atlantic-temperature-weekly-oa_v3_rt.temperature59648.subset /home/
>>>> tloubrie/dods/dodsG/OA_V302/coriolis-atlantic-temperature-weekly-
>>>> oa_v3_rt.ctl temperature 101 301 260 401 7 7 242 242
>>>>
>>>> output  = /export/home/logiciels/jakarta-tomcat-4.1.30/temp/gds/
>>>> grads/dods/subset/coriolis-atlantic-temperature-weekly-
>>>> oa_v3_rt.temperature59648.subset
>>>>
>>>> dataset = /home/tloubrie/dods/dodsG/OA_V302/coriolis-atlantic-
>>>> temperature-weekly-oa_v3_rt.ctl
>>>>
>>>> varname = temperature
>>>>
>>>> 101 -> 301 Varying
>>>>
>>>> 260 -> 401 Varying
>>>>
>>>> 7 -> 7
>>>>
>>>> 242 -> 242
>>>>
>>>> 2 Varying dimensions
>>>>
>>>>  Xsize = 421  Ysize = 597  Zsize = 59  Tsize = 255
>>>>
>>>> varname=temperature
>>>>
>>>> Output from q dims:
>>>>
>>>> Default file number is: 1
>>>>
>>>> X is varying   Lon = -100 to 40   X = 1 to 421
>>>>
>>>> Y is varying   Lat = -69.9662 to 69.9662   Y = 1 to 597
>>>>
>>>> Z is fixed     Lev = 5  Z = 1
>>>>
>>>> T is fixed     Time = 00Z05JAN2000  T = 1
>>>>
>>>>
>>>>
>>>> Output from q fwrite:
>>>>
>>>> FWrite file is closed
>>>>
>>>> FWrite file name is: /export/home/logiciels/jakarta-tomcat-4.1.30/
>>>> temp/gds/grads/dods/subset/coriolis-atlantic-temperature-weekly-
>>>> oa_v3_rt.temperature59648.subset
>>>>
>>>> FWrite byte order is big_endian ; machine byte order is
>>>> little_endian
>>>>
>>>>
>>>>
>>>> Wrote 28542 of 28542 elements to /export/home/logiciels/jakarta-
>>>> tomcat-4.1.30/temp/gds/grads/dods/subset/coriolis-atlantic-
>>>> temperature-weekly-oa_v3_rt.temperature59648.subset as Stream
>>>> Big_Endian
>>>>
>>>>
>>>>
>>>> No hardcopy metafile open
>>>>
>>>> GX package terminated
>>>>
>>>>
>>>> The size of the output file is : 114168 octets
>>>>
>>>> Thanks for your help,
>>>>
>>>> Thomas
>>>>
>>>>
>>>> Jennifer Adams wrote:
>>>>
>>>>
>>>>
>>>>> Hi, Thomas --
>>>>> Thanks for providing so much helpful information. At the moment,
>>>>> everything you've done looks right. The "ran out of input while
>>>>> spooling..." message means that the size of the subset that got
>>>>> written out did not match the requested grid. To find out what
>>>>> happened, you have to put your server in debug mode <log
>>>>> level=debug>
>>>>> and then make the request. You'll see in the log file a long
>>>>> command
>>>>> (with "nice" in front of it) that invokes GrADS, runs the subset.gs
>>>>> script and outputs to a file under the $GDSHOME/temp directory. If
>>>>> you could please run this command manually, outside the GDS, and
>>>>> send
>>>>> along the complete text of the GrADS session as well as the size of
>>>>> the output file it created. Please do the same with the request
>>>>> that
>>>>> works too.
>>>>> Jennifer
>>>>>
>>>>>
>>>>> On Jun 9, 2005, at 4:25 AM, Thomas LOUBRIEU wrote:
>>>>>
>>>>>
>>>>>
>>>>>> Hello,
>>>>>>
>>>>>> I am now testing the new version of gds and grads with my netcdf
>>>>>> files
>>>>>> and I have some problems.
>>>>>>
>>>>>> I got various errors depending on my dods request :
>>>>>> I cannot really find out why some request fail and other work.
>>>>>> The gds.log shows :
>>>>>>
>>>>>> Jun 09 09:11:07 error: [ http8080-Processor4 br144-162 GET
>>>>>> /coriolis-atlantic-temperature-weekly-oa_v3_rt.dods?temperature
>>>>>> [241:1:241][6:1:6][259:1:425][150:1:267]
>>>>>> ]  io error during data send: class java.io.IOException: ran
>>>>>> out of
>>>>>> input while spooling 78824 bytes
>>>>>>
>>>>>> or
>>>>>>
>>>>>> Jun 09 10:02:00 error: [ http8080-Processor3 br144-162 GET
>>>>>> /coriolis-atlantic-temperature-weekly-oa_v3_rt.asc?temperature
>>>>>> [241:1:241][6:1:6][259:1:425][151:1:267]
>>>>>> ]  io error during data send: class java.io.EOFException: null
>>>>>>
>>>>>> while, this one works :
>>>>>>
>>>>>> Jun 09 09:59:49 [ http8080-Processor2 br144-162 GET
>>>>>> /coriolis-atlantic-temperature-weekly-oa_v3_rt.dods?temperature
>>>>>> [241:1:241][6:1:6][259:1:400][100:1:300]
>>>>>> ] FINISHED (success in 114ms)
>>>>>>
>>>>>> Unfortunatly, my test server is not public.
>>>>>>
>>>>>> I uses :
>>>>>> grads-1.9b4
>>>>>> gds-1.3
>>>>>> java version "1.4.2"
>>>>>> Tomcat/4.1.30
>>>>>> on linux red hat  9
>>>>>>
>>>>>> The gds.xml file contains :
>>>>>>
>>>>>>     <dataset name = "coriolis-atlantic-temperature-weekly-
>>>>>> oa_v3_rt"
>>>>>>
>>>>>>        file="/home/tloubrie/dods/dodsG/OA_V302/coriolis-atlantic-
>>>>>> temperature-weekly-oa_v3_rt.ctl"
>>>>>>
>>>>>>        format="nc">
>>>>>>
>>>>>>      </dataset>
>>>>>>
>>>>>>
>>>>>> I have configured my dataset that way (ctl file) :
>>>>>>
>>>>>> DSET  /home/coriolis_exp/spool/co04/co0401/co040102/oa/field/TEMP/
>>>>>> OA_%y4%m2%d2.nc
>>>>>>
>>>>>> DTYPE netcdf
>>>>>>
>>>>>> TITLE CORIOLIS Analysed temperature
>>>>>>
>>>>>> UNDEF 32767
>>>>>>
>>>>>> UNPACK scale_factor add_offset
>>>>>>
>>>>>> OPTIONS template
>>>>>>
>>>>>> XDEF 421 LEVELS
>>>>>>
>>>>>> -100 -99.66666 -99.33334 -99 -98.66666 -98.33334 -98
>>>>>>
>>>>>>    -97.66666 -97.33334 -97 -96.66666 -96.33334 -96 -95.66666
>>>>>>
>>>>>>    -95.33334 -95 -94.66666 -94.33334 -94 -93.66666 -93.33334 -93
>>>>>>
>>>>>>    -92.66666 -92.33334 -92 -91.66666 -91.33334 -91 -90.66666
>>>>>>
>>>>>>    -90.33334 -90 -89.66666 -89.33334 -89 -88.66666 -88.33334 -88
>>>>>>
>>>>>> ....
>>>>>>
>>>>>>  36.33333 36.66667 37 37.33333 37.66667 38 38.33333 38.66667 39
>>>>>>
>>>>>>    39.33333 39.66667 40
>>>>>>
>>>>>> YDEF 597 LEVELS
>>>>>>
>>>>>> -69.96624 -69.85174 -69.7366 -69.62084 -69.50445 -69.38742
>>>>>>
>>>>>>    -69.26975 -69.15144 -69.03249 -68.91288 -68.79263 -68.67171
>>>>>>
>>>>>> ...
>>>>>>
>>>>>>   68.42793 68.55015 68.67171 68.79263 68.91288 69.03249 69.15144
>>>>>>
>>>>>>    69.26975 69.38742 69.50445 69.62084 69.7366 69.85174 69.96624
>>>>>>
>>>>>> ZDEF 59 LEVELS
>>>>>>
>>>>>> 5 10 20 30 40 50 60 80 100 120 140 160 180 200 220
>>>>>>
>>>>>>    240 260 280 300 320 360 400 440 480 520 560 600 640 680
>>>>>>
>>>>>>    720 760 800 840 880 920 960 1000 1040 1080 1120 1160 1200
>>>>>>
>>>>>>    1240 1280 1320 1360 1400 1440 1480 1520 1560 1600 1650 1700
>>>>>>
>>>>>>    1750 1800 1850 1900 1950
>>>>>>
>>>>>> TDEF 255 LINEAR 05jan2000 7dy
>>>>>>
>>>>>> VARS 2
>>>>>>
>>>>>> temperature 59 t,z,x,y sea temperature
>>>>>>
>>>>>> pct_variance 59 t,z,x,y variance a priori
>>>>>>
>>>>>> ENDVARS
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> And the ncdump -h output for one of my netcdf file is :
>>>>>>
>>>>>> netcdf OA_20040317 {
>>>>>>
>>>>>> dimensions:
>>>>>>
>>>>>>        time = 1 ;
>>>>>>
>>>>>>        depth = 59 ;
>>>>>>
>>>>>>        latitude = 597 ;
>>>>>>
>>>>>>        longitude = 421 ;
>>>>>>
>>>>>> variables:
>>>>>>
>>>>>>        float time(time) ;
>>>>>>
>>>>>>                time:units = "days since 1950-01-01 00:00:00" ;
>>>>>>
>>>>>>        float latitude(latitude) ;
>>>>>>
>>>>>>                latitude:units = "degree_north" ;
>>>>>>
>>>>>>                latitude:valid_min = -90.f ;
>>>>>>
>>>>>>                latitude:valid_max = 90.f ;
>>>>>>
>>>>>>        float longitude(longitude) ;
>>>>>>
>>>>>>                longitude:units = "degree_east" ;
>>>>>>
>>>>>>                longitude:valid_min = -180.f ;
>>>>>>
>>>>>>                longitude:valid_max = 180.f ;
>>>>>>
>>>>>>        short depth(depth) ;
>>>>>>
>>>>>>                depth:units = "m" ;
>>>>>>
>>>>>>                depth:positive = "down" ;
>>>>>>
>>>>>>                depth:valid_min = 0s ;
>>>>>>
>>>>>>                depth:valid_max = 2000s ;
>>>>>>
>>>>>>        short temperature(time, depth, latitude, longitude) ;
>>>>>>
>>>>>>                temperature:long_name = "Temperature" ;
>>>>>>
>>>>>>                temperature:missing_value = 32767s ;
>>>>>>
>>>>>>                temperature:_FillValue = 32767s ;
>>>>>>
>>>>>>                temperature:units = "degree_Celsius" ;
>>>>>>
>>>>>>                temperature:valid_min = -3s ;
>>>>>>
>>>>>>                temperature:valid_max = 40s ;
>>>>>>
>>>>>>                temperature:add_offset = 20. ;
>>>>>>
>>>>>>                temperature:scale_factor = 0.001 ;
>>>>>>
>>>>>>                temperature:comment = "Temperature estimate (by
>>>>>> objective analysis)" ;
>>>>>>
>>>>>>        short pct_variance(time, depth, latitude, longitude) ;
>>>>>>
>>>>>>                pct_variance:long_name = "Estimate variance / a
>>>>>> priori variance (percent)" ;
>>>>>>
>>>>>>                pct_variance:missing_value = 32767s ;
>>>>>>
>>>>>>                pct_variance:_FillValue = 32767s ;
>>>>>>
>>>>>>                pct_variance:units = "percent" ;
>>>>>>
>>>>>>                pct_variance:valid_min = 0s ;
>>>>>>
>>>>>>                pct_variance:valid_max = 100s ;
>>>>>>
>>>>>>                pct_variance:scale_factor = 0.01 ;
>>>>>>
>>>>>>
>>>>>>
>>>>>> // global attributes:
>>>>>>
>>>>>>                :CONVENTIONS = "COARDS" ;
>>>>>>
>>>>>>                :producer_agengy = "IFREMER" ;
>>>>>>
>>>>>>                :project_name = "CORIOLIS" ;
>>>>>>
>>>>>>                :creation_time = "14-Jan-2005 19:19:52" ;
>>>>>>
>>>>>>                :software_version = "OA_V3.02" ;
>>>>>>
>>>>>>                :product_version = "OA3_02" ;
>>>>>>
>>>>>>                :data_set = "STD3_02" ;
>>>>>>
>>>>>>                :data_manager = "Emmanuelle Autret, Thomas
>>>>>> Loubrieu" ;
>>>>>>
>>>>>>                :estimate_date = "17-Mar-2004" ;
>>>>>>
>>>>>>                :south_latitude = "-70" ;
>>>>>>
>>>>>>                :north_latitude = "70" ;
>>>>>>
>>>>>>                :west_longitude = "-100" ;
>>>>>>
>>>>>>                :east_longitude = "40" ;
>>>>>>
>>>>>> }
>>>>>>
>>>>>>
>>>>>>
>>>>>> Does anyone knows what happens ?
>>>>>>
>>>>>> Thanks,
>>>>>>
>>>>>> Thomas
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>
>>>>
>>>>
>> netcdf OA_20040811 {
>> dimensions:
>>         time = 1 ;
>>         depth = 59 ;
>>         latitude = 597 ;
>>         longitude = 421 ;
>> variables:
>>         float time(time) ;
>>                 time:units = "days since 1950-01-01 00:00:00" ;
>>         float latitude(latitude) ;
>>                 latitude:units = "degree_north" ;
>>                 latitude:valid_min = -90.f ;
>>                 latitude:valid_max = 90.f ;
>>         float longitude(longitude) ;
>>                 longitude:units = "degree_east" ;
>>                 longitude:valid_min = -180.f ;
>>                 longitude:valid_max = 180.f ;
>>         short depth(depth) ;
>>                 depth:units = "m" ;
>>                 depth:positive = "down" ;
>>                 depth:valid_min = 0s ;
>>                 depth:valid_max = 2000s ;
>>         short temperature(time, depth, latitude, longitude) ;
>>                 temperature:long_name = "Temperature" ;
>>                 temperature:missing_value = 32767s ;
>>                 temperature:_FillValue = 32767s ;
>>                 temperature:units = "degree_Celsius" ;
>>                 temperature:valid_min = -3s ;
>>                 temperature:valid_max = 40s ;
>>                 temperature:add_offset = 20. ;
>>                 temperature:scale_factor = 0.001 ;
>>                 temperature:comment = "Temperature estimate (by
>> objective analysis)" ;
>>         short pct_variance(time, depth, latitude, longitude) ;
>>                 pct_variance:long_name = "Estimate variance / a
>> priori variance (percent)" ;
>>                 pct_variance:missing_value = 32767s ;
>>                 pct_variance:_FillValue = 32767s ;
>>                 pct_variance:units = "percent" ;
>>                 pct_variance:valid_min = 0s ;
>>                 pct_variance:valid_max = 100s ;
>>                 pct_variance:scale_factor = 0.01 ;
>>
>> // global attributes:
>>                 :CONVENTIONS = "COARDS" ;
>>                 :producer_agengy = "IFREMER" ;
>>                 :project_name = "CORIOLIS" ;
>>                 :creation_time = "14-Jan-2005 23:36:22" ;
>>                 :software_version = "OA_V3.02" ;
>>                 :product_version = "OA3_02" ;
>>                 :data_set = "STD3_02" ;
>>                 :data_manager = "Emmanuelle Autret, Thomas Loubrieu" ;
>>                 :estimate_date = "11-Aug-2004" ;
>>                 :south_latitude = "-70" ;
>>                 :north_latitude = "70" ;
>>                 :west_longitude = "-100" ;
>>                 :east_longitude = "40" ;
>> data:
>>
>>  time = 19946 ;
>>
>>  latitude = -69.96624, -69.85174, -69.7366, -69.62084, -69.50445,
>> -69.38742,
>>     -69.26975, -69.15144, -69.03249, -68.91288, -68.79263, -68.67171,
>>     -68.55015, -68.42793, -68.30503, -68.18148, -68.05725, -67.93236,
>>     -67.80679, -67.68053, -67.5536, -67.42599, -67.29768,
>> -67.16869, -67.039,
>>     -66.90862, -66.77753, -66.64575, -66.51326, -66.38006, -66.24615,
>>     -66.11152, -65.97617, -65.84011, -65.70332, -65.5658, -65.42755,
>>     -65.28857, -65.14885, -65.00839, -64.8672, -64.72525, -64.58255,
>>     -64.4391, -64.29491, -64.14995, -64.00423, -63.85774, -63.71049,
>>     -63.56247, -63.41367, -63.2641, -63.11375, -62.96262,
>> -62.81071, -62.658,
>>     -62.50451, -62.35021, -62.19513, -62.03924, -61.88255, -61.72505,
>>     -61.56674, -61.40763, -61.24769, -61.08694, -60.92537, -60.76297,
>>     -60.59975, -60.4357, -60.27082, -60.1051, -59.93855, -59.77115,
>>     -59.60291, -59.43382, -59.26389, -59.0931, -58.92146, -58.74895,
>>     -58.57559, -58.40137, -58.22628, -58.05033, -57.8735, -57.6958,
>>     -57.51722, -57.33776, -57.15743, -56.97621, -56.79411, -56.61111,
>>     -56.42722, -56.24245, -56.05677, -55.8702, -55.68272, -55.49435,
>>     -55.30507, -55.11488, -54.92377, -54.73176, -54.53884, -54.34499,
>>     -54.15023, -53.95455, -53.75795, -53.56042, -53.36197, -53.16258,
>>     -52.96227, -52.76102, -52.55884, -52.35572, -52.15167, -51.94667,
>>     -51.74074, -51.53386, -51.32603, -51.11727, -50.90755, -50.69688,
>>     -50.48527, -50.2727, -50.05918, -49.8447, -49.62927, -49.41288,
>>     -49.19553, -48.97723, -48.75796, -48.53773, -48.31654, -48.09439,
>>     -47.87127, -47.64719, -47.42214, -47.19613, -46.96915, -46.7412,
>>     -46.51228, -46.2824, -46.05155, -45.81973, -45.58693, -45.35318,
>>     -45.11845, -44.88274, -44.64608, -44.40844, -44.16983, -43.93025,
>>     -43.68971, -43.44819, -43.20571, -42.96226, -42.71784, -42.47246,
>>     -42.22611, -41.97879, -41.73051, -41.48127, -41.23106, -40.9799,
>>     -40.72777, -40.47469, -40.22064, -39.96564, -39.70969, -39.45279,
>>     -39.19493, -38.93612, -38.67636, -38.41566, -38.15401, -37.89143,
>>     -37.6279, -37.36343, -37.09803, -36.8317, -36.56443, -36.29624,
>>     -36.02712, -35.75708, -35.48612, -35.21424, -34.94145, -34.66775,
>>     -34.39314, -34.11763, -33.84122, -33.56391, -33.28571, -33.00661,
>>     -32.72664, -32.44578, -32.16404, -31.88143, -31.59795, -31.3136,
>>     -31.02839, -30.74232, -30.45541, -30.16764, -29.87904, -29.58959,
>>     -29.29932, -29.00821, -28.71628, -28.42354, -28.12998, -27.83562,
>>     -27.54046, -27.2445, -26.94775, -26.65022, -26.35192, -26.05284,
>>     -25.75299, -25.45239, -25.15103, -24.84893, -24.54609, -24.24251,
>>     -23.93821, -23.63319, -23.32746, -23.02102, -22.71389, -22.40606,
>>     -22.09755, -21.78837, -21.47852, -21.168, -20.85683, -20.54502,
>>     -20.23257, -19.91949, -19.60579, -19.29148, -18.97656, -18.66105,
>>     -18.34495, -18.02827, -17.71101, -17.3932, -17.07483, -16.75592,
>>     -16.43648, -16.1165, -15.79601, -15.47501, -15.15352, -14.83153,
>>     -14.50907, -14.18613, -13.86273, -13.53888, -13.21459, -12.88987,
>>     -12.56473, -12.23918, -11.91322, -11.58687, -11.26014, -10.93304,
>>     -10.60558, -10.27777, -9.949614, -9.62113, -9.292328, -8.963216,
>>     -8.633805, -8.304108, -7.974132, -7.64389, -7.313393, -6.982651,
>>     -6.651675, -6.320475, -5.989064, -5.657451, -5.325648, -4.993666,
>>     -4.661516, -4.329208, -3.996755, -3.664166, -3.331455, -2.99863,
>>     -2.665704, -2.332689, -1.999594, -1.666432, -1.333213, -0.9999492,
>>     -0.6666516, -0.3333315, 0, 0.3333315, 0.6666516, 0.9999492,
>> 1.333213,
>>     1.666432, 1.999594, 2.332689, 2.665704, 2.99863, 3.331455,
>> 3.664166,
>>     3.996755, 4.329208, 4.661516, 4.993666, 5.325648, 5.657451,
>> 5.989064,
>>     6.320475, 6.651675, 6.982651, 7.313393, 7.64389, 7.974132,
>> 8.304108,
>>     8.633805, 8.963216, 9.292328, 9.62113, 9.949614, 10.27777,
>> 10.60558,
>>     10.93304, 11.26014, 11.58687, 11.91322, 12.23918, 12.56473,
>> 12.88987,
>>     13.21459, 13.53888, 13.86273, 14.18613, 14.50907, 14.83153,
>> 15.15352,
>>     15.47501, 15.79601, 16.1165, 16.43648, 16.75592, 17.07483,
>> 17.3932,
>>     17.71101, 18.02827, 18.34495, 18.66105, 18.97656, 19.29148,
>> 19.60579,
>>     19.91949, 20.23257, 20.54502, 20.85683, 21.168, 21.47852,
>> 21.78837,
>>     22.09755, 22.40606, 22.71389, 23.02102, 23.32746, 23.63319,
>> 23.93821,
>>     24.24251, 24.54609, 24.84893, 25.15103, 25.45239, 25.75299,
>> 26.05284,
>>     26.35192, 26.65022, 26.94775, 27.2445, 27.54046, 27.83562,
>> 28.12998,
>>     28.42354, 28.71628, 29.00821, 29.29932, 29.58959, 29.87904,
>> 30.16764,
>>     30.45541, 30.74232, 31.02839, 31.3136, 31.59795, 31.88143,
>> 32.16404,
>>     32.44578, 32.72664, 33.00661, 33.28571, 33.56391, 33.84122,
>> 34.11763,
>>     34.39314, 34.66775, 34.94145, 35.21424, 35.48612, 35.75708,
>> 36.02712,
>>     36.29624, 36.56443, 36.8317, 37.09803, 37.36343, 37.6279,
>> 37.89143,
>>     38.15401, 38.41566, 38.67636, 38.93612, 39.19493, 39.45279,
>> 39.70969,
>>     39.96564, 40.22064, 40.47469, 40.72777, 40.9799, 41.23106,
>> 41.48127,
>>     41.73051, 41.97879, 42.22611, 42.47246, 42.71784, 42.96226,
>> 43.20571,
>>     43.44819, 43.68971, 43.93025, 44.16983, 44.40844, 44.64608,
>> 44.88274,
>>     45.11845, 45.35318, 45.58693, 45.81973, 46.05155, 46.2824,
>> 46.51228,
>>     46.7412, 46.96915, 47.19613, 47.42214, 47.64719, 47.87127,
>> 48.09439,
>>     48.31654, 48.53773, 48.75796, 48.97723, 49.19553, 49.41288,
>> 49.62927,
>>     49.8447, 50.05918, 50.2727, 50.48527, 50.69688, 50.90755,
>> 51.11727,
>>     51.32603, 51.53386, 51.74074, 51.94667, 52.15167, 52.35572,
>> 52.55884,
>>     52.76102, 52.96227, 53.16258, 53.36197, 53.56042, 53.75795,
>> 53.95455,
>>     54.15023, 54.34499, 54.53884, 54.73176, 54.92377, 55.11488,
>> 55.30507,
>>     55.49435, 55.68272, 55.8702, 56.05677, 56.24245, 56.42722,
>> 56.61111,
>>     56.79411, 56.97621, 57.15743, 57.33776, 57.51722, 57.6958,
>> 57.8735,
>>     58.05033, 58.22628, 58.40137, 58.57559, 58.74895, 58.92146,
>> 59.0931,
>>     59.26389, 59.43382, 59.60291, 59.77115, 59.93855, 60.1051,
>> 60.27082,
>>     60.4357, 60.59975, 60.76297, 60.92537, 61.08694, 61.24769,
>> 61.40763,
>>     61.56674, 61.72505, 61.88255, 62.03924, 62.19513, 62.35021,
>> 62.50451,
>>     62.658, 62.81071, 62.96262, 63.11375, 63.2641, 63.41367, 63.56247,
>>     63.71049, 63.85774, 64.00423, 64.14995, 64.29491, 64.4391,
>> 64.58255,
>>     64.72525, 64.8672, 65.00839, 65.14885, 65.28857, 65.42755,
>> 65.5658,
>>     65.70332, 65.84011, 65.97617, 66.11152, 66.24615, 66.38006,
>> 66.51326,
>>     66.64575, 66.77753, 66.90862, 67.039, 67.16869, 67.29768,
>> 67.42599,
>>     67.5536, 67.68053, 67.80679, 67.93236, 68.05725, 68.18148,
>> 68.30503,
>>     68.42793, 68.55015, 68.67171, 68.79263, 68.91288, 69.03249,
>> 69.15144,
>>     69.26975, 69.38742, 69.50445, 69.62084, 69.7366, 69.85174,
>> 69.96624 ;
>>
>>  longitude = -100, -99.66666, -99.33334, -99, -98.66666, -98.33334,
>> -98,
>>     -97.66666, -97.33334, -97, -96.66666, -96.33334, -96, -95.66666,
>>     -95.33334, -95, -94.66666, -94.33334, -94, -93.66666,
>> -93.33334, -93,
>>     -92.66666, -92.33334, -92, -91.66666, -91.33334, -91, -90.66666,
>>     -90.33334, -90, -89.66666, -89.33334, -89, -88.66666,
>> -88.33334, -88,
>>     -87.66666, -87.33334, -87, -86.66666, -86.33334, -86, -85.66666,
>>     -85.33334, -85, -84.66666, -84.33334, -84, -83.66666,
>> -83.33334, -83,
>>     -82.66666, -82.33334, -82, -81.66666, -81.33334, -81, -80.66666,
>>     -80.33334, -80, -79.66666, -79.33334, -79, -78.66666,
>> -78.33334, -78,
>>     -77.66666, -77.33334, -77, -76.66666, -76.33334, -76, -75.66666,
>>     -75.33334, -75, -74.66666, -74.33334, -74, -73.66666,
>> -73.33334, -73,
>>     -72.66666, -72.33334, -72, -71.66666, -71.33334, -71, -70.66666,
>>     -70.33334, -70, -69.66666, -69.33334, -69, -68.66666,
>> -68.33334, -68,
>>     -67.66666, -67.33334, -67, -66.66666, -66.33334, -66, -65.66666,
>>     -65.33334, -65, -64.66666, -64.33334, -64, -63.66667,
>> -63.33333, -63,
>>     -62.66667, -62.33333, -62, -61.66667, -61.33333, -61, -60.66667,
>>     -60.33333, -60, -59.66667, -59.33333, -59, -58.66667,
>> -58.33333, -58,
>>     -57.66667, -57.33333, -57, -56.66667, -56.33333, -56, -55.66667,
>>     -55.33333, -55, -54.66667, -54.33333, -54, -53.66667,
>> -53.33333, -53,
>>     -52.66667, -52.33333, -52, -51.66667, -51.33333, -51, -50.66667,
>>     -50.33333, -50, -49.66667, -49.33333, -49, -48.66667,
>> -48.33333, -48,
>>     -47.66667, -47.33333, -47, -46.66667, -46.33333, -46, -45.66667,
>>     -45.33333, -45, -44.66667, -44.33333, -44, -43.66667,
>> -43.33333, -43,
>>     -42.66667, -42.33333, -42, -41.66667, -41.33333, -41, -40.66667,
>>     -40.33333, -40, -39.66667, -39.33333, -39, -38.66667,
>> -38.33333, -38,
>>     -37.66667, -37.33333, -37, -36.66667, -36.33333, -36, -35.66667,
>>     -35.33333, -35, -34.66667, -34.33333, -34, -33.66667,
>> -33.33333, -33,
>>     -32.66667, -32.33333, -32, -31.66667, -31.33333, -31, -30.66667,
>>     -30.33333, -30, -29.66667, -29.33333, -29, -28.66667,
>> -28.33333, -28,
>>     -27.66667, -27.33333, -27, -26.66667, -26.33333, -26, -25.66667,
>>     -25.33333, -25, -24.66667, -24.33333, -24, -23.66667,
>> -23.33333, -23,
>>     -22.66667, -22.33333, -22, -21.66667, -21.33333, -21, -20.66667,
>>     -20.33333, -20, -19.66667, -19.33333, -19, -18.66667,
>> -18.33333, -18,
>>     -17.66667, -17.33333, -17, -16.66667, -16.33333, -16, -15.66667,
>>     -15.33333, -15, -14.66667, -14.33333, -14, -13.66667,
>> -13.33333, -13,
>>     -12.66667, -12.33333, -12, -11.66667, -11.33333, -11, -10.66667,
>>     -10.33333, -10, -9.666667, -9.333333, -9, -8.666667, -8.333333,
>> -8,
>>     -7.666667, -7.333333, -7, -6.666667, -6.333333, -6, -5.666667,
>> -5.333333,
>>     -5, -4.666667, -4.333333, -4, -3.666667, -3.333333, -3, -2.666667,
>>     -2.333333, -2, -1.666667, -1.333333, -1, -0.6666667,
>> -0.3333333, 0,
>>     0.3333333, 0.6666667, 1, 1.333333, 1.666667, 2, 2.333333,
>> 2.666667, 3,
>>     3.333333, 3.666667, 4, 4.333333, 4.666667, 5, 5.333333,
>> 5.666667, 6,
>>     6.333333, 6.666667, 7, 7.333333, 7.666667, 8, 8.333333,
>> 8.666667, 9,
>>     9.333333, 9.666667, 10, 10.33333, 10.66667, 11, 11.33333,
>> 11.66667, 12,
>>     12.33333, 12.66667, 13, 13.33333, 13.66667, 14, 14.33333,
>> 14.66667, 15,
>>     15.33333, 15.66667, 16, 16.33333, 16.66667, 17, 17.33333,
>> 17.66667, 18,
>>     18.33333, 18.66667, 19, 19.33333, 19.66667, 20, 20.33333,
>> 20.66667, 21,
>>     21.33333, 21.66667, 22, 22.33333, 22.66667, 23, 23.33333,
>> 23.66667, 24,
>>     24.33333, 24.66667, 25, 25.33333, 25.66667, 26, 26.33333,
>> 26.66667, 27,
>>     27.33333, 27.66667, 28, 28.33333, 28.66667, 29, 29.33333,
>> 29.66667, 30,
>>     30.33333, 30.66667, 31, 31.33333, 31.66667, 32, 32.33333,
>> 32.66667, 33,
>>     33.33333, 33.66667, 34, 34.33333, 34.66667, 35, 35.33333,
>> 35.66667, 36,
>>     36.33333, 36.66667, 37, 37.33333, 37.66667, 38, 38.33333,
>> 38.66667, 39,
>>     39.33333, 39.66667, 40 ;
>>
>>  depth = 5, 10, 20, 30, 40, 50, 60, 80, 100, 120, 140, 160, 180,
>> 200, 220,
>>     240, 260, 280, 300, 320, 360, 400, 440, 480, 520, 560, 600,
>> 640, 680,
>>     720, 760, 800, 840, 880, 920, 960, 1000, 1040, 1080, 1120,
>> 1160, 1200,
>>     1240, 1280, 1320, 1360, 1400, 1440, 1480, 1520, 1560, 1600,
>> 1650, 1700,
>>     1750, 1800, 1850, 1900, 1950 ;
>> }
>>



More information about the gradsusr mailing list