Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

[hwsd2-wrf] ERROR: Error in ext_pkg_write_field #1

Open
SurAurora opened this issue Sep 16, 2024 · 0 comments
Open

[hwsd2-wrf] ERROR: Error in ext_pkg_write_field #1

SurAurora opened this issue Sep 16, 2024 · 0 comments

Comments

@SurAurora
Copy link

Hi jobonaf,

I hope you're doing well. My name is Aurora, and I'm using your data from the hwsd2-wrf repository. I wanted to thank you for sharing this project.

However, I'm having some trouble understanding how to convert the TIFF files to WPS binary data. In particular, I am using the following command:

SoilType_i1 $ convert_geotiff -c 16 -b 3 -w 2 /the-path/SoilType_depth000to030cm_i1.tif

This process seems to finish successfully, creating many files with names like "00601-00700.00501-00600" along with an "index" file. Here is the content of the index file:

projection = regular_ll
known_x = 1
known_y = 504
known_lat = -31.691668
known_lon = -71.283333
dx = 8.332824e-03
dy = 8.332443e-03
type = categorical
signed = yes
units = "NO UNITS"
description = "NO DESCRIPTION"
wordsize = 2
tile_x = 100
tile_y = 100
tile_z = 1
category_min = 1
category_max = 17
tile_bdr = 3
missing_value = 17.000000
scale_factor = 1.000000
row_order = bottom_top
endian = little

I modified the GEOGRID.TBL to point to this directory and updated the 'geog_data_res' field in the namelist.wps to use this data option. However, when running geogrid, it doesn't finish successfully and exits with the following error:

ERROR: Error in ext_pkg_write_field

I suspect I might be using the convert_geotiff command incorrectly, and I was hoping you could help me figure out what's going wrong. Below is my namelist.wps file:

&share
wrf_core = 'ARW'
max_dom = 1
start_date = '2024-08-00_03:00:00'
end_date = '2024-08-00_18:00:00'
interval_seconds = 21600
io_form_geogrid = 2
debug_level = 0
/

&geogrid
parent_id = 1
parent_grid_ratio = 1
i_parent_start = 1
j_parent_start = 1
e_we = 32
e_sn = 28
geog_data_res = 'CLsoil'
dx = 5000
dy = 5000
map_proj = 'lambert'
ref_lat = -33.652
ref_lon = -69.399
truelat1 = -33.720
truelat2 = -33.720
pole_lat = 90
pole_lon = 0
stand_lon = -69.291
geog_data_path = '/home/proyecto_mma/DATA/ARCHIVOS_MMA/WRF/geog/'
opt_geogrid_tbl_path = './geogrid/'
/

&ungrib
out_format = 'WPS'
prefix = 'FILE'
/

&metgrid
fg_name = 'FILE'
io_form_metgrid = 2
opt_metgrid_tbl_path = './metgrid'
/
Additionally, here is a portion of my GEOGRID.TBL:

name=SOILCTOP
priority=1
dest_type=categorical
z_dim_name=soil_cat
dominant=SCT_DOM
interp_option = CLsoil:nearest_neighbor
rel_path= CLsoil:SoilType_i1/

Any advice or guidance on how to resolve this issue would be greatly appreciated. Thank you in advance for your time and help.

Best regards,
Aurora

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant