[gradsusr] Control File for Mask
Justin Hicks
jhicks2014 at gmail.com
Mon Jun 27 13:19:24 EDT 2016
Jeff,
This is the code:
ga-> open GLDAS_MOS10_M.ctl
Scanning description file: GLDAS_MOS10_M.ctl
Data file
/data2/gds/GLDAS_V1/GLDAS_MOS10_M/%y4/GLDAS_MOS10_M.A%y4%m2.001.grb is open
as file 1
LON set to 0 360
LAT set to -59.5 89.5
LEV set to 1 1
Time values set: 1979:1:1:0 1979:1:1:0
E set to 1 1
ga-> open Region15_AFR_maskv3.nc
Scanning description file: Region15_AFR_maskv3.nc
Warning: OPTIONS keyword "template" is used, but the
DSET entry contains no substitution templates.
Data file ~/Downloads/Region15_AFR_maskv3.nc is open as file 2
ga-> d maskout(evap,data.2(time=00Z01JUN1990)-0.5)
Data Request Warning: Request is completely outside file limits
Entire grid contents are set to missing data
Grid limits of file: X = 1 360 Y = 1 180 Z = 1 1 T = 1 1 E = 1 1
Grid limits of request: X = 180 541 Y = 31 180 Z = 1 1 T = 138 138 E
= 1 1
Warning issued for variable = data.2
Cannot contour grid - all undefined values
I sent you the control file for the mask before, so here is the control
file for the GLDAS data:
dset /data2/gds/GLDAS_V1/GLDAS_MOS10_M/%y4/GLDAS_MOS10_M.A%y4%m2.001.grb
index /data2/control/GLDAS_MOS10_M.idx
options template
undef 9.999E+20
title /data2/gds/GLDAS_V1/GLDAS_MOS10_M/1979/GLDAS_MOS10_M.A197901.001.grb
* produced by grib2ctl v0.9.13
* command line options:
/data2/gds/GLDAS_V1/GLDAS_MOS10_M/1979/GLDAS_MOS10_M.A197901.001.grb
dtype grib 0
ydef 150 linear -59.500000 1
xdef 360 linear -179.500000 1.000000
tdef 436 linear 00Z01jan1979 1mo
zdef 1 linear 1 1
vars 24
AvgSurfT 0 138,1,0 ** Average Surface Temperature K
Canopint 0 71,1,0 ** Total canopy water storage kg/m^2
Evap 0 57,1,0 ** Total Evapotranspiration kg/m^2/s
LWdown 0 205,1,0 ** Surface incident longwave radiation W/m^2
LWnet 0 112,1,0 ** Net Longwave Radiation W/m^2
PSurf 0 1,1,0 ** Surface pressure Pa
Qair 0 51,1,0 ** Near surface specific humidity kg/kg
Qg 0 155,1,0 ** Ground Heat Flux W/m^2
Qh 0 122,1,0 ** Sensible Heat Flux W/m^2
Qle 0 121,1,0 ** Latent Heat Flux W/m^2
Qs 0 235,1,0 ** Surface Runoff kg/m^2/s
Qsb 0 234,1,0 ** Subsurface Runoff kg/m^2/s
Qsm 0 99,1,0 ** Snowmelt kg/m^2/s
Rainf 0 132,1,0 ** Rainfall rate kg/m^2/s
SWE 0 65,1,0 ** Snow Water Equivalent kg/m^2
SWdown 0 204,1,0 ** Surface incident shortwave radiation W/m^2
SWnet 0 111,1,0 ** Net Shortwave Radiation W/m^2
Snowf 0 131,1,0 ** Snowfall rate kg/m^2/s
TSoil 0 85,112,50001 ** Deep Soil Temperature K
SoilM1 0 86,112,50001 ** 0-2 cm Average layer 1 soil moisture kg/m^2
SoilM2 0 86,112,40002 ** 2-150 cm Average layer 2 soil moisture kg/m^2
SoilM3 0 86,112,30003 ** 150-350 cm Average layer 3 soil moisture
kg/m^2
Tair 0 11,1,0 ** Near surface air temperature K
Wind 0 32,1,0 ** Near surface wind magnitude m/s
ENDVARS
Now that I look at it, would it be the different 'ydef' between the two
files that is causing this problem? And if so, how could I change this?
-Justin
On Mon, Jun 27, 2016 at 1:03 PM, Jeff Duda <jeffduda319 at gmail.com> wrote:
> Justin,
> Show me the exact code you're running, including any and all control files
> opened by the code.
>
> Jeff
>
> On Mon, Jun 27, 2016 at 12:00 PM, Justin Hicks <jhicks2014 at gmail.com>
> wrote:
>
>> Jeff,
>>
>> I tried to use manipulation strings using the following command: 'd
>> maskout(evap,data.2(time=00Z01JUN1990)-0.5)'.
>> Note how I changed the variable from 'land' to 'data' just because of the
>> ncdump contents. I received this error after entering this command: Data
>> Request Warning:
>>
>> Request is completely outside file limits
>> Entire grid contents are set to missing data
>> Grid limits of file: X = 1 360 Y = 1 180 Z = 1 1 T = 1 1 E = 1
>> 1
>> Grid limits of request: X = 180 541 Y = 31 180 Z = 1 1 T = 138 138
>> E = 1 1
>> Warning issued for variable = data.2
>> Cannot contour grid - all undefined values
>>
>>
>> I'll post the contents of ncdump -c below:
>>
>> ncdump -c Region15_AFR_maskv3.nc
>> netcdf Region15_AFR_maskv3 {
>> dimensions:
>> Longitude = 360 ;
>> Latitude = 180 ;
>> variables:
>> float Longitude(Longitude) ;
>> float Latitude(Latitude) ;
>> float data(Latitude, Longitude) ;
>> data:
>>
>> Longitude = -179.5, -178.5, -177.5, -176.5, -175.5, -174.5, -173.5,
>> -172.5,
>> -171.5, -170.5, -169.5, -168.5, -167.5, -166.5, -165.5, -164.5,
>> -163.5,
>> -162.5, -161.5, -160.5, -159.5, -158.5, -157.5, -156.5, -155.5,
>> -154.5,
>> -153.5, -152.5, -151.5, -150.5, -149.5, -148.5, -147.5, -146.5,
>> -145.5,
>> -144.5, -143.5, -142.5, -141.5, -140.5, -139.5, -138.5, -137.5,
>> -136.5,
>> -135.5, -134.5, -133.5, -132.5, -131.5, -130.5, -129.5, -128.5,
>> -127.5,
>> -126.5, -125.5, -124.5, -123.5, -122.5, -121.5, -120.5, -119.5,
>> -118.5,
>> -117.5, -116.5, -115.5, -114.5, -113.5, -112.5, -111.5, -110.5,
>> -109.5,
>> -108.5, -107.5, -106.5, -105.5, -104.5, -103.5, -102.5, -101.5,
>> -100.5,
>> -99.5, -98.5, -97.5, -96.5, -95.5, -94.5, -93.5, -92.5, -91.5, -90.5,
>> -89.5, -88.5, -87.5, -86.5, -85.5, -84.5, -83.5, -82.5, -81.5, -80.5,
>> -79.5, -78.5, -77.5, -76.5, -75.5, -74.5, -73.5, -72.5, -71.5, -70.5,
>> -69.5, -68.5, -67.5, -66.5, -65.5, -64.5, -63.5, -62.5, -61.5, -60.5,
>> -59.5, -58.5, -57.5, -56.5, -55.5, -54.5, -53.5, -52.5, -51.5, -50.5,
>> -49.5, -48.5, -47.5, -46.5, -45.5, -44.5, -43.5, -42.5, -41.5, -40.5,
>> -39.5, -38.5, -37.5, -36.5, -35.5, -34.5, -33.5, -32.5, -31.5, -30.5,
>> -29.5, -28.5, -27.5, -26.5, -25.5, -24.5, -23.5, -22.5, -21.5, -20.5,
>> -19.5, -18.5, -17.5, -16.5, -15.5, -14.5, -13.5, -12.5, -11.5, -10.5,
>> -9.5, -8.5, -7.5, -6.5, -5.5, -4.5, -3.5, -2.5, -1.5, -0.5, 0.5, 1.5,
>> 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5, 10.5, 11.5, 12.5, 13.5, 14.5,
>> 15.5, 16.5, 17.5, 18.5, 19.5, 20.5, 21.5, 22.5, 23.5, 24.5, 25.5,
>> 26.5,
>> 27.5, 28.5, 29.5, 30.5, 31.5, 32.5, 33.5, 34.5, 35.5, 36.5, 37.5,
>> 38.5,
>> 39.5, 40.5, 41.5, 42.5, 43.5, 44.5, 45.5, 46.5, 47.5, 48.5, 49.5,
>> 50.5,
>> 51.5, 52.5, 53.5, 54.5, 55.5, 56.5, 57.5, 58.5, 59.5, 60.5, 61.5,
>> 62.5,
>> 63.5, 64.5, 65.5, 66.5, 67.5, 68.5, 69.5, 70.5, 71.5, 72.5, 73.5,
>> 74.5,
>> 75.5, 76.5, 77.5, 78.5, 79.5, 80.5, 81.5, 82.5, 83.5, 84.5, 85.5,
>> 86.5,
>> 87.5, 88.5, 89.5, 90.5, 91.5, 92.5, 93.5, 94.5, 95.5, 96.5, 97.5,
>> 98.5,
>> 99.5, 100.5, 101.5, 102.5, 103.5, 104.5, 105.5, 106.5, 107.5, 108.5,
>> 109.5, 110.5, 111.5, 112.5, 113.5, 114.5, 115.5, 116.5, 117.5, 118.5,
>> 119.5, 120.5, 121.5, 122.5, 123.5, 124.5, 125.5, 126.5, 127.5, 128.5,
>> 129.5, 130.5, 131.5, 132.5, 133.5, 134.5, 135.5, 136.5, 137.5, 138.5,
>> 139.5, 140.5, 141.5, 142.5, 143.5, 144.5, 145.5, 146.5, 147.5, 148.5,
>> 149.5, 150.5, 151.5, 152.5, 153.5, 154.5, 155.5, 156.5, 157.5, 158.5,
>> 159.5, 160.5, 161.5, 162.5, 163.5, 164.5, 165.5, 166.5, 167.5, 168.5,
>> 169.5, 170.5, 171.5, 172.5, 173.5, 174.5, 175.5, 176.5, 177.5, 178.5,
>> 179.5 ;
>>
>> Latitude = -89.5, -88.5, -87.5, -86.5, -85.5, -84.5, -83.5, -82.5,
>> -81.5,
>> -80.5, -79.5, -78.5, -77.5, -76.5, -75.5, -74.5, -73.5, -72.5, -71.5,
>> -70.5, -69.5, -68.5, -67.5, -66.5, -65.5, -64.5, -63.5, -62.5, -61.5,
>> -60.5, -59.5, -58.5, -57.5, -56.5, -55.5, -54.5, -53.5, -52.5, -51.5,
>> -50.5, -49.5, -48.5, -47.5, -46.5, -45.5, -44.5, -43.5, -42.5, -41.5,
>> -40.5, -39.5, -38.5, -37.5, -36.5, -35.5, -34.5, -33.5, -32.5, -31.5,
>> -30.5, -29.5, -28.5, -27.5, -26.5, -25.5, -24.5, -23.5, -22.5, -21.5,
>> -20.5, -19.5, -18.5, -17.5, -16.5, -15.5, -14.5, -13.5, -12.5, -11.5,
>> -10.5, -9.5, -8.5, -7.5, -6.5, -5.5, -4.5, -3.5, -2.5, -1.5, -0.5,
>> 0.5,
>> 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5, 10.5, 11.5, 12.5, 13.5,
>> 14.5, 15.5, 16.5, 17.5, 18.5, 19.5, 20.5, 21.5, 22.5, 23.5, 24.5,
>> 25.5,
>> 26.5, 27.5, 28.5, 29.5, 30.5, 31.5, 32.5, 33.5, 34.5, 35.5, 36.5,
>> 37.5,
>> 38.5, 39.5, 40.5, 41.5, 42.5, 43.5, 44.5, 45.5, 46.5, 47.5, 48.5,
>> 49.5,
>> 50.5, 51.5, 52.5, 53.5, 54.5, 55.5, 56.5, 57.5, 58.5, 59.5, 60.5,
>> 61.5,
>> 62.5, 63.5, 64.5, 65.5, 66.5, 67.5, 68.5, 69.5, 70.5, 71.5, 72.5,
>> 73.5,
>> 74.5, 75.5, 76.5, 77.5, 78.5, 79.5, 80.5, 81.5, 82.5, 83.5, 84.5,
>> 85.5,
>> 86.5, 87.5, 88.5, 89.5 ;
>> }
>>
>> -Justin
>>
>> On Mon, Jun 27, 2016 at 12:00 PM, Jeff Duda <jeffduda319 at gmail.com>
>> wrote:
>>
>>> Justin,
>>> It very well could be the TDEF line. GrADS is expecting the data in all
>>> files to be valid at the same time. If it isn't, you have to use
>>> manipulation strings to get things to plot right. For example, if the time
>>> of your data is today, 12Z 27 June 2016, then you would have to enter 'd
>>> maskout(evap,land.2(time=12Z27JUN2016)-0.5)'
>>>
>>> You might also want to consider using xdfopen instead of sdfopen. While
>>> you still have to write a control file to use xdfopen, you don't have to
>>> put nearly as many details in the control file as you would for just a
>>> regular open, which requires a full blown control file like you created.
>>>
>>> Jeff
>>>
>>> On Mon, Jun 27, 2016 at 10:21 AM, Justin Hicks <jhicks2014 at gmail.com>
>>> wrote:
>>>
>>>> GrADS Users,
>>>>
>>>> I am trying to open up a netcdf mask file using 'sdfopen', but GrADS is
>>>> not opening it up. The mask should show values ONLY for Africa. I decided
>>>> to make a control file instead, but I am confused as how to write one for a
>>>> mask.
>>>>
>>>> For this mask, land is given a value of 1 and ocean value are given a
>>>> value of 0.
>>>>
>>>> The control file I have is below:
>>>>
>>>> DSET ~/Downloads/Region15_AFR_maskv3.nc
>>>> options template
>>>> undef 0
>>>> Xdef 360 LINEAR -179.5 1.0
>>>> Ydef 180 LINEAR 89.5 1.0
>>>> Zdef 1 LINEAR 1.0 1.0
>>>> Tdef 1 LINEAR 00Z01jan1979 1mo
>>>> Vars 1
>>>> land=>land 0 t,y,x land_mask
>>>> endvars
>>>>
>>>>
>>>> I have tried changing the undef value, but I get the same error each
>>>> time I run the command 'd maskout(evap,land.2-0.5)':
>>>>
>>>> "Cannot contour grid - all undefined values"
>>>>
>>>> I am not sure why all values are coming up undefined (the land area I
>>>> want isn't showing up either, instead it's the default world map). Perhaps
>>>> it is the Tdef? I'm not sure why a mask would need a tdef.
>>>>
>>>> -Justin
>>>>
>>>>
>>>> _______________________________________________
>>>> gradsusr mailing list
>>>> gradsusr at gradsusr.org
>>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>>
>>>>
>>>
>>>
>>> --
>>> Jeff Duda
>>> Post-doctoral research associate
>>> University of Oklahoma School of Meteorology
>>>
>>> _______________________________________________
>>> gradsusr mailing list
>>> gradsusr at gradsusr.org
>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>
>>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>
>
> --
> Jeff Duda
> Post-doctoral research associate
> University of Oklahoma School of Meteorology
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20160627/f475733a/attachment-0001.html
More information about the gradsusr
mailing list