[ncl-talk] select sub_region
Karin Meier-Fleischer
meier-fleischer at dkrz.de
Thu Dec 16 06:40:45 MST 2021
1. The longitudes of tauuo_ssp are in range -180. to 180. degrees which
means that you have to flip them to the range 0 to 359 degrees.
longitude (degrees_east) : min=-180 max=180 (=> 361 elements)
2. And the longitude array has one element too much (first element ==
last element) that's why you have to skip the last element.
3. If you want to keep the variable name tauuo_ssp don't forget to use
the reassignment operator := when extracting a sub-region.
diri = "./"
fname = "tauuo_ssp.nc"
f = addfile(diri+fname, "r")
tauuo_ssp = f->tauuo_ssp(:,0:359)
tauuo_ssp = lonFlip(tauuo_ssp)
tauuo_ssp := tauuo_ssp({10:75},{120:240})
wks = gsn_open_wks("pdf", "plot_tauuo")
res = True
res at gsnAddCyclic = False
res at mpCenterLonF = 180.
res at cnFillOn = True
plot = gsn_csm_contour_map(wks, tauuo_ssp, res)
Am 16.12.21 um 13:21 schrieb Chathurika:
>
> Dear Karin,
>
>
> Thank you so much for your kind answer. Actually It is my fault that I
> have not updated the script after creating the one time nc file to
> send you. The data that I used has the variable nametauuo. I created
> the one time file (variable name was tauuo_ssp instead of
> tauuo),unfortunately forgot to update the information in the script
> before sent it to you. I did not get any error while running the
> script. For some reason it just does not select the right region. You
> too can check it. Even though you give the latitude and longitude
> values associated with the Pacific Ocean, You won't get the defined
> region. That's the only problem I have.
>
>
> Anyway thanks for your valuable time and kind response.
>
>
> Best regards,
>
> Chathu
>
>
>
> -----Original Messages-----
> *From:*"Karin Meier-Fleischer" <meier-fleischer at dkrz.de>
> *Sent Time:*2021-12-16 19:51:40 (Thursday)
> *To:* Chathurika <chatu at scsio.ac.cn>
> *Cc:* "ncl-talk at ucar.edu" <ncl-talk at ucar.edu>
> *Subject:* Re: [ncl-talk] select sub_region
>
> Running your script causes the following error which tells you
> exactly what is going wrong.
>
> ncl NCL_Talk.ncl
> Copyright (C) 1995-2019 - All Rights Reserved
> University Corporation for Atmospheric Research
> NCAR Command Language Version 6.6.2
> The use of this software is governed by a License Agreement.
> See http://www.ncl.ucar.edu/ for more details.
> fatal:["Execute.c":133]:variable *(tauuo)* is not in file
> (./tauuo_ssp.nc)
> fatal:["Execute.c":6378]:variable *(tauuo)* is not in file (f)
> fatal:["Execute.c":8637]:Execute: *Error occurred at or near line
> 9 in file NCL_Talk.ncl*
>
> fatal:Variable (tauuo_ssp) is undefined
> fatal:["Execute.c":8637]:Execute: Error occurred at or near line
> 10 in file NCL_Talk.ncl
>
> fatal:Variable (tauuo_ssp) is undefined
> fatal:["Execute.c":8637]:Execute: Error occurred at or near line
> 11 in file NCL_Talk.ncl
>
>
> The variable you try to read has the name tauuo_ssp and not tauuo!
> And NCL is so kind to tell you at which line you have to look first.
> Change the variable name from tauuo to tauuo_ssp and it works fine.
>
> If you do not know the contents of a netCDF file you should have a
> look at it with the command line programs /*ncl_filedump*/ or
> /*ncdump*/.
>
>
> Am 15.12.21 um 13:43 schrieb Chathurika:
>>
>> Dear Karin,
>>
>>
>> Thank you very much for your reply. Please find the attached
>> script and data file (one time step). I want to select the
>> pacific ocean (120E-120W).
>>
>>
>> Thanks again and best regards,
>>
>> Chathu
>>
>>
>> -----Original Messages-----
>> *From:*"Karin Meier-Fleischer" <meier-fleischer at dkrz.de>
>> *Sent Time:*2021-12-15 19:00:13 (Wednesday)
>> *To:* Chathurika <chatu at scsio.ac.cn>
>> *Cc:* "ncl-talk at ucar.edu" <ncl-talk at ucar.edu>
>> *Subject:* Re: [ncl-talk] select sub_region
>>
>> Hi Chatchurika,
>>
>> it is very hard to say what is going wrong without the data
>> and only the information 'It does not select the Pacifik
>> Ocean'. Can you send one timestep of your data and best would
>> be to provide your NCL script.
>>
>> -Karin
>>
>> Am 14.12.21 um 18:17 schrieb Chathurika:
>>>
>>> Dear Karin,
>>>
>>>
>>> Many many thanks for your email. I am sorry It is a typo. I
>>> actually did try this
>>>
>>>> = f->tauuo(:,{10:75},{-120:120}) and vice versa.
>>>> = f->tauuo(:,{10:75},300:60}) and vice versa.
>>>>
>>> However, It does not select the Pacific Ocean.
>>>
>>> Summary of longitude is;
>>>
>>>
>>>
>>> Variable: lon (coordinate)
>>>
>>> Type: double
>>>
>>> Total Size: 2888 bytes
>>>
>>> 361 values
>>>
>>> Number of Dimensions: 1
>>>
>>> Dimensions and sizes: [lon | 361]
>>>
>>> Coordinates:
>>>
>>> Number Of Attributes: 4
>>>
>>> standard_name : longitude
>>>
>>> long_name : longitude
>>>
>>> units : degrees_east
>>>
>>> axis : X
>>>
>>> (0) -180
>>>
>>> (1) -179
>>>
>>> (2) -178
>>>
>>> (3) -177
>>>
>>> (4) -176
>>>
>>> (5) -175
>>>
>>> (6) -174
>>>
>>> (7) -173
>>>
>>> (8) -172
>>>
>>> (9) -171
>>>
>>> (10) -170
>>>
>>> (11) -169
>>>
>>> (12) -168
>>>
>>> (13) -167
>>>
>>> (14) -166
>>>
>>> (15) -165
>>>
>>> (16) -164
>>>
>>> (17) -163
>>>
>>> (18) -162
>>>
>>> (19) -161
>>>
>>> (20) -160
>>>
>>> (21) -159
>>>
>>> (22) -158
>>>
>>> (23) -157
>>>
>>> (24) -156
>>>
>>> (25) -155
>>>
>>> (26) -154
>>>
>>> (27) -153
>>>
>>> (28) -152
>>>
>>> (29) -151
>>>
>>> (30) -150
>>>
>>> (31) -149
>>>
>>> (32) -148
>>>
>>> (33) -147
>>>
>>> (34) -146
>>>
>>> (35) -145
>>>
>>> (36) -144
>>>
>>> (37) -143
>>>
>>> (38) -142
>>>
>>> (39) -141
>>>
>>> (40) -140
>>>
>>> (41) -139
>>>
>>> (42) -138
>>>
>>> (43) -137
>>>
>>> (44) -136
>>>
>>> (45) -135
>>>
>>> (46) -134
>>>
>>> (47) -133
>>>
>>> (48) -132
>>>
>>> (49) -131
>>>
>>> (50) -130
>>>
>>> (51) -129
>>>
>>> (52) -128
>>>
>>> (53) -127
>>>
>>> (54) -126
>>>
>>> (55) -125
>>>
>>> (56) -124
>>>
>>> (57) -123
>>>
>>> (58) -122
>>>
>>> (59) -121
>>>
>>> (60) -120
>>>
>>> (61) -119
>>>
>>> (62) -118
>>>
>>> (63) -117
>>>
>>> (64) -116
>>>
>>> (65) -115
>>>
>>> (66) -114
>>>
>>> (67) -113
>>>
>>> (68) -112
>>>
>>> (69) -111
>>>
>>> (70) -110
>>>
>>> (71) -109
>>>
>>> (72) -108
>>>
>>> (73) -107
>>>
>>> (74) -106
>>>
>>> (75) -105
>>>
>>> (76) -104
>>>
>>> (77) -103
>>>
>>> (78) -102
>>>
>>> (79) -101
>>>
>>> (80) -100
>>>
>>> (81) -99
>>>
>>> (82) -98
>>>
>>> (83) -97
>>>
>>> (84) -96
>>>
>>> (85) -95
>>>
>>> (86) -94
>>>
>>> (87) -93
>>>
>>> (88) -92
>>>
>>> (89) -91
>>>
>>> (90) -90
>>>
>>> (91) -89
>>>
>>> (92) -88
>>>
>>> (93) -87
>>>
>>> (94) -86
>>>
>>> (95) -85
>>>
>>> (96) -84
>>>
>>> (97) -83
>>>
>>> (98) -82
>>>
>>> (99) -81
>>>
>>> (100) -80
>>>
>>> (101) -79
>>>
>>> (102) -78
>>>
>>> (103) -77
>>>
>>> (104) -76
>>>
>>> (105) -75
>>>
>>> (106) -74
>>>
>>> (107) -73
>>>
>>> (108) -72
>>>
>>> (109) -71
>>>
>>> (110) -70
>>>
>>> (111) -69
>>>
>>> (112) -68
>>>
>>> (113) -67
>>>
>>> (114) -66
>>>
>>> (115) -65
>>>
>>> (116) -64
>>>
>>> (117) -63
>>>
>>> (118) -62
>>>
>>> (119) -61
>>>
>>> (120) -60
>>>
>>> (121) -59
>>>
>>> (122) -58
>>>
>>> (123) -57
>>>
>>> (124) -56
>>>
>>> (125) -55
>>>
>>> (126) -54
>>>
>>> (127) -53
>>>
>>> (128) -52
>>>
>>> (129) -51
>>>
>>> (130) -50
>>>
>>> (131) -49
>>>
>>> (132) -48
>>>
>>> (133) -47
>>>
>>> (134) -46
>>>
>>> (135) -45
>>>
>>> (136) -44
>>>
>>> (137) -43
>>>
>>> (138) -42
>>>
>>> (139) -41
>>>
>>> (140) -40
>>>
>>> (141) -39
>>>
>>> (142) -38
>>>
>>> (143) -37
>>>
>>> (144) -36
>>>
>>> (145) -35
>>>
>>> (146) -34
>>>
>>> (147) -33
>>>
>>> (148) -32
>>>
>>> (149) -31
>>>
>>> (150) -30
>>>
>>> (151) -29
>>>
>>> (152) -28
>>>
>>> (153) -27
>>>
>>> (154) -26
>>>
>>> (155) -25
>>>
>>> (156) -24
>>>
>>> (157) -23
>>>
>>> (158) -22
>>>
>>> (159) -21
>>>
>>> (160) -20
>>>
>>> (161) -19
>>>
>>> (162) -18
>>>
>>> (163) -17
>>>
>>> (164) -16
>>>
>>> (165) -15
>>>
>>> (166) -14
>>>
>>> (167) -13
>>>
>>> (168) -12
>>>
>>> (169) -11
>>>
>>> (170) -10
>>>
>>> (171) -9
>>>
>>> (172) -8
>>>
>>> (173) -7
>>>
>>> (174) -6
>>>
>>> (175) -5
>>>
>>> (176) -4
>>>
>>> (177) -3
>>>
>>> (178) -2
>>>
>>> (179) -1
>>>
>>> (180) 0
>>>
>>> (181) 1
>>>
>>> (182) 2
>>>
>>> (183) 3
>>>
>>> (184) 4
>>>
>>> (185) 5
>>>
>>> (186) 6
>>>
>>> (187) 7
>>>
>>> (188) 8
>>>
>>> (189) 9
>>>
>>> (190) 10
>>>
>>> (191) 11
>>>
>>> (192) 12
>>>
>>> (193) 13
>>>
>>> (194) 14
>>>
>>> (195) 15
>>>
>>> (196) 16
>>>
>>> (197) 17
>>>
>>> (198) 18
>>>
>>> (199) 19
>>>
>>> (200) 20
>>>
>>> (201) 21
>>>
>>> (202) 22
>>>
>>> (203) 23
>>>
>>> (204) 24
>>>
>>> (205) 25
>>>
>>> (206) 26
>>>
>>> (207) 27
>>>
>>> (208) 28
>>>
>>> (209) 29
>>>
>>> (210) 30
>>>
>>> (211) 31
>>>
>>> (212) 32
>>>
>>> (213) 33
>>>
>>> (214) 34
>>>
>>> (215) 35
>>>
>>> (216) 36
>>>
>>> (217) 37
>>>
>>> (218) 38
>>>
>>> (219) 39
>>>
>>> (220) 40
>>>
>>> (221) 41
>>>
>>> (222) 42
>>>
>>> (223) 43
>>>
>>> (224) 44
>>>
>>> (225) 45
>>>
>>> (226) 46
>>>
>>> (227) 47
>>>
>>> (228) 48
>>>
>>> (229) 49
>>>
>>> (230) 50
>>>
>>> (231) 51
>>>
>>> (232) 52
>>>
>>> (233) 53
>>>
>>> (234) 54
>>>
>>> (235) 55
>>>
>>> (236) 56
>>>
>>> (237) 57
>>>
>>> (238) 58
>>>
>>> (239) 59
>>>
>>> (240) 60
>>>
>>> (241) 61
>>>
>>> (242) 62
>>>
>>> (243) 63
>>>
>>> (244) 64
>>>
>>> (245) 65
>>>
>>> (246) 66
>>>
>>> (247) 67
>>>
>>> (248) 68
>>>
>>> (249) 69
>>>
>>> (250) 70
>>>
>>> (251) 71
>>>
>>> (252) 72
>>>
>>> (253) 73
>>>
>>> (254) 74
>>>
>>> (255) 75
>>>
>>> (256) 76
>>>
>>> (257) 77
>>>
>>> (258) 78
>>>
>>> (259) 79
>>>
>>> (260) 80
>>>
>>> (261) 81
>>>
>>> (262) 82
>>>
>>> (263) 83
>>>
>>> (264) 84
>>>
>>> (265) 85
>>>
>>> (266) 86
>>>
>>> (267) 87
>>>
>>> (268) 88
>>>
>>> (269) 89
>>>
>>> (270) 90
>>>
>>> (271) 91
>>>
>>> (272) 92
>>>
>>> (273) 93
>>>
>>> (274) 94
>>>
>>> (275) 95
>>>
>>> (276) 96
>>>
>>> (277) 97
>>>
>>> (278) 98
>>>
>>> (279) 99
>>>
>>> (280) 100
>>>
>>> (281) 101
>>>
>>> (282) 102
>>>
>>> (283) 103
>>>
>>> (284) 104
>>>
>>> (285) 105
>>>
>>> (286) 106
>>>
>>> (287) 107
>>>
>>> (288) 108
>>>
>>> (289) 109
>>>
>>> (290) 110
>>>
>>> (291) 111
>>>
>>> (292) 112
>>>
>>> (293) 113
>>>
>>> (294) 114
>>>
>>> (295) 115
>>>
>>> (296) 116
>>>
>>> (297) 117
>>>
>>> (298) 118
>>>
>>> (299) 119
>>>
>>> (300) 120
>>>
>>> (301) 121
>>>
>>> (302) 122
>>>
>>> (303) 123
>>>
>>> (304) 124
>>>
>>> (305) 125
>>>
>>> (306) 126
>>>
>>> (307) 127
>>>
>>> (308) 128
>>>
>>> (309) 129
>>>
>>> (310) 130
>>>
>>> (311) 131
>>>
>>> (312) 132
>>>
>>> (313) 133
>>>
>>> (314) 134
>>>
>>> (315) 135
>>>
>>> (316) 136
>>>
>>> (317) 137
>>>
>>> (318) 138
>>>
>>> (319) 139
>>>
>>> (320) 140
>>>
>>> (321) 141
>>>
>>> (322) 142
>>>
>>> (323) 143
>>>
>>> (324) 144
>>>
>>> (325) 145
>>>
>>> (326) 146
>>>
>>> (327) 147
>>>
>>> (328) 148
>>>
>>> (329) 149
>>>
>>> (330) 150
>>>
>>> (331) 151
>>>
>>> (332) 152
>>>
>>> (333) 153
>>>
>>> (334) 154
>>>
>>> (335) 155
>>>
>>> (336) 156
>>>
>>> (337) 157
>>>
>>> (338) 158
>>>
>>> (339) 159
>>>
>>> (340) 160
>>>
>>> (341) 161
>>>
>>> (342) 162
>>>
>>> (343) 163
>>>
>>> (344) 164
>>>
>>> (345) 165
>>>
>>> (346) 166
>>>
>>> (347) 167
>>>
>>> (348) 168
>>>
>>> (349) 169
>>>
>>> (350) 170
>>>
>>> (351) 171
>>>
>>> (352) 172
>>>
>>> (353) 173
>>>
>>> (354) 174
>>>
>>> (355) 175
>>>
>>> (356) 176
>>>
>>> (357) 177
>>>
>>> (358) 178
>>>
>>> (359) 179
>>>
>>> (360) 180
>>>
>>>
>>> It always select the middle part instead of Pacific.
>>>
>>>
>>> Thank you very very much and best regards,
>>>
>>> Chathu
>>>
>>> -----Original Messages-----
>>> *From:*"Karin Meier-Fleischer via ncl-talk"
>>> <ncl-talk at mailman.ucar.edu>
>>> *Sent Time:*2021-12-14 20:33:47 (Tuesday)
>>> *To:* ncl-talk at mailman.ucar.edu
>>> *Cc:*
>>> *Subject:* Re: [ncl-talk] select sub_region
>>>
>>> Hi Chatchu,
>>>
>>> you should have received an error message because you
>>> try to extract data from a 4D variable tauuo which seems
>>> to be a 3D variable (with a different name tauuo_ssp ??).
>>>
>>> Variable: tauuo_ssp
>>> ...
>>> Dimensions and sizes: [time | 720] x [lat | 66] x [lon |
>>> 360]
>>> Coordinates:
>>> time: [569784..744360]
>>> lat: [ 10.. 75]
>>> lon: [ 180..-180]
>>>
>>> Try
>>>
>>> variable = f->tauuo_ssp(:,:,{-120.:120.})
>>>
>>> -Karin
>>>
>>>
>>> Am 14.12.21 um 13:12 schrieb Chathurika via ncl-talk:
>>>>
>>>> Hello Everyone,
>>>>
>>>>
>>>> I have a dataset and the variable summary as follows;
>>>>
>>>>
>>>> Variable: tauuo_ssp
>>>> Type: float
>>>> Total Size: 45809280 bytes
>>>> 11452320 values
>>>> Number of Dimensions: 3
>>>> Dimensions and sizes: [time | 720] x [lat | 66] x [lon
>>>> | 360]
>>>> Coordinates:
>>>> time: [569784..744360]
>>>> lat: [ 10.. 75]
>>>> lon: [ 180..-180]
>>>> Number Of Attributes: 7
>>>> depth : 0
>>>> standard_name : surface_downward_x_stress
>>>> long_name : Surface downward x stress
>>>> units : N m-2
>>>> code : 548
>>>> _FillValue : -9e+33
>>>> missing_value : -9e+33
>>>>
>>>> I want to select the pacific ocean (-120:120). However,
>>>> this does not work. I tried out several ways to choose
>>>> this subsection, but I always got the part that I
>>>> wanted to subtract between -120:120 (means Atlantic
>>>> plus the Indian Ocean instead of Pacific Ocean).
>>>>
>>>>
>>>> I tried this
>>>>
>>>>
>>>> = f->tauuo(:,0,{10:75},{-120:120}) and vice versa.
>>>> = f->tauuo(:,0,{10:75},300:60}) and vice versa.
>>>>
>>>> Nothing selected the Pacific Ocean from 120E to 120W.
>>>>
>>>>
>>>> Could you please give me a hint on how I could do this?
>>>> That would be a great help.
>>>>
>>>>
>>>> Thank you very very much and best regards,
>>>> Chathu
>>>>
>>>> ------------------------------------------------------------------------
>>>>
>>>>
>>>>
>>>>
>>>> //
>>>> //
>>>> //
>>>> //
>>>>
>>>>
>>>> _______________________________________________
>>>> ncl-talk mailing listncl-talk at mailman.ucar.edu List instructions, subscriber options, unsubscribe:https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>>
>>>
>>> ------------------------------------------------------------------------
>>>
>>>
>>>
>>> Wickramage Chathurika Hemamali
>>> Msc in Physical Oceanography
>>> State Key Laboratory of Tropical Oceanography
>>> South China Sea Institute of Oceanology
>>> University of Chinese Academy of Science
>>> China
>>>
>>> /Specialized in Oceanography and Marine Geology (Bachelor)/
>>> /University of Ruhuna/
>>> /Matara/
>>> /Sri Lanka/
>>>
>>> Email : wickramagechathurika at rocketmail.com
>>> chatu at scsio.ac.cn
>>>
>>
>>
>>
>>
>> ------------------------------------------------------------------------
>>
>>
>>
>> Wickramage Chathurika Hemamali
>> Msc in Physical Oceanography
>> State Key Laboratory of Tropical Oceanography
>> South China Sea Institute of Oceanology
>> University of Chinese Academy of Science
>> China
>>
>> /Specialized in Oceanography and Marine Geology (Bachelor)/
>> /University of Ruhuna/
>> /Matara/
>> /Sri Lanka/
>>
>> Email : wickramagechathurika at rocketmail.com
>> chatu at scsio.ac.cn
>>
>
> --
> Dipl. Geophys. Karin Meier-Fleischer
> Visualization Group - NCL, CDO, Python
> Application Support
>
> Deutsches Klimarechenzentrum GmbH (DKRZ)
> Bundesstrasse 45a - D-20146 Hamburg - Germany
>
> Phone: +49 (0)40 460094 126
> Fax: +49 (0)40 460094 270
> E-Mail:meier-fleischer at dkrz.de URL:www.dkrz.de Geschäftsführer: Prof. Dr. Thomas Ludwig
> Sitz der Gesellschaft: Hamburg
> Amtsgericht Hamburg HRB 39784
>
>
>
>
> ------------------------------------------------------------------------
>
>
>
> Wickramage Chathurika Hemamali
> Msc in Physical Oceanography
> State Key Laboratory of Tropical Oceanography
> South China Sea Institute of Oceanology
> University of Chinese Academy of Science
> China
>
> /Specialized in Oceanography and Marine Geology (Bachelor)/
> /University of Ruhuna/
> /Matara/
> /Sri Lanka/
>
> Email : wickramagechathurika at rocketmail.com
> chatu at scsio.ac.cn
>
--
Dipl. Geophys. Karin Meier-Fleischer
Visualization Group - NCL, CDO, Python
Application Support
Deutsches Klimarechenzentrum GmbH (DKRZ)
Bundesstrasse 45a - D-20146 Hamburg - Germany
Phone: +49 (0)40 460094 126
Fax: +49 (0)40 460094 270
E-Mail:meier-fleischer at dkrz.de
URL:www.dkrz.de
Geschäftsführer: Prof. Dr. Thomas Ludwig
Sitz der Gesellschaft: Hamburg
Amtsgericht Hamburg HRB 39784
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20211216/334ff45d/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: plot_tauuo.png
Type: image/png
Size: 101300 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20211216/334ff45d/attachment.png>
More information about the ncl-talk
mailing list