[ncl-talk] reordering variable attributes
Tabish Ansari
tabishumaransari at gmail.com
Thu Jan 5 11:51:16 MST 2023
Hi Rick,
Thanks very much for your reply.
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.
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.
Thanks again
and kind regards,
Tabish
-------------------------------
Dr Tabish Ansari
Research Associate
Air Quality Modelling Group
IASS-Potsdam
Germany
On Thu, 5 Jan 2023 at 17:44, Rick Brownrigg <brownrig at ucar.edu> wrote:
> Hi Tabish,
>
> 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?
>
> Here's a thread from a user forum that seems relevant to what you are
> doing -- perhaps there may be a clue in there?
>
>
> https://bb.cgd.ucar.edu/cesm/threads/pio_support-pio_die-myrank-1-error-pionfget_mod-f90-in-some-f-case.5663/
>
> I hope that helps...
> Rick
>
>
> On Thu, Jan 5, 2023 at 9:22 AM Tabish Ansari via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> Hi
>>
>> 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.
>>
>> These are the details of the old file that works:
>> netcdf edgar_HTAPv3_BIGALK_repaired2 {
>> dimensions:
>> time = UNLIMITED ; // (228 currently)
>> lon = 144 ;
>> lat = 96 ;
>> variables:
>> float time(time) ;
>> time:standard_name = "time" ;
>> time:long_name = "time" ;
>> time:units = "days since 1997-01-01 00:00:00+0:00" ;
>> time:calendar = "gregorian" ;
>> time:axis = "T" ;
>> time:_FillValue = 9.96921e+36f ;
>> int date(time) ;
>> date:units = "YYYYMMDD" ;
>> date:long_name = "Date" ;
>> double lon(lon) ;
>> lon:units = "degrees_east" ;
>> lon:long_name = "longitude" ;
>> lon:axis = "X" ;
>> lon:standard_name = "longitude" ;
>> double lat(lat) ;
>> lat:units = "degrees_north" ;
>> lat:long_name = "latitude" ;
>> lat:axis = "Y" ;
>> lat:standard_name = "latitude" ;
>> float HTAPv3_1_International_Shipping(time, lat, lon) ;
>> HTAPv3_1_International_Shipping:units = "molecules/cm2/s"
>> ;
>> HTAPv3_1_International_Shipping:description =
>> "International_Shipping" ;
>> HTAPv3_1_International_Shipping:long_name =
>> "International Shipping" ;
>> HTAPv3_1_International_Shipping:remap = "remapped via
>> ESMF_regrid_with_weights: First-order Conservative" ;
>> HTAPv3_1_International_Shipping:missing_value =
>> 9.96921e+36f ;
>> HTAPv3_1_International_Shipping:_FillValue = 9.96921e+36f
>> ;
>> float HTAPv3_2_1_Domestic_Aviation(time, lat, lon) ;
>> HTAPv3_2_1_Domestic_Aviation:units = "molecules/cm2/s" ;
>> HTAPv3_2_1_Domestic_Aviation:description =
>> "Domestic_Aviation" ;
>> HTAPv3_2_1_Domestic_Aviation:long_name = "Domestic
>> Aviation" ;
>> HTAPv3_2_1_Domestic_Aviation:remap = "remapped via
>> ESMF_regrid_with_weights: First-order Conservative" ;
>> HTAPv3_2_1_Domestic_Aviation:missing_value = 9.96921e+36f
>> ;
>> HTAPv3_2_1_Domestic_Aviation:_FillValue = 9.96921e+36f ;
>> }
>>
>> and here are the details of the new file that doesn't work:
>>
>> dimensions:
>> time = UNLIMITED ; // (228 currently)
>> lon = 144 ;
>> lat = 96 ;
>> variables:
>> float time(time) ;
>> time:standard_name = "time" ;
>> time:long_name = "time" ;
>> time:units = "days since 1997-01-01 00:00:00+0:00" ;
>> time:calendar = "gregorian" ;
>> time:axis = "T" ;
>> time:_FillValue = 9.96921e+36f ;
>> int date(time) ;
>> date:units = "YYYYMMDD" ;
>> date:long_name = "Date" ;
>> double lon(lon) ;
>> lon:units = "degrees_east" ;
>> lon:long_name = "longitude" ;
>> lon:axis = "X" ;
>> lon:standard_name = "longitude" ;
>> double lat(lat) ;
>> lat:units = "degrees_north" ;
>> lat:long_name = "latitude" ;
>> lat:axis = "Y" ;
>> lat:standard_name = "latitude" ;
>> float HTAPv3_1_International_Shipping(time, lat, lon) ;
>> HTAPv3_1_International_Shipping:remap = "remapped via
>> ESMF_regrid_with_weights: First-order Conservative" ;
>> HTAPv3_1_International_Shipping:long_name =
>> "International Shipping" ;
>> HTAPv3_1_International_Shipping:description =
>> "International_Shipping" ;
>> HTAPv3_1_International_Shipping:units = "molecules/cm2/s"
>> ;
>> HTAPv3_1_International_Shipping:missing_value =
>> 9.96921e+36f ;
>> HTAPv3_1_International_Shipping:_FillValue = 9.96921e+36f
>> ;
>> float HTAPv3_2_1_Domestic_Aviation(time, lat, lon) ;
>> HTAPv3_2_1_Domestic_Aviation:remap = "remapped via
>> ESMF_regrid_with_weights: First-order Conservative" ;
>> HTAPv3_2_1_Domestic_Aviation:long_name = "Domestic
>> Aviation" ;
>> HTAPv3_2_1_Domestic_Aviation:description =
>> "Domestic_Aviation" ;
>> HTAPv3_2_1_Domestic_Aviation:units = "molecules/cm2/s" ;
>> HTAPv3_2_1_Domestic_Aviation:missing_value = 9.96921e+36f
>> ;
>> HTAPv3_2_1_Domestic_Aviation:_FillValue = 9.96921e+36f ;
>> float HTAPv3_2_2_International_Aviation(time, lat, lon) ;
>> HTAPv3_2_2_International_Aviation:remap = "remapped via
>> ESMF_regrid_with_weights: First-order Conservative" ;
>> HTAPv3_2_2_International_Aviation:long_name =
>> "International Aviation" ;
>> HTAPv3_2_2_International_Aviation:description =
>> "International_Aviation" ;
>> HTAPv3_2_2_International_Aviation:units =
>> "molecules/cm2/s" ;
>> HTAPv3_2_2_International_Aviation:missing_value =
>> 9.96921e+36f ;
>> HTAPv3_2_2_International_Aviation:_FillValue =
>> 9.96921e+36f ;
>> }
>>
>>
>> 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?
>>
>> Thank you
>>
>> with kind regards,
>>
>> Tabish
>>
>> -------------------------------
>> Dr Tabish Ansari
>> Research Associate
>> Air Quality Modelling Group
>> IASS-Potsdam
>> Germany
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at mailman.ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20230105/4209cd41/attachment.htm>
More information about the ncl-talk
mailing list