Skip to content

Commit

Permalink
Files impacted:
Browse files Browse the repository at this point in the history 
D    test/em_fire/input_sounding
D    test/em_fire/namelist.input
M    phys/module_physics_init.F
M    phys/module_fddaobs_driver.F
M    phys/module_fddaobs_rtfdda.F
M    phys/module_fr_sfire_driver.F
M    phys/module_fr_sfire_model.F
M    phys/module_fdda_psufddagd.F
M    share/wrf_fddaobs_in.F
M    Registry/Registry.EM_SST
M    Registry/Registry.EM_CHEM
M    Registry/Registry.EM
M    Registry/registry.dimspec
M    Registry/Registry.wrfvar
M    dyn_em/module_initialize_real.F
M    dyn_em/start_em.F


1) send a note out to users concerning max map factors in d01 in the real program
ARW only, no changes in results
M      dyn_em/module_initialize_real.F

2) allow users to get the dominant category for vegetables and soil from the fields in geogrid
ARW only, a new nml option setting, no changes in results with orig values (mods to ARW's GEOGRID.TBL files)
M      dyn_em/module_initialize_real.F

3a) obs nudging fixes from Al Bourgeois
Fix default values for Regime-4 vertical scaling, Add stnidprt, needed for proper obs diagnostic report, Add obsstid for dimensioning stnidprt
M      Registry/Registry.EM_SST
M      Registry/Registry.EM_CHEM
M      Registry/Registry.EM
M      Registry/registry.dimspec
3b) Overhaul diagnostics to fix potential for bad indices into print arrays when handling soundings. Clarify a few user messages.
M      share/wrf_fddaobs_in.F
3c) Add stnidprt, needed for proper obs diagnostic report,
M      phys/module_fddaobs_driver.F
3d) Add individual nudge flags to fdaobs_init, needed to correct/clarify user messages. Fix/improve printouts for vertical influence function settings. Add stnid_prt.
M      phys/module_fddaobs_rtfdda.F
3e) Add individual nudge flags to pass to fddaobs_init
M      phys/module_physics_init.F

4a) Fire updates to allow restarts
ARW only, no changes for results
M      phys/module_fr_sfire_driver.F
M      phys/module_fr_sfire_model.F
4b) remove the wrong namelist.input and sounding for ideal
M      test/em_fire/input_sounding
M      test/em_fire/namelist.input

5) vertically interpolate RH, not Qv
ARW only, changes results, need to use diagnosed Qv (with dry pressure) for initialization, Qv from RH with full pressure added at end of routine
M      dyn_em/module_initialize_real.F

6) if surface values are missing, set use_surface flag to false
ARW only, changes results, changes to ARW METGRID.TBL files
M      dyn_em/module_initialize_real.F

7) initialize half-level hydrostatic pressure for 0-h output, if no restart
ARW only, diagnostic only
M      dyn_em/start_em.F

8) analysis nudging smoother goes out of bounds, reduce indexing bounds by 1 on all sides
ARW only, changes results
M      phys/module_fdda_psufddagd.F

9) insure vertical interpolation array has top pressure value with close vertical spacing
ARW only, NaNs if you hit this, big zap_close_levels and small delta press levels cause trouble
M      dyn_em/module_initialize_real.F

10) input RH to real constrained: if( (p<100 mb) and (rh>0.1%)), then assign rh = 1.e-5%
ARW only, changes results
M      dyn_em/module_initialize_real.F




git-svn-id: https://svn-wrf-model.cgd.ucar.edu/trunk@4254 b0b5d27b-6f0f-0410-a2a3-cb1e977edc3d
  • Loading branch information
@davegill
davegill committed Mar 22, 2010
1 parent 654f844 commit 36e0c26
Show file tree
Hide file tree
Showing 16 changed files with 545 additions and 629 deletions.
19 changes: 11 additions & 8 deletions Registry/Registry.EM
View file
Expand Up @@ -95,6 +95,8 @@ state real tsk_gc ij dyn_em 1 - i1 "SKINTEMP
state real tavgsfc ij dyn_em 1 - i1 "TAVGSFC" "daily mean of surface air temperature" "K"
state real tmn_gc ij dyn_em 1 - i1 "SOILTEMP" "annual mean deep soil temperature" "K"
state real pslv_gc ij dyn_em 1 - i1 "PMSL" "sea level pressure" "Pa"
state real sct_dom_gc ij dyn_em 1 - i1 "SCT_DOM" "Dominant soil (top) category from GEOGRID" "cat"
state real scb_dom_gc ij dyn_em 1 - i1 "SCB_DOM" "Dominant soil (bottom) category from GEOGRID" "cat"
state real greenfrac imj dyn_em 1 Z i1 "GREENFRAC" "monthly greenness fraction" "0 - 1 fraction"
state real albedo12m imj dyn_em 1 Z i1 "ALBEDO12M" "background albedo" "0 - 1 fraction"
state real pd_gc igj dyn_em 1 Z - "PD" "dry pressure" "Pa"
Expand Down Expand Up @@ -1507,20 +1509,20 @@ rconfig real obs_nudgezfullr1_uv namelist,fdda 1
rconfig real obs_nudgezrampr1_uv namelist,fdda 1 50 rh "obs_nudgezrampr1_uv" "Vert infl ramp-to-zero height for LML obs, regime 1, winds" ""
rconfig real obs_nudgezfullr2_uv namelist,fdda 1 50 rh "obs_nudgezfullr2_uv" "Vert infl full weight height for LML obs, regime 2, winds" ""
rconfig real obs_nudgezrampr2_uv namelist,fdda 1 50 rh "obs_nudgezrampr2_uv" "Vert infl ramp-to-zero height for LML obs, regime 2, winds" ""
rconfig real obs_nudgezfullr4_uv namelist,fdda 1 50 rh "obs_nudgezfullr4_uv" "Vert infl full weight height for LML obs, regime 4, winds" ""
rconfig real obs_nudgezrampr4_uv namelist,fdda 1 -5000 rh "obs_nudgezrampr4_uv" "Vert infl ramp-to-zero height for LML obs, regime 4, winds" ""
rconfig real obs_nudgezfullr4_uv namelist,fdda 1 -5000 rh "obs_nudgezfullr4_uv" "Vert infl full weight height for LML obs, regime 4, winds" ""
rconfig real obs_nudgezrampr4_uv namelist,fdda 1 50 rh "obs_nudgezrampr4_uv" "Vert infl ramp-to-zero height for LML obs, regime 4, winds" ""
rconfig real obs_nudgezfullr1_t namelist,fdda 1 50 rh "obs_nudgezfullr1_t" "Vert infl full weight height for LML obs, regime 1, temperature" ""
rconfig real obs_nudgezrampr1_t namelist,fdda 1 50 rh "obs_nudgezrampr1_t" "Vert infl ramp-to-zero height for LML obs, regime 1, temperature" ""
rconfig real obs_nudgezfullr2_t namelist,fdda 1 50 rh "obs_nudgezfullr2_t" "Vert infl full weight height for LML obs, regime 2, temperature" ""
rconfig real obs_nudgezrampr2_t namelist,fdda 1 50 rh "obs_nudgezrampr2_t" "Vert infl ramp-to-zero height for LML obs, regime 2, temperature" ""
rconfig real obs_nudgezfullr4_t namelist,fdda 1 50 rh "obs_nudgezfullr4_t" "Vert infl full weight height for LML obs, regime 4, temperature" ""
rconfig real obs_nudgezrampr4_t namelist,fdda 1 -5000 rh "obs_nudgezrampr4_t" "Vert infl ramp-to-zero height for LML obs, regime 4, temperature" ""
rconfig real obs_nudgezfullr4_t namelist,fdda 1 -5000 rh "obs_nudgezfullr4_t" "Vert infl full weight height for LML obs, regime 4, temperature" ""
rconfig real obs_nudgezrampr4_t namelist,fdda 1 50 rh "obs_nudgezrampr4_t" "Vert infl ramp-to-zero height for LML obs, regime 4, temperature" ""
rconfig real obs_nudgezfullr1_q namelist,fdda 1 50 rh "obs_nudgezfullr1_q" "Vert infl full weight height for LML obs, regime 1, moisture" ""
rconfig real obs_nudgezrampr1_q namelist,fdda 1 50 rh "obs_nudgezrampr1_q" "Vert infl ramp-to-zero height for LML obs, regime 1, moisture" ""
rconfig real obs_nudgezfullr2_q namelist,fdda 1 50 rh "obs_nudgezfullr2_q" "Vert infl full weight height for LML obs, regime 2, moisture" ""
rconfig real obs_nudgezrampr2_q namelist,fdda 1 50 rh "obs_nudgezrampr2_q" "Vert infl ramp-to-zero height for LML obs, regime 2, moisture" ""
rconfig real obs_nudgezfullr4_q namelist,fdda 1 50 rh "obs_nudgezfullr4_q" "Vert infl full weight height for LML obs, regime 4, moisture" ""
rconfig real obs_nudgezrampr4_q namelist,fdda 1 -5000 rh "obs_nudgezrampr4_q" "Vert infl ramp-to-zero height for LML obs, regime 4, moisture" ""
rconfig real obs_nudgezfullr4_q namelist,fdda 1 -5000 rh "obs_nudgezfullr4_q" "Vert infl full weight height for LML obs, regime 4, moisture" ""
rconfig real obs_nudgezrampr4_q namelist,fdda 1 50 rh "obs_nudgezrampr4_q" "Vert infl ramp-to-zero height for LML obs, regime 4, moisture" ""
rconfig real obs_nudgezfullmin namelist,fdda 1 50 rh "obs_nudgezfullmin" "Minimum depth through which vertical influence fcn remains 1.0" "m"
rconfig real obs_nudgezrampmin namelist,fdda 1 50 rh "obs_nudgezrampmin" "Minimum depth through which vertical influence fcn decreases from 1.0 to 0.0" "m"
rconfig real obs_nudgezmax namelist,fdda 1 3000 rh "obs_nudgezmax" "Maximum depth in which vertical influence function is nonzero" "m"
Expand Down Expand Up @@ -1800,10 +1802,10 @@ package dfi_dfl dfi_opt==1 - state:dfi_u,df
package dfi_ddfi dfi_opt==2 - state:dfi_u,dfi_v,dfi_w,dfi_ph,dfi_phb,dfi_ph0,dfi_php,dfi_t,dfi_p,dfi_ww,dfi_mu,dfi_tke,dfi_pb,dfi_al,dfi_alt,dfi_TSLB,dfi_SMOIS,dfi_SNOW,dfi_SNOWH,dfi_CANWAT,dfi_SMFR3D,dfi_KEEPFR3DFLAG,dfi_TSK,dfi_SOILT1,dfi_TSNAV,dfi_SNOWC,dfi_QVG,dfi_rh,dfi_tten_rad
package dfi_tdfi dfi_opt==3 - state:dfi_u,dfi_v,dfi_w,dfi_ph,dfi_phb,dfi_ph0,dfi_php,dfi_t,dfi_p,dfi_ww,dfi_mu,dfi_tke,dfi_pb,dfi_al,dfi_alt,dfi_TSLB,dfi_SMOIS,dfi_SNOW,dfi_SNOWH,dfi_CANWAT,dfi_SMFR3D,dfi_KEEPFR3DFLAG,dfi_TSK,dfi_SOILT1,dfi_TSNAV,dfi_SNOWC,dfi_QVG,dfi_rh,dfi_tten_rad
ifdef RUC_CLOUD
package realonly use_wps_input==1 - state:u_gc,v_gc,t_gc,rh_gc,ght_gc,p_gc,xlat_gc,xlong_gc,ht_gc,tsk_gc,tavgsfc,tmn_gc,pslv_gc,greenfrac,albedo12m,pd_gc,psfc_gc,intq_gc,pdhs,sh_gc,qv_gc,qr_gc,qc_gc,qs_gc,qi_gc,qg_gc,qni_gc
package realonly use_wps_input==1 - state:u_gc,v_gc,t_gc,rh_gc,ght_gc,p_gc,xlat_gc,xlong_gc,ht_gc,tsk_gc,tavgsfc,tmn_gc,pslv_gc,sct_dom_gc,scb_dom_gc,greenfrac,albedo12m,pd_gc,psfc_gc,intq_gc,pdhs,sh_gc,qv_gc,qr_gc,qc_gc,qs_gc,qi_gc,qg_gc,qni_gc
endif
ifndef RUC_CLOUD
package realonly use_wps_input==1 - state:u_gc,v_gc,t_gc,rh_gc,ght_gc,p_gc,xlat_gc,xlong_gc,ht_gc,tsk_gc,tavgsfc,tmn_gc,pslv_gc,greenfrac,albedo12m,pd_gc,psfc_gc,intq_gc,pdhs,sh_gc,qv_gc
package realonly use_wps_input==1 - state:u_gc,v_gc,t_gc,rh_gc,ght_gc,p_gc,xlat_gc,xlong_gc,ht_gc,tsk_gc,tavgsfc,tmn_gc,pslv_gc,sct_dom_gc,scb_dom_gc,greenfrac,albedo12m,pd_gc,psfc_gc,intq_gc,pdhs,sh_gc,qv_gc
endif


Expand Down Expand Up @@ -2006,6 +2008,7 @@ typedef fdob_type real latprt [ - 1 -
typedef fdob_type real lonprt [ - 1 - - "lonprt" "obs longitude for diagnostic printout"
typedef fdob_type real mlatprt [ - 1 - - "mlatprt" "model latitude at obs location"
typedef fdob_type real mlonprt [ - 1 - - "mlonprt" "model longitude at obs location"
typedef fdob_type integer stnidprt {obsstid}{[} - 1 - - "stnidprt" "obs station id for diagnostic printout"
typedef fdob_type real base_state k - 1 - - "base_state" "base-state height on half (mass) levels" "meters"
state fdob_type fdob - -
Expand Down
19 changes: 11 additions & 8 deletions Registry/Registry.EM_CHEM
View file
Expand Up @@ -88,6 +88,8 @@ state real tsk_gc ij dyn_em 1 - i1 "SKINTEMP
state real tavgsfc ij dyn_em 1 - i1 "TAVGSFC" "daily mean of surface air temperature" "K"
state real tmn_gc ij dyn_em 1 - i1 "SOILTEMP" "annual mean deep soil temperature" "K"
state real pslv_gc ij dyn_em 1 - i1 "PMSL" "sea level pressure" "Pa"
state real sct_dom_gc ij dyn_em 1 - i1 "SCT_DOM" "Dominant soil (top) category from GEOGRID" "cat"
state real scb_dom_gc ij dyn_em 1 - i1 "SCB_DOM" "Dominant soil (bottom) category from GEOGRID" "cat"
state real greenfrac imj dyn_em 1 Z i1 "GREENFRAC" "monthly greenness fraction" "0 - 1 fraction"
state real albedo12m imj dyn_em 1 Z i1 "ALBEDO12M" "background albedo" "0 - 1 fraction"
state real pd_gc igj dyn_em 1 Z - "PD" "dry pressure" "Pa"
Expand Down Expand Up @@ -1481,20 +1483,20 @@ rconfig real obs_nudgezfullr1_uv namelist,fdda 1
rconfig real obs_nudgezrampr1_uv namelist,fdda 1 50 rh "obs_nudgezrampr1_uv" "Vert infl ramp-to-zero height for LML obs, regime 1, winds" ""
rconfig real obs_nudgezfullr2_uv namelist,fdda 1 50 rh "obs_nudgezfullr2_uv" "Vert infl full weight height for LML obs, regime 2, winds" ""
rconfig real obs_nudgezrampr2_uv namelist,fdda 1 50 rh "obs_nudgezrampr2_uv" "Vert infl ramp-to-zero height for LML obs, regime 2, winds" ""
rconfig real obs_nudgezfullr4_uv namelist,fdda 1 50 rh "obs_nudgezfullr4_uv" "Vert infl full weight height for LML obs, regime 4, winds" ""
rconfig real obs_nudgezrampr4_uv namelist,fdda 1 -5000 rh "obs_nudgezrampr4_uv" "Vert infl ramp-to-zero height for LML obs, regime 4, winds" ""
rconfig real obs_nudgezfullr4_uv namelist,fdda 1 -5000 rh "obs_nudgezfullr4_uv" "Vert infl full weight height for LML obs, regime 4, winds" ""
rconfig real obs_nudgezrampr4_uv namelist,fdda 1 50 rh "obs_nudgezrampr4_uv" "Vert infl ramp-to-zero height for LML obs, regime 4, winds" ""
rconfig real obs_nudgezfullr1_t namelist,fdda 1 50 rh "obs_nudgezfullr1_t" "Vert infl full weight height for LML obs, regime 1, temperature" ""
rconfig real obs_nudgezrampr1_t namelist,fdda 1 50 rh "obs_nudgezrampr1_t" "Vert infl ramp-to-zero height for LML obs, regime 1, temperature" ""
rconfig real obs_nudgezfullr2_t namelist,fdda 1 50 rh "obs_nudgezfullr2_t" "Vert infl full weight height for LML obs, regime 2, temperature" ""
rconfig real obs_nudgezrampr2_t namelist,fdda 1 50 rh "obs_nudgezrampr2_t" "Vert infl ramp-to-zero height for LML obs, regime 2, temperature" ""
rconfig real obs_nudgezfullr4_t namelist,fdda 1 50 rh "obs_nudgezfullr4_t" "Vert infl full weight height for LML obs, regime 4, temperature" ""
rconfig real obs_nudgezrampr4_t namelist,fdda 1 -5000 rh "obs_nudgezrampr4_t" "Vert infl ramp-to-zero height for LML obs, regime 4, temperature" ""
rconfig real obs_nudgezfullr4_t namelist,fdda 1 -5000 rh "obs_nudgezfullr4_t" "Vert infl full weight height for LML obs, regime 4, temperature" ""
rconfig real obs_nudgezrampr4_t namelist,fdda 1 50 rh "obs_nudgezrampr4_t" "Vert infl ramp-to-zero height for LML obs, regime 4, temperature" ""
rconfig real obs_nudgezfullr1_q namelist,fdda 1 50 rh "obs_nudgezfullr1_q" "Vert infl full weight height for LML obs, regime 1, moisture" ""
rconfig real obs_nudgezrampr1_q namelist,fdda 1 50 rh "obs_nudgezrampr1_q" "Vert infl ramp-to-zero height for LML obs, regime 1, moisture" ""
rconfig real obs_nudgezfullr2_q namelist,fdda 1 50 rh "obs_nudgezfullr2_q" "Vert infl full weight height for LML obs, regime 2, moisture" ""
rconfig real obs_nudgezrampr2_q namelist,fdda 1 50 rh "obs_nudgezrampr2_q" "Vert infl ramp-to-zero height for LML obs, regime 2, moisture" ""
rconfig real obs_nudgezfullr4_q namelist,fdda 1 50 rh "obs_nudgezfullr4_q" "Vert infl full weight height for LML obs, regime 4, moisture" ""
rconfig real obs_nudgezrampr4_q namelist,fdda 1 -5000 rh "obs_nudgezrampr4_q" "Vert infl ramp-to-zero height for LML obs, regime 4, moisture" ""
rconfig real obs_nudgezfullr4_q namelist,fdda 1 -5000 rh "obs_nudgezfullr4_q" "Vert infl full weight height for LML obs, regime 4, moisture" ""
rconfig real obs_nudgezrampr4_q namelist,fdda 1 50 rh "obs_nudgezrampr4_q" "Vert infl ramp-to-zero height for LML obs, regime 4, moisture" ""
rconfig real obs_nudgezfullmin namelist,fdda 1 50 rh "obs_nudgezfullmin" "Minimum depth through which vertical influence fcn remains 1.0" "m"
rconfig real obs_nudgezrampmin namelist,fdda 1 50 rh "obs_nudgezrampmin" "Minimum depth through which vertical influence fcn decreases from 1.0 to 0.0" "m"
rconfig real obs_nudgezmax namelist,fdda 1 3000 rh "obs_nudgezmax" "Maximum depth in which vertical influence function is nonzero" "m"
Expand Down Expand Up @@ -1772,10 +1774,10 @@ package dfi_dfl dfi_opt==1 - state:dfi_u,df
package dfi_ddfi dfi_opt==2 - state:dfi_u,dfi_v,dfi_w,dfi_ph,dfi_phb,dfi_ph0,dfi_php,dfi_t,dfi_p,dfi_ww,dfi_mu,dfi_tke,dfi_pb,dfi_al,dfi_alt,dfi_TSLB,dfi_SMOIS,dfi_SNOW,dfi_SNOWH,dfi_CANWAT,dfi_SMFR3D,dfi_KEEPFR3DFLAG,dfi_TSK,dfi_SOILT1,dfi_TSNAV,dfi_SNOWC,dfi_QVG,dfi_rh,dfi_tten_rad
package dfi_tdfi dfi_opt==3 - state:dfi_u,dfi_v,dfi_w,dfi_ph,dfi_phb,dfi_ph0,dfi_php,dfi_t,dfi_p,dfi_ww,dfi_mu,dfi_tke,dfi_pb,dfi_al,dfi_alt,dfi_TSLB,dfi_SMOIS,dfi_SNOW,dfi_SNOWH,dfi_CANWAT,dfi_SMFR3D,dfi_KEEPFR3DFLAG,dfi_TSK,dfi_SOILT1,dfi_TSNAV,dfi_SNOWC,dfi_QVG,dfi_rh,dfi_tten_rad
ifdef RUC_CLOUD
package realonly use_wps_input==1 - state:u_gc,v_gc,t_gc,rh_gc,ght_gc,p_gc,xlat_gc,xlong_gc,ht_gc,tsk_gc,tavgsfc,tmn_gc,pslv_gc,greenfrac,albedo12m,pd_gc,psfc_gc,intq_gc,pdhs,sh_gc,qv_gc,qr_gc,qc_gc,qs_gc,qi_gc,qg_gc,qni_gc
package realonly use_wps_input==1 - state:u_gc,v_gc,t_gc,rh_gc,ght_gc,p_gc,xlat_gc,xlong_gc,ht_gc,tsk_gc,tavgsfc,tmn_gc,pslv_gc,sct_dom_gc,scb_dom_gc,greenfrac,albedo12m,pd_gc,psfc_gc,intq_gc,pdhs,sh_gc,qv_gc,qr_gc,qc_gc,qs_gc,qi_gc,qg_gc,qni_gc
endif
ifndef RUC_CLOUD
package realonly use_wps_input==1 - state:u_gc,v_gc,t_gc,rh_gc,ght_gc,p_gc,xlat_gc,xlong_gc,ht_gc,tsk_gc,tavgsfc,tmn_gc,pslv_gc,greenfrac,albedo12m,pd_gc,psfc_gc,intq_gc,pdhs,sh_gc,qv_gc
package realonly use_wps_input==1 - state:u_gc,v_gc,t_gc,rh_gc,ght_gc,p_gc,xlat_gc,xlong_gc,ht_gc,tsk_gc,tavgsfc,tmn_gc,pslv_gc,sct_dom_gc,scb_dom_gc,greenfrac,albedo12m,pd_gc,psfc_gc,intq_gc,pdhs,sh_gc,qv_gc
endif

# only need to specify these once; not for every io_form* variable
Expand Down Expand Up @@ -1977,6 +1979,7 @@ typedef fdob_type real latprt [ - 1 -
typedef fdob_type real lonprt [ - 1 - - "lonprt" "obs longitude for diagnostic printout"
typedef fdob_type real mlatprt [ - 1 - - "mlatprt" "model latitude at obs location"
typedef fdob_type real mlonprt [ - 1 - - "mlonprt" "model longitude at obs location"
typedef fdob_type integer stnidprt {obsstid}{[} - 1 - - "stnidprt" "obs station id for diagnostic printout"
typedef fdob_type real base_state k - 1 - - "base_state" "base-state height on half (mass) levels" "meters"

state fdob_type fdob - -
Expand Down

0 comments on commit 36e0c26

Please sign in to comment.