留言板
- 侠骨柔情
- 前辈，你好，想问下，wrfchem输出的PM2.5和PM10均为零，其他变量比如no2、O3都有数值，这是怎么回事？我用的WRF3.9,namelist.input 设置如下:
&time_control
run_days = 3,
run_hours = 13,
run_minutes = 0,
run_seconds = 0,
start_year = 2013,
start_month = 06,
start_day = 25,
start_hour = 12,
start_minute = 00,
start_second = 00,
end_year = 2013,
end_month = 06,
end_day = 29,
end_hour = 00,
end_minute = 00,
end_second = 00,
interval_seconds = 21600,
input_from_file = .true.,
history_interval = 60,
frames_per_outfile = 1,
restart = .false.,
restart_interval = 43200,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
auxinput5_interval_m = 60,
io_form_auxinput5 = 2
auxinput6_interval_m = 43200,
auxinput6_inname ='wrfbiochemi_d01',
io_form_auxinput6 = 2,
debug_level = 0,
auxinput11_interval_s = 120, 120, 120, 120,
auxinput11_end_h = 192, 30, 30, 30,
/
&domains
time_step = 60,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
e_we = 301,
e_sn = 301,
e_vert = 35,
p_top_requested = 5000,
num_metgrid_levels = 27,
num_metgrid_soil_levels = 4,
dx = 6000,
dy = 6000,
grid_id = 1,
parent_id = 0,
i_parent_start = 1,
j_parent_start = 1,
parent_grid_ratio = 1,
parent_time_step_ratio = 1,
feedback = 1,
smooth_option = 0
/
&physics
mp_physics = 10,
progn = 1,
ra_lw_physics = 1,
ra_sw_physics = 1,
radt = 10,
sf_sfclay_physics = 1,
sf_surface_physics = 1,
bl_pbl_physics = 1,
bldt = 1,
cu_physics = 1,
cudt = 0,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
num_land_cat = 21,
sf_urban_physics = 0,
num_urban_layers = 1200,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
cu_rad_feedback = .false.
/
&fdda
obs_nudge_opt = 0,0,0,0,0
max_obs = 150000,
fdda_start = 0.,0., 0., 0., 0.
fdda_end = 14400.
obs_nudge_wind = 1,1,1,0,0
obs_coef_wind = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
obs_nudge_temp = 1,0,0,0,0
obs_coef_temp = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
obs_nudge_mois = 1,1,1,1,1
obs_coef_mois = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
obs_rinxy = 1080.,360.,90.,30.,180
obs_rinsig = 0.1,
obs_twindo = 0.6666667,0.6666667,0.6666667,0.6666667,0.6666667,
obs_npfi = 10,
obs_ionf = 2,2, 2, 2, 2,
obs_idynin = 0,
obs_dtramp = 40.,
obs_prt_freq = 10, 10, 10, 10, 10,
obs_prt_max = 10
obs_ipf_errob = .true.
obs_ipf_nudob = .true.
obs_ipf_in4dob = .true.
obs_ipf_init = .true.
/
&dynamics
w_damping = 1,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 2,
base_temp = 290.,
damp_opt = 0,
zdamp = 5000.,
dampcoef = 0.2,
khdif = 0,
kvdif = 0,
non_hydrostatic = .true.,
moist_adv_opt = 1,
scalar_adv_opt = 1,
chem_adv_opt = 1,
tracer_adv_opt = 1,
tke_adv_opt = 1,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true.,
nested = .false.,
/
&grib2
/
&chem
kemit = 1,
chem_opt = 5,
bioemdt = 30,
photdt = 30.,
chemdt = 10.0,
io_style_emissions = 2,
emiss_inpt_opt = 1,
emiss_opt = 13,
chem_in_opt = 1,
phot_opt = 0,
gas_drydep_opt = 1,
aer_drydep_opt = 1,
bio_emiss_opt = 0,
gas_bc_opt = 1,
gas_ic_opt = 1,
aer_bc_opt = 1,
aer_ic_opt = 1,
gaschem_onoff = 1,
aerchem_onoff = 1,
wetscav_onoff = 0,
cldchem_onoff = 0,
vertmix_onoff = 1,
chem_conv_tr = 0,
conv_tr_wetscav= 1,
conv_tr_aqchem = 1,
biomass_burn_opt = 0,
aer_ra_feedback = 0,
ne_area = 300,
opt_pars_out = 1,
have_bcs_chem = .true.,
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
- 小肥虾619
- 言深深前辈您好，我想请教下，想用Fortran一行一行判断，如果遇到条件成立的，则把这一样的数据都写到一个文件里，program tiqushuju
! parameter(date=12,val=6)
! dimension a(i),b(i,j)
!character*10 a(date)
real b(3,6)
!character*8 stid(sta)
open(11,file='D:\vf\visit-1.txt')
do i=1,3
read(11,*) (b(i,j),j=1,6)
write(*,*) (b(i,j),j=1,6)
enddo
!pause
open(3,file='D:\vf\visit-NEW-1.txt')
do i=1,3
do j=1,6
if(b(i,j)<=500) THEN
write(3,*) (b(i,k),k=1,6)
endif
enddo
enddo
end
我这样写就会出现，一行中多个符合条件的，这行就会被重复写。
129.0 106.0 124.0 2339.0 2339.0 2339.0
14041.0 200 11167.0 9608.0 14413.0 10776.0
6232.0 6771.0 2214.0 2025.0 7216.0 8735.0
输出结果
129.0000 106.0000 124.0000 2339.000 2339.000 2339.000
129.0000 106.0000 124.0000 2339.000 2339.000 2339.000
129.0000 106.0000 124.0000 2339.000 2339.000 2339.000
14041.00 200.0000 11167.00 9608.000 14413.00 10776.00
- 云一林
- 言深深前辈您好，我参考了论坛的一个eof程序，但是fortran没报错的情况下运行不出来，于是我在Fortran板块发了个eof的帖子把问题表述了一遍，可以麻烦前辈有空的话能帮我看看是哪里出问题了吗？谢谢前辈
- 江淮梅雨天
- 闫师兄现在还在兰州吗？
- 江淮梅雨天
- 闫师兄现在还在兰州吗？
- ruyufe123
- 您好，求助一下，在用2.0 的时候，365 day calendars are no longer supported by sdfopen，这得怎么解决呀？着急用啊
- rosedream22
- 您好，言深深前辈，想问一下就是Fortran批量读取文件时，如何跳过空文件啊？谢谢了！
- 宋斌
- 你好，言深深！感谢你们对气象家园的付出，给我们提供一个非常好的交流平台！我有一个建议，你们看是否合适。是否可以开一个关于各种学术会议的通知的版块？因为各个单位都会办很多的学术会议，但除了一些比较大型的例如气象学年会会议，很多大家都不知道。如果有一个平台让大家把各自高校或单位的学术会议（或者自己知道的国内外会议）通知发出来，相信会更利于学习和交流！你们看这样是否合适？谢谢！
- 郑智勇
- 为啥删我帖子
- green_tea789
- 谢谢
言深深: 我没有操作权限，不过问题已经帮你反馈了
- green_tea789
- 免积分提问
言深深: 哪一个版块？我的版块吗？
- green_tea789
- 抱歉，管理员设置了本版块最后回复于 30 天以前的主题自动关闭，不再接受新回复
言深深: 什么问题不能回复？论坛似乎没有这种设定啊，我好几年以前的帖子还是可以回复的。
- green_tea789
- 论坛已经有几个帖子问这个问题了，版主把过去的帖子设为不能回复，所以重新开了一个。我觉得这样不好，因为即使短时间没人回复，以后也还是有人会想回复的。
- yanan_1018
- 深深，我新写了篇sci，专家审稿意见：
There is a shift in the oscillation period in many of the parameters after 1970 as seen in the wavelet analysis. The discussion should include that this may be associated with climate change occurred in the decade 1970-1980 in many of the global parameters (i.e. sst).
我理解的是对比1970时期的变化，一些影响全球气候变化的关联性参量，比如sst数据的相关性应加入到讨论中。
由于我不是气象专业的，不知道该怎么分析，能否给予指点？
- 林夕
- 常来玩啊~
言深深: 来了···
- zjgi756
- 言大神，你对于小波分析的理解深感佩服，我能问问小波分析之后，怎么对周期分解和类似的逆变呢?求推荐链接，拜谢了！
- dragonloon
- 版主您好！真不好意思要过年了还打扰您！想请教个问题。REOF各特征向量的方差贡献该如何计算？望您不吝赐教！非常感谢！
- arthurmaple
- 你好，深深姐，我是一个在读博士，最近开始接触WRF，请问可以向您讨教一些基础的问题吗？谢谢
- Jasonluo
- 您好，深深，在吗？能加下您QQ号吗？
- Mc.Fish
- 您好。看了您对
http://bbs.06climate.com/forum.php?mod=viewthread&tid=22537
的回复，这个MK检验比较乱，在50年代末、60年代、80年代和20世纪UF和UB都有交点，且交点都在检验线以内，为什么您认为八十年代末存在突变呢？
初学者，
非常感谢。