Skip to content
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion var/da/da_buoy/da_get_innov_vector_buoy.inc
Original file line number Diff line number Diff line change
Expand Up @@ -134,7 +134,7 @@ subroutine da_get_innov_vector_buoy( it,num_qcstat_conv, grid, ob, iv)
! Surface data assimilation approach 2
! -----------------------------------

! 1.2.1 Surface assmiilation approach 2(10-m u, v, 2-m t, q,
! 1.2.1 Surface assimilation approach 2(10-m u, v, 2-m t, q,
! and sfc_p)

call da_interp_lin_2d (grid%xb%u10, iv%info(buoy), 1,model_u)
Expand Down
2 changes: 1 addition & 1 deletion var/da/da_buoy/da_transform_xtoy_buoy.inc
Original file line number Diff line number Diff line change
Expand Up @@ -71,7 +71,7 @@ subroutine da_transform_xtoy_buoy (grid, iv, y)
deallocate (vb)

else if (sfc_assi_options == sfc_assi_options_2) then
! [2.0] Surface assmiilation approach 2
! [2.0] Surface assimilation approach 2
call da_transform_xtopsfc(grid,iv,buoy,iv%buoy(:),y%buoy(:))
end if

Expand Down
2 changes: 1 addition & 1 deletion var/da/da_main/da_solve.inc
Original file line number Diff line number Diff line change
Expand Up @@ -531,7 +531,7 @@
read(cvt_unit) cvt
close(cvt_unit)
else
write(unit=message(1),fmt='(a)') "cvt file '"//trim(cvtfile)//"' does not exists, initiallizing cvt."
write(unit=message(1),fmt='(a)') "cvt file '"//trim(cvtfile)//"' does not exists, initializing cvt."
call da_message(message(1:1))
call da_initialize_cv (cv_size, cvt)
end if
Expand Down
2 changes: 1 addition & 1 deletion var/da/da_main/da_update_firstguess.inc
Original file line number Diff line number Diff line change
@@ -1,7 +1,7 @@
subroutine da_update_firstguess(grid, out_filename)

!---------------------------------------------------------------------------
! Purpose: Only replce the fields touched by WRFDA, other than re-generate
! Purpose: Only replace the fields touched by WRFDA, rather than re-generate
! whole wrfvar_output from the scratch.
!
! Update : jliu@ucar.edu Apr 1, 2013
Expand Down
2 changes: 1 addition & 1 deletion var/da/da_main/da_wrfvar_main.f90
Original file line number Diff line number Diff line change
Expand Up @@ -23,7 +23,7 @@ program da_wrfvar_main

implicit none

! Split initialisation into 2 parts so we can start and stop trace here
! Split initialization into 2 parts so we can start and stop trace here
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

There are zillions of British spellings of initialisation within DA. I'd say leave these alone. Here's a grep in the var directory ...

find . -name *.inc -exec grep -i initialis {} ; -print > foo.txt

! Initialise columns k+1:n to columns of the unit matrix
./external/lapack/dorg2r.inc
! Initialise columns 1:n-k to columns of the unit matrix
./external/lapack/dorg2l.inc
! [1.0]: Initialise:
./da/da_tools/da_1d_eigendecomposition.inc
! [1.0] Initialise components of innovation vector:
./da/da_obs/da_setup_pseudo_obs.inc
! Initialise obs error factors (which will be overwritten in use_obs_errfac)
./da/da_obs/da_fill_obs_structures.inc
! Initialise obs error factors (which will be overwritten in use_obs_errfac)
./da/da_obs/da_fill_obs_structures_rain.inc
! [1.0] Initialise random number sequence:
./da/da_obs/da_random_omb_all.inc
! Initialise obs error factors (which will be overwritten in use_obs_errfac)
./da/da_obs/da_fill_obs_structures_radar.inc
! [2.0] Initialise components of innovation vector:
./da/da_ssmi/da_get_innov_vector_ssmt1.inc
! [3.0] Calculate rv, rf at OBS location and initialise components of &
./da/da_radar/da_get_innov_vector_radar.inc
! [2.0] Initialise components of innovation vector:
./da/da_rain/da_get_innov_vector_rain.inc
! [1.0] Initialise bins:
./da/da_statistics/da_data_distribution.inc
! [1.0] Initialise:
./da/da_define_structures/da_allocate_background_errors.inc
! Purpose: Initialises the Y-array
./da/da_define_structures/da_zero_y.inc
! [2.0] Initialise components of innovation vector:
./da/da_sound/da_get_innov_vector_sound.inc
! [1.0] Initialise:
./da/da_recursive_filter/da_recursive_filter_1d_adj.inc
! [1.0]: Initialise:
! [1.0]: Initialise:
./da/da_recursive_filter/da_transform_through_rf_adj_dual_res.inc
! [1.0]: Initialise:
! [1.0]: Initialise:
./da/da_recursive_filter/da_transform_through_rf_adj.inc
! [1.0]: Initialise:
./da/da_recursive_filter/da_transform_through_rf_dual_res.inc
! [1.0]: Initialise:
./da/da_recursive_filter/da_transform_through_rf.inc
! [1.0]: Initialise:
./da/da_recursive_filter/da_calculate_rf_factors.inc
! [1.0] Initialise:
./da/da_recursive_filter/da_apply_rf_1v.inc
! [1.0] Initialise:
./da/da_recursive_filter/da_recursive_filter_1d.inc
! [2.0] Initialise components of innovation vector:
./da/da_tamdar/da_get_innov_vector_tamdar.inc
! [1.0] Initialise:
./da/da_setup_structures/da_get_vertical_truncation.inc
! 5.0 Initialise control variable space components of header:
./da/da_setup_structures/da_setup_be_regional.inc
! [1.0] Initialise:
./da/da_setup_structures/da_rescale_background_errors.inc
! 4.0 Initialise control variable space components of header:
./da/da_setup_structures/da_setup_be_nmm_regional.inc
! [1.0] Initialise:
./da/da_test/da_check_cvtovv_adjoint.inc
! [1.0] Initialise:
./da/da_test/da_check_vvtovp_adjoint.inc
! [1.0] Initialise:
./da/da_test/da_check_vtox_adjoint.inc
! [1.0]: Initialise:
./da/da_test/da_check_balance.inc
! [1.0] Initialise:
./da/da_test/da_check_xtoy_adjoint.inc
! [1.0]: initialise:
./da/da_test/da_setup_testfield.inc
! [1.0] Initialise:
./da/da_test/da_check_vptox_adjoint.inc
! initialise index:
./da/da_tracing/da_trace_int_sort.inc
! Purpose: Initialise tracing
./da/da_tracing/da_trace_init.inc
! initialise index:
./da/da_tracing/da_trace_real_sort.inc
! [2.0] Initialise components of innovation vector:
./da/da_mtgirs/da_get_innov_vector_mtgirs.inc
! [2.0] Initialise wrfvar parameters:
./da/da_main/da_solve.inc
! [2.0] Initialise components of innovation vector:
./da/da_profiler/da_get_innov_vector_profiler.inc
! [2.0] Initialise components of innovation vector:
./da/da_pilot/da_get_innov_vector_pilot.inc
! [2.0] Initialise components of innovation vector:
./da/da_bogus/da_get_innov_vector_bogus.inc
! [2.0] Initialise components of innovation vector:
./da/da_airsr/da_get_innov_vector_airsr.inc
! [2.0] Initialise components of innovation vector:
./da/da_airep/da_get_innov_vector_airep.inc
! [1.0] Initialise:
./da/da_ffts/da_solve_poissoneqn_fct.inc
! [1.0] Initialise:
./da/da_ffts/da_solve_poissoneqn_fct_adj.inc
! [1.0] Initialise:
./da/da_ffts/da_solve_poissoneqn_fst_adj.inc
! [1.0] Initialise:
./da/da_ffts/da_solve_poissoneqn_fst.inc
! 0.0 initialise QC by flags assuming good obs
./da/da_radiance/da_qc_iasi.inc
! 0.0 initialise QC by flags assuming good obs
./da/da_radiance/da_qc_airs.inc
! 0.0 initialise QC by flags assuming good obs
./da/da_radiance/da_qc_hirs.inc
! 0.0 initialise QC by flags assuming good obs
./da/da_radiance/da_qc_atms.inc
! 0.0 initialise QC by flags assuming good obs
./da/da_radiance/da_qc_mwts.inc
! 0.0 initialise QC by flags assuming good obs
./da/da_radiance/da_qc_seviri.inc
! 0.0 initialise QC by flags assuming good obs
./da/da_radiance/da_qc_amsua.inc
! 0.0 initialise QC by flags assuming good obs
./da/da_radiance/da_qc_amsr2.inc
! 0.0 initialise QC flags by assuming good obs
./da/da_radiance/da_qc_mwhs.inc
! 0.0 initialise QC flags by assuming good obs
./da/da_radiance/da_qc_mhs.inc
gooddata_map(i,j) = .TRUE. ! Initialised for step 6)
./da/da_radiance/ATMS_Spatial_Average.inc
! 0.0 initialise QC by flags assuming good obs
./da/da_radiance/da_qc_ssmis.inc
! 0.0 initialise QC flags by assuming good obs
./da/da_radiance/da_qc_amsub.inc
x = 0.5*(x1 + x2) ! Initialise the guess for the root
./da/da_spectral/da_get_gausslats.inc
! Initialise Power spectrum
./da/da_spectral/da_test_spectral.inc
! Initialise the varibles
./da/da_gpspw/da_transform_xtoy_gpsztd_adj.inc
! Initialise the varibles
./da/da_gpspw/da_transform_xtoy_gpspw_adj.inc
! [1.0]: Initialise:
./da/da_vtox_transforms/da_check_eof_decomposition.inc
! [1.0] Initialise:
./da/da_vtox_transforms/da_transform_bal.inc
! [1.0] Initialise:
./da/da_vtox_transforms/da_transform_bal_adj.inc
! [2.0] Initialise components of innovation vector:
./da/da_metar/da_get_innov_vector_metar.inc
! [1.0] Initialise:
./da/da_dynamics/da_wpec_constraint_adj.inc
! [1.0] Initialise i and set multiplicative constants
./da/da_dynamics/da_wpec_constraint.inc
! [1.0] Initialise i and set multiplicative constants
./da/da_dynamics/da_wpec_constraint_lin.inc
! [1.0] Initialise:
./da/da_dynamics/da_uv_to_divergence.inc
! [1.0] Initialise iand set Multipilactive constants
./da/da_dynamics/da_balance_equation_lin.inc
! [1.0] Initialise:
./da/da_dynamics/da_wpec_constraint_geoterm.inc
! [1.0] Initialise:
./da/da_dynamics/da_uv_to_divergence_adj.inc
! [1.0] Initialise:
./da/da_dynamics/da_balance_cycloterm.inc
! [1.0] Initialise:
./da/da_dynamics/da_wpec_constraint_cycloterm.inc
! [1.0] Initialise:
./da/da_dynamics/da_uv_to_vorticity.inc
! [1.0] Initialise:
./da/da_dynamics/da_balance_cycloterm_lin.inc
! [1.0] Initialise:
./da/da_dynamics/da_balance_equation_adj.inc
! [1.0] Initialise:
./da/da_dynamics/da_balance_geoterm.inc
! [1.0] Initialise:
./da/da_dynamics/da_balance_cycloterm_adj.inc
! [1.0] initialise:
./da/da_physics/da_tp_to_qs.inc
! [1.0] initialise:
./da/da_physics/da_tp_to_qs_lin1.inc
! [1.0] initialise:
./da/da_physics/da_tprh_to_q_adj.inc
! [1.0] initialise:
! [1.0] initialise:
./da/da_physics/da_tp_to_qs1.inc
! [1.0] initialise:
./da/da_physics/da_prho_to_t_adj.inc
! [1.0] initialise:
./da/da_physics/da_tp_to_qs_adj.inc
! Purpose: Initialises the Y-array
./da/da_minimisation/da_calculate_j.inc
! [2.0] Initialise components of innovation vector:
./da/da_qscat/da_get_innov_vector_qscat.inc
! [4.0] Allocate and initialise residual structure:
./obsproc/MAP_plot/Dir_map/DA_Setup_Obs_Structures.inc

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@davegill I changed that one because in the beginning of that same subroutine (a few lines earlier) the American spelling is used. But I can leave this one alone I guess.


call disable_quilting

Expand Down
2 changes: 1 addition & 1 deletion var/da/da_metar/da_transform_xtoy_metar.inc
Original file line number Diff line number Diff line change
Expand Up @@ -76,7 +76,7 @@ subroutine da_transform_xtoy_metar (grid, iv, y)
deallocate (vb)

else if (sfc_assi_options == sfc_assi_options_2) then
! [2.0] Surface assmiilation approach 2
! [2.0] Surface assimilation approach 2
call da_transform_xtopsfc(grid,iv,metar,iv%metar(:),y%metar(:))
end if

Expand Down
2 changes: 1 addition & 1 deletion var/da/da_minimisation/da_write_diagnostics.inc
Original file line number Diff line number Diff line change
@@ -1,7 +1,7 @@
subroutine da_write_diagnostics(it, grid,num_qcstat_conv, ob, iv, re, y, j)

!---------------------------------------------------------------------------
! Purpose: Output data assmiilation diagnostics
! Purpose: Output data assimilation diagnostics
!---------------------------------------------------------------------------

implicit none
Expand Down
2 changes: 1 addition & 1 deletion var/da/da_physics/da_transform_xtopsfc.inc
Original file line number Diff line number Diff line change
Expand Up @@ -36,7 +36,7 @@ subroutine da_transform_xtopsfc(grid, iv, obs_index, synop, y_synop)
allocate (t(1,iv%info(obs_index)%n1:iv%info(obs_index)%n2))
allocate (q(1,iv%info(obs_index)%n1:iv%info(obs_index)%n2))

! 2.0 Surface assmiilation approach 2 (2-m t and q, 10-m u and v)
! 2.0 Surface assimilation approach 2 (2-m t and q, 10-m u and v)

call da_interp_lin_2d(grid%xa%u10, iv%info(obs_index), 1, u)
call da_interp_lin_2d(grid%xa%v10, iv%info(obs_index), 1, v)
Expand Down
10 changes: 5 additions & 5 deletions var/da/da_setup_structures/da_setup_be_ncep_gfs.inc
Original file line number Diff line number Diff line change
Expand Up @@ -182,7 +182,7 @@ subroutine da_setup_be_ncep_gfs( grid, be )
ALLOCATE ( be % wgvz(its:ite,jts:jte,kts:kte) )
!

! 2, load the WFR model latitude and map factor:
! 2, load the WRF model latitude and map factor:

#ifdef DM_PARALLEL
ij = ( ide- ids+1)*( jde- jds+1)
Expand Down Expand Up @@ -420,7 +420,7 @@ subroutine da_setup_be_ncep_gfs( grid, be )
end do
end do

! potentail velocity:
! potential velocity:
do k=kts,kte
do j=jts,jte
do i=its,ite
Expand Down Expand Up @@ -578,7 +578,7 @@ if ( max_ext_its > 1 ) then
endif
enddo
!
! Horizotal scales:
! Horizontal scales:
if (rootproc) then
write (be_rf_unit) hwll
xsum = sum (hwll*hwll)
Expand Down Expand Up @@ -630,7 +630,7 @@ if ( max_ext_its > 1 ) then
! write (*,'(5x,i3,1x,a,2x,"sum^2=",e20.12)') ii, vname, xsum
enddo
!
! Horizotal scales:
! Horizontal scales:
! write (*,'(3x,a)') 'HORIZONTAL SCALES:'
write (be_rf_unit) hwll
xsum = sum (hwll*hwll)
Expand Down Expand Up @@ -680,7 +680,7 @@ if ( max_ext_its > 1 ) then
write (be_print_unit,'(5x,i3,1x,a,2x,"sum^2=",e20.12)') ii, vname, xxsum(1)
enddo

! Horizotal scales:
! Horizontal scales:
if (rootproc) &
write (be_print_unit,'(3x,a)') 'READ IN THE HORIZONTAL SCALES:'
if (rootproc) then
Expand Down
2 changes: 1 addition & 1 deletion var/da/da_ships/da_get_innov_vector_ships.inc
Original file line number Diff line number Diff line change
Expand Up @@ -131,7 +131,7 @@ subroutine da_get_innov_vector_ships( it,num_qcstat_conv, grid, ob, iv)
else if (sfc_assi_options == sfc_assi_options_2) then
! Surface data assimilation approach 2

! 1.2.1 Surface assmiilation approach 2(10-m u, v, 2-m t, q, and
! 1.2.1 Surface assimilation approach 2(10-m u, v, 2-m t, q, and
! sfc_p)

call da_interp_lin_2d (grid%xb%u10, iv%info(ships), 1, model_u)
Expand Down
2 changes: 1 addition & 1 deletion var/da/da_ships/da_transform_xtoy_ships.inc
Original file line number Diff line number Diff line change
Expand Up @@ -67,7 +67,7 @@ subroutine da_transform_xtoy_ships (grid, iv, y)
deallocate (ub)
deallocate (vb)
else if (sfc_assi_options == sfc_assi_options_2) then
! [2.0] Surface assmiilation approach 2
! [2.0] Surface assimilation approach 2
call da_transform_xtopsfc(grid,iv,ships,iv%ships(:),y%ships(:))
end if

Expand Down
2 changes: 1 addition & 1 deletion var/da/da_synop/da_get_innov_vector_synop.inc
Original file line number Diff line number Diff line change
Expand Up @@ -132,7 +132,7 @@ subroutine da_get_innov_vector_synop( it,num_qcstat_conv, grid, ob, iv)
! Surface data assimilation approach 2
!------------------------------------

! 1.2.1 Surface assmiilation approach 2(10-m u, v, 2-m t, q, and sfc_p)
! 1.2.1 Surface assimilation approach 2(10-m u, v, 2-m t, q, and sfc_p)

call da_interp_lin_2d (grid%xb%u10, iv%info(synop), 1, model_u(1,:))
call da_interp_lin_2d (grid%xb%v10, iv%info(synop), 1, model_v(1,:))
Expand Down
2 changes: 1 addition & 1 deletion var/da/da_tools/da_intpsfc_tem.inc
Original file line number Diff line number Diff line change
Expand Up @@ -5,7 +5,7 @@ subroutine da_intpsfc_tem (val, ho, po, to, height, pressure, temp, kts, kte)
!
! Reference:
! ---------
! The use of surface observations in four dimensional data assmiilation
! The use of surface observations in four dimensional data assimilation
! Using a mesoscale model.
! Ruggiero et al., 1996, Monthly Weather Review, Volume 124, 1018-1033
!
Expand Down
4 changes: 2 additions & 2 deletions var/da/da_transfer_model/da_transfer_wrftltoxa_adj.inc
Original file line number Diff line number Diff line change
@@ -1,7 +1,7 @@
subroutine da_transfer_wrftltoxa_adj(grid, config_flags, filnam, timestr)

!---------------------------------------------------------------------------
! Purpose: Convert analysis increments into WRFAD increments
! Purpose: Convert analysis increments into WRFDA increments
! (following xatowrf, but only keep the increments)
!---------------------------------------------------------------------------

Expand Down Expand Up @@ -176,7 +176,7 @@ subroutine da_transfer_wrftltoxa_adj(grid, config_flags, filnam, timestr)
grid%xa%v(is:ie, js:je, ks:ke) = 0.0

!---------------------------------------------------------------------------
! [4.0] CONVERT THETA inCREMENTS TO T inCREMENTS
! [4.0] CONVERT THETA INCREMENTS TO T INCREMENTS
!---------------------------------------------------------------------------

! In the inverse, g_ph information is lost. This should be investigated later!!!
Expand Down
6 changes: 3 additions & 3 deletions var/da/da_transfer_model/da_transfer_wrftoxb.inc
Original file line number Diff line number Diff line change
Expand Up @@ -68,7 +68,7 @@ subroutine da_transfer_wrftoxb(xbx, grid, config_flags)
!rizvi xbx%big_header%bhi(8,5) = grid%start_hour

!---------------------------------------------------------------------------
! WRF-specific fitelds:
! WRF-specific fields:
!---------------------------------------------------------------------------

ptop = grid%p_top
Expand Down Expand Up @@ -387,7 +387,7 @@ subroutine da_transfer_wrftoxb(xbx, grid, config_flags)
(base_pres*(aln+albn)))**cpovcv
! total density
grid%xb%rho(i,j,k)= 1.0 / (albn+aln)
! pressure purtubation:
! pressure perturbation:
grid%p(i,j,k) = grid%xb%p(i,j,k) - ppb
end do
end do
Expand Down Expand Up @@ -510,7 +510,7 @@ subroutine da_transfer_wrftoxb(xbx, grid, config_flags)
! Convert to specific humidity from mixing ratio of water vapor:
grid%xb%q(i,j,k)=grid%xb%q(i,j,k)/(1.0+grid%xb%q(i,j,k))

! Background qrn needed for radar radial velocity assmiilation:
! Background qrn needed for radar radial velocity assimilation:

if (size(grid%moist,dim=4) >= 4) then
grid%xb%qcw(i,j,k) = max(grid%moist(i,j,k,p_qc), 0.0)
Expand Down
2 changes: 1 addition & 1 deletion var/da/da_vtox_transforms/da_transform_xtoxa.inc
Original file line number Diff line number Diff line change
@@ -1,7 +1,7 @@
subroutine da_transform_xtoxa(grid)

!--------------------------------------------------------------------------------------------------
! Purpose: Transfers fields from WRF TL fields to dignostic fields needed by observational operators
! Purpose: Transfers fields from WRF TL fields to diagnostic fields needed by observational operators
! Updated for Analysis on Arakawa-C grid
! Author: Syed RH Rizvi, MMM/ESSL/NCAR, Date: 10/22/2008
!--------------------------------------------------------------------------------------------------
Expand Down
2 changes: 1 addition & 1 deletion var/da/da_vtox_transforms/da_transform_xtoxa_adj.inc
Original file line number Diff line number Diff line change
@@ -1,7 +1,7 @@
subroutine da_transform_xtoxa_adj(grid)

!--------------------------------------------------------------------------------------------------
! Purpose: Transfers fields from WRF TL fields to dignostic fields needed by observational operators
! Purpose: Transfers fields from WRF TL fields to diagnostic fields needed by observational operators
! Adjoint
!--------------------------------------------------------------------------------------------------

Expand Down
2 changes: 1 addition & 1 deletion var/da/da_vtox_transforms/da_transform_xtoxa_adj_all.inc
Original file line number Diff line number Diff line change
@@ -1,7 +1,7 @@
subroutine da_transform_xtoxa_adj_all(grid)

!--------------------------------------------------------------------------------------------------
! Purpose: Transfers fields from WRF TL fields to dignostic fields needed by observational operators
! Purpose: Transfers fields from WRF TL fields to diagnostic fields needed by observational operators
! Adjoint
!--------------------------------------------------------------------------------------------------

Expand Down
4 changes: 2 additions & 2 deletions var/da/da_vtox_transforms/da_transform_xtoxa_all.inc
Original file line number Diff line number Diff line change
@@ -1,7 +1,7 @@
subroutine da_transform_xtoxa_all(grid)

!--------------------------------------------------------------------------------------------------
! Purpose: Transfers fields from WRF TL fields to dignostic fields needed by observational operators
! Purpose: Transfers fields from WRF TL fields to diagnostic fields needed by observational operators
! Updated for Analysis on Arakawa-C grid
! Author: Syed RH Rizvi, MMM/ESSL/NCAR, Date: 10/22/2008
!--------------------------------------------------------------------------------------------------
Expand Down Expand Up @@ -71,7 +71,7 @@ subroutine da_transform_xtoxa_all(grid)

call da_pt_to_rho_lin(grid)

!by lixin update perturbated geoh
!by lixin update perturbed geoh

do j=jts,jte
do i=its,ite
Expand Down