add_grid

General questions and postings pertaining to the use of HYSPLIT regardless of the platform. For platform specific questions, use the HYSPLIT Platform forums.
Post Reply
federica
Posts: 31
Joined: November 12th, 2014, 5:17 pm
Registered HYSPLIT User: Yes

add_grid

Post by federica »

Hello,

I'm trying to use the add_grid routine to increase the resolution of my meteo data. I'm using a 1deg resolution file and I'd like to obtain a 0.5 deg file.

This is what I'm doing:

Enter meteorological directory ...
./
Enter meteorological file name ...
extract_26568.bin
Enter grid resolution divisor: 1
2
Add one i-coordinate (Y or N) ?
Y
NOTICE pakset:
Number of index records = 1
Number of records /time = 164
NOTICE pakini: start initialization rec = 1
NOTICE pakini: end initialization rec = 164
NOTICE pakini: start initialization rec = 165
NOTICE pakini: end initialization rec = 328
NOTICE pakini: start initialization rec = 329
NOTICE pakini: end initialization rec = 492
NOTICE pakini: start initialization rec = 493
NOTICE pakini: end initialization rec = 656
NOTICE pakini: start initialization rec = 657
NOTICE pakini: end initialization rec = 820
NOTICE pakini: start initialization rec = 821
NOTICE pakini: end initialization rec = 984
NOTICE pakini: start initialization rec = 985
NOTICE pakini: end initialization rec = 1148
NOTICE pakini: start initialization rec = 1149
NOTICE pakini: end initialization rec = 1312
NOTICE pakini: start initialization rec = 1313
NOTICE pakini: end initialization rec = 1476
NOTICE pakini: start initialization rec = 1477
NOTICE pakini: end initialization rec = 1640
NOTICE pakini: start initialization rec = 1641
NOTICE pakini: end initialization rec = 1804
NOTICE pakini: start initialization rec = 1805
NOTICE pakini: end initialization rec = 1968
NOTICE pakini: start initialization rec = 1969
NOTICE pakini: end initialization rec = 2132
NOTICE pakini: start initialization rec = 2133
NOTICE pakini: end initialization rec = 2296
NOTICE pakini: start initialization rec = 2297
NOTICE pakini: end initialization rec = 2460
NOTICE pakini: start initialization rec = 2461
NOTICE pakini: end initialization rec = 2624
NOTICE pakini: start initialization rec = 2625
NOTICE pakini: end initialization rec = 2788
NOTICE pakini: start initialization rec = 2789
NOTICE pakini: end initialization rec = 2952
NOTICE pakini: start initialization rec = 2953
NOTICE pakini: end initialization rec = 3116
NOTICE pakini: start initialization rec = 3117
NOTICE pakini: end initialization rec = 3280
NOTICE pakini: start initialization rec = 3281
NOTICE pakini: end initialization rec = 3444
NOTICE pakini: start initialization rec = 3445
NOTICE pakini: end initialization rec = 3608
NOTICE pakini: start initialization rec = 3609
NOTICE pakini: end initialization rec = 3772
NOTICE pakini: start initialization rec = 3773
NOTICE pakini: end initialization rec = 3936
NOTICE pakini: start initialization rec = 3937
NOTICE pakini: end initialization rec = 4100
NOTICE pakini: start initialization rec = 4101
NOTICE pakini: end initialization rec = 4264
NOTICE pakini: start initialization rec = 4265
NOTICE pakini: end initialization rec = 4428
NOTICE pakini: start initialization rec = 4429
NOTICE pakini: end initialization rec = 4592
NOTICE pakini: start initialization rec = 4593
NOTICE pakini: end initialization rec = 4756
NOTICE pakini: start initialization rec = 4757
NOTICE pakini: end initialization rec = 4920
NOTICE pakini: start initialization rec = 4921
NOTICE pakini: end initialization rec = 5084
NOTICE pakini: start initialization rec = 5085
NOTICE pakini: end initialization rec = 5248
NOTICE pakini: start initialization rec = 5249
NOTICE pakini: end initialization rec = 5412
NOTICE pakini: start initialization rec = 5413
NOTICE pakini: end initialization rec = 5576
NOTICE pakini: start initialization rec = 5577
NOTICE pakini: end initialization rec = 5740
NOTICE pakini: start initialization rec = 5741
NOTICE pakini: end initialization rec = 5904
NOTICE pakini: start initialization rec = 5905
NOTICE pakini: end initialization rec = 6068
NOTICE pakini: start initialization rec = 6069
NOTICE pakini: end initialization rec = 6232
NOTICE pakini: start initialization rec = 6233
NOTICE pakini: end initialization rec = 6396
NOTICE pakini: start initialization rec = 6397
NOTICE pakini: end initialization rec = 6560
NOTICE pakini: start initialization rec = 6561
NOTICE pakini: end initialization rec = 6724
NOTICE pakini: start initialization rec = 6725
NOTICE pakini: end initialization rec = 6888
NOTICE pakini: start initialization rec = 6889
NOTICE pakini: end initialization rec = 7052
NOTICE pakini: start initialization rec = 7053
NOTICE pakini: end initialization rec = 7216
NOTICE pakini: start initialization rec = 7217
NOTICE pakini: end initialization rec = 7380
NOTICE pakini: start initialization rec = 7381
NOTICE pakini: end initialization rec = 7544
NOTICE pakini: start initialization rec = 7545
NOTICE pakini: end initialization rec = 7708
NOTICE pakini: start initialization rec = 7709
NOTICE pakini: end initialization rec = 7872
NOTICE pakini: start initialization rec = 7873
NOTICE pakini: end initialization rec = 8036
NOTICE pakini: start initialization rec = 8037
NOTICE pakini: end initialization rec = 8200
NOTICE pakini: start initialization rec = 8201
NOTICE pakini: end initialization rec = 8364
NOTICE pakini: start initialization rec = 8365
NOTICE pakini: end initialization rec = 8528
NOTICE pakini: start initialization rec = 8529
NOTICE pakini: end initialization rec = 8692
NOTICE pakini: start initialization rec = 8693
NOTICE pakini: end initialization rec = 8856
NOTICE pakini: start initialization rec = 8857
NOTICE pakini: end initialization rec = 9020
NOTICE pakini: start initialization rec = 9021
NOTICE pakini: end initialization rec = 9184

The output files are:

ADDGRID.CFG (attached)
MESSAGE (attached)
addgrid.bin

However, when I check the addgrid.bin file, all the variables are equal to 0 (attached a profile.txt obtained from addgrid.bin)

Thanks for your help
Regards
Federica
Attachments
MESSAGE.txt
(146 Bytes) Downloaded 65 times
ADDGRID.txt
(2.01 KiB) Downloaded 65 times
profile.txt
(2.77 KiB) Downloaded 65 times
barbara.stunder
Posts: 425
Joined: November 9th, 2012, 4:23 pm
Registered HYSPLIT User: Yes

Re: add_grid

Post by barbara.stunder »

It should have worked. Try running chk_file on addgrid.bin to see if the (1,1) values are non-zero given in the ASCII header are non-zero.

If they are zero, have you successfully run HYSPLIT with your 1-degree file? Or run chk_file on your 1-degree file to check it. Clearly your 1-degree file must be correct.

If the (1,1) values of addgrid.bin are not zero and your 1-deg file is alright, then can you put your 1-deg file on an ftpserver where I can get it and do a test run?
federica
Posts: 31
Joined: November 12th, 2014, 5:17 pm
Registered HYSPLIT User: Yes

Re: add_grid

Post by federica »

I've tried again and it worked. I set 2 as resolution divisor and N in Add one i-coordinate.

Is there a threshold on the resolution divisor? I've tried with 10 but it did not work.

Attached the link to the gdas file that I'm using.

https://drive.google.com/open?id=1Mekup ... dJD2bJvdIs

Many thanks for your help
Federica
barbara.stunder
Posts: 425
Joined: November 9th, 2012, 4:23 pm
Registered HYSPLIT User: Yes

Re: add_grid

Post by barbara.stunder »

I tried 10 on your dataset and it similarly failed. In looking at the code at a first glance it seems like a higher divisor number should work, but I do not have the time right now to spend on this. If you revise the code to handle what you are trying to do, we can put it in the repository.

Why are you doing this? The output finer resolution grid is not necessarily better because it is linearly interpolated, as is done in HYSPLIT.
federica
Posts: 31
Joined: November 12th, 2014, 5:17 pm
Registered HYSPLIT User: Yes

Re: add_grid

Post by federica »

Thanks for the reply.

I'm trying to use the add_grid routine because I'd like to improve the spatial resolution of my meteo data.
Alternatively, I'm doing WRF simulations and they seems to work.

Regards
Federica
barbara.stunder
Posts: 425
Joined: November 9th, 2012, 4:23 pm
Registered HYSPLIT User: Yes

Re: add_grid

Post by barbara.stunder »

Taking a 1-degree grid and changing it to be a 0.5 degree grid by interpolation will not show the small scale features that a WRF run will show. The small scale features in the atmosphere do not usually vary linearly in space (and time) as given by the output from add_grid.
Post Reply