Search found 77 matches

by bartbrashers
Wed Mar 24, 2021 12:50 am
Forum: Special Running Options
Topic: MYNN2.5 with CFL errors at the edge of nested domains
Replies: 1
Views: 413

MYNN2.5 with CFL errors at the edge of nested domains

I have been testing the MYNN2.5 PBL scheme in stable conditions (wintertime in Alaska) for several projects, some successfully, some not. Using WRF-4.2.2, with nested 12, 4, and 1.33 km domains, I see lots of CFL errors in the "margins", but never anywhere else. My d02 is 202 x 202 grid points, and ...
by bartbrashers
Wed Mar 24, 2021 12:31 am
Forum: Special Running Options
Topic: Adaptive time_step in very stable conditions
Replies: 3
Views: 480

Re: Adaptive time_step in very stable conditions

For a different test using WRF-4.2.2, I used this block: use_adaptive_time_step = .true., step_to_output_time = .true., adaptation_domain = 3, target_cfl = 0.6, 0.6, 0.6, target_hcfl = 0.6, 0.6, 0.6, max_step_increase_pct = 3, 51, 51, starting_time_step = 6, 2, 1, max_time_step = 36, 12, 4, min_time...
by bartbrashers
Wed Mar 24, 2021 12:23 am
Forum: Special Running Options
Topic: Adaptive time_step in very stable conditions
Replies: 3
Views: 480

Adaptive time_step in very stable conditions

I'm running WRF-4.1.5 using 12/4/1.33 km nested domains to simulate central Alaska in the wintertime. There are strong surface-based inversions observed from time to time, and I would like to use adaptive time stepping so it can speed up when possible. From my namelist.input: use_adaptive_time_step ...
by bartbrashers
Tue Mar 23, 2021 11:06 pm
Forum: Special Running Options
Topic: Array out of bounds in WRF-4.2.2 using OBS nudging
Replies: 7
Views: 652

Re: Array out of bounds in WRF-4.2.2 using OBS nudging

I have traced this to what appears to be a bad sounding. I have debug_level = 500, so it's quite verbose. Many "pages" back, it was giving me a hint: # grep -B5 problem rsl.error.0000 rsl.error.0000-d03 2019-12-09_12:00:00 input_wrf: end, fid = 9 rsl.error.0000-d03 2019-12-09_12:00:00 Input data pro...
by bartbrashers
Wed Feb 17, 2021 3:51 pm
Forum: metgrid
Topic: (RESOLVED) WPS using OI SST NetCDF data
Replies: 10
Views: 1498

Re: (RESOLVED) WPS using OI SST NetCDF data

If you are using https://github.com/bbrashers/WPS-ghrsst-to-intermediate , then I wrote that program, not Ming Chen. Is that the program you mean? Where did you download the data from? Useful so other users can download and try themselves. What other steps to debug have you tried? Generally speaking...
by bartbrashers
Sat Feb 13, 2021 8:03 pm
Forum: metgrid
Topic: metgrid crashes: get_min(): No items left in the heap
Replies: 23
Views: 2411

Re: metgrid crashes: get_min(): No items left in the heap

I never run with constant SST, so I don't have much advice for you.

Perhaps your SST:2014-07-05_00 file has a field named LANDSEA in it? You don't want it to. I suggest you use ghrsst-to-intermediate to create a static SST_MASK file, and give that to metgrid via the constants_name path as well.
by bartbrashers
Thu Feb 11, 2021 12:49 am
Forum: Special Running Options
Topic: Array out of bounds in WRF-4.2.2 using OBS nudging
Replies: 7
Views: 652

Re: Array out of bounds in WRF-4.2.2 using OBS nudging

Ping! Any other suggestions anyone has as to how to debug this?

Has anyone else hit an out-of-bounds error in phys/module_fddaobs_rtfdda.F?
by bartbrashers
Wed Feb 10, 2021 12:20 am
Forum: metgrid
Topic: metgrid crashes: get_min(): No items left in the heap
Replies: 23
Views: 2411

Re: metgrid crashes: get_min(): No items left in the heap

Try plotting just your SST_MASK variable. ghtsst-to-intermediate will make one in a separate file. If you append the mask field to an SST file, then you can't see it with int2nc - there's a bug, so only the first one is shown. Or there's a bug in my program.
by bartbrashers
Tue Feb 09, 2021 5:26 pm
Forum: metgrid
Topic: metgrid crashes: get_min(): No items left in the heap
Replies: 23
Views: 2411

Re: metgrid crashes: get_min(): No items left in the heap

All that warning means is that when creating SST_MASK, it can't find a mask to apply.

I think you can safely ignore it. Do the SKINTEMP fields look OK, especially near the land-sea boundaries?
by bartbrashers
Sat Jan 30, 2021 10:13 pm
Forum: Special Running Options
Topic: Array out of bounds in WRF-4.2.2 using OBS nudging
Replies: 7
Views: 652

Re: Array out of bounds in WRF-4.2.2 using OBS nudging

I am using WPS-4.2, the force_use_old_data = true is just a habit.

I have confirmed the out-of-bounds error happens exactly the same (same thread, same time, etc.) with that option removed from the namelist.

Any more hints about how to solve this?
by bartbrashers
Thu Jan 28, 2021 6:24 pm
Forum: Special Running Options
Topic: Array out of bounds in WRF-4.2.2 using OBS nudging
Replies: 7
Views: 652

Re: Array out of bounds in WRF-4.2.2 using OBS nudging

I'm using both of these settings, so it should be backward compatible with the changes to the defaults in WRF-4.0. &time_control force_use_old_data = .true., &dynamics hybrid_opt = 0, Good thought, though. I will try with WRF-3.8.1 (configured with the -D debug switch) and let you know.
by bartbrashers
Wed Jan 27, 2021 2:28 am
Forum: Special Running Options
Topic: Array out of bounds in WRF-4.2.2 using OBS nudging
Replies: 7
Views: 652

Array out of bounds in WRF-4.2.2 using OBS nudging

I'm using WRF-4.2.2 configured with the -D option (debug compilation), along with output from OBSGRID-3.8 (the OBSDOMAIN* files) for OBS nudging. I'm also using the metoa_em* files, FWIW. I've used a similar setup for years (without the debugging part). I am currently using debug_level = 300, and ru...
by bartbrashers
Tue Jan 26, 2021 6:47 pm
Forum: metgrid
Topic: (RESOLVED) WPS using OI SST NetCDF data
Replies: 10
Views: 1498

Re: (RESOLVED) WPS using OI SST NetCDF data

The CSH I was referring to is https://github.com/bbrashers/WPS-ghrsst ... ghrsst.csh, so already posted.
by bartbrashers
Tue Jan 26, 2021 6:44 pm
Forum: metgrid
Topic: (RESOLVED) WPS using OI SST NetCDF data
Replies: 10
Views: 1498

Re: (RESOLVED) WPS using OI SST NetCDF data

The data from JPL's PODAAC (the GHRSST project) are all once-daily. Because UNGRIB normally interpolates in time, you need to interpolate this to your planned interval (6h or 3h, or whatever) as a separate step. You can use https://github.com/bbrashers/WPS-interp-intermediate to interpolate the file...
by bartbrashers
Sun Jan 24, 2021 5:42 pm
Forum: wrf.exe
Topic: SegFault in MYNNSFC
Replies: 4
Views: 235

Re: SegFault in MYNNSFC

Adding a check for (fx2-fx1) being too small (if smaller than 1.e-6, set to 1.e-6) made WRF continue running past the point it failed in 2 previous test runs, so I think that's the culprit.

What's next? Do you want me to file a bug report on the Github site?
by bartbrashers
Fri Jan 22, 2021 6:12 pm
Forum: wrf.exe
Topic: SegFault in MYNNSFC
Replies: 4
Views: 235

Re: SegFault in MYNNSFC

Thanks for the reply. I made a "debug" version as instructed, and ran it. The tail of rsl.error.0006 (running on compute node c08) shows: d03 2019-11-29_12:10:04+07/25 DEBUG wrf_timetoa(): returning with str = [2019-11-29_12:10:04] d03 2019-11-29_12:10:04+07/25 call radiation_driver d03 2019-11-29_1...
by bartbrashers
Sun Jan 17, 2021 10:36 pm
Forum: wrf.exe
Topic: SegFault in MYNNSFC
Replies: 4
Views: 235

SegFault in MYNNSFC

I'm running a 12/4/1.33 km simulation in central Alaska in the wintertime (T2 ranges from -40 degC to -5 degC). For a different WRF domain covering the North Slope, the MYNN2.5 PBL scheme worked very well, so I'd like to try it in this case too. With both WRF-4.2 and WRF-4.2.2, I get a Segmentation ...
by bartbrashers
Wed Jan 13, 2021 7:01 pm
Forum: ungrib
Topic: SST data converted from netCDF to grib2 format - metadata scrambled
Replies: 3
Views: 903

Re: SST data converted from netCDF to grib2 format - metadata scrambled

I wrote a program that will read GHRSST (the Global High Resolution Sea Surface Temperature project) files that come in a self-consistent netCDF format, and convert them to WPS Intermediate file format version 5. It does not convert them to GRIB. The idea here is that we write a program to replace U...
by bartbrashers
Tue Jul 07, 2020 1:47 pm
Forum: wrf.exe
Topic: ACM2 PBL: frequent RIBX never exceeds RIC
Replies: 2
Views: 629

Re: ACM2 PBL: frequent RIBX never exceeds RIC

(1) For the the very beginning of the WRF run one should set pxlsm_smois_init =0 and pxlsm_soil_nudge =1. Subsequent 5-day runs that re-start from the end of the first 5-day init should have pxlsm_smois_init =1 and pxlsm_soil_nudge =1. But my first init only runs for part of a day before crashing wi...
by bartbrashers
Sun Jul 05, 2020 9:30 pm
Forum: Meteorological Input Data
Topic: smooth_cg_topo and ERA-5 / ERA-Interim data
Replies: 3
Views: 712

Re: smooth_cg_topo and ERA-5 / ERA-Interim data

As soon as I clicked Submit, I thought of a better answer - the system probably uses the output from geogrid.exe, HGT_M, to construct SOILGEO. That's why it still works without downloading 'geopotential' for ERA5.
by bartbrashers
Sun Jul 05, 2020 9:23 pm
Forum: Meteorological Input Data
Topic: smooth_cg_topo and ERA-5 / ERA-Interim data
Replies: 3
Views: 712

Re: smooth_cg_topo and ERA-5 / ERA-Interim data

For ECMWF and ERA, the Vtable says it uses grib code 129 level 0 for both SOILGEO and SOILHGT. I built some invariant files for ERA-I years ago, but my ERA-5 download script doesn't include anything like that. I have been following https://dreambooker.site/2019/10/03/Initializing-the-WRF-model-with-...
by bartbrashers
Sun Jul 05, 2020 8:23 pm
Forum: wrf.exe
Topic: ACM2 PBL: frequent RIBX never exceeds RIC
Replies: 2
Views: 629

ACM2 PBL: frequent RIBX never exceeds RIC

I'm trying to test the effect of using soil moisture nudging in Australia using the PX LSM and ACM2 PBL scheme. I have two runs that are identical expect for the values of pxlsm_smois_init and pxlsm_soil_nudge: # grep -e bl_pbl -e sf_sfclay -e sf_surface SM_nudge/2017-10-18/namelist.input.wrf sf_sur...
by bartbrashers
Fri Jun 26, 2020 10:48 pm
Forum: Miscellaneous DA Topics
Topic: WRFDA reanalyses for analysis nudging?
Replies: 0
Views: 608

WRFDA reanalyses for analysis nudging?

This is a very obvious use-case to me, but I can't find any evidence that someone has done this yet. For hindcasting (simulating the past, for multi-year periods) we often use analysis nudging. For domains with finer grid spacing than the input dataset (which is most domains for me) we objectively r...
by bartbrashers
Sat Jun 20, 2020 12:28 am
Forum: metgrid
Topic: (RESOLVED) LANDSEA mask best practices
Replies: 12
Views: 1763

Re: LANDSEA mask best practices

A summary of LAND-SEA Mask best practices, as I understand them, just to close out this thread. When using an SST dataset that's finer resolution than the rest of your input data (e.g. ERA5 or GFS or ...) it's best to supply a land-sea mask. Do NOT name the field "LANDSEA" or it will conflict with t...
by bartbrashers
Fri Jun 12, 2020 7:28 pm
Forum: metgrid
Topic: Cleaning up examples
Replies: 2
Views: 461

Re: Cleaning up examples

There are a number of links (mostly from old "how-to's") which point to https://www2.mmm.ucar.edu/wrf/src/, from where it's easy to find the outdated netcdf-to-intermediate.f file. https://publish.illinois.edu/mesomodel/2016/01/20/lecture-1-getting-started-compiling-gordon-infrastructure/ https://st...

Go to advanced search