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.

Sign up for GitHub

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

Jump to bottom

Create an initial wind field for mass-consistent wind downscaling #36

Closed
janmandel opened this issue Feb 1, 2021 · 6 comments
Closed

Create an initial wind field for mass-consistent wind downscaling #36

janmandel opened this issue Feb 1, 2021 · 6 comments

Comments

@janmandel
Copy link
Member

janmandel commented Feb 1, 2021 •
edited

Interpolate U and V to a 3D boundary layer at fire mesh resolution and given heights.
janmandel added a commit that referenced this issue Feb 1, 2021
@janmandel
starting branch develop-36 from develop #36
57bf778
@janmandel
Copy link
Member Author

WRF Tutorial 2004 https://www.powershow.com/viewfl/536f0a-YTQ2N/Tutorial_Notes_WRF_Software_powerpoint_ppt_presentation

Slides from Winter 2021 Basic WRF Tutorial https://www2.mmm.ucar.edu/wrf/users/tutorial/tutorial_presentation_winter_2021.htm

janmandel added a commit that referenced this issue Feb 8, 2021
@janmandel
adding registry entries for fmw packages #36
ebbadfc
janmandel added a commit that referenced this issue Feb 8, 2021
@janmandel
adding fmw variables to check_grid_alloc #36
8566a13
janmandel added a commit that referenced this issue Feb 12, 2021
@janmandel
fixing fmw variables in registry #36
605661f
janmandel added a commit that referenced this issue Feb 12, 2021
@janmandel
check_size -> check_shape #36
3162ea3
@janmandel
Copy link
Member Author

janmandel commented Feb 12, 2021 •
edited

If I define arrays in registry like this:
dimspec num_fmw - namelist=nfmw z "number of fire wind vertical layers" "1"
rconfig integer nfmw namelist,fire 1 20 - "NFMW" "number of fire wind vertical layers" "1"
state real v_fmw i{num_fmw}j fire 1 - h "V_FMW" "fire wind v" "m/s"
state real u0_fmw i{num_fmw}j fire 1 - h "U0_FMW" "initial fire wind u" "m/s"
state real v0_fmw i{num_fmw}j fire 1 - h "V0_FMW" "initial fire wind v" "m/s"
state real ht_fmw i{num_fmw}j fire 1 z h "HT_FMW" "heights of the bottom of fire wind layers above the terrain" "m/s"
state real lambda_fmw *i{num_fmw}*j fire 1 - hr "LAMBDA_FMW" "langrange multiplier" "m^3/s"

I get NetCDF errors and only the first variable V_FMW is present in wrfout:
d01 0001-01-01_00:00:00 open_hist_w : opening wrfout_d01_0001-01-01_00:00:00 for writing.
d01 0001-01-01_00:00:00 NetCDF error: NetCDF: String match to name in use
d01 0001-01-01_00:00:00 NetCDF error in wrf_io.F90, line 2601
d01 0001-01-01_00:00:00 NetCDF error: NetCDF: String match to name in use
d01 0001-01-01_00:00:00 NetCDF error in wrf_io.F90, line 2601
d01 0001-01-01_00:00:00 NetCDF error: NetCDF: String match to name in use
d01 0001-01-01_00:00:00 NetCDF error in wrf_io.F90, line 2601
d01 0001-01-01_00:00:00 NetCDF error: NetCDF: String match to name in use
d01 0001-01-01_00:00:00 NetCDF error in wrf_io.F90, line 2601
d01 0001-01-01_00:00:00 NetCDF error: NetCDF: String match to name in use
d01 0001-01-01_00:00:00 NetCDF error in wrf_io.F90, line 2601
d01 0001-01-01_00:00:00 Information: NOFILL being set for writing to wrfout_d01_0001-01-01_00:00:00
d01 0001-01-01_00:00:00 Warning LENGTH != DRY RUN LENGTH for |U0_FMW | dim 3 in wrf_io.F90, line 2725
d01 0001-01-01_00:00:00 LENGTH 5 DRY RUN LENGTH 0
d01 0001-01-01_00:00:00 Warning LENGTH != DRY RUN LENGTH for |V0_FMW | dim 3 in wrf_io.F90, line 2725
d01 0001-01-01_00:00:00 LENGTH 5 DRY RUN LENGTH 0
d01 0001-01-01_00:00:00 Warning LENGTH != DRY RUN LENGTH for |HT_FMW | dim 3 in wrf_io.F90, line 2725
d01 0001-01-01_00:00:00 LENGTH 5 DRY RUN LENGTH 0
d01 0001-01-01_00:00:00 Warning LENGTH != DRY RUN LENGTH for |LAMBDA_FMW | dim 3 in wrf_io.F90, line 2725
d01 0001-01-01_00:00:00 LENGTH 5 DRY RUN LENGTH 0
d01 0001-01-01_00:00:00 Warning LENGTH != DRY RUN LENGTH for |U_FMW | dim 3 in wrf_io.F90, line 2725
d01 0001-01-01_00:00:00 LENGTH 5 DRY RUN LENGTH 0
d01 0001-01-01_00:00:00 med_hist_out: opened wrfout_d01_0001-01-01_00:00:00 as DATASET=HISTORY
Timing for Writing wrfout_d01_0001-01-01_00:00:00 for domain 1: 0.46010 elapsed seconds

I need this on the fire mesh ("submesh") but it happens with the usual atmosphere mesh also as seen here.

Running in test/em_fire/firewind 8440103

janmandel added a commit that referenced this issue Feb 12, 2021
@janmandel
testing 3D arrays in registry #36
8440103 
it seems arrays with defined k dimension give the io system trouble for some reason
@janmandel
Copy link
Member Author

janmandel commented Feb 12, 2021 •
edited

How to create 4d arrays not needed for output: wrf-model#528
Examples in WRF/Registry/registry.dimspec
Latest registry tutorial https://www2.mmm.ucar.edu/wrf/users/tutorial/202101/gill_registry.pdf

janmandel added a commit that referenced this issue Feb 14, 2021
@janmandel
making _fmw arrays in wrfout though at atm resolution #36
c3019c8
janmandel added a commit that referenced this issue Feb 14, 2021
@janmandel
writing 3D fire arrays to wrfout #36
bdb80cd
janmandel added a commit that referenced this issue Feb 18, 2021
@janmandel
interpolating to u0_fmw v0_fmw #36 #5 openwfm/wrf-fire-matlab#2
681005e
janmandel added a commit that referenced this issue Feb 21, 2021
@janmandel
generating u0_fmw v0_fmw at midpoints #36
0579746 
layers heights same at all locations, stored in 1d
janmandel added a commit that referenced this issue Feb 22, 2021
@janmandel
again netcdf errors on first wrfout, cleaning up and debugging #36
879d9c0
@janmandel
Copy link
Member Author

janmandel commented Feb 24, 2021 •
edited

HEAD 879d9c0 w/WRF-SFIRE-femwind: zeros, again netcdf errors, ht_fmw is 1D
HEAD^ 0579746 w/WRF-SFIRE-femwind3: zero arrays errors, ht_fmw is 1D, bad dimspec
HEAD^^ 162f1ad w/WRF-SFIRE-femwind2 makes 3D u0_fmw v0_fmw ht_fmw, but interpolation not finished yet
HEAD^^^ 681005e w/WRF-SFIRE-femwind4 makes 3D u0_fmw v0_fmw ht_fmw, but interpolation not finished yet

@janmandel
Copy link
Member Author

janmandel commented Feb 25, 2021 •
edited

Branch develop-36-2 b778b8d: generating 3D ht_fmw, memory fix of 879d9c0

janmandel added a commit that referenced this issue Feb 25, 2021
@janmandel
interpolated u0_fmw, v0_fmw in wrfout, ht_fmw is 1D #36 #5 openwfm/wr…
8c31830 
…f-fire-matlab#2
janmandel added a commit that referenced this issue Feb 26, 2021
@janmandel
align femwind initial wind with terrain gradient #36 #5
8964aa4 
cleanup debug print in log interpolation
janmandel added a commit that referenced this issue Mar 14, 2021
@janmandel
prints in interpolation of initial wind, testing fmw run flags on ==1…
9da405d 
… per the registry #36
janmandel added a commit that referenced this issue Jun 3, 2021
@janmandel
Merge branch 'develop-36' into develop-18
1a939c8 
Conflicts:
	phys/module_fr_sfire_driver.F

Clean merge, different lines added. #36 #18
janmandel added a commit that referenced this issue Jun 3, 2021
@janmandel
adding write u0_fmw v0_fmw w0_fmw by pnetcdf; #ifdef PNETCDF all pnet…
21ca893 
…cdf code #18 #36

NOTE: does not support openmp, operates on whole patch inside the parallel do master loop
janmandel added a commit that referenced this issue Jun 23, 2021
@janmandel
bypassing pnetcdf if ifmw_qes=0 #18 #36
4bb3d0b
janmandel added a commit that referenced this issue Jul 2, 2021
@janmandel
adding computation of chsum0_fmw #18 #36
547ddc7 
commenting out scalar pnetcdf write that does not wor:wk
janmandel added a commit that referenced this issue Jul 8, 2021
@janmandel
adding izb in interpolate_wind2fire_height to go all the way to boundary
c5b77d2 
from continued wind and geopotential to prevent zero strip in u0_fwm and v0_fmw at the boundary #36
janmandel added a commit that referenced this issue Jul 23, 2021
@janmandel
adding subroutine pnetcdf_read_arr #36 #18
6b01ed3 
also fixing end subroutine pnetcdf_read_int
@janmandel
Copy link
Member Author

Wrapped into #18

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Assignees
No one assigned
Labels
None yet
Projects
None yet
Milestone
No milestone
Development

No branches or pull requests
1 participant
@janmandel