[ncl-talk] select sub_region

Thu Dec 16 05:21:59 MST 2021

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 name tauuo. 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 list ncl-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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20211216/15323f51/attachment-0001.html>