From 52d6059185a8d6d86071e9df350fa1fa7b647ea4 Mon Sep 17 00:00:00 2001 From: kdraeder Date: Wed, 10 Jul 2024 09:17:29 -0600 Subject: [PATCH 1/2] Moved inflation file name description from netcdf_inflation_files.rst to inflation.rst --- guide/inflation.rst | 21 +++++++------ guide/netcdf_inflation_files.rst | 52 -------------------------------- 2 files changed, 12 insertions(+), 61 deletions(-) delete mode 100644 guide/netcdf_inflation_files.rst diff --git a/guide/inflation.rst b/guide/inflation.rst index 2ca23137a8..0f8a5c36f4 100644 --- a/guide/inflation.rst +++ b/guide/inflation.rst @@ -90,9 +90,11 @@ In the namelist each entry has two values. The first is for Prior inflation and ``&filter_nml :: inf_initial_from_restart`` *valid values:* .true. or .false. - If true, read the inflation values from an inflation restart file named ``input_{prior,post}inf_mean.nc.`` An initial - run could be done to let spatially-varying inflation values evolve in a spinup phase, and then the saved values can - be read back in and used as fixed values in further runs. Or if time-varying inflation is used, then the restart file + If true, read the inflation values from an inflation restart file named ``input_{prior,post}inf_mean{_d##}.nc``, + (d\#\# is the domain number, if needed). + An initial run could be done to let spatially-varying inflation values evolve in a spinup phase, + and then the saved values can be read back in and used as fixed values in further runs. + Or if time-varying inflation is used, then the restart file from the previous job step must be supplied as an input file for the next step. ``&filter_nml :: inf_initial`` @@ -127,7 +129,7 @@ In the namelist each entry has two values. The first is for Prior inflation and ``&filter_nml :: inf_sd_initial_from_restart`` *valid values:* .true. or .false. - If true, read the inflation standard deviation values from an restart file named ``input_{prior,post}inf_sd.nc.`` See + If true, read the inflation standard deviation values from an restart file named ``input_{prior,post}inf_sd{_d##.nc.`` See the comments above about ``inflation_initial_from_restart``. ``&filter_nml :: inf_sd_initial`` @@ -188,7 +190,8 @@ The suggested procedure for testing inflation options is to start without any (b ``inf_damping`` > 0.). Then enable Prior state space, spatially-varying inflation, with no Posterior inflation (set ``inf_flavor`` to [2, 0]). Then try damped inflation (set ``inf_damping`` to 0.9 and set ``inf_sd_initial`` and ``inf_sd_lower_bound`` to 0.6). The inflation values and standard deviation are written out to files with -``_{prior,post}inf_{mean,sd}`` in their names. These NetCDF files can be viewed with common tools (we often use +``_{prior,post}inf_{mean,sd}{_d##}.nc``, in their names +(d\#\# is the domain number, if needed). These NetCDF files can be viewed with common tools (we often use `ncview `__ ). Expected inflation values are generally in the 1 to 30 range; if values grow much larger than this it usually indicates a problem with the assimilation. @@ -207,11 +210,11 @@ Here's an example of using ncap2 to set the T,U and V inf values: :: - ncap2 -s 'T=1.0;U=1.0;V=1.0' wrfinput_d01 input_priorinf_mean.nc - ncap2 -s 'T=0.6;U=0.6;V=0.6' wrfinput_d01 input_priorinf_sd.nc + ncap2 -s 'T=1.0;U=1.0;V=1.0' wrfinput_d01 input_priorinf_mean_d01.nc + ncap2 -s 'T=0.6;U=0.6;V=0.6' wrfinput_d01 input_priorinf_sd_d01.nc -or- - ncap2 -s 'T(:,:,:)=1.0;U(:,:,:)=1.0;V(:,:,:)=1.0' wrfinput_d01 input_priorinf_mean.nc - ncap2 -s 'T(:,:,:)=0.6;U(:,:,:)=0.6;V(:,:,:)=0.6' wrfinput_d01 input_priorinf_sd.nc + ncap2 -s 'T(:,:,:)=1.0;U(:,:,:)=1.0;V(:,:,:)=1.0' wrfinput_d01 input_priorinf_mean_d01.nc + ncap2 -s 'T(:,:,:)=0.6;U(:,:,:)=0.6;V(:,:,:)=0.6' wrfinput_d01 input_priorinf_sd_d01.nc Some versions of the NCO utilities change the full 3D arrays into a single scalar. If that's your result (check your output with ``ncdump -h``) use the alternate syntax or a more recent version of the NCO tools. diff --git a/guide/netcdf_inflation_files.rst b/guide/netcdf_inflation_files.rst deleted file mode 100644 index b813018771..0000000000 --- a/guide/netcdf_inflation_files.rst +++ /dev/null @@ -1,52 +0,0 @@ -Netcdf Inflation Files -====================== - -The filter_nml now read restart and inflation files directly from NetCDF files - -Netcdf inflation files are no longer special files. DART format inflation files were always 2 copies in one file (mean -and standard devation). Taking away this special status of inflation files has the advantage that all copies (restarts, -ensemble mean, ensemble standard deviation, inflation mean, inflation sd, etc.) can all be treated the same for IO -purposes. Since there are two inflation files when reading/writing netcdf the filenames are different to DART format -restart files. - -The names of the netcdf inflation files are now fixed. - -**Input inflation file names** - -The filter_nml option: - -``inf_in_file_name = prior_inflation_ics, post_inflation_ics`` - -has been **deprecated** and for 1 domain filter is expecting to read: - -| input_{priorinf,postinf}_mean.nc -| input_{priorinf,postinf}_sd.nc - -For multiple domains filter is expecting to read: - -| input_{priorinf,postinf}_mean_d01.nc -| input_{priorinf,postinf}_sd_d01.nc -| input_{priorinf,postinf}_mean_d02.nc -| input_{priorinf,postinf}_sd_d02.nc - -where d0\* is the domain number. - -**Output inflation file names** - -The filter_nml option: - -``inf_out_file_name = prior_inflation_restart, post_inflation_restart`` - -has been **deprecated** and for 1 domain filter is expecting to read: - -| output_{priorinf,postinf}_mean.nc -| output_{priorinf,postinf}_sd.nc - -For multiple domains filter is expecting to write: - -| prior_inflation_restart_mean_d01 -| prior_inflation_restart_sd_d01 -| prior_inflation_restart_mean_d02 -| prior_inflation_restart_sd_d02 - -where d0\* is the domain number. From e3a4f826e2be423684601b672c28d3b8142c764d Mon Sep 17 00:00:00 2001 From: kdraeder Date: Wed, 10 Jul 2024 09:30:12 -0600 Subject: [PATCH 2/2] typo --- guide/inflation.rst | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/guide/inflation.rst b/guide/inflation.rst index 0f8a5c36f4..e26e09672d 100644 --- a/guide/inflation.rst +++ b/guide/inflation.rst @@ -129,8 +129,8 @@ In the namelist each entry has two values. The first is for Prior inflation and ``&filter_nml :: inf_sd_initial_from_restart`` *valid values:* .true. or .false. - If true, read the inflation standard deviation values from an restart file named ``input_{prior,post}inf_sd{_d##.nc.`` See - the comments above about ``inflation_initial_from_restart``. + If true, read the inflation standard deviation values from an restart file named ``input_{prior,post}inf_sd{_d##}.nc.`` + See the comments above about ``inflation_initial_from_restart``. ``&filter_nml :: inf_sd_initial`` *valid values:* ≤ 0.0 to disable evolution of inflation, > 0.0 otherwise