<div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Hi Rick,</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Thanks very much for your reply.</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">I have indeed used the ESMF regrid function in NCL to generate these files. The raw files were at a different resolution. They look OK when viewed through ncdump -h or ncview. But the model somehow doesn't accept them. </div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">I does accept the file which only contains 2 variables, whose header is shown in blue in the previous email. But when I add a third variable to it and write it out as a new file, I get the error again. Between the two files, the only difference I could spot was the order of attributes. But I will look further into it. </div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Thanks again</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">and kind regards,</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Tabish</div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><font face="trebuchet ms, sans-serif" color="#666666">-------------------------------<br><span style="background-color:rgb(255,255,255)">Dr Tabish Ansari</span></font></div><div><font face="trebuchet ms, sans-serif" color="#666666">Research Associate </font></div><div><font face="trebuchet ms, sans-serif" color="#666666">Air Quality Modelling Group</font></div><div><span style="background-color:rgb(255,255,255)"><span style="font-weight:normal"><font face="trebuchet ms, sans-serif" color="#666666">IASS-Potsdam </font></span></span></div><div><font face="trebuchet ms, sans-serif" color="#666666">Germany</font></div></div></div></div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, 5 Jan 2023 at 17:44, Rick Brownrigg <<a href="mailto:brownrig@ucar.edu">brownrig@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi Tabish,</div><div><br></div><div>I would be very surprised if the attribute ordering was the issue. The only reference to those symbols PIO_SUPPORT/PIO_DIE in NCL that I can find are coming from the ESMF library. That library in turn integrates a PIO package (parallel I/O for NetCDF). Are you attempting to regrid the file?<br></div><div><br></div><div>Here's a thread from a user forum that seems relevant to what you are doing -- perhaps there may be a clue in there?</div><div><br></div><div><a href="https://bb.cgd.ucar.edu/cesm/threads/pio_support-pio_die-myrank-1-error-pionfget_mod-f90-in-some-f-case.5663/" target="_blank">https://bb.cgd.ucar.edu/cesm/threads/pio_support-pio_die-myrank-1-error-pionfget_mod-f90-in-some-f-case.5663/</a></div><div><br></div><div>I hope that helps...</div><div>Rick</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jan 5, 2023 at 9:22 AM Tabish Ansari via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif">Hi</div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif"><br></div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif">I'm getting an error "pio_support::pio_die:: myrank= -1 : ERROR" in the CAM-Chem model when using new emission files. The only difference between the new files and the old ones that work smoothly seems to be the order of variable attributes.</div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif"><br></div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif">These are the details of the old file that works:</div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif"><font color="#0000ff">netcdf edgar_HTAPv3_BIGALK_repaired2 {<br>dimensions:<br>        time = UNLIMITED ; // (228 currently)<br>        lon = 144 ;<br>        lat = 96 ;<br>variables:<br>        float time(time) ;<br>                time:standard_name = "time" ;<br>                time:long_name = "time" ;<br>                time:units = "days since 1997-01-01 00:00:00+0:00" ;<br>                time:calendar = "gregorian" ;<br>                time:axis = "T" ;<br>                time:_FillValue = 9.96921e+36f ;<br>        int date(time) ;<br>                date:units = "YYYYMMDD" ;<br>                date:long_name = "Date" ;<br>        double lon(lon) ;<br>                lon:units = "degrees_east" ;<br>                lon:long_name = "longitude" ;<br>                lon:axis = "X" ;<br>                lon:standard_name = "longitude" ;<br>        double lat(lat) ;<br>                lat:units = "degrees_north" ;<br>                lat:long_name = "latitude" ;<br>                lat:axis = "Y" ;<br>                lat:standard_name = "latitude" ;<br>        float HTAPv3_1_International_Shipping(time, lat, lon) ;<br>                HTAPv3_1_International_Shipping:units = "molecules/cm2/s" ;<br>                HTAPv3_1_International_Shipping:description = "International_Shipping" ;<br>                HTAPv3_1_International_Shipping:long_name = "International Shipping" ;<br>                HTAPv3_1_International_Shipping:remap = "remapped via ESMF_regrid_with_weights: First-order Conservative" ;<br>                HTAPv3_1_International_Shipping:missing_value = 9.96921e+36f ;<br>                HTAPv3_1_International_Shipping:_FillValue = 9.96921e+36f ;<br>        float HTAPv3_2_1_Domestic_Aviation(time, lat, lon) ;<br>                HTAPv3_2_1_Domestic_Aviation:units = "molecules/cm2/s" ;<br>                HTAPv3_2_1_Domestic_Aviation:description = "Domestic_Aviation" ;<br>                HTAPv3_2_1_Domestic_Aviation:long_name = "Domestic Aviation" ;<br>                HTAPv3_2_1_Domestic_Aviation:remap = "remapped via ESMF_regrid_with_weights: First-order Conservative" ;<br>                HTAPv3_2_1_Domestic_Aviation:missing_value = 9.96921e+36f ;<br>                HTAPv3_2_1_Domestic_Aviation:_FillValue = 9.96921e+36f ;<br>}<br></font></div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif"><br></div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif">and here are the details of the new file that doesn't work:</div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif"><br></div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif"><font color="#ff0000">dimensions:<br>        time = UNLIMITED ; // (228 currently)<br>        lon = 144 ;<br>        lat = 96 ;<br>variables:<br>        float time(time) ;<br>                time:standard_name = "time" ;<br>                time:long_name = "time" ;<br>                time:units = "days since 1997-01-01 00:00:00+0:00" ;<br>                time:calendar = "gregorian" ;<br>                time:axis = "T" ;<br>                time:_FillValue = 9.96921e+36f ;<br>        int date(time) ;<br>                date:units = "YYYYMMDD" ;<br>                date:long_name = "Date" ;<br>        double lon(lon) ;<br>                lon:units = "degrees_east" ;<br>                lon:long_name = "longitude" ;<br>                lon:axis = "X" ;<br>                lon:standard_name = "longitude" ;<br>        double lat(lat) ;<br>                lat:units = "degrees_north" ;<br>                lat:long_name = "latitude" ;<br>                lat:axis = "Y" ;<br>                lat:standard_name = "latitude" ;<br>        float HTAPv3_1_International_Shipping(time, lat, lon) ;<br>                HTAPv3_1_International_Shipping:remap = "remapped via ESMF_regrid_with_weights: First-order Conservative" ;<br>                HTAPv3_1_International_Shipping:long_name = "International Shipping" ;<br>                HTAPv3_1_International_Shipping:description = "International_Shipping" ;<br>                HTAPv3_1_International_Shipping:units = "molecules/cm2/s" ;<br>                HTAPv3_1_International_Shipping:missing_value = 9.96921e+36f ;<br>                HTAPv3_1_International_Shipping:_FillValue = 9.96921e+36f ;<br>        float HTAPv3_2_1_Domestic_Aviation(time, lat, lon) ;<br>                HTAPv3_2_1_Domestic_Aviation:remap = "remapped via ESMF_regrid_with_weights: First-order Conservative" ;<br>                HTAPv3_2_1_Domestic_Aviation:long_name = "Domestic Aviation" ;<br>                HTAPv3_2_1_Domestic_Aviation:description = "Domestic_Aviation" ;<br>                HTAPv3_2_1_Domestic_Aviation:units = "molecules/cm2/s" ;<br>                HTAPv3_2_1_Domestic_Aviation:missing_value = 9.96921e+36f ;<br>                HTAPv3_2_1_Domestic_Aviation:_FillValue = 9.96921e+36f ;<br>        float HTAPv3_2_2_International_Aviation(time, lat, lon) ;<br>                HTAPv3_2_2_International_Aviation:remap = "remapped via ESMF_regrid_with_weights: First-order Conservative" ;<br>                HTAPv3_2_2_International_Aviation:long_name = "International Aviation" ;<br>                HTAPv3_2_2_International_Aviation:description = "International_Aviation" ;<br>                HTAPv3_2_2_International_Aviation:units = "molecules/cm2/s" ;<br>                HTAPv3_2_2_International_Aviation:missing_value = 9.96921e+36f ;<br>                HTAPv3_2_2_International_Aviation:_FillValue = 9.96921e+36f ;<br>}<br></font></div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif"><br></div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif"><br></div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif">Therefore, I want to change the order of the variable attributes to that in the old file. Is there a way I can achieve this in NCL?</div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif"><br></div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif">Thank you</div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif"><br></div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif">with kind regards,</div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif"><br></div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif">Tabish</div><div class="gmail_default" style="font-family:"trebuchet ms",sans-serif"><br></div><div><div dir="ltr"><div dir="ltr"><div><font face="trebuchet ms, sans-serif" color="#666666">-------------------------------<br><span style="background-color:rgb(255,255,255)">Dr Tabish Ansari</span></font></div><div><font face="trebuchet ms, sans-serif" color="#666666">Research Associate </font></div><div><font face="trebuchet ms, sans-serif" color="#666666">Air Quality Modelling Group</font></div><div><span style="background-color:rgb(255,255,255)"><span style="font-weight:normal"><font face="trebuchet ms, sans-serif" color="#666666">IASS-Potsdam </font></span></span></div><div><font face="trebuchet ms, sans-serif" color="#666666">Germany</font></div></div></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div>
</blockquote></div>