<div dir="ltr">Hello!<div><br></div><div>I wanted to plot a contour map of lat v/s time for Precipitation anomaly over certain longitude. The time coordinates were given in "Hours since 1900-01-01 00:00:00". I converted them to YYYYMMDD using 'cd_calendar' and saved it as a new time coordinates. When i tried plotting, i got the following error message</div><div>------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------</div><div>warning:IrTransInitialize: trXCoordPoints contains invalid coordinate array: defaulting trXAxisType to LinearAxis</div><div>warning:_NhlCreateSplineCoordApprox: Attempt to create spline approximation for X axis failed: consider adjusting trXTensionF value</div><div>warning:_NhlCreateSplineCoordApprox: Attempt to create spline approximation for X axis failed: consider adjusting trXTensionF value</div><div>warning:_NhlCreateSplineCoordApprox: Attempt to create spline approximation for X axis failed: consider adjusting trXTensionF value</div><div>fatal:IrTransInitialize: spline coordinate approximation not possible:</div><div>fatal:Unable to initialize layer-Can't Create</div><div>fatal:ContourPlotInitialize: Error creating transformation object</div><div>fatal:ContourPlotInitialize: error setting up transformation</div><div>fatal:Unable to initialize layer-Can't Create</div><div>fatal:Unable to access object with id:-4</div><div>fatal:PID #-4 can't be found in NhlSetValues</div><div>fatal:NhlGetValues:PID #-4 is invalid</div><div>fatal:["Execute.c":8640]:Execute: Error occurred at or near line 2222 in file /apps/codes/gnu/ncl_stable/lib/ncarg/nclscripts/csm/gsn_code.ncl</div><div><br></div><div>fatal:["Execute.c":8640]:Execute: Error occurred at or near line 9571 in file /apps/codes/gnu/ncl_stable/lib/ncarg/nclscripts/csm/gsn_code.ncl</div><div><br></div><div>fatal:["Execute.c":8640]:Execute: Error occurred at or near line 13393 in file /apps/codes/gnu/ncl_stable/lib/ncarg/nclscripts/csm/gsn_csm.ncl</div><div><br></div><div>fatal:["Execute.c":8640]:Execute: Error occurred at or near line 224 in file pract8.ncl </div><div>-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------</div><div><br></div><div>My data includes only one year of daily precipitation data with specifications as below.</div><div><br></div><div>----------------------------------------------------------------------------------------------------------------------------</div><div><div>Variable: p</div><div>Type: float</div><div>Total Size: 52560 bytes</div><div> 13140 values</div><div>Number of Dimensions: 3</div><div>Dimensions and sizes:<span style="white-space:pre"> </span>[time | 365] x [lat | 6] x [lon | 6]</div><div>Coordinates: </div><div> time: [990552..999288]</div><div> lat: [28.75..31.25]</div><div> lon: [78.25..80.75]</div><div>Number Of Attributes: 14</div><div> _FillValue :<span style="white-space:pre"> </span>-9.96921e+36</div><div> var_desc :<span style="white-space:pre"> </span>Precipitation</div><div> valid_range :<span style="white-space:pre"> </span>( 0, 1000 )</div><div> units :<span style="white-space:pre"> </span>mm</div><div> statistic :<span style="white-space:pre"> </span>Total</div><div> parent_stat :<span style="white-space:pre"> </span>Other</div><div> missing_value :<span style="white-space:pre"> </span>-9.96921e+36</div><div> long_name :<span style="white-space:pre"> </span>Daily total of precipitation</div><div> level_desc :<span style="white-space:pre"> </span>Surface</div><div> dataset :<span style="white-space:pre"> </span>CPC Global Precipitation</div><div> cell_methods :<span style="white-space:pre"> </span>time: sum</div><div> avg_period :<span style="white-space:pre"> </span>0000-00-01 00:00:00</div><div> actual_range :<span style="white-space:pre"> </span>( 0, 473.6278 )</div><div> DIMENSION_LIST :<span style="white-space:pre"> </span>-2147483647</div></div><div>----------------------------------------------------------------------------------------------------------------------------------</div><div><br></div><div><b>How do i plot contour map of lat v/s time with time coordinates in YYYYMMDD format for a certain time period? </b></div><div><br></div><div>My code for the above is as follows:</div><div><br></div><div>---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------</div><div><div>read = 1</div><div>if (read.ne.0) then</div><div>f = addfile("<a href="http://precip_2013.nc">precip_2013.nc</a>","r")</div><div>p = f->precip(:,{28.75:31.25},{78.25:80.75})</div><div>pr = p(lat|:,lon|:,time|:) </div><div>prm = dim_avg_Wrap(pr)</div><div>prnew = new((/6,6,365/),float)</div><div>t = f->time</div><div>tim = cd_calendar(t,2)</div><div><br></div><div>;computing anomoly at all points</div><div>do i = 0,5</div><div> do j = 0,5</div><div> do k = 0,364</div><div> prnew(i,j,k) = (pr(i,j,k)-prm(i,j))</div><div> end do </div><div> end do</div><div>end do</div><div>prnew!0 = "lat"</div><div>prnew!1 = "lon"</div><div>prnew!2 = "time"</div><div>prn = prnew(time|:,lat|:,lon|:)</div><div><br></div><div>;applying fft over anomoly at all points</div><div>ft = new((/2,365,6,6/),float)</div><div>do h = 0,5</div><div> do k = 0,5</div><div> gt = cfftf(prn(:,h,k),0.0,0)</div><div> ft(0,:,h,k) = gt(0,:)</div><div> ft(1,:,h,k) = gt(1,:)</div><div> end do</div><div>end do</div><div><br></div><div>;filtering out unwanted frequencies</div><div>do h = 0,5</div><div> do k = 0,5</div><div> ft(0,13:352,h,k) = 0.0</div><div> ft(1,13:352,h,k) = 0.0</div><div> ft(0,360:,h,k) = 0.0</div><div> ft(1,360:,h,k) = 0.0</div><div> ft(0,0:5,h,k) = 0.0</div><div> ft(1,0:5,h,k) = 0.0</div><div> end do</div><div>end do</div><div><br></div><div>;ifft of required frequencies</div><div>ift = new((/365,6,6/),float)</div><div>do h = 0,5</div><div> do k = 0,5</div><div> kt = cfftb(ft(:,:,h,k),1)</div><div> ift(:,h,k) = kt(:)</div><div> end do</div><div>end do</div><div>ift!0 = "time"</div><div>ift!1 = "lat"</div><div>ift!2 = "lon"</div><div>iftnew = ift(lat|:,lon|:,time|:)</div><div><br></div><div><br></div><div>copy_VarAtts(pr,iftnew)</div><div>copy_VarCoords(pr,iftnew)</div><div>oxn = iftnew(lon|:,lat|:,time|:)</div><div>oxn@long_name = "30-60 day anomoly of precipitation"</div><div>oxn&time = tim</div><div>wks = gsn_open_wks ("png", "lat_time" ) ; send graphics to PNG file</div><div><br></div><div> res = True ; plot mods desired</div><div> res@cnFillOn = True ; turn on color fill</div><div> res@cnFillPalette = "BlWhRe" ; set color map</div><div> res@cnLinesOn = False</div><div> res@tiMainString = "Uttarkhand at 79.25E" ; title</div><div> res@tiYAxisString = "Latitude"</div><div> res@cnLevelSelectionMode = "AutomaticLevels" ; manual contour levels</div><div> res@lbOrientation = "Vertical"</div><div> </div><div> plot = gsn_csm_contour(wks, oxn({79.25},:,160:170), res) </div><div> </div><div>end if</div></div><div>----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------</div><div><div><div class="gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><b>By,</b></div><div><b>Dilip.V</b></div><div><b>Master's Student,</b></div><div><b>Indian Institute of Science Education and Research,</b></div><div><b>Pune, India</b></div></div></div></div></div></div>
</div></div>