苹果tokenpocket官网app下载安卓|cctm
CCTV节目官网_央视网
>CCTV节目官网_央视网
热榜
探索·发现发现黄帝城探索·发现《秘境追踪第五辑》探索·发现(亚洲版)发现之旅《探索发现》考古晋国——发现晋都《告台湾同胞书》发表40周年纪念会发现之路发现国宝·发现
历史记录清除记录
热门搜索
1今日说法
2新闻联播
3百家讲坛
4天下足球
5NBA
6今日关注
7宗师列传
8天网
9一路向前
10焦点访谈
正在热播
精彩推荐
渔舟上的协奏曲
梅家荡——浙北水乡渔歌悠扬
穿越千年 沉浸式体验美食之旅
哈里-凯恩:冠军的渴望
从“Sora”的横空出世看人工智能
看百年码头如何重焕生机
《每个人都了不起》
《一路向前》
《瓣瓣同心——京津冀协同发展十年》
《千秋诗颂》
《守护明天》
《荣誉解码》
《2024年总台元宵晚会》
《2024年总台春节联欢晚会》
《2024总台春节戏曲晚会》
《山海经奇》
《吾家吾国》
《2024新春喜剧之夜》
《翠云廊》
《跃龙门 春晚幕后纪实》
广告
广告
5关闭
直 播
电视剧
电视剧 >
高山清渠
穆桂英挂帅
武林外传
节目单
动画片
动画片 >
变形金刚
葫芦兄弟
小马宝莉
频道大全
纪 录 片
纪录片 >
寻古中国
山海经奇
奇妙中国
栏目大全
特别节目
特别节目 >
诗画中国
从长安到罗马
博物馆乐游记
主 持 人
4K专区
听音
扫一扫 沉浸听
热榜
焦点访谈
高端访谈
方圆剧阵
健康中国
新闻调查
军事科技
角儿来了
篮球公园
视界
更多节目+
ai妙招儿
健康中国
千方百剂
味道中国
ai美食
微视频
手机看
今日热门
海外博主看中国
从《年锦》看中式美学 探索中国传统纹样魅力
“丝路海运”航线开行突破1.4万艘次
黄金价格上涨 消费者购金热情不减
春运盘点 客流总量创新高 春运出行变化大
赴港澳“个人游” 市民游客出行意愿高涨
金火伴月天象将于3月8日黎明上演
江苏:春菜抢“鲜”上市 品种丰富价格实惠
安徽黄山:梅花绽放枝头俏 风景如画收益好
北极地区出现“无冰日” 或比此前预测提前
巴黎奥运会开幕式门票方案敲定
日本东京都政府大楼上演灯光秀
英国:超百座乐高积木雕塑齐亮相
国家二级保护动物马鹿现身湿地公园
直播导视
节目单
点击刷新
综合
新闻
财经
综艺
亚洲
体育
电影
国防军事
电视剧
纪录
科教
戏曲
社会与法
少儿
音乐
奥林匹克
农业农村
欧洲
美洲
体育赛事
本周 CCTV-1
人口
本周 CCTV-2
消费主张
本周 CCTV-3
我爱满堂彩
本周 CCTV-4
今日亚洲
本周 CCTV-5
冠军欧洲
本周 CCTV-10
地理·中国
本周 CCTV-12
法律讲堂
本周 CCTV-13
新闻1+1
揭示生命奥秘
传递生命科学
向公众普及健康知识
节目反映卫生健康与人口发展领域民生话题与社会热点,为“健康中国”建设传递正能量。
人口
有价值的消费资讯
关注消费生活动态变化
中国消费市场的守护者
关注消费变化,引领消费时尚潮流,为广大消费者提供有价值的资讯、观点和判断依据。
消费主张
绝技绝活 惊叹全场
曲艺创新 融合之美
用文化的力量去启迪受众
展现最新、最好的语言节目,开发精彩、新颖、观众喜闻乐见的曲艺表演。
我爱满堂彩
离你最近的热点新闻
给你最快的新闻现场
亚太热点事件深度剖析解读
独家资源、权威学者、专业制作,倾力打造亚太时事述评新闻高地。
今日亚洲
欧冠赛事精彩集锦
足坛豪门传奇故事
一起享受最高端的足球赛事
全面的展示欧冠联赛每个比赛日的精彩内容,令球迷全方面享受世界上最成功的足球俱乐部赛事。
冠军欧洲
游历祖国大好山河
体验人与自然和谐共生
新奇的地质地貌
珍惜自然,带您感受大自然的神奇魅力。
地理·中国
传播实用法律知识
服务实际生活需求
解答生活中的法律疑问
传播实用法律知识,告诉你遇到棘手事情该怎么办,解决的途径有哪些?
法律讲堂
还原新闻的全貌
最新最热的新闻话题
解读事件的真相
从时事政策、公共话题、突发事件中选取新闻话题,力求以精度、纯度和锐度为新闻导向,呈现最质朴的新闻。
新闻1+1
微视频
搜片库
找栏目
职场健康课
今日说法
听音
天网
片库
电视剧
动画片
纪录片
特别节目
《青春之城》
《青春之城》
一群父辈“拓荒牛”,在风云激荡的时代,以深沉的家国情怀再次创业,让原创研发的生物医药走向了更广阔的舞台。
都市
都市
《此心安处是吾乡》
《此心安处是吾乡》
平安中国建设者们,直面纷繁复杂的社会矛盾,与各利益集团展开较量、守护人民群众的平安幸福的故事。
现代
现代
《大宅门》
《大宅门》
讲述了医药世家白府随着国家、民族的历史发展而几度兴衰的过程,展现一幅波澜壮阔的社会生活画卷。
家庭
家庭
《觉醒年代》
《觉醒年代》
该剧展现从新文化运动、五四运动到中国共产党建立这段历史画卷,讲述觉醒年代社会风情和百态人生。
革命
革命
《一路向前》
《一路向前》
以修建新中国成立后的第一条铁路——成渝铁路为背景,展示了筑路民工自力更生、艰苦奋斗的崇高精神。
年代
年代
《问苍茫》
《问苍茫》
讲述了青年毛泽东如何从一个追随者成长为先行者、开拓者,并探索出中国革命正确道路的热血故事。
革命
革命
《正阳门下》
《正阳门下》
讲述韩春明和他的朋友乘着改革开放的东风,用辛勤劳动和不懈坚持,积累财富,收获爱情的故事。
家庭
家庭
《今生有你》
《今生有你》
钟汉良、李小冉再度携手,讲述了一对久别重逢的恋人,相互治愈,重拾爱情的故事。
都市
都市
《鹿精灵》
《鹿精灵》
一部融合奇幻、喜剧、冒险等元素的动画片,讲述主人公小鹿迪尔的冒险之旅。
搞笑
搞笑
《熊出没·变形记》
《熊出没·变形记》
光头强和两熊被“缩小机”意外缩小,踏上了追赶强爸和机器的“微观世界”冒险之旅。
冒险
冒险
《百变校巴11》
《百变校巴11》
同学们到户外写生,画下可爱的动物,神奇的自然现象等,遇到困难时,百变校巴变身超级形态解决难题。
冒险
冒险
《大山里的梦之队》
《大山里的梦之队》
以贵州乡村篮球赛事“村BA”为载体,从青少年成长、体育运动、美丽乡村、文化自信四个纬度展开。
励志
励志
《少年英雄小哪吒》
《少年英雄小哪吒》
讲述了哪吒神奇降世,而后从文习武、惩恶扬善、带领人们兴茶马、走商贸、富民强兵,心怀天下的故事。
冒险
冒险
《猫头鹰布布》
《猫头鹰布布》
围绕布布、比尔、邦尼的冒险故事。在森林中他们会认识很多的小动物,与他们的朋友一起探索大自然的奥秘。
益智
益智
《敦煌的故事》
《敦煌的故事》
通过生动多样的动画形式,再现千年敦煌的盛世风华,带领观众感受敦煌不同历史时期的独特魅力。
历史
历史
《班班和莉莉的小王国》
《班班和莉莉的小王国》
讲述了仙女公主莉莉和她最好的朋友小精灵班班在魔法小王国的冒险故事。
冒险
冒险
《石窟中国》
《石窟中国》
通过讲述几代石窟艺术的保护者、研究者和传播者的故事,体现中国石窟是中外文化交流的重要载体。
人文历史
人文历史
《超级装备》第二季
《超级装备》第二季
聚焦当今中国的顶尖装备,通过创新、优化的叙事方式,展现出令人惊叹的工业奇观。
军事
军事
《大敦煌》
《大敦煌》
以敦煌为史实例证,讲述东西方交流历史,传达中国一以贯之的“和而不同”、“美美与共”的文明价值观。
人文历史
人文历史
《越山河》
《越山河》
通过以小见大的手法,让观众深切地感受和理解一座座桥梁将高山深谷“变平”,反复突破世界纪录背后的精神伟力。
人文历史
人文历史
《红色密档》
《红色密档》
回望建党百年的光辉历程,以一系列极具代表意义的档案为切入点,聚焦于大的时代背景下档案背后的真实故事。
人文历史
人文历史
《长城之歌》
《长城之歌》
带您穿越2000多年历史长河,以全球视野、全新维度领略长城代表的人类文明奇迹,探寻长城隐含的中华文明密码。
人文历史
人文历史
《如果国宝会说话》第四季
《如果国宝会说话》第四季
通过讲述一件文物,介绍国宝背后的中国精神、中国审美和中国价值观,以全新视角带您读懂中华文化。
人文历史
人文历史
《大国小鲜 石榴说》
《大国小鲜 石榴说》
以第一人称视角,讲述石榴前世今生的精彩故事;解读中国人源于食物灵感,而汇聚的独特审美。
人文历史
人文历史
《诗画中国》
《诗画中国》
以祖国的江河文脉为线索,以诗画为载体,绘就“现代文明”的崭新画卷。
文化
文化
《大国工匠 2023年度人物》
《大国工匠 2023年度人物》
10位“大国工匠年度人物”来自装备制造、水利、环境、人工智能、电力、文物保护、冶炼、气象等行业。
新闻
新闻
《美美与共》
《美美与共》
品位多元文化之韵,欣赏交流互鉴之美,讲述共建“一带一路”的动人故事。
文化
文化
《宗师列传·唐宋八大家》
《宗师列传·唐宋八大家》
生动还原八位宗师的生活场景,再现了宗师磅礴跌宕的人生故事。
文化
文化
《非遗里的中国》第二季
《非遗里的中国》第二季
在全国各省份打造非遗街区,以互动体验方式呈现当地多项非遗项目。
文化
文化
《零容忍》
《零容忍》
选取典型案例,采访纪检监察干部、调查对象及涉案人员,讲述全面从严治党、推进反腐败斗争的故事。
其他
其他
《军工巡礼》
《军工巡礼》
走进科研一线、生产一线,零距离感受国防科技工业践行“强国、强军、富民”目标过程中取得的辉煌成就。
军事
军事
《2024总台春节联欢晚会》
《2024总台春节联欢晚会》
“思想+艺术+技术”融合传播,用创新精品节目和互动方式传递真情实感,为全球华人献上一台精彩的除夕文化大餐。
综艺
综艺
广告
看点
新闻30分今日说法第一时间军迷行天下九州大戏台冠军欧洲中国文艺田园帮帮团
朱迅海上工作初体验 遇见海上“最美奋斗者”
每个人都了不起
美高超音速导弹现身关岛 渲染“威胁”谋求“优势”
防务新观察
入睡难总失眠?您这些地方火旺了,得降!
健康之路
乳腺结节=定时炸弹?别慌,听听医生怎么说
人口
折耳根奶茶 老干妈披萨!极致混搭你敢尝尝吗?
越战越勇
一手拿枪一手拿锹 另一个“战场”上的英雄往事
星推荐
川渝人的麻辣清单 麻辣不同但快乐相似
探索·发现
探访“云峰八寨” 看600年前“调北征南”军屯变迁
军迷行天下
文学借光影映刻中国精神、中国价值、中国力量
今日影评
“牛”县养成记 阳信的牛肉到底“牛”在哪?
谁知盘中餐
高燃又硬核!南昌舰指挥员祁伟光带你长“舰”识
军事纪实
王泽山院士为什么被称为中国“火炸药之王”?
国防故事
氛围感拉满!澎湃国潮成为消费新风尚
消费主张
穿越千年领略唐风宋雅 来一场活色生香的美食之旅
走遍中国
隐秘于乡村的铁匠如何打造出名扬天下的章丘铁锅?
超级农人秀
寻觅散落在滕州的鲁班遗迹 重新认识这位发明家
跟着书本去旅行
记住这几点 让自制盒饭像刚出锅一样
回家吃饭
春季荨麻疹高发 中医妙招巧应对
生活圈
收视TOP榜
人形机器人开始进厂“打工”了
脑机接口何时可以 “人机合一”
睡眠不好可能是这里出了问题
国货崛起 年轻人带“飞”新国潮
川渝人三餐四季 一口麻辣滚烫
春季常起的“风团”究竟是个啥?
神准!有一种传球叫德布劳内式
自制盒饭有妙招 万能公式要记好
信义文化对法律制度有何影响?
“叮当”声里见证传奇铁锅传承路
全部榜单
广告
王忠心
导弹兵王养成记
广告
刘清伟
伶仃洋里的“美猴王”
广告
苏炳添
超越自我 全力奔跑
闫卓
冰场上的追梦人
周嘉鹰
中国女冰“五连扑”有多难
今年金价是涨还是跌?投资黄金还是好时机吗?
消费主张
“领养”变“强卖”?揭开宠物救助领养的骗局
法律讲堂
春雷动 万物生!惊蛰时节吃这些美味更应景
回家吃饭
有了这本书,在家就能看熊猫界的顶流大明星啦!
读书
他把空心贡面做成“网红”,有什么特别之处?
超级农人秀
来岭南独有的水上花市 感受千年水上花墟
记住乡愁 第十季
6G时代加速来临 “鸿蒙”迎新“盘古”发力
今日亚洲
揭秘土耳其无人机航母 雄心勃勃欲打造“全球海军”
兵器面面观
19岁欧冠捧杯 皇马中场核心卡马文加的励志故事
欧冠开场哨
演习数量翻倍 美韩大规模军演“点燃”半岛
防务新观察
燃爆味蕾!军营大厨的拿手菜来袭
军事纪实
万千毛中拣一毫!湖笔为什么能成为“笔中之冠”?
探索·发现
风雅工具还是艺术载体 唐人为何钟情于屏风?
百家讲坛
探寻豆类的百变吃法 解锁百变“豆”趣
谁知盘中餐
关注涉侨文物和文化遗产保护 百年码头重焕生机
今日说法
跑步锻炼看似简单,但你真的跑对了吗?
健康之路
北纬20度的海上 究竟藏着什么惊喜?
生财有道
让科技不再是压力 这款一体机老年人也能轻松掌控
时尚科技秀
听音
大宋梦华录
宋人的生活有着哪些别样风貌
长城之歌
何以中国 何以文明 长城作答
千古中医故事
走近一代名医扁鹊、李时珍、张仲景
康震说李清照
她为什么拥有“一代词宗”的美誉?
诸葛亮之道
他的“道”对后世产生什么样的影响
广告
两会汇聚全球目光 世界对中国充满期待
环球视线
冯远征:话剧也需要国家级最高奖项
中国电影报道
开新篇!中国新能源高质量发展跑出“加速度”
今日亚洲
为和平出发!他们坚守使命诠释中国担当
军武零距离
古老文物遇上鲜活创意 博物馆如何花式“整活”
走遍中国
我的观众我来宠!江西艺术中心宠粉程度令人惊讶
艺览天下
探秘AI的奇妙世界 它究竟能带给我们什么?
创新进行时
硬核伙食 部队野战“移动厨房”超出你想象
军事纪实
因爱相聚!他们是一群热爱音乐的普通人
开门大吉
走进朱村 来看那些适合春天的非遗级时尚单品
记住乡愁
百工祖师!七行八作的匠人为何都要拜鲁班?
跟着书本去旅行
德布劳内的传球:让对手望尘莫及的技艺!
天下足球
噱头还是真实用?助睡眠保健食品真能改善睡眠?
生活圈
别具匠心 唐代宫廷中都有哪些各具华彩的家具?
百家讲坛
一看二摸三拍照 教你轻松排查脊柱侧弯风险
健康之路
“万千毛中拣一毫”的湖笔是如何制成的?
超级农人秀
2024消费新观察:今年你最想去哪儿旅游?
消费主张
《她们:魏晋南北朝女子图鉴》为8位女性立传
读书
经典热播剧场
《远方的山楂树》书写知青的平凡与伟大
《霞光》反套路轻喜剧演绎谍战剧新类型
《和平之舟》中国海军舰行万里守卫和平
《幸福二重奏》“高婚商”夫妻的幸福之道
《汉武大帝》再现一代雄主刘彻传奇一生
经典热播剧场
广告
广告
广告
广告
广告
广告
广告
广告
广告
广告
广告
CCTV-1综 合
CCTV-2财 经
CCTV-3综 艺
CCTV-4中文国际
CCTV-5体 育
CCTV-5+体育赛事
CCTV-6电 影
CCTV-7国防军事
CCTV-8电视剧
CCTV-9纪 录
CCTV-10科 教
CCTV-11戏 曲
CCTV-12社会与法
CCTV-13新 闻
CCTV-14少 儿
CCTV-15音 乐
CCTV-16奥林匹克
CCTV-17农业农村
新闻联播
焦点访谈
等着我
晚间新闻
经典咏流传
挑战不可能
故事里的中国
今日说法
时代楷模发布厅
开讲啦
我有传家宝
第一动画乐园
动物世界
人与自然
正大综艺
人口
国际艺苑
生活提示
中华民族
秘境之眼
典籍里的中国
中国诗词大会
生活圈
对话
经济半小时
经济信息联播
第一时间
消费主张
欢乐大猜想
创业英雄汇
一槌定音
生财有道
职场健康课
央视财经评论
回家吃饭
中国经济大讲堂
正点财经
天下财经
生活家
魅力中国城
金牌喜剧班
我的艺术清单
走在回家的路上
你好生活
舞蹈世界
文化十分
越战越勇
黄金100秒
星光大道
开门大吉
非常6+1
我要上春晚
向幸福出发
艺览天下
幸福账单
天天把歌唱
回声嘹亮
喜上加喜
中国文艺报道
走遍中国
环球综艺秀
美食中国
深度国际
中国新闻
今日关注
今日亚洲
中国文艺
海峡两岸
远方的家
中国舆论场
国家记忆
记住乡愁
国宝·发现
今日环球
华人故事
中国缘
鲁健访谈
中华医药
欧冠开场哨
天下足球
体育晨报
体坛快讯
体育新闻
体育世界
体谈
足球之夜
篮球公园
棋牌乐
冰球冰球
ATP周刊
冰天雪地
今日影评
中国电影报道
谁是终极英雄
军事科技
军营的味道
军迷行天下
五星剧场
国防军事早报
防务新观察
军事报道
军歌嘹亮
砺剑
讲武堂
军事纪实
老兵你好
正午国防军事
军事制高点
军事纪录
第二战场
兵器面面观
世界战史
国防故事
军武零距离
剧说很好看
星推荐
寰宇视野
发现
人文地理
时代
特别呈现
真相
自然
万象
9视频
故事·中国
活力·源
微9
百家讲坛
地理·中国
自然传奇
味道
健康之路
探索·发现
透视新科技
实验现场
科幻地带
科学动物园
时尚科技秀
人物·故事
考古公开课
跟着书本去旅行
百家说故事
读书
解码科技史
创新进行时
中国影像方志
戏曲青年说
一鸣惊人
九州大戏台
CCTV空中剧院
中国京剧音配像精粹
梨园闯关我挂帅
梨园周刊
角儿来了
过把瘾
名段欣赏
青春戏苑
名家书场
宝贝亮相吧
典藏
中国戏曲像音像集萃
了不起的戏曲
警察特训营
小区大事
法律讲堂生活版
一线
法律讲堂文史版
夜线
热线12
方圆剧阵
夕阳红
生命线
见证
天网
从心开始
现场
道德观察
心理访谈
热心话
法治深壹度
全网追踪
律师来了
决不掉队
焦点访谈
面对面
国际时讯
环球视线
每周质量报告
法治在线
新闻直播间
共同关注
午夜新闻
新闻30分
朝闻天下
24小时
新闻调查
新闻1+1
世界周刊
东方时空
新闻周刊
军情时间到
高端访谈
周末动画片
快乐体验
小小智慧树
智力快车
音乐快递
英雄出少年
快乐大巴
动漫世界
新闻袋袋裤
大风车
七巧板
大手牵小手
智慧树
动感特区
动画大放映
看我72变
快乐童行
风车剧场
乐享汇
全球中文音乐榜上榜
一起音乐吧
影视留声机
CCTV音乐厅
风华国乐
民歌·中国
精彩音乐汇
中国音乐电视
乐游天下
音乐人生
音乐公开课
聆听时刻
中国节拍
合唱先锋
音乐周刊
童声唱
北京2022
艺术里的奥林匹克
逐冰追雪
五环纪事
奥林匹克传奇
奥秘无穷
奥林匹克人
我的奥林匹克
奥林匹克新鲜说
乡土中国
共富经
三农长短说
乐游新乡村
超级农人秀
田园帮帮团
三农群英汇
乡村大舞台
乡理乡亲
谁知盘中餐
大地讲堂
农业气象
直播
节目单
频道大全
栏目大全
综合
财经
综艺
中文国际
体育
电影
国防军事
电视剧
纪录
科教
戏曲
社会与法
新闻
少儿
音乐
体育赛事
农业农村
新闻联播
焦点访谈
等着我
晚间新闻
经典咏流传
挑战不可能
故事里的中国
今日说法
时代楷模发布厅
开讲啦
我有传家宝
第一动画乐园
动物世界
人与自然
正大综艺
人口
国际艺苑
生活提示
中华民族
秘境之眼
典籍里的中国
中国诗词大会
生活圈
对话
经济半小时
经济信息联播
第一时间
消费主张
欢乐大猜想
创业英雄汇
一槌定音
生财有道
职场健康课
央视财经评论
回家吃饭
中国经济大讲堂
正点财经
天下财经
生活家
魅力中国城
金牌喜剧班
我的艺术清单
走在回家的路上
你好生活
舞蹈世界
文化十分
越战越勇
黄金100秒
星光大道
开门大吉
非常6+1
我要上春晚
向幸福出发
艺览天下
幸福账单
天天把歌唱
回声嘹亮
喜上加喜
中国文艺报道
走遍中国
环球综艺秀
美食中国
深度国际
中国新闻
今日关注
今日亚洲
中国文艺
海峡两岸
远方的家
中国舆论场
国家记忆
记住乡愁
国宝·发现
今日环球
华人故事
中国缘
鲁健访谈
中华医药
欧冠开场哨
天下足球
体育晨报
体坛快讯
体育新闻
体育世界
体谈
足球之夜
篮球公园
棋牌乐
今日影评
中国电影报道
谁是终极英雄
军事科技
军营的味道
军迷行天下
五星剧场
国防军事早报
防务新观察
军事报道
军歌嘹亮
砺剑
讲武堂
军事纪实
老兵你好
正午国防军事
军事制高点
军事纪录
第二战场
兵器面面观
世界战史
国防故事
军武零距离
剧说很好看
星推荐
寰宇视野
发现
人文地理
时代
特别呈现
真相
自然
万象
9视频
故事·中国
活力·源
微9
百家讲坛
地理·中国
自然传奇
味道
健康之路
探索·发现
透视新科技
实验现场
科幻地带
科学动物园
时尚科技秀
人物·故事
考古公开课
跟着书本去旅行
百家说故事
读书
解码科技史
创新进行时
中国影像方志
戏曲青年说
一鸣惊人
九州大戏台
CCTV空中剧院
中国京剧音配像精粹
梨园闯关我挂帅
梨园周刊
角儿来了
过把瘾
名段欣赏
青春戏苑
名家书场
宝贝亮相吧
典藏
中国戏曲像音像集萃
了不起的戏曲
警察特训营
小区大事
法律讲堂生活版
一线
法律讲堂文史版
夜线
热线12
方圆剧阵
夕阳红
生命线
见证
天网
从心开始
现场
道德观察
心理访谈
热心话
法治深壹度
全网追踪
律师来了
决不掉队
焦点访谈
面对面
国际时讯
环球视线
每周质量报告
法治在线
新闻直播间
共同关注
午夜新闻
新闻30分
朝闻天下
24小时
新闻调查
新闻1+1
世界周刊
东方时空
新闻周刊
军情时间到
高端访谈
周末动画片
快乐体验
小小智慧树
智力快车
音乐快递
英雄出少年
快乐大巴
动漫世界
新闻袋袋裤
大风车
七巧板
大手牵小手
智慧树
动感特区
动画大放映
看我72变
快乐童行
风车剧场
乐享汇
全球中文音乐榜上榜
一起音乐吧
影视留声机
CCTV音乐厅
风华国乐
民歌·中国
精彩音乐汇
中国音乐电视
乐游天下
音乐人生
音乐公开课
聆听时刻
中国节拍
合唱先锋
音乐周刊
童声唱
冰球冰球
ATP周刊
冰天雪地
乡土中国
共富经
三农长短说
乐游新乡村
超级农人秀
田园帮帮团
三农群英汇
乡村大舞台
乡理乡亲
谁知盘中餐
大地讲堂
农业气象
片库
热榜
看点
微视频
视频百科
微故事
ai美食
主持人
全部
全部
直播
节目单
热榜
主持人
听音
频道大全
栏目大全
新闻联播
焦点访谈
等着我
今日说法
海峡两岸
今日关注
今日亚洲
更多
片库
电视剧
动画片
纪录片
特别节目
更多
看点
微视频
视频百科
微故事
ai美食
热榜
探索·发现发现黄帝城探索·发现《秘境追踪第五辑》探索·发现(亚洲版)发现之旅《探索发现》考古晋国——发现晋都《告台湾同胞书》发表40周年纪念会发现之路发现国宝·发现
历史记录清除记录
热门搜索
1今日说法
2新闻联播
3百家讲坛
4天下足球
5NBA
6今日关注
7宗师列传
8天网
9一路向前
10焦点访谈
广告
今日热门
直播导视
精品
片库
看点
央视大全
手机访问
扫一扫手机继续看
返回顶部
广告
空气质量模型_百度百科
模型_百度百科 网页新闻贴吧知道网盘图片视频地图文库资讯采购百科百度首页登录注册进入词条全站搜索帮助首页秒懂百科特色百科知识专题加入百科百科团队权威合作下载百科APP个人中心收藏查看我的收藏0有用+10空气质量模型播报讨论上传视频数学工具本词条由“科普中国”科学百科词条编写与应用工作项目 审核 。空气质量模型是基于人类对大气物理和化学过程科学认识的基础上,运用气象学原理及数学方法,从水平和垂直方向在大尺度范围内对空气质量进行仿真模拟,再现污染物在大气中输送、反应、清除等过程的数学工具,是分析大气污染时空演变规律、内在机理、成因来源、建立“污染减排”与“质量改善”间定量关系及推进我国环境规划和管理向定量化、精细化过渡的重要技术方法。 [1]中文名空气质量模型外文名air quality model分 类城市模型、区域模型和全球模型基本理论梯度输送理论、统计理论等典型模型法规化中小尺度模型等作 用分析大气污染时空演变规律等目录1简介2基本理论3发展历程4分类5质量模型6存在的问题7相关建议简介播报编辑空气质量模型是用数学方法来模拟影响大气污染物的扩散和反应的物理和化学过程。基于输入的气象数据和污染源信息如排放率,烟囱高度等,这些模型可以模拟直接排入大气的一次污染物和由于复杂的化学反应形成的二次污染物。这些模型对空气质量管理是非常重要的,因为他们被许多机构用来测算源分担率,同时帮助制定有效的削减污染物排放的政策。例如空气质量模型可以用来预测一个新的污染源会不会达标排放,如果超标的话,还可以给出适当的控制措施。此外,空气质量模型还可以预测未来新的政策法规实施后的污染物的浓度。可以估计政策法规的有效性以及减少人类和环境暴露。 [2]基本理论播报编辑空气质量模型一般考虑了以下大气过程:排放(人为和自然源排放)、输送(水平平流和垂直对流)、扩散(水平和垂直扩散)、化学转化(气、液、固相化学反应)、清除机制(干湿沉降)等。其理论研究一直是沿着湍流扩散三个理论体系发展起来的,即梯度输送理论(K理论),统计理论和相似理论。(1)梯度输送理论(K理论)是在湍流半经验理论的基础上发展起来的。其缺陷体现在:一方面,它把无规则的湍涡看成分子热运动,假定湍涡是流体微团,与分子输送模型具有相同属性,由此得到的梯度与通量之间的线性关系,实质上这只是一种假定。另一方面,近地层流场情况十分复杂,湍流输送的性质远非简单的线性关系,尤其是湍流交换系数,它随大气湍流场的性质及空间尺度而改变,其形式难以确定。因此梯度输送理论在小尺度预测上缺陷很突出,但它在处理大尺度污染扩散问题上具有一定优越性,能够利用观测的风速廓线资料,不需假定某种分布形式,即可得到污染物的浓度分布。(2)统计理论是从湍流场的统计特征量出发,描述流场中扩散物质的散布规律。泰勒把扩散系数和湍流脉动场的统计特征量联系起来,用气象参数来表达这些统计特征量,找出扩散参数和气象条件的联系,导出了适用于连续运动扩散过程的泰勒公式。该理论的核心是扩散粒子关于时间和空间的概率分布,通过概率分布函数描述扩散粒子浓度的空间分布和时间变化。泰勒公式是在均匀、定常的假设条件下导出,而实际大气并不符合这种条件,只有在下垫面开阔平坦、气流稳定的小尺度扩散处理中,才近似满足这样的条件。(3)相似理论是在量纲分析基础上发展起来的,是研究近地层大气湍流的一种有效理论方法。其基本原理是关于拉格朗日相似性的假设,假定流场的拉格朗日性质仅决定于表征流场欧拉性质的已知参数,粒子扩散的特征与流场的拉格朗日性质相联系。在上述假定下,可以把大气扩散和风速及温度的空间分布联系起来,但由于量纲分析的复杂性和不确定性,目前主要在小尺度的铅直扩散问题中比较成功。 [1]发展历程播报编辑鉴于空气质量模型在大气污染控制中的重要地位,开发和推广新型的空气质量模型显得尤为重要。自1970年到现在,USEPA或其他机构共资助开发了三代空气质量模型:70年代到80年代,EPA推出了第一代空气质量模型,这些模型又分为箱式模型、高斯扩散模型和拉格朗日轨迹模型,其中高斯扩散模型主要有ISC、AERMOD、ADMS等,拉格朗日模型如OZIP/EKMA、CALPUFF等;80年代到90年代的第二代空气质量模型主要包括UAM、ROM、RADM在内的欧拉网格模型;90年代以后出现的第三代空气质量模型是以CMAQ、CAMx、WRF-CHEM、NAQPMS为代表的综合空气质量模型,即“一个大气”的模拟系统。(1)第一代空气质量模型主要包括了基于质量守恒定律的箱式模型、基于湍流扩散统计理论的高斯模型和拉格朗日轨迹模式。当时的模型一般以Pasquill和Gifford等研究者得出的离散不同稳定度条件下的大气扩散参数曲线和Pasquill方法确定的扩散参数为基础,采用简单的、参数化的线性机制描述复杂的大气物理过程,适用于模拟惰性污染物的长期平均浓度。高斯模式(如ISC、AERMOD、ADMS)由于其结构简单,对输入数据的要求不高以及计算简便,20世纪60年代以后,在大气环境问题中得到了最为广泛的应用。但近年来城市及区域环境问题如细粒子、光化学烟雾等往往与污染物在大气中的化学反应紧密相关,而第一代模型没有或仅有简单的化学反应模块,这使它们的应用受到了很大限制。但是这些模型结构简单、运算速度快、长期浓度模拟的准确度高,至今仍在常规污染物模拟方面被广泛使用。值得注意的是第一代空气质量模型的划分并不是非常明确,例如ADMS、AERMOD、CALPUFF模型应用了90年代以来大气研究的最新成果,与传统的第一代模型已有很大不同。(2)20世纪70年代末80年代初,随着对大气边界层湍流特征的研究,研究者开展了大量室内试验、数值试验和现场野外观测等工作,发现高斯模型对许多问题都无法解答,这逐渐推动了第二代空气质量模型的发展。第二代欧拉数值空气质量模型中加入了比较复杂的气象模式和非线性反应机制,并将被模拟的区域分成许多三维网格单元。模型将模拟每个单元格大气层中的化学变化过程、云雾过程,以及位于该网格周边的其他单元格内的大气状况,这包括污染源对网格区域内的影响以及所产生的干、湿沉降作用等。这类模型在1980-1990年期间被广泛应用。这一时期一些三维城市尺度光化学污染模式(如CIT、UAM等模式)、区域尺度光化学模式ROM以及酸沉降模式(RADM、ADOM、STEM等模式)开始得到研究。我国第二代空气质量模型主要有中国科学院雷孝恩基RADM模型建立的高分辨率对流层化学模式HRCM,中国科学院大气物理所等研发的区域空气质量模式RAQM和三维时变欧拉型区域酸沉降模式RegADM等。(3)第二代空气质量模式在设计上仅考虑了单一的大气污染问题,对于各污染物间的相互转化和相互影响考虑不全面,而实际大气中各种污染物之间存在着复杂的物理、化学反应过程。因此,20世纪90年代末美国环保局基于“一个大气”理念,设计研发了第三代空气质量模式系统Medels-3/CMAQ,CMAQ是一个多模块集成、多尺度网格嵌套的三维欧拉模型,突破了传统模式针对单一物种或单相物种的模拟,考虑了实际大气中不同物种之间的相互转换和互相影响,开创了模式发展的新理念。当前主流的第三代空气质量模式还包括CAMx、WRF-CHEM等。特别是美国大气研究中心NCAR开发的WRF-CHEM模式考虑了气象和大气污染的双向反馈过程,在一定程度上代表了区域大气模式未来发展的主流方向。中国的第三代空气质量模式以中国科学院大气物理所自主研发的嵌套网格空气质量预报模式NAQPMS为代表,目前已在北京、上海、深圳、郑州等城市空气质量实时预报业务中得以应用。 [1]分类播报编辑无论是一代、二代或三代空气质量模型按照尺度划分,大致可以分为城市模型、区域模型和全球模型(如GEOS-Chem);按机理划分,可分为统计模型和数值模型,前者是以现有的大量数据为基础做统计分析建立的模型,后者则是对污染物在大气中发生的物理化学过程(如传输、扩散、化学反应等)进行数学抽象所建立;从流体力学的角度看,空气质量模型又分为拉格朗日模型和欧拉模型,前者由跟随流体移动的空气微团来描述污染物浓度的变化,后者则相对于固定坐标系研究污染物的运动,以空间内固定的微元为研究对象;从模型研究对象来看,空气质量模型又分为扩散模式、光化学氧化模式、酸沉降模式、气溶胶细粒子模式和综合性空气质量模式。 [1]质量模型播报编辑空气质量数值模式已经有数十年的发展历史,在世界范围内产生了数十个不同的模式,当前国际上典型的空气质量模式主要包括ISC3、AERMOD、ADMS、CALPUFF等法规化中小尺度模型,NAQPMS、CAMX、WRF-CHEM、CMAQ等综合型区域尺度模型和GEOS-CHEM等全球尺度空气质量模型。法规化中小尺度模型ISC3、AREMOD、ADMS、CALPUFF均属于第一代空气质量模型,是最典型的法规化中小尺度模型。按照模型法规化进程划分,ISC3属于第一代法规性模型,而AREMOD、ADMS、CALPUFF为第二代法规性模型。四个模型的优点均在于结构简单、计算速度快、基础数据要求低等。其简单易用的优点奠定了其成为法规化模型的基础。不足之处体现在适用尺度相对较小、没有化学过程或化学过程较为简化,基本理论假设过于理想,不能很好的模拟O3、PM2.5、酸雨等区域性复合型大气污染过程。从实践应用来看,ADMS、AREMOD、CALPUFF模型多用于环境影响评价和城市尺度一次污染物的模拟,尤其是在国内外环境影响评价领域发挥了主力军作用,已被多个国家定为法规化模型。国家环保部发布的《环境影响评价技术导则大气环境》(HJ2.2-2008)奠定了ADMS、AERMOD、CALPUFF三个模型的在我国环境影响评价领域的法规地位。1、ISC3模型ISC3(Industrial Source Complex 3)模式属于第一代法规模式,是美国环保局开发的一个复合工业源空气质量扩散模式,其公式利用稳态封闭型高斯扩散方程。ISC3模式的适用范围一般小于50km,模拟物质一般为一次污染物。模式可处理各种烟气抬升和扩散过程,如静风、风廓线指数、烟囱顶端尾流、城市建筑下洗、污染物转化、沉积和沉降等。可对点源、面源、线源、体源等多种污染源进行模拟;可输出多种污染物浓度以及颗粒物的沉积和干、湿沉降量等计算结果;污染物可选取SO2、TSP、PM10、NOx等;可选择逐时、数小时、日、月及年等多种平均模拟时段。ISC3与AERMOD、ADMS对比,它的最大的优势是其操作简单,ISC3需要的输入数据相对较少,而且可以利用NWS(美国国家气象局)航空数据。当污染物质为惰性物质,气象条件单一时,除污染源排放参数以外,ISC3要求气象数据为风向、风向角、大气稳定度、混合层高度、接受点地形高度、建筑物的维度。ISC3的主要劣势是大气边界层结构知识已经发展进步到新的阶段,而模型对湍流扩散过程的模拟并没有跟上时代发展。ISC3的局限性如下:(1)没有考虑建筑物对周围点源扩散的影响;(2)没有考虑流线反射对烟羽轨迹的影响;(3)没有考虑烟羽抬升过程中尾流速度缺失的影响;(4)没有解决近处尾流截获远处尾流中物质的问题;(5)两个下洗方程的接口不连续;(6)没有考虑低矮建筑物对周围风向影响;(7)小风稳定条件下,污染浓度估计过大。2、AERMOD模型AERMOD由美国国家环保局联合美国气象学会组建的法规模式改善委员会开发。其目标是开发一个能完全替代ISC3的法规模型,新模型将采用ISC3的输入与输出结构、应用最新的扩散理论和计算机技术。20世纪90年代中后期,法规模式改善委员会在ISC3模型框架的基础上成功开发出AERMOD扩散模型,AERMOD系统包括AERMET气象、AREMAP地形、AERMOD扩散三个模块,适用范围一般小于50km。该系统以高斯统计扩散理论为出发点,假设污染物的浓度分布在一定程度上服从高斯分布。可用于乡村环境和城市环境、平坦地形和复杂地形、低矮面源和高架点源等多种排放扩散情形的模拟和预测。AERMOD是一种稳态烟羽模型。在稳定边界层(SBL),将垂直和水平方向的浓度分布看作高斯分布。在对流边界层(CBL),将水平分布也看作是高斯分布,但是垂直分布考虑用概率密度函数来描述。另外,在对流边界层中,AERMOD考虑“烟羽抬举”(plumelofting)现象:从浮力源出来的部分烟羽物质,先是升到边界层顶部附近并在那里停留一段时间,然后混合入对流边界层内部中。AERMOD计算穿透进入稳定层的部分烟羽,允许它在某些情况下重新返回边界层内。无论在稳定边界层还是在对流边界层中,AERMOD均考虑了弯曲烟羽导致的水平扩散加强现象。AERMOD具有以下特点:(1)以行星边界层(PBL)湍流结构及理论为基础。按空气湍流结构和尺度概念,湍流扩散由参数化方程给出,稳定度用连续参数表示;(2)中等浮力通量对流条件采用非正态的PDF模式;(3)考虑了对流条件下浮力烟羽和混合层顶的相互作用;(4)对简单地形和复杂地形进行了一体化的处理;(5)可以计算城市边界层,建筑物下洗,以及干、湿沉降等清除过程。3、ADMS模型ADMS模型是由英国剑桥环境研究中(CERC)开发的一套先进的三维高斯型大气扩散模型,属新一代大气扩散模型,适用范围一般小于50km。ADMS可模拟点源、面源、线源和体源排放出的污染物在短期(小时平均、日平均)和长期(年平均)的浓度分布,还包括一个街道窄谷模型,适用于简单和复杂地形,同时也可考虑建筑物下洗、湿沉降、重力沉降和干沉降以及化学反应等功能。ADMS模型耦合了大气边界层研究的最新进展,利用常规气象要素来定义边界层结构。ADMS模型与其它大气扩散模型的一个显著区别是:使用了最小莫宁-奥布霍夫(Monin-Obukhov)长度和边界结构的最新理论,精确定义边界层特征参数;另外ADMS模型在不稳定条件下摒弃了高斯模式体系,采用高斯概率密度函数(PDF)及小风对流模式。ADMS利用莫宁-奥布霍夫(Monin-Obukhov)长度表示大气稳定程度,定义用Lo表示。在白天,由于地表受热,大气处于不稳定状态,这时Lo是负值;而在夜间,由于地表辐射冷却,大气处于稳定状态,这时Lo是正值。如果Lo绝对值接近于零,表明大气非常不稳定(负值时)或非常稳定(正值时)。在城市区域由于地表障碍物(如建筑物)产生的机械扰动会使得边界层趋向中性。因此,在城市区域的稳定时间段(夜间),估算的莫宁长度值可能比实际情况要偏小,即偏向稳定。为了解决这个问题,在模式中稳定时间段里设置一个最小的Lo值。最小Lo值根据障碍物高度对区域流场影响的大小确定。4、CALPUFF模型CALPUFF是三维非稳态拉格朗日扩散模式系统,与传统的稳态高斯扩散模式相比,它能更好的处理长距离污染物传输(50km以上的距离范围)。它由西格玛研究公司(Sigma Research Corporation)开发,是USEPA长期支持开发的首选法规化模型。CALPUFF模型系统包括三部分:CALMET、CALPUFF、CALPOST,以及一系列对常规气象、地理数据进行预处理的程序。CALMET气象模型用于在三维网格模型区域上生成小时风场和温度场。CALPUFF非稳态三维拉格朗日烟团输送模型利用CALMET生成的风场和温度场文件,输送污染源排放的污染物烟团,模拟扩散和转化过程。CALPOST通过处理CALPUFF输出文件,生成所需浓度文件用于后处理及可视化。CALPUFF具有以下自身的优势和特点:(1)能用于模拟从几十到几百公里中等尺度的环境问题;(2)能模拟一些非稳态的情况(静小风、熏烟、环流、地形和海岸效应);(3)气象模型包括了陆上和水上边界层模型,可以利用MM5或WRF中尺度气象模式输出的网格风场作为观测数据,或者作为初始猜测风场;(4)采用地形动力学、坡面流参数方法对初始猜测风场进行分析,适合于粗糙、复杂地形条件下的模拟;(5)加入了处理针对面源浮力抬升和扩散的功能模块。综合型区域尺度模型事实上大气污染过程异常复杂,各种污染物之间存在极为复杂的物理、化学反应及气固两相转化过程。尽管我国已确定了部分法规化模型,但这些过于简单的空气质量模型很难再现真实的大气污染过程。由于科学研究与环境决策的目的在于追求大气模拟的真实性、内生原理性和污染过程的系统性,因此在科学研究与环境决策领域较少使用ADMS、AREMOD、CALPUFF等已有法规化模型,应用最多的均为第三代综合性空气质量模型,如:NAQPMS、CAMx、WRF-CHEM及CMAQ等。这些模型均具有以下共同优点:(1)充分考虑了各种大气物理过程和各污染物间的化学反应及气固两相转化过程,可模拟多污染物间的协同效应;(2)基于嵌套网格设计,可用于模拟局地、区域等多种尺度的大气环境问题;(3)基于“一个大气”的设计理念,通过一次工作可以同时模拟各种大气环境问题,特别适用于模拟O3、PM2.5、酸雨等区域性复合型大气污染过程。存在的不足是:(1)对气象、污染源等基础数据要求过于苛刻。尤其是对排放清单中的污染物排放量要求具体到每一个化学物种、每一个网格及每小时。由于排放清单的复杂性,排放清单编制已成为一个新的研究领域;(2)功能灵活多样,但可操作性降低。为增加模型开发及应用的灵活性,第三代模型均无可视化操作界面,采用模块化集成设计方式,使用者必须熟悉模型的架构、基本物理化学原理及模型程序代码;(3)计算机专业知识要求大幅度提高。第三代空气质量模型计算量极大,多运行在基于LINUX操作系统的高性能集群计算机平台,需较高硬件资源及专门人员负责平台的日常管理和维护;(4)海量输入输出数据需要分析及可视化。第三代空气质量模型输入输出数据少则数百G,多则数千G,海量数据的管理、分析及可视化大幅增加了工作成本。1、NAQPMS模型嵌套网格空气质量预报系统(NAQPMS)由中国科学院大气物理研究所自主开发研制。该模式系统经历了近20年的发展,通过集成自主开发的一系列城市、区域尺度空气质量模式发展而成。NAQPMS为三维欧拉输送模式,垂直坐标采用地形追随坐标,垂直方向不等距分为18层;水平结构为多重嵌套网格,采用单向和双向嵌套技术,水平分辨率一般为3-81公里。NAQPMS由四个子系统组成,分别为基础数据系统、中尺度天气预报系统、空气污染预报系统和预报结果分析系统。NAQPMS可用于多尺度污染问题的研究,不但可以研究区域尺度的空气污染问题(如臭氧、细颗粒物、酸雨、沙尘等污染物的跨界跨国输送等),还可以研究城市尺度空气污染的发生机理及其变化规律,以及不同尺度之间的相互影响过程。NAQPMS模式目前主要应用于空气质量预报领域,已在北京、上海、深圳、郑州等城市的空气质量预报业务中得到大量应用。2、WRF-CHEM模型WRF-Chem模式是美国最新发展的区域大气动力-化学耦合模式,是在NCAR开发的中尺度数值预报气象模式(WRF)中加入大气化学模块集成而成。中尺度数值预报模式(WRF)是一个完全可压非静力模式,对湍流交换、大气辐射、积云降水、云微物理及陆面等多种物理过程均有不同的参数化方案,可以为化学模式在线提供大气流场,模拟污染物输送(包括平流、扩散和对流过程)、干湿沉降、气相化学、气溶胶形成、辐射和光分解率、生物所产生的放射、气溶胶参数化和光解频率等过程。WRF-CHEM的最大优点是气象模式与化学传输模式在时间和空间分辨率上完全耦合,实现真正的在线反馈。该模式尚在我国尚处于探索研究阶段,应用案例相对较少。3、CAMx模型CAMx模式是美国ENVIRON公司在UAM-V模式基础上开发的综合空气质量模式,它将“科学级”的空气质量模型所需要的所有技术特征合成为单一系统,可用来对气态和颗粒物态的大气污染物在城市和区域的多种尺度上进行综合性评估。CAMx除具有第三代空气质量模型的典型特征之外,CAMx最著名的特点包括:双向嵌套及弹性嵌套、网格烟羽(PiG)模块、臭氧源分配技术(OSAT)、颗粒物源分配技术(PSAT)等。CAMx可以在三种笛卡儿地图投影体系中进行模拟:通用的横截墨卡托圆柱投影(Universal Transverse Mercator)、旋转的极地立体投影(Rotated Polar Stereographic)和兰伯特圆锥正形投影(Lambert Conic Conformal)。CAMx也提供在弯曲的线性测量经纬度网格体系中运算的选项。此外,垂直分层结构是从外部定义的,所以各层高度可以定义为任意的空间或时间的函数。这种在定义水平和垂直网格结构方面的灵活性,使CAMx能适应任何用来为环境模型提供输入场的气象模型。4、Models-3/CMAQ模型CMAQ模型是我国应用最广泛、最为成熟的第三代空气质量模型,由USEPA于1998年第一次正式发布。CMAQ最初设计的目的在于将复杂的空气污染问题如对流层的臭氧、PM、毒化物、酸沉降及能见度等问题综合处理,为此Models-3/CMAQ模式最大的特色即采用了“One-Atmosphere”的设计理念,能对多种尺度、各种复杂的大气环境污染问题进行系统模拟,CMAQ模型目前已成为美国EPA应用于环境规划、管理及决策的准法规化模型。该模型的特点在于:(1)可以同时模拟多种大气污染物,包括臭氧、PM、酸沉降以及能见度等各种环境污染问题在不同空间尺度范围内的行为;(2)充分利用了最新的计算机硬件和软件技术,如高性能计算、模块化设计、可视化技术等,使空气质量模拟技术更高效、更精确,且应用领域趋于多元化。Models-3/CMAQ系统由排放清单处理模型(SMOKE)、中尺度气象模型(mm5模型或WRF模型等)和通用多尺度空气质量模型(CMAQ)三部分组成,其中CMAQ是整个系统的核心。CMAQ模型主要由边界条件模块BCON、初始条件模块ICON、光解速率模块JPROC、气象-化学预处理模块MCIP和化学输送模块CCTM构成。CMAQ模型的关键部分是化学输送模块CCTM,污染物在大气中的扩散和输送过程、气相化学过程、气溶胶化学过程、液相化学过程、云化学过程以及动力学过程都由该模块模拟完成。其它模块的主要功能主要是为CCTM提供输入数据和相关参数。CCTM模块提供了多种气相化学机制和气溶胶化学机制供使用者选择,输出结果包括各种气态污染物和气溶胶组分在内的污染物逐时浓度,以及逐时的能见度和干湿沉降。CMAQ模式需要MM5或WRF气象模式提供模拟所需的气象资料。最新发布的CMAQ5.0版本已实现气象模式与化学传输模式在线耦合,吸收了WRF-CHEM模型优点。2.3全球性大尺度模型比较著名的全球尺度空气质量模式包括MOZART和GEOS-CHEM等模型,其主要用于模拟全球尺度大气污染物的长距离传输及化学反应过程,此外全球尺度空气质量模型还可以与卫星遥感资料结合反演近地面空气污染物浓度,如AaronvanDonkelaar等基于卫星观测的气溶胶光学厚度和GEOS-CHEM模型成功反演了2001-2006年全球PM2.5浓度,首次揭示了全球PM2.5浓度的空间分布特征。除全球性空气质量模型外,全球性气候模型及气候与空气质量双向耦合模型近几年也有了迅速发展。 [1]存在的问题播报编辑1、不能准确把握模型的适用条件,盲目使用模型开展空气质量模拟工作首先应选取合适的空气质量模型,空气质量模型的选取应结合研究问题的特点、模拟的时空尺度、数据可获性、技术复杂性、软硬件平台要求等多种因素,不是所有研究问题都要选择最先进的空气质量模型。先进的空气质量模型一般对基础数据、硬件平台、专业门槛等要求非常苛刻,工作成本相对普通空气质量模型会增加数倍。因此选择合适的空气质量模型是以最小成本实现研究目标的关键。从空气质量模型应用实践来看,使用者常难以准确把握模型的适用条件,或将简单问题复杂化,或将复杂问题简单化,以至于得出不符合客观事实的决策结论。2、空气质量模型愈来愈复杂,专业门槛不断提高第三代空气质量模型系统一般都有三个基本模块:气象过程模拟模块(mm5、WRF等中尺度气象模型)、污染源排放清单模块(SMOKE等排放清单模型)、空气质量模块(CMAQ等模型),因此开展一次空气质量模拟工作必须依次完成气象场模拟、排放清单编制及空气质量模拟三个步骤。而mm5、WRF等中尺度气象学模型,SMOKE等排放清单模型,CMAQ等空气质量模型均较复杂,各种模型都以源代码形式发布,需在LINUX高性能集群计算机平台编译、运行。空气质量模型模拟技术跨越了环境、物理、化学、数学及计算机等多个学科,因此建设一个好的空气质量模拟平台需要多学科专业人员协作完成。除空气质量模拟系统本身的复杂性之外,空气质量模型对气象资料、排放清单等基础输入数据的内容及格式要求也非常苛刻,模型输出的海量数据同样需花费较大人力进行维护、管理、分析及可视化。第三代空气质量模型的“科学性”与“易用性”之间的矛盾非常突出,已严重影响到其推广及应用,这也成为建立新一代法规化空气质量模型的难点。3、排放清单多样化,导致模拟结果无可比性排放清单的准确性是影响空气质量模型模拟结果的最重要因素。我国污染源排放数据一直数出多门,官方掌握的污染源普查、环境统计、总量减排核查、排污申报及收费等数据仅限于环保系统内部使用,详细数据未向社会公开。为此,国内外科研院所纷纷展开了对中国排放清单的估算及研究工作,排放清单呈现“多样化,各自为战”的研究局面,其中空间尺度最大、涵盖物种最多、应用最广的排放清单为清华大学参与编制的MEIC(2010年)、INTEX-B(2006年)、TRACE-P(2000年)排放清单。在这种格局下,我国逐步形成了“官方数据”与“科研数据”两套相对独立的排放清单体系。除官方数据外,不同研究单位建立排放清单的方法学差异同样较大,相互之间无法对比。排放数据来源的多样化直接导致模型模拟结果缺乏可比性,针对同一个问题采用不同的排放清单数据或采用相同的排放清单数据采取不同的技术处理规则,模型模拟结果差异可能非常大。因此,当前我国空气质量模型研究工作大多就事论事,不同研究工作,不同时空范围的研究成果很难相互印证、对照。4、空气质量观测手段单一、观测指标少,关键性参数难以率定尽管空气质量模型是基于科学理论及大量野外观测的基础之上开发而成,但开展空气质量模拟工作的一个重要环节即利用实际观测数据对空气质量模型模拟结果进行验证,将模型中的一些关键参数本地化,使空气质量模型能够真实的反映所研究区域的大气污染过程。我国城市空气质量常规监测主要以近地面观测为主,观测点位多集中在城市中心,且点位密度小,高空探测、卫星遥测等先进的观测手段应用较少,难以获得高时空分辨率的空气质量观测数据。此外,常规空气质量监测指标仅限于SO2、PM10、NO2三项污染物,对VOCs、O3、PM2.5等涉及大气化学反应过程的关键性指标尚未大范围开展监测。因此,对于复合型大气污染过程进行模拟,现有观测数据较难对模型中的关键性参数进行有效率定。 [1]相关建议播报编辑1、建立适用于环境规划与决策的法规化模型从国外空气质量管理的模式来看,大致经历了“排放标准控制—总量控制—质量控制”三个阶段,而污染物总量控制是我国目前大气污染防治最主要的手段之一。在总量控制管理实践中,建立“总量减排”与“质量改善”之间的关系,评估总量减排的环境效果,制定环境效果最大化的污染减排策略已成为决策者无法回避的科学问题。此外,随着大气污染问题的复杂化,我国环境管理模式必然由“总量控制模式”逐步向“质量控制模式”转变,基于质量改善目标制定污染减排策略将成为一个决策难题。而空气质量模型恰好能用于建立“污染减排”与“质量改善”之间的定量化关系,解决上述环境管理中所面临的科学难题,因此建立适用于环境规划、管理及决策的法规化模型是提高大气环境管理与决策水平的必然要求。虽然目前我国已将ADMS、AERMOD、CALPUFF三个模型作为环境影响评价领域的行业性法规化模型,但由于第一代空气质量模型自身设计的局限,无法用于模拟区域性复合型大气污染过程。因此,以CMAQ为代表的第三代空气质量模型应作为我国建立适用于环境规划、管理及决策法规化模型的首选对象。2、法规化模型建设的核心在于制定空气质量模型使用规范我国不同地区的空气污染问题及拥有的软硬件技术条件差异较大,一刀切的制定法规化模型将难以推广应用。因此,建立法规化模型的目的不是强制要求模型应用的单一化,核心在于制定空气质量模型使用的技术规范。模型使用技术规范指在综合考虑不同环境问题、研究尺度、数据可得性、技术难度等因素基础上,建立一套模型选取、方案设定、参数率定、结果验证及结果分析的技术方法和基本流程,针对不同的环境问题推荐出可选的模型系列,以此规范模型的使用规则。从模型适用性来看,每一个空气质量模型都有其优缺点,适用性各有差异,在实际应用中只有最合适的模型没有最好的模型。比如:模拟一次污染物扩散及简易的轨迹可采用第一代空气质量模型,推荐的模式主要包括:ADMS、AERMOD及CALPUFF;模拟大气中各污染物间复杂的反应行为及气固相间的转化作用,应采用第三代空气质量模式,推荐的模式主要包括:CMAQ、CAMx、NAQPMS、WRF-CHEM等。3、排放清单编制技术的标准化是建立第三代法规化模型的前提空气质量模型的法规化不仅仅是模型工具及其技术方法的法规化,排放清单作为空气质量模型最重要的输入数据,将直接影响空气质量模型的模拟结果。因此,空气质量模型法规化的前提是排放清单编制技术的标准化。我国应借鉴欧美等国家国家排放清单的编制经验,改进污染源排放的监测及统计方法,增加日常监测的污染物指标,从监测技术、统计方法、质量控制等多方面提升我国环境统计水平,环境统计应与国际排放清单编制的技术规范逐步接轨。同时应开展重点行业大气污染物的化学成份谱监测,基于污染物化学成份谱建立涵盖多种化学物种(特别是VOCs、PM的关键组份)、能适用于第三代空气质量模型模拟的国家排放清单及其方法学体系。4、建立国家环境基础信息数据库系统,夯实法规模型的应用基础数据是所有模型的基础。空气质量模型所需基础数据涉及地形高程、土地利用、气象观测资料、空气质量监测数据及污染源排放数据等大量环境基础信息数据,但这些数据掌握在不同管理部门,目前我国尚未建立起国家级的环境基础数据库,基础信息缺乏互通、共享机制。因此,我国应尽快建立国家级的环境基础信息数据库,对各种数据进行集约化管理,加大对数据归类、分析、处理、融合及二次挖掘的投入,保障数据来源的畅通及规范。此外应将数据信息通过互联网对外共享,以开放、交流的模式与国内外科研机构加大合作,接受第三方监督、质疑,通过信息公开和大量实践应用提高数据质量。5、在法规化模型的基础之上推进空气质量模型的工程化建设大气环境管理及决策贯穿于“经济发展—能源消费—污染物排放控制—环境效应—健康效益”等多个环节。所谓工程模型既是将经济预测模型、能源消费模型、排放清单模型、空气质量模型及健康效益模型等工具高度集成,发挥综合性决策功能,国际上知名的工程化模型有IIASA开发的GAINS、RAINS等模型。此外仅从空气质量模型自身来看,空气质量模拟系统至少涉及排放清单、气象模拟及空气质量模拟三个模块,建立工程化空气质量模型同样是将排放清单、气象模拟及空气质量模拟三者进行高度集成,使空气质量管理及决策过程简单化、快速化、标准化。 [1]新手上路成长任务编辑入门编辑规则本人编辑我有疑问内容质疑在线客服官方贴吧意见反馈投诉建议举报不良信息未通过词条申诉投诉侵权信息封禁查询与解封©2024 Baidu 使用百度前必读 | 百科协议 | 隐私政策 | 百度百科合作平台 | 京ICP证030173号 京公网安备110000020000cmaq实践学习-飓风案例运行bcon icon mcip cctm脚本修改运行_cmaq运行-CSDN博客
>cmaq实践学习-飓风案例运行bcon icon mcip cctm脚本修改运行_cmaq运行-CSDN博客
cmaq实践学习-飓风案例运行bcon icon mcip cctm脚本修改运行
全服霸主
已于 2022-10-25 13:32:14 修改
阅读量1.1k
收藏
14
点赞数
文章标签:
学习
于 2022-10-25 13:29:16 首次发布
版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/a3637474883/article/details/127511179
版权
没有排放清单也可以使用CMAQ和wrf模型
(48条消息) WRF模式案例运行初体验--飓风示例全过程记录_oceanography-Rookie的博客-CSDN博客_wrf
把这个飓风案例生成的wrf文件导入mcip处理后,接着处理icon与bcon部分最后使用ccmt模块
cctm修改
时间
时间,这个按照mcip设置的修改,可以看到是20161006日3点到20161008日0点
模拟配置
VRSN是intel编译器的版本 选择并行计算mpi MECH是气溶胶id APPL是我设置的文件名字
输入输出工作文件放置
根据 脚本默认的放到对应位置
这些是cctm用到的文件
bcon 和icon 编译后的文件放在这里,文件格式修改为.nc
这些是wrf处理后的文件放在这里哦
运行设置线程
4*4设置为电脑的16线程 mpi并行计算
修改地理信息log名字
来自于
里面的
2016_12SE1
设置层数
使用默认Z轴35层
关闭海洋污染物排放(因为我这次没有排放数据)
第一个Ocean选择N
输入文件位置
根据默认脚本设定好的文件夹变量以及路径,放在对应位置,上面已经提到了
BCON文件名字修改
关闭线面排放
setenv 是设置变量 都为0
mcip脚本修改
cmaq配置脚本绝对路径
设置Gridname
设置对应文件所在位置
添加met后的nc文件
由于wrf输出时间间隔为3小时,添加很烦人,下次改为输出时间间隔12小时或者24小时
文件少一点
设置地理静态数据
这里也是绝对路径,d01是单层模拟
修改时间
时间间隔要大于等于wrf的,wrf设置是3小时,这里180分也可以
开始时间要晚于wrf的开始时间,我为了取整数,开始时间设置为3点
接下来运行脚本ok
脚本都要copy一份,后面修改错了可以再来
bcon脚本修改
修改profile
因为wrf运行时候没有嵌套所以为profile,这里卡了半天
修改GRIDDESC
手动改为绝对路径
接下来就可以运行脚本了icon修改和它一样的
我都csdn名字修改好了,终于是全服霸主了哈哈,前面没修改名字所以图片@一串数字
脚本还在run中
优惠劵
全服霸主
关注
关注
0
点赞
踩
14
收藏
觉得还不错?
一键收藏
知道了
0
评论
cmaq实践学习-飓风案例运行bcon icon mcip cctm脚本修改运行
cmaq实践学习-飓风案例运行bcon icon mcip cctm脚本修改运行
复制链接
扫一扫
CMAQ依赖库.zip
06-16
CMAQ-5.2,curl,hdf5,ioapi,netcdf-c,netcdf-fortran,openmpi,zlib
【Linux】CMAQ实践部署
qq_51228157的博客
08-19
448
UOS Redhat Centos CMAQ部署实践
参与评论
您还未登录,请先
登录
后发表或查看评论
01zzCMAQ
早睡身体好zzZZ的博客
07-27
2942
CMAQ入门回顾
CMAQ在CMAQ各个模块的个人理解
CMAQv5.1的个人说明
CMAQ4.6使用手册, 英文
06-22
cmaq英文教程,不错, 主要讲了CMAQ组件的功能,产生的文件,官方手册,4.6版本
WRF-CMAQ library.zip
10-10
WRF-CMAQ模型安装所需要的库
CMAQ空气质量模式
weixin_55561616的博客
03-01
386
更多资讯,关注公众号:Ai科研学术社;(联系方式见文章底部)空气质量预报模式系统(WRF-CMAQ)和污染源处理技术是目前大范围灰霾天气预警及综合治理的重要手段,其在全国的推广应用将有利于提高实际的业务预报水平、增强防灾减灾能力、取得显著的社会经济效益。【专家】:程老师(研究员),长期从事空气质量模拟改进科研工作,自主研发了多种污染源反演方法和基于三维变分的气溶胶激光雷达资料同化技术,拥有十几年使用CMAQ模式的实际操作经验,在多所高校讲授过CMAQ模式。
SMOKE-CMAQ丨SMOKE多模式排放清单处理及EDGAR/MEIC清单制作方法、CMAQ空气质量模式运行实践及改进技术
最新发布
WangYan2022的博客
09-26
142
从实时的数据监测公布到空气质量数值预报及预报产品的发布,我国在空气质量监测和预报方面取得了一定进展。随着计算机技术的高速发展、空气污染监测手段的提高和人们对大气物理化学过程认识的深入,开发并利用先进的大气化学模式进行我国空气质量的预测预报、对于减少大气污染灾害、提高人民生活质量都具有积极的意义。空气质量预报模式系统(WRF-CMAQ)和污染源处理技术是目前大范围灰霾天气预警及综合治理的重要手段,其在全国的推广应用将有利于提高实际的业务预报水平、增强防灾减灾能力、取得显著的社会经济效益。
CMAQ-5.2安装
点滴
06-16
1679
CMAQ安装教程
前置安装
依赖库安装
1、zlib
tar -zxvf zlib-1.2.8.tar.gz
cd zlib-1.2.8
./configure --prefix=/$home/netcdf_install
make
make check
make install
2、curl
tar -zxvf curl-7.55.1.tar.gz
cd -zxvf curl-7.55.1
./configure --prefix=/$home/netcdf_install.
利用ioapi工具制作CMAQ-ISAM 区域文件(注意:文件生成之后的地区变量是int格式,需要将其转为float格式才可以正常使用)
hy19990903的博客
04-18
378
以河北地区为例,结果如下图所示,标记的区域为红色,数值为1,其余区域为蓝色,数值为0。CMAQ运行所需的OCEAN文件也可以用此方法制作。最近在制作CMAQ-MASK文件,使用了python制作但是没有办法识别,在论坛里寻找一番得知可以通过ioapi的m3mask程序制作,故写下此文记录。运行完毕之后即可得到名为mask.nc的区域划分文件,可使用ncview进行查看。
那些年CMAQ报过的错
六出的博客
04-01
3253
CMAQ是比较常用的多尺度空气质量建模系统,学习CMAQ也有小一年了,期间遇到了很多挑战和困难:安装困难,跑起来各种报错,误删文件等等。下面分享一下下经验,希望有机会帮到同样与CMAQ较劲的你。
一、安装部分:
CMAQ版本很多,本人安装过5.1,5.2,5.2.1,5.3四个版本,大体安装过程没有较大区别,但是花式报错各不相同,有时候真怀疑安装软件也是个玄学。(由于现在多用新版本,所以下面以5.3和5.2版本报错为主进行分享,安装过程具体参考CMAQ-5.2安装及测试数据运行和CMAQ-5.3安装及测
CMAQ-5.3用户手册.pdf
06-16
CMAQ用户手册
CMAQ-Install:Ubuntu 20.04的CMAQ安装
04-04
Ubuntu 20.05的CMAQ安装
CMAQ5.2、数据集CMAQv5.2-Benchmark-SingleDay-Input-09-12-2017.tar.gz
02-23
CMAQ5.2的运行程序,以及测试数据集
基于MM5-Models-3/CMAQ的中国地区大气污染模拟
07-03
使用气象模拟系统,模拟大气污染情况。介绍MM5、Models-3的用法和示例。
CMAQ64KMtif.zip_IDl_cmaq_model
07-14
采用IDL语言对cmaq模式数据进行处理,转换成tiff图像格式
CMAQ空气质量模型(5.2版本)中文操作说明书.pdf
02-22
伯鑫 CMAQ空气质量模型(5.2版本)中文操作说明书.pdf 空气质量模型是基于对大气物理和化学过程科学认识的基础上,运用气象学原理及数学方法,从水平和垂直方向在一定尺度范围内对空气质量进行仿真模拟,再现污染物...
基于CMAQ模型的邯郸市霾污染来源的模拟研究
06-17
本文应用CMAQ模型进行了河北及周边地区区域尺度的模拟计算,通过情景分析估算了周边各地区对邯郸市PM2.5及其主要成分、消光系数的贡献。结果表明河北省对邯郸市PM2.5的贡献最大,达到67.9%。其次为山西省,其贡献率为...
CMAQ:美国EPA的社区多尺度空气质量模型(CMAQ)的代码库。 有关CMAQ的其他背景,请访问
05-04
CMAQv5.3.2 社区多尺度空气质量模型US EPA CMAQ网站:( )。 CMAQ是美国EPA研究与开发办公室的一个活跃的开源开发项目,包括用于进行空气质量模型模拟的一... 教程-有关WRF-CMAQ,ISAM和修改化学机理的新教程。 -
论文研究-64位Linux并行计算大气模型效率优化研究.pdf
07-22
研究了CMAQ大气模型在64位Linux操作系统上不同CPU核心数目并行计算模拟耗时以及结果的差异情况。研究结果表明,并行计算能大幅缩短CMAQ模拟耗时,以16个CPU核心并行处理为性价比最佳值;此时连续模拟中国区域37天...
MEIAT-CMAQ
09-12
MEIAT-CMAQ是一个用于模拟大气污染的工具,其全称为"Model of Emissions of Gases and Aerosols from Interactive Tutorials - Community Multiscale Air Quality"。根据引用,MEIAT-CMAQ提供了多个版本供下载,包括内测版本。在使用MEIAT-CMAQ之前,需要进行一些初始化和配置工作。根据引用和引用,你需要先下载并解压MEIAT-CMAQ的压缩包,然后进入相关目录,安装必要的库文件,配置Python GIS环境。接下来,你可以通过修改namelist.input文件来调整参数,并通过命令行运行相应的脚本来进行模拟计算。引用中提到了一个粗网格排放进行空间分配的脚本,你可以使用python指令运行该脚本。如果在使用过程中遇到任何问题或有任何交流需求,引用中提到你可以通过邮件或者CSDN与开发者取得联系。
所以,MEIAT-CMAQ是一个用于模拟大气污染的工具,它提供了多个版本供下载,包括内测版本。你需要下载并解压相关文件,进行初始化和配置工作,然后根据需要调整参数并运行相应的脚本进行模拟计算。如果有任何问题或交流需求,你可以通过邮件或CSDN与开发者取得联系。
“相关推荐”对你有帮助么?
非常没帮助
没帮助
一般
有帮助
非常有帮助
提交
全服霸主
CSDN认证博客专家
CSDN认证企业博客
码龄4年
长安大学
3
原创
163万+
周排名
114万+
总排名
1374
访问
等级
38
积分
1
粉丝
0
获赞
1
评论
16
收藏
私信
关注
热门文章
cmaq实践学习-飓风案例运行bcon icon mcip cctm脚本修改运行
1172
国控点污染物根据小时求均值代码
111
国控站点污染物求月均值
89
您愿意向朋友推荐“博客详情页”吗?
强烈不推荐
不推荐
一般般
推荐
强烈推荐
提交
最新文章
国控站点污染物求月均值
国控点污染物根据小时求均值代码
2022年3篇
目录
目录
最新文章
国控站点污染物求月均值
国控点污染物根据小时求均值代码
2022年3篇
目录
评论
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
查看更多评论
添加红包
祝福语
请填写红包祝福语或标题
红包数量
个
红包个数最小为10个
红包总金额
元
红包金额最低5元
余额支付
当前余额3.43元
前往充值 >
需支付:10.00元
取消
确定
下一步
知道了
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝
规则
hope_wisdom 发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额
0
抵扣说明:
1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。 2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。
余额充值
CMAQ - 上海交大超算平台用户手册 Documentation
CMAQ - 上海交大超算平台用户手册 Documentation
Skip to content
上海交大超算平台用户手册 Documentation
CMAQ
Initializing search
SJTU HPC Docs
上海交大超算平台用户手册 Documentation
SJTU HPC Docs
Getting Started
快速上手
平台硬件资源
账号及充值
可视化平台
登录交我算集群
文件系统与数据访问
提交作业
容器
软件
软件
软件模块使用方法
编译器、MPI库、数学库、应用工具
基准测试
AI 计算
原子分子工程计算
原子分子工程计算
ABINIT
Amber
Basilisk
Calculix
Cantera
CASM
CESM
CP2k
cplex
Gerris
gnuplot
GAUSSIAN
GROMACS
LAMMPS
LAMMPS-RBE
Nektar++
Neper-FEPX
NWChem
Octave
OpenRadioss
WhiskyTHC
OpenFOAM
Osiris
OVITO
PSI4
PySPH
Paraview
Quantum ESPRESSO
QuickPIC
SIESTA
STAR-CCM+
SUMO
VASP
VMD
WRF
Elphbolt
gmx_MMPBSA
CMAQ
CMAQ
Table of contents
简介
可用版本
算例获取路径
集群上的CMAQ
思源一号上的CMAQ
运行CMAQ的流程
首先拷贝CMAQ_Project到本地
然后将算例解压到CMAQ_Project下的data目录
接下来根据运行核数修改可执行文件的内部参数
执行脚本设置如下
Pi2.0上的CMAQ
ARM上的CMAQ
运行结果(单位为:秒,越低越好)
思源一号上CMAQ的运行时间
Pi2.0上CMAQ的运行时间
ARM上CMAQ的运行时间
参考资料
thirdorder
Gift-BTE
DAMSK
ROMS
PROJ
GAMESS
DeePMD-kit
CHROMA
FLASH
colmap
pymultinest
LORENE
BoltzTraP
Croco
Openseespy
WannierTools
NAMD
Xerus
生物信息计算
常见问题
文档下载
Table of contents
简介
可用版本
算例获取路径
集群上的CMAQ
思源一号上的CMAQ
运行CMAQ的流程
首先拷贝CMAQ_Project到本地
然后将算例解压到CMAQ_Project下的data目录
接下来根据运行核数修改可执行文件的内部参数
执行脚本设置如下
Pi2.0上的CMAQ
ARM上的CMAQ
运行结果(单位为:秒,越低越好)
思源一号上CMAQ的运行时间
Pi2.0上CMAQ的运行时间
ARM上CMAQ的运行时间
参考资料
CMAQ¶
简介¶
CMAQ(The Community Multiscale Air Quality Modeling System)是美国环境保护局开发的可用于空气质量模型模拟的开源项目,该软件结合了大气科学和空气质量建模方面的现有知识,同时与多处理器计算技术相结合,可有效预测臭氧、颗粒物、有毒物质和酸沉降。
可用版本¶
版本
平台
构建方式
模块名
5.3.2
源码
cmaq/5.3.2-oneapi-2021.4.0 思源一号
5.4
源码
cmaq/5.4-oneapi-2021.4.0 思源一号
5.3.3
源码
cmaq/5.3.3-intel-2021.4.0 Pi2.0
5.4
源码
cmaq/5.4-intel-2021.4.0 Pi2.0
5.4
源码
cmaq/5.4-intel-2021.4.0 ARM
算例获取路径¶
/dssg/share/sample/cmaq/CMAQv5.3.2_Benchmark_2Day_Input.tar.gz
/dssg/share/sample/cmaq/CMAQv5.4_2018_12NE3_Benchmark_2Day_Input.tar.gz
/lustre/share/samples/cmaq/CMAQv5.3.2_Benchmark_2Day_Input.tar.gz
/lustre/share/samples/cmaq/CMAQv5.4_2018_12NE3_Benchmark_2Day_Input.tar.gz
集群上的CMAQ¶
思源一号上的CMAQ
Pi2.0上的CMAQ
ARM上的CMAQ
思源一号上的CMAQ¶
可执行文件所在的目录
CCTM_v532.exe
/dssg/opt/icelake/linux-centos8-icelake/oneapi-2021.4.0/wrf_cmaq/cmaq/CMAQ_Project/CCTM
BCON_v532.exe、ICON_v532.exe、mcip.exe
/dssg/opt/icelake/linux-centos8-icelake/oneapi-2021.4.0/wrf_cmaq/cmaq/CMAQ_Project/PREP
您还可以 module show cmaq获取文件安装目录
module avail cmaq
module show cmaq/5.3.2-oneapi-2021.4.0
运行CMAQ的流程¶
此处使用CMAQ的核心可执行文件CCTM_v532.exe。
首先拷贝CMAQ_Project到本地¶
mkdir ~/cmaq
cd ~/cmaq
cp -r /dssg/opt/icelake/linux-centos8-icelake/oneapi-2021.4.0/wrf_cmaq/cmaq/CMAQ_Project ./
然后将算例解压到CMAQ_Project下的data目录¶
cd ~/cmaq/CMAQ_Project/data
cp -r /dssg/share/sample/cmaq/CMAQv5.3.2_Benchmark_2Day_Input.tar.gz ./
tar xf CMAQv5.3.2_Benchmark_2Day_Input.tar.gz
mv CMAQv5.3.2_Benchmark_2Day_Input/* ./
接下来根据运行核数修改可执行文件的内部参数¶
修改的可执行文件为:~/cmaq/CMAQ_Project/CCTM/scripts/run_cctm_Bench_2016_12SE1.csh
@ NPCOL = 16; @ NPROW = 8
### NPCOL*NPROW=运行数据的总核数。比如思源上使用2个节点共128核运行数据,参数配置如上
set END_DATE = "2016-07-02"
### 运行时间为:7月1日-7月2日
执行脚本设置如下¶
脚本的位置应在: ~/cmaq/CMAQ_Project/CCTM/scripts/
#!/bin/csh
#SBATCH --job-name=cmaq
#SBATCH --partition=64c512g
#SBATCH -N 2
#SBATCH --ntasks-per-node=64
#SBATCH --output=%j.out
#SBATCH --error=%j.err
module load cmaq
setenv INPDIR /dssg/home/acct-hpc/hpchgc/data/cmaq/cmaq_test2/CMAQ_Project/data/2016_12SE1
csh run_cctm_Bench_2016_12SE1.csh
Pi2.0上的CMAQ¶
可执行文件所在的目录
CCTM_v54.exe
/lustre/opt/cascadelake/linux-rhel8-skylake_avx512/intel-2021.4.0/wrf_cmaq/cmaq-5.4/CMAQ-CMAQv5.4_14Oct2022/CCTM/scripts/BLD_CCTM_v54_intel
BCON_v54.exe、ICON_v54.exe、mcip.exe
/lustre/opt/cascadelake/linux-rhel8-skylake_avx512/intel-2021.4.0/wrf_cmaq/cmaq-5.4/CMAQ-CMAQv5.4_14Oct2022/PREP
ARM上的CMAQ¶
可执行文件所在的目录
CCTM_v54.exe
/lustre/opt/linux-openeuler22-aarch64/linux-centos8-aarch64/gcc-10.3.1/cmaq/CMAQ-CMAQv5.4_14Oct2022/CCTM/scripts/BLD_CCTM_v54_gcc
BCON_v54.exe、ICON_v54.exe、mcip.exe
/lustre/opt/linux-openeuler22-aarch64/linux-centos8-aarch64/gcc-10.3.1/cmaq/CMAQ-CMAQv5.4_14Oct2022/PREP
运行结果(单位为:秒,越低越好)¶
思源一号上CMAQ的运行时间¶
cmaq/5.3.2-oneapi-2021.4.0
核数
64
128
256
Exec time
0:06:41
0:05:18
0:04:26
Pi2.0上CMAQ的运行时间¶
cmaq/5.4-intel-2021.4.0
核数
64
128
256
Exec time
0:11:37
0:07:57
0:04:55
ARM上CMAQ的运行时间¶
cmaq/5.4-gcc-10.3.1
核数
32
64
128
Exec time
0:38:44
0:23:06
0:17:00
参考资料¶
CMAQ 官网
Last update:
Mar 03, 2024
Previous
gmx_MMPBSA
Next
thirdorder
© Copyright 2024, 上海交通大学网络信息中心.
Created using
Sphinx
7.1.2.
and
Sphinx-Immaterial
基于教研共同体的CCTM网络集体备课模式探索 - 知乎
基于教研共同体的CCTM网络集体备课模式探索 - 知乎切换模式写文章登录/注册基于教研共同体的CCTM网络集体备课模式探索前进https://www.zhihu.com/video/1461387460927860736一、时代背景 (一)国家层面: 2018年,教育部在教育信息化2.0行动计划及实施全国中小学教师信息技术应用能力提升工程2.0中指出:要“探索基于互联网+的教研组织形式”,“将教育信息化作为教育系统性变革的内生变量,支撑引领教育现代化发展,推动教育理念更新、模式变革、体系重构。”2021年,中央网络安全和信息化委员会印发《教育信息化“十四五”规划》提出:加快建设中国教育专用网络和“互联网+教育”大平台,构建泛在的网络学习空间,支撑各类创新型教学的常态化应用,推动优质教育资源开放共享,缩小区域、城乡、校际之间的差距,实现更加公平更有质量的教育。(二)在市级层面:1. 在我市,莞式慕课改革推动教育均衡发展,网络教研成为常态,2. 自2020年疫情开始,我市目前已构建了初中信息技术网络教研集体备课共同体,网络备课已初成规模。(三)学校层面:教研共同体中,厚街湖景中学、石排中学均为东莞市慕课实验学校,市教育技术实验学校、信息化2.0试点学校,它们在为推动学校优质化发展,两所学校均在不断探索基于互联网+与信息化2.0的教学、教研新模式并取得了良好的成绩。二、问题的提出随着信息技术逐渐深入,网络教研由于很好地突破“时空限制、资源共享”,逐渐得到大家的青睐,网络集体备课也应运而生。以网络为平台,打破传统集体备课在时间、地域等方面的局限,促进教师跨时空的合作与共享,提高区域集体备课质量,对实现教育均衡发展起着重要作用。然而,在网络集体备课中,常存在以下情况:(一)网络备课没有明确的组织和管理、缺乏统一的备课计划,监督、评价机制,后期难以持续开展;(二)网络备课、研修无备而来;(三)由于各学校学情不同而导致学情分析不到位或缺少学情分析;(四)网络备课变成学科带头人、教学能手或骨干教师备课,其他老师参与度不高,没有充分调动全体老师的积极性、参与性;(五)个别老师奉行“拿来主义”,只索取不分享;(六)教学设计、课程开发不规范,网络备课质量不高,没有进行二次、三次备课,也没有及时修订,形成数字资源存储;(七)部分信息技术兼职维护或其他事务,集体备课没有定时开展;如何加强对网络集体备课机制与文化建设,营造网络集体备课的氛围,充分发挥教师的集体智慧,提升区域集体备课的质量,对实现区域教育均衡发展有着举足轻重的作用。三、理论依据 实践共同体作为一个完整的概念,最早由lave和Wenger在《情境认知中的合法边际参与》一书中提出,它表达了一种“以知识为基础的社会结构”。Wenger等在《实践社团学习型组织知识管理指南》中对实践社团实践共同体做了一个明确的定义,即“实践共同体有着共同关注点、同样问题或者对同一个话题的热情,通过在不断发展的基础上互相影响,加深在这一领域的知识和专业技能。”网络教研集体备课共同体:由具有共同愿景,教师自发组织,基于互联网环境,依托信息技术与手段,以提高集体备课质量和专业发展为根本宗旨,实现区域教师共同体进行集体备课及进行相关教研活动并注重成员之间的社会确认、经验资源共享的一种开放、动态、交互、生态的新型教师学习型组织。四、CCTM网络集体备课 CCTM集体备课指的是以教研共同体(Community)为中心,集体备课(Collective lesson preparation)为核心任务,以技术支持(Technology)为支撑、以体制机制(Mechanism)为保障的CCTM网络集体备课模式,如下图1所示:图1其中,集体备课活动主要指以主备课人备课、网络研讨、二次备课、反思修改、课程包的发布与分享、个性化修改使用为流程的主题集体备课活动;实践共同体指参与集体备课教研活动的所有人,包括主体、协作者、实验学校以及教研共同体内的角色分工;技术支撑主要是指支持教研共同体的技术工具,如钉钉平台、Ittools信息技术教学辅助平台、幕布、教师协同办公等;保障机制主要指区域教研共同体中的规则,主要包括政策激励、积分制度、领导机构、制度保障等。在该区域网络教研模式中,为解决集体备课实际问题,学习共同体以技术平台为支撑,以体制机制为保障,通过参与基于网络的主题教研活动,最终达到提高教师专业能力、提升集体备课质量、实现区域教育均衡发展的目的。(一)前期调研,组建团队在组建网络集体备课共同体前,首先利用问卷星对全市初中信息技术教师展开问卷调查,了解教师所在学校、年级、科组或个人集体备课的现状、面临的问题、自身的特长,是否愿意加入网络集体备课团队并享受相关的权利及义务等。通过对问卷数据调查,了解一线教师的需求,并对有意向参加的老师合理进行分组(课程设计组、技术组、专家组、资源收集组等),组成网络集体备课组,确定通过线上、线下结合的方式召开会议,进一步明确网络备课的目标、意义及分工等。为了更好地做好统筹安排,在分工上,我们可以根据教师的专业特长在团队里面建立活动计划组、导师组、课程开发组、课程审核组、技术支撑组,从而形成一支执行力强的团队。活动计划组一般由各片区的片长、学科带头人一起构成;导师组一般由经验丰富的名师、学科带头人或教学能手构成;课程开发组则由广大一线老师参与;课程审核组主要有各备课组的组长及专家导师组构成;技术组主要协助平台开发与完善、服务器搭建、微视频制作及多媒体处理等。通常,集体备课都是在学校规定的时间、指定的地点、有明确的主题而开展的常态的集体教研活动。网络集体备课是否能够有效、有序、持续的开展,活动也必须定期按计划开展。活动计划组在新学期开学前的假期,就必须做好下一个学期的备课计划,落实相关的人员与任务,确定每一章节的备课团队与主备课人。在集体备课中,各组均以“工作项目化、项目责任化、任务具体化”为导向,将各项工作任务分解成具体目标,责任到人,人人有事做。(二)平台的测试和选择在线集体备课依赖网络平台与网络技术,“平台”、“技术”将“社团',、“实践”和“领域”牢牢结合在一起。选择适合学科的、合理的集体备课、教研平台显得尤为重要。通过前期调研和测试,在网络集体备课方面,我们有针对性地选择了Ittools与钉钉平台。Ittools信息技术教学辅助平台在我市中小学已得到广泛应用,利用Ittools平台我们可以快速地添加新课程、导入课程、课程打包、课程分类、课程查询、课程共享、课程管理、教案管理、课程表管理、管理等。这一切,在技术层面上为集体备课提供了有效的支撑。目前,我市基于Ittools的联盟学校已有29间,联盟学校的老师们都非常熟练利用Ittools平台备课和分享课程资源。钉钉:一款免费的沟通和协同的多端平台 ,提供PC版,Web版和手机版,支持手机和电脑间文件互传。另外,钉钉提供了强大的直播、视频会议、项目任务、在线文档修订等功能,方便开展在线教研与集体备课。(三)保障激励机制不以规矩不能成方圆。有效的机制是确保网络集体备课能持续、有效开展的重要保障。这些机制包括制定实践共同体联盟学校协议、备课的流程规范、教学设计的规范(可以根据常态课、公开课、理论课、作品课等不同课程性质制定相关的规范)、课程的审核制度、网络备课教研活动的计划、还应该建立相关的积分制度,如资源开发、设计、上传、下载及承当公开课、讲座的积分奖励政策等。保障制度进一步明确了联盟体每个学校、每个成员的职责、任务和公约、确保每一个参与的个体、共同体学校服从大局,积极配合、主动作为。另外,为了网络集体备课更有效执行,需要各年级集体备课团队首先通过在线研讨,制定集体备课的计划,并确定各年级组的组长,并做成详细的计划表,通过钉钉项目管理功能,把计划落实到人与具体的时间、具体的事项,并通过钉钉项目管理及时进行查看和落实。五、技术支撑下的网络教研共同体集体备课实施过程如图2所示:图2(一)确定主题,合理分工活动筹划组根据学期授课内容,确定备课的章节内容、时间安排,备课组老师通过任务认领方式对备课内容进行认领,并组建相关的备课团队。(二)线上交流,打磨设计主备课人组织课题组研讨课题,收集建议,负责撰写教学设计,各成员在线阅读教学设计,利用钉钉在线修订功能提出修改建议,完善设计。(三)网络研讨,完善设计主备利用钉钉平台进行说课,其他老师在线聆听、研讨,提出修改意见,专家指导协同设计与修改。(四)二次备课,建设资源主备二次备课,其他成员协助提供相关教学资源和微课制作等,形成课程包,审核组对教学设计进行审核,确认后再进行共享;(五)课例展示,精雕课例授课人用钉钉群直播授课,共同体老师利用钉钉等工具进行在线听课、评课,修改教学设计,重新修订课程包,发布与分享。 (六)课例分享,个性教学实验学校教师下载课程包,并根据实际学情进行微调,彰显教学个性。温格( Wenger,E)在《实践共同体:学习、意义和身份》中提出了实践共同体的三要素,即相互介入、共同事业和共享的技艺库。如图3所示:在集体备课共同体活动中,老师们对学科、课堂教学都怀着情怀、有着共同的兴趣和追求,把集体备课、课堂教学作为共同事业去经营。他们都能以一种平等、信任、合作、理解、尊重和相互促进的方式介入到集体备课中,通过钉钉、Ittools教学辅助平台及各种网络技术工具进行沟通、交流、在线备课、研讨等,把自己多年扎根课堂、对教学进行反思、积累、提炼的教学经验、对教育教学的思考、见解与研究成果与共同体的老师进行分享、交流,分享技能及技术,相互影响、取长补短,最终形成领域和实践的共同理解,从而促进专业化成长。在这个合作的社团里,教师们有着共同的愿景,主动参与、持续学习、相互承认、共同探究、无私分享,每位老师都奉献着自己的一份力量,从而为共同体带来新的信息、新的资源与能量,并从中体验、获得认同感、存在感、责任感、成就感,收获职业的幸福感,这参与的过程是生态的、继承性的、更是动态的、活力的。六、展望基于共同体的CTMM网络集体备课模式作为“互联网+”背景下传统集体备课活动的有效补充,通过钉钉、Ittools信息技术辅助平台及慧教育协同办公等技术赋能,突破地域限制,跨越时空界限,有效地把学校、片区、全市的师资进行有机整合,对提升区域集体备课质量、实现资源共享、促进教师的专业化成长及学校之间的教育均衡发展有着积极的意义,以基于Ittols的校际联盟体及“莞式慕课”实验基地集体备课活动的开展为实践案例,很好地验证了该模式的有效性和可操作性。如何吸引更多学校的老师积极参与,做好合理的分组、分工,提高集体备课的质量与效率,建立有效的保障与激励机制,并结合线上、线下定期开展相关的备课、教研活动,并对成果进行推广应用等值得我们继续进一步探索与实践。编辑于 2022-01-03 17:28教师教研工具慕课推荐赞同 11 条评论分享喜欢收藏申请
中国诚通金属集团有限公司
中国诚通金属集团有限公司
CCTM之設定 | Focus on Air Quality
CCTM之設定 | Focus on Air Quality Link Search Menu Expand Document Focus on Air Quality HomeAQ Data AnalysisGlobal AQ Data Analysis MOZART CAM-chem WACCM NCF to IOAPI Converter CAMS Global Forecast ECMWF ReAnalysis ECMWF Near Realtime MCD19A2.006 NASA EARTHDATARegional AQ Data 大氣污染情報網站圖面之下載 LGHAPTaiwan AQ Analysis 依據Voronoi圖進行空間外插 日均值計算 歷線篩選與繪圖程式 鄉鎮區平均值計算wind modelsCWB Obs. Data Inquire Sys. CODiS日報表下載整併 批次下載CODiS數據 CODiS數據轉成little_r格式 計算距離反比加權值 由CWB數據計算軌跡 地面氣象要素之鄉鎮區平均cwb WRF_3Km 中央氣象局WRF_3Km數值預報產品之下載 WRF_3Km空間內插 grb2wrfout_d04轉檔 相同網格系統之grb2轉檔 將CWB數據填入WRF客觀分析場NCEP FNL檔案之自動化下載 高空觀測數據之自動下載 地面觀測數據之自動下載海溫的讀取全球模式GFSEast Asia Regional Reanalysis EARR 數據批次下載 EARR2csv gfs2csv LGHAP2csvFDSWRFWPS namelist.wps 地形、土地使用之定義與處理 dowps.sh dowpsYMD.shOBSGRID 程式修改及編譯 namelist.oa obsYYMM_run.shREAL & WRF namelist.input doreal_4Nests.sh ndown dowrfWRF-chem 配置及編譯 數據之準備 namelist.input WRF-chem的後處理 輸出揚沙排放量WRFOUT 對時間軸整併任意批次WRF之結果 轉換wrfout格式Global/Regional EmissionRegional Emission inventory in ASia REAS數據之下載 地面排放檔之讀取及轉換 電廠排放之處理Ship Traffic Emission Assessment Model 船隻排放之處理_CAMx 船隻排放之處理_CMAQEmis. Database for Global Atm. Res. EDGAR to CMAQ Redistribution of Ship EmissionsTEDS Pythondbf2csv.pyArea Sources prep_areagridLL.py prep_TimVar DEF's used area_YYMM cnty/town processingBiogenic Sources Monthly Breakdown Daily Breakdown and Fill ncMobile Sources Prepare for MobileS Construct the MobileS Taipei Traffic DataPoint Sources Prepare for PointS TimVar for Elevated PTS Elev. PTse for CAMx Writing Elev PTse TimVar for Ground PTS Grnd PTse for CAMx Elev PTse DowngradeMarine Sources Elev. PointS for HarborsPlume ModelsCO Pathways and Compilation ISC/AERMOD程式之執行SO Pathways 臺灣地區既有點源排放資訊 BPIP程式 BPIP遠端執行 Rotate KML iscParserRE & TG Pathways Prepare of AERMAP Exec. and Effects of AERMAP 地形前處理文字解析與執行控制程式 全臺AERMAP之執行ME Pathways MMIF MMIF Remote Excution System 風花圖之繪製_wdrose.pyOU Pathways PLOTFILE to KML Ozone Limiting MethodCALINE CALINE3 I/O caline3遠端計算服務 CALINE4 I/O Emission FactorsRemote Processing AERMODs 遠端執行系統 ISCST/AERMOD 主程式Trajectory ModelsCALMET calwrf calmet.inp calmet batch processingCALPUFF 軌跡模式的法律定位與風波 calpuff.inp calpuff.inp from V7.0 eg 濃度預報系統之實現 calpuff遠端計算CALPOST calpost.inp CALPUFF模擬結果之轉檔 CALPUFF結果轉nc檔HYSPLIT地面二維軌跡分析 ftuv10.py addVI.py daily_traj.cs 公司版本與新增功能WRF三維軌跡分析 三維軌跡分析 三維反軌跡線之計算 choose10.py程式說明 km.py程式說明 acc_prob.py程式說明 局部反軌跡成果分析 正軌跡分析CMAQ Model SystemMet. Chem. Interface Proc. mcip程式的編譯 wrfout檔案之連結 mcip程式執行腳本 mcip結果初始小時值的延伸 整併wrfout成為全月檔案 解讀GFS之MCIP版本Boundary Condition 全球模擬結果的垂直內插 初始濃度檔案序列之產生 水平內插與污染項目對照 EAC4再分析濃度之填入 CMAQ邊界檔案之產生 大型網格系統切割邊界濃度 逐日循序執行bcon.exeInitial ConditionGeography and Land Data 亞洲土地使用檔案 農作物 土壤 gridded land cover/land use CTCI landuseGround Level Emissions 地面排放檔案轉換 揚塵排放之準備Elev. Point Sources CMAQ點源常數檔案之準備 pt_const程式說明 CMAQ點源變數檔案之準備 pt_timvar程式說明CCTM Main Program CCTM之編譯 CCTM之設定 CCTM之連續執行與監控 CCTM之科學設定Post Processing 綜合空品項目之計算 跨日結果之篩選整併 境外O3 8小時值佔比之計算 境外PM2.5佔台灣平均值比例之計算 增量濃度分析程序與檢討 測站模擬值之讀取 datesIn.py程式說明Abundant NoG Runs CWBWRF_15k HUADON_3kISAM Analysis 執行CMAQ-ISAM ISAM結果檔案之讀取(PM25_IONS) d01地理分區檔案之準備Recommend System 快速啟動 輸出入檔案及傳輸 背景及增量排放量 執行檔與程式庫 後製工具 國網環境與Slurm 費用估算 全年執行批次Forecast System 台灣地區CMAQ模擬預報 逐日mcip之執行腳本 逐日排放檔案之切割程式 10 days of leading time 跨批次串連與均勻化 逐3小時啟動CMAQ 初始時段之均勻化 執行預報腳本之分段說明 WACCM全球預報之應用 預報系統圖檔之拮取 CMAQ分工執行腳本 單/雙工作站版本之切換 地面濃度動畫批次製作 將預報系統移轉到國網 臺灣地區CMAQ預報比較CAMx Model Systemmet for CAMx WRF模擬結果的轉接 MM5模擬結果的轉接 CAMx氣象檔案的修補 氣象檔案的調整迴圈Compilation and Running CAMx之編譯 CAM.in之內容與產生ICBC for CAMx CAMx 模式的ICBC CAMx 模式的BC 從空品檔案切割邊界濃度 BNDEXTR BCON轉.bc檔ozone columns and TUVs 臭氧垂直濃度之處理 紫外線數據與光解係數 化學機制之設定emis for CAMx CMAQ/CAMx排放量檔案之轉換 CAMx地面排放量之調整 REASv3.1之面源 CMAQ/CAMx排放量檔案之轉換 D4範圍地面排放檔案之轉換 CAMx2IOAPI檔案之轉換 CMAQ2CAMx之單向轉換ptsrce for CAMx CAMx點源格式說明 pt2emNest.f程式說明 hsinda3G.py程式說明 ptsEnHRBE.py程式說明 v6_v7格式之轉換 REASv3.1之點源Postprocess of CAMx CAMx模擬結果之壓縮 從空品檔案挑選特定物種時間之內容 CAMx模擬結果之比對(aok) 濃度差異轉成排放量 按日切割avrg檔案 avrg檔案乘上常數 就時間軸合併2個avrg檔案 nc檔版本之shk By_DIST.pyForecast Systems5 days version 台灣地區CMAQ模擬預報 逐日mcip之執行腳本 逐日排放檔案之切割程式 執行預報腳本之分段說明10 days version 預報時距之延長 CMAQ分工執行腳本 單/雙工作站版本之切換 電力、網路及儲存空間 初始/邊界/氣象檔案之檢核IC related problems 跨批次串連與均勻化 逐3小時啟動CMAQ 初始時段之均勻化BC related problems bcon檔案之視覺化 WACCM全球預報之應用Emission ProcessingWorks on NCHC 東亞地區初始及邊界條件 地面及高空排放量 國網上執行wrf專案 國網上執行MCIP5.1 cctm相關腳本之修改 後處理方式的差異FCST Post Processing 整合與切割 預報系統圖檔之拮取 地面濃度動畫批次製作 臺灣地區CMAQ預報比較 空品預報網站 NCDR網站示範習作 D03模擬結果並列播放 D03預報結果比較討論 以github備份模擬圖檔 time-bar html etcPaper ReviewsModels and Comparisons CAMx vs CMAQ CMAQ SO4 issues CALPUFF應用相關論文Disease and Air Pollution 空氣污染與疾病模型框架 呼吸道疾病(肺癌、氣喘、COPD) 心臟血管疾病 腦血管病變與腦中風 阿茲海默症、巴金森氏症等腦部、神經病變 糖尿病 臭氧與疾病發生率的關係 疾病發生率與PM2.5年均值之相關性分析 期初報告 HRA for PM & O3Large Seaside PtSrc EIA 高斯模式及穩定度 煙流下洗現象 煙流受地形效應之影響 林園地區臭氧事件分析 工業園區總量管制上限模式 地形、環流、大氣擴散與空氣品質管理 污染源在空品測站上的表現 AERMOD在臺灣應用之實務探討Terr. EffectsWho's WhoUtilitiesCGI-pythons AERMOD.py CALINE.py CALPUFF.py MMIF計算服務 OLM bpipprim.py drawings CGI-pythons isc_parser rd_kmlFull rotate_kml save_chrpleth save_trjMarble save_wrose surf_trajLL2 煙流模式地形前處理 煙流模式地形前處理CGI主程式 terrain.py說明 空品測站正/反軌跡kml檔案產生系統 網頁計算服務軟體系統的建置Crawlers Download Hiway Forecasting Data Download Taipei Traffic Data PDF檔案之下載與整理 Regulation Reader 環保專案報告之下載DB_servers GPU supported OLAP 線上分析程序——OLAPDates and Times EAC的時間標籤 Julian Day to Calendar Day(j2c) WRF的時間標籤 dt2jul and jul2dt grib2的時間標籤 m3nc的時間標籤 月曆日轉儒略日(c2j)GIS Relatives GDAL GeoTiff檔之讀取 Noise Removal of a Raster Data OpenMaps arcgis收費方案 g0v/twgeojson m3nc檔轉成GeoTiff檔 python解析GML檔 python解析KML檔 shape_to_raster 地圖數位板 空品測站之Delaunay圖 空品測站之Voronoi圖 等值圖KML檔之撰寫 網格濃度在行政區範圍內之平均 行政區範圍格點化 鄉鎮行政區多邊形轉1公里tiff檔 鏈結資訊之地圖展現(uMap) 點狀資訊KML檔之撰寫Graphics CaaS to Graphs HTML MeteoInfo NCL Programs Surfer IO Programs VERDI d3js earth excalidraw imageMagicks leaflet matplotlib Programs wrf-python 由圖面顏色讀取濃度值Markdown VS Code提供Markdown語言的功能 local Jekyll Webblog math in Jekyll md2ipynb note in JekyllNetCDF Relatives 2個nc檔案間的差值 NCKS空品模式應用 NCf之8小時移動平均 NC檔案多維度批次篩選 add_lastHr add_ncatt compilation of IOAPI json converters ncdump應用範例 ncgen & pncgen nc檔案時空最大值 sum of NCFs unMask a NCF 列印m3.nc的時間標籤 列印wrf檔案的時間標籤 按日拆分m3.nc檔案 矩陣階層numpy.newaxis(None)的用法 程式庫之編譯Operation System API之設定 CPU使用率之個人加總 Linux Administration Mac Fortran/C compiler Mac http Ports Mac網路設定及服務 Mac讀寫ntfs格式外接硬碟 Mac鍵盤快速鍵 TMUX 終端機多工器 VS Code安裝使用 git and github grep history linux entry linux之啟動 sed unix環境中的自動排程cron unix系統小工具 數位筆記系統使用心得 自動添加並啟用IP黑名單Parallel Computation mpich編譯及應用 工作站防火牆管理grib Relatives CWB_WRF grib數據檔的時間內插 ncl_convert2nc的應用 eccodes程式編譯與應用Monday, December 25, 2023Monday, November 27, 2023Saturday, February 17, 2024Saturday, November 25, 2023Thursday, December 21, 2023Thursday, December 28, 2023Thursday, February 22, 2024Thursday, January 4, 2024Thursday, September 14, 2023Tuesday, August 29, 2023 NLPTuesday, January 16, 2024Tuesday, June 20, 2023Tuesday, November 28, 2023Wednesday, December 20, 2023view point of education This site uses Just the Docs, a documentation theme for Jekyll. Focus on Air Quality at GitHub CMAQ Model SystemCCTM Main ProgramCCTM之設定 CCTM之設定 Table of contents 前言/背景架構範例位置CCTM編譯相關設定總體設定時間、空間及反應機制設定,執行程式位置設定主要輸出/入檔案目錄位置設定日期時間及mpi控制水平與垂直網格設定時間步階、容錯範圍與科學設定診斷報表輸入檔案之目錄逐日之設定逐日迴圈之準備及啟動逐日之IC/BC/臭氧/光學/氣象檔案逐日排放相關檔案之設定在線排放之設定程序分析之設定輸出檔案名稱、層數執行有關之環境設定控制程式之環境變數設定記憶體與路徑設定、執行程式批次時間報表Reference 前言/背景 架構CCTM的執行過程採用c-shell的環境變數進行設定,因此必須使用批次檔來控制。主要段落包括:編譯相關設定(視機器而異)總體設定時間、空間及反應機制設定(不建議修改),執行程式位置設定(視機器而異)主要輸出/入檔案目錄位置設定(不建議修改)日期時間及mpi控制(不建議修改)時間步階、容錯範圍與科學設定(不建議修改)診斷報表(不建議修改)輸入檔案之目錄(不建議修改)逐日進行設定逐日迴圈之準備及啟動(不建議修改)逐日之IC/BC/臭氧/光學/氣象檔案(不建議修改)逐日排放相關檔案之設定在線排放之設定(未測試)程序分析之設定(未測試)輸出檔案名稱(不建議修改)執行有關之環境設定控制程式之環境變數設定(不建議修改)執行環境資源與路徑設定(視機器而異) 範例位置https://github.com/sinotec2/cmaq_relatives/blob/master/run_cctmMM_RR_DM.csh CCTM編譯相關設定除非測試不同編譯的差異、或移轉到不同電腦平台上,否則這一段應無修改的必要。 設定內容 包括:編譯程式(範例為gcc)。由於netCDF、IOAPI等都必須同一編譯軟體進行編譯,因此其選項為系統性須配套。CCTM編譯環境(範例為/opt/CMAQ_Project/config_cmaq.csh $compiler $compilerVrsn)連結程式庫位置(範例為/opt/netcdf/netcdf4_gcc/lib)偵錯階段可以打開DIAG_LVL,否則設為0,有較少的輸出(範例為setenv CTM_DIAG_LVL 0)範例 1 #!/bin/csh -f2
3 # ===================== CCTMv5.3 Run Script =========================
4 # Usage: run.cctm >&! cctm_v53.log &
5 #
6 # To report problems or request help with this script/program:
7 # http://www.epa.gov/cmaq (EPA CMAQ Website)
8 # http://www.cmascenter.org (CMAS Website)
9 # ===================================================================
10 # ===================================================================
11 # ===================================================================
12 #> Runtime Environment Options
13 # ===================================================================
14
15 echo 'Start Model Run At ' `date`
16
17 #> Toggle Diagnostic Mode which will print verbose information to
18 #> standard output
19 setenv CTM_DIAG_LVL 0
20 setenv CMAQ_HOME $cwd
21 setenv compiler gcc
22 #> Choose compiler and set up CMAQ environment with correct
23 #> libraries using config.cmaq. Options: intel | gcc | pgi
24 #setenv compiler $argv[1]
25 if ( ! $?compiler ) then
26 setenv compiler intel
27 endif
28 if ( ! $?compilerVrsn ) then
29 setenv compilerVrsn Empty
30 endif
31
32 #> Source the config.cmaq file to set the build environment
33 cd /opt/CMAQ_Project
34 source ./config_cmaq.csh $compiler $compilerVrsn
35 setenv LD_LIBRARY_PATH /opt/netcdf/netcdf4_gcc/lib
36 cd $CMAQ_HOME
總體設定 時間、空間及反應機制設定,執行程式位置設定CCTM的執行是由粗網格執行後再進行細網格、時間上也是分月、分批次(與wrf相同)進行,此處以引數方式輸入。 反應機制與執行程式除非偵錯或敏感性測試階段、或移轉到不同電腦平台上,否則其位置應無修改的必要。 設定內容包括:年份(APYR):由所在目錄名稱讀取引數(順序及意義如檔名所示run_cctmMM_RR_DM.csh)月份(MO):第1引數,MO=01, 02,…12批次(RUN):第2引數,RUN=1, 2, … 12 (1~4為前月範圍並不重複執行)層數(DM):第3引數,DM=d02、d04層數與巢狀網格名稱具有對照關係。如要新增或修改,須注意${GRID_NAME}配套層數與ICTYPE/BCTYPE亦有搭配反應機制:與污染物項目名稱有搭配執行程式位置與名稱:程式修改及編譯:詳見另處筆記(CMAQ compilations)倘若修改,建議配套修改輸出目錄位置,以避免覆蓋過原有結果 - 範例 37 set APYR = `echo $CMAQ_HOME|cut -d'/' -f4|cut -c3-4`
38 set MO = $argv[1]
39 set RUN = $argv[2]
40 set DM = $argv[3]
41 set APYM = $APYR$MO
42 set APPL = ${APYM}_run${argv[2]}
43 if ( $DM == 'd02' ) then
44 setenv GRID_NAME sChina_27k # 16-character maximum
45 # set ICTYPE = profile
46 set ICTYPE = regrid
47 else if( $DM == 'd04' ) then
48 setenv GRID_NAME TWN_3X3 # 16-character maximum
49 set ICTYPE = regrid
50 else
51 echo "Error input d02/d04"
52 exit 1
53 endif
54
55 #> Set General Parameters for Configuring the Simulation
56 set VRSN = v53 #> Code Version
57 set PROC = mpi #> serial or mpi
58 set MECH = cb6r3_ae7_aq #> Mechanism ID
59
60
61 #> Define RUNID as any combination of parameters above or others. By default,
62 #> this information will be collected into this one string, $RUNID, for easy
63 #> referencing in output binaries and log files as well as in other scripts.
64 setenv RUNID ${VRSN}_${compilerString}_${APPL}
65
66 #> Set the build directory (this is where the CMAQ executable
67 #> is located by default).
68
69 set BLD = /opt/CMAQ_Project/CCTM/scripts/BLD_CCTM_${VRSN}_${compilerString}ENKF
70 set EXEC = CCTM_${VRSN}.exe
主要輸出/入檔案目錄位置設定由於輸出/入檔案不小,建議設到足夠大的磁碟機為宜。 傳輸的速度也要注意, 設定內容包括:WRKDIR:工作目錄,家目錄。CTM會暫存在此,監控執行進度。OUTDIR:輸出目錄,會在其下繼續開POST,進行後處理。INPDIR:輸入目錄LOGDIR:CTM檔案之目錄,文字檔,偵錯用。NMLpath:反應機制設定檔目錄。由於檔案不大,可以設在與程式相同路徑。範例 71
72 #> Output Each line of Runscript to Log File
73 if ( $CTM_DIAG_LVL != 0 ) set echo
74
75 #> Set Working, Input, and Output Directories
76 setenv WORKDIR ${CMAQ_HOME} #> Working Directory. Where the runscript is.
77 setenv OUTDIR /nas1/cmaqruns/2016base/data/output_CCTM_${RUNID} #> Output Directory
78 setenv INPDIR ${CMAQ_HOME}/data #Input Directory
79 setenv LOGDIR ${OUTDIR}/LOGS #> Log Directory Location
80 setenv NMLpath ${BLD} #> Location of Namelists. Common places are:
81 #> ${WORKDIR} | ${CCTM_SRC}/MECHS/${MECH} | ${BLD}
82
83 echo ""
84 echo "Working Directory is $WORKDIR"
85 echo "Build Directory is $BLD"
86 echo "Output Directory is $OUTDIR"
87 echo "Log Directory is $LOGDIR"
88 echo "Executable Name is $EXEC"
日期時間及mpi控制除非測試不同編譯的差異、或移轉到不同電腦平台上,否則這一段應無修改的必要。 設定內容 包括:如果$RUN是5,表示此模擬批次包括當月的第一天,NEW_START為真,必須特別準備初始值檔案,其餘批次則連結到前一批次執行結果即可BEGD:和其他前處理一樣,此值為前一月的15日,為各批次的起始日期,由該日期計算各批次的起迄日期。批次的時間長度為4天STIME:起始時間,須涵蓋在mcip氣象檔案、BC範圍,並與ICON檔案完全一樣。NSTEPS:執行CCTM的時間範圍。由於CCTM5.31.exe並不會跳日(詳見CMAQ compilations ),必須一天天執行,此值只能設為240000序列或平行運作之選擇:測試除錯階段可能選擇序列運轉,正常作業時必須為平行運作,選項有serial 及mpi。CPU分配以均勻為原則,以充分運用計算資源。如不能除盡網格數,程式仍然能夠繼續進行(nodes=20、supermicro=100、mac=6)。CPU總數如超過電腦實際核心數:仍可執行、但不穩定、故不建議電腦計算速度除CPU設定之外、硬碟IO亦會影響CLOBBER_DATA:將會刪除並覆蓋既有的成果。範例 92 # =====================================================================
93
94 #> Set Start and End Days for looping
95 if ( $RUN == 5 ) then
96 setenv NEW_START TRUE #> Set to FALSE for model restart
97 else
98 setenv NEW_START FALSE #> Set to FALSE for model restart
99 endif
100 set BEGD = `date -ud "20${APYR}-${MO}-15 - 1months" +%Y-%m-%d`
101 @ A = $RUN - 1; @ DD = $A * 4 ; @ ED = $A * 4 + 3
102 setenv EMLAYS_MX 13
103 set START_DATE = `date -ud "$BEGD +${DD}days" +%Y-%m-%d` #> beginning date (July 1, 2016)
104 set END_DATE = `date -ud "$BEGD +${ED}days" +%Y-%m-%d` #> ending date (July 14, 2016)
105 set STD_YMD = `date -ud "${START_DATE}" +%Y%m%d` #> Convert YYYY-MM-DD to YYYYMMDD
106
107 #> Set Timestepping Parameters
108 set STTIME = 0 #> beginning GMT time (HHMMSS)
109 set NSTEPS = 240000 #> time duration (HHMMSS) for this run
110 set TSTEP = 10000 #> output time step interval (HHMMSS)
111
112 #> Horizontal domain decomposition
113 if ( $PROC == serial ) then
114 setenv NPCOL_NPROW "1 1"; set NPROCS = 1 # single processor setting
115 else
116 @ NPCOL = 8; @ NPROW = 12
117 @ NPROCS = $NPCOL * $NPROW
118 setenv NPCOL_NPROW "$NPCOL $NPROW";
119 endif
120
121 #> Define Execution ID: e.g. [CMAQ-Version-Info]_[User]_[Date]_[Time]
122 setenv EXECUTION_ID "CMAQ_CCTM${VRSN}_`id -u -n`_`date -u +%Y%m%d_%H%M%S_%N`" #> Inform IO/API of the Execution ID
123 echo ""
124 echo "---CMAQ EXECUTION ID: $EXECUTION_ID ---"
125
126 #> Keep or Delete Existing Output Files
127 set CLOBBER_DATA = TRUE
128
129 #> Logfile Options
130 #> Master Log File Name; uncomment to write standard output to a log, otherwise write to screen
131 #setenv LOGFILE $CMAQ_HOME/$RUNID.log
132 if (! -e $LOGDIR ) then
133 mkdir -p $LOGDIR
134 endif
135 setenv PRINT_PROC_TIME Y #> Print timing for all science subprocesses to Logfile
水平與垂直網格設定除非網格座標系統改變(增減、平移),否則這一段應無修改的必要。 設定內容 包括:GRIDDESC必須和mcip一致、配套。排放量、ICBC、土地使用、土壤等都須符合統一的座標系統。NZ(無作用),模式垂直網格數內設與WRF、MCIP一致。NCELLS(僅檢視用)輸出檔案的成分與垂直網格數,如不檢討高空煙流的擴散情況,可以不必輸出高空結果以節省空間。CONC(每步階)ACONC(逐時平均)成分與層數成分須以double quote括住,或為“ALL”(如需計算所有粒狀物,建議設為ALL)起訖層數如不指定,則為所有層數。範例137 setenv STDOUT T #> Override I/O-API trying to write information to both the processor
138 #> logs and STDOUT [ options: T | F ]
139
140 setenv GRIDDESC $INPDIR/mcip/${APPL}/${GRID_NAME}/GRIDDESC #> grid description file
141
142 #> Retrieve the number of columns, rows, and layers in this simulation
143 set NZ = 35
144 set NX = `grep -A 1 ${GRID_NAME} ${GRIDDESC} | tail -1 | sed 's/ */ /g' | cut -d' ' -f6`
145 set NY = `grep -A 1 ${GRID_NAME} ${GRIDDESC} | tail -1 | sed 's/ */ /g' | cut -d' ' -f7`
146 set NCELLS = `echo "${NX} * ${NY} * ${NZ}" | bc -l`
147
148 #> Output Species and Layer Options
149 #> CONC file species; comment or set to "ALL" to write all species to CONC
150 setenv CONC_SPCS "O3 NO ANO3I ANO3J NO2 FORM ISOP NH3 ANH4I ANH4J ASO4I ASO4J"
151 setenv CONC_BLEV_ELEV #" 1 1" #> CONC file layer range; comment to write all layers to CONC
152
153 #> ACONC file species; comment or set to "ALL" to write all species to ACONC
154 #setenv AVG_CONC_SPCS "O3 NO CO NO2 ASO4I ASO4J NH3"
155 setenv AVG_CONC_SPCS "ALL"
156 setenv ACONC_BLEV_ELEV #" 1 1" #> ACONC file layer range; comment to write all layers to ACONC
157 setenv AVG_FILE_ENDTIME N #> override default beginning ACONC timestamp [ default: N ]
158
時間步階、容錯範圍與科學設定為使模式有比較性,這段盡量使用內設值,應無修改的必要。 設定內容詳見CCTM之科學設定,包括:同步時間步階:縮短時間步階可提高正確性。加長步階秒數可有效使用CPU。CTM_OCEAN_CHEM(海洋飛沫在風大時會增加PM10濃度、尚未測試)CTM_WB_DUST(風蝕揚塵、可能會與TEDS重複,未測試)CTM_BIDI_FERT_NH3:如取消程式無法進行範例159 #> Synchronization Time Step and Tolerance Options
160 setenv CTM_MAXSYNC 300 #> max sync time step (sec) [ default: 720 ]
161 setenv CTM_MINSYNC 60 #> min sync time step (sec) [ default: 60 ]
162 setenv SIGMA_SYNC_TOP 0.7 #> top sigma level thru which sync step determined [ default: 0.7 ]
163 #setenv ADV_HDIV_LIM 0.95 #> maximum horiz. div. limit for adv step adjust [ default: 0.9 ]
164 setenv CTM_ADV_CFL 0.95 #> max CFL [ default: 0.75]
165 #setenv RB_ATOL 1.0E-09 #> global ROS3 solver absolute tolerance [ default: 1.0E-07 ]
166
167 #> Science Options
168 setenv CTM_OCEAN_CHEM N #> Flag for ocean halgoen chemistry and sea spray aerosol emissions [ default: Y ]
169 setenv CTM_WB_DUST N #> use inline windblown dust emissions [ default: Y ]
170 setenv CTM_WBDUST_BELD BELD3 #> landuse database for identifying dust source regions
171 #> [ default: UNKNOWN ]; ignore if CTM_WB_DUST = N
172 setenv CTM_LTNG_NO N #> turn on lightning NOx [ default: N ]
173 setenv CTM_WVEL Y #> save derived vertical velocity component to conc
174 #> file [ default: N ]
175 setenv KZMIN Y #> use Min Kz option in edyintb [ default: Y ],
176 #> otherwise revert to Kz0UT
177 setenv CTM_MOSAIC N #> landuse specific deposition velocities [ default: N ]
178 setenv CTM_FST N #> mosaic method to get land-use specific stomatal flux
179 #> [ default: N ]
180 setenv PX_VERSION Y #> WRF PX LSM
181 setenv CLM_VERSION N #> WRF CLM LSM
182 setenv NOAH_VERSION N #> WRF NOAH LSM
183 setenv CTM_ABFLUX Y #> ammonia bi-directional flux for in-line deposition
184 #> velocities [ default: N ]
185 setenv CTM_BIDI_FERT_NH3 T #> subtract fertilizer NH3 from emissions because it will be handled
186 #> by the BiDi calculation [ default: Y ]
187 setenv CTM_HGBIDI N #> mercury bi-directional flux for in-line deposition
188 #> velocities [ default: N ]
189 setenv CTM_SFC_HONO Y #> surface HONO interaction [ default: Y ]
190 setenv CTM_GRAV_SETL Y #> vdiff aerosol gravitational sedimentation [ default: Y ]
191 setenv CTM_BIOGEMIS N #> calculate in-line biogenic emissions [ default: N ]
192
193 #> Vertical Extraction Options
194 setenv VERTEXT N
195 setenv VERTEXT_COORD_PATH ${WORKDIR}/lonlat.csv
診斷報表為使模式有比較性,這段盡量使用內設值,應無修改的必要。 設定內容 包括:IOAPI_OFFSET_64 :support大檔案APMDIAG_BLEV_ELEV: 診斷各層之主要點源排放量範例197 #> I/O Controls
198 setenv IOAPI_LOG_WRITE F #> turn on excess WRITE3 logging [ options: T | F ]
199 setenv FL_ERR_STOP N #> stop on inconsistent input files
200 setenv PROMPTFLAG F #> turn on I/O-API PROMPT*FILE interactive mode [ options: T | F ]
201 setenv IOAPI_OFFSET_64 YES #> support large timestep records (>2GB/timestep record) [ options: YES | NO ]
202 setenv IOAPI_CHECK_HEADERS N #> check file headers [ options: Y | N ]
203 setenv CTM_EMISCHK N #> Abort CMAQ if missing surrogates from emissions Input files
204 setenv EMISDIAG F #> Print Emission Rates at the output time step after they have been
205 #> scaled and modified by the user Rules [options: F | T or 2D | 3D | 2DSUM ]
206 #> Individual streams can be modified using the variables:
207 #> GR_EMIS_DIAG_## | STK_EMIS_DIAG_## | BIOG_EMIS_DIAG
208 #> MG_EMIS_DIAG | LTNG_EMIS_DIAG | DUST_EMIS_DIAG
209 #> SEASPRAY_EMIS_DIAG
210 #> Note that these diagnostics are different than other emissions diagnostic
211 #> output because they occur after scaling.
212
213 #> Diagnostic Output Flags
214 setenv CTM_CKSUM Y #> checksum report [ default: Y ]
215 setenv CLD_DIAG Y #> cloud diagnostic file [ default: N ]
216
217 setenv CTM_PHOTDIAG Y #> photolysis diagnostic file [ default: N ]
218 setenv NLAYS_PHOTDIAG #1" #> Number of layers for PHOTDIAG2 and PHOTDIAG3 from
219 #> Layer 1 to NLAYS_PHOTDIAG [ default: all layers ]
220 #setenv NWAVE_PHOTDIAG "294 303 310 316 333 381 607" #> Wavelengths written for variables
221 #> in PHOTDIAG2 and PHOTDIAG3
222 #> [ default: all wavelengths ]
223
224 setenv CTM_PMDIAG N #> Instantaneous Aerosol Diagnostic File [ default: Y ]
225 setenv CTM_APMDIAG Y #> Hourly-Average Aerosol Diagnostic File [ default: Y ]
226 setenv APMDIAG_BLEV_ELEV #"1 1" #> layer range for average pmdiag = NLAYS
227
228 setenv CTM_SSEMDIAG N #> sea-spray emissions diagnostic file [ default: N ]
229 setenv CTM_DUSTEM_DIAG N #> windblown dust emissions diagnostic file [ default: N ];
230 #> Ignore if CTM_WB_DUST = N
231 setenv CTM_DEPV_FILE N #> deposition velocities diagnostic file [ default: N ]
232 setenv VDIFF_DIAG_FILE N #> vdiff & possibly aero grav. sedimentation diagnostic file [ default: N ]
233 setenv LTNGDIAG N #> lightning diagnostic file [ default: N ]
234 setenv B3GTS_DIAG N #> BEIS mass emissions diagnostic file [ default:
輸入檔案之目錄前述為主要檔案的目錄(上層),也包括輸出檔案,此處定義輸入檔部分,直到檔案名稱之前。 設定內容 包括:如果不需要輸入,可以在等號後保持空白OMIpath :採內設值範例236 # =====================================================================
237 #> Input Directories and Filenames
238 # =====================================================================
239
240 set ICpath = $INPDIR/icon #> initial conditions input directory
241 set BCpath = $INPDIR/bcon #> boundary conditions input directory
242 set EMISpath = $INPDIR/emis #> gridded emissions input directory
243 set EMISpath2 = #> gridded surface residential wood combustion emissions directory
244 set IN_PTpath = $INPDIR/ptse #> point source emissions input directory
245 set IN_LTpath = #> lightning NOx input directory
246 set METpath = $INPDIR/mcip/$APPL/$GRID_NAME #> meteorology input directory
247 #set JVALpath = $INPDIR/jproc #> offline photolysis rate table directory
248 set OMIpath = $BLD #> ozone column data for the photolysis model
249 set LUpath = $INPDIR/land #> BELD landuse data for windblown dust model
250 set SZpath = $INPDIR/land #> surf zone file for in-line seaspray emissions
逐日之設定 逐日迴圈之準備及啟動CCTM.exe為每24小時模擬的方式來進行,在同一批次內的4天中,依序進行4次的CCTM.exe。批次檔內以當日(TODAY)及昨日(YESTERDAY)2日輪轉方式來簡化換日的檔名問題。設定內容包括:YYYYMMDD、YYYYMM、YYMMDD、YYYYJJJ等多種不同的日期表示方式,以因應檔名的方便性STKCASEE為排放清冊標記,如有排放量版本之差異,應在此設定。此標記將會出現在輸出檔名之尾端。範例252 # =====================================================================
253 #> Begin Loop Through Simulation Days
254 # =====================================================================
255 set rtarray = ""
256
257 set TODAYG = ${START_DATE}
258 set TODAYJ = `date -ud "${START_DATE}" +%Y%j` #> Convert YYYY-MM-DD to YYYYJJJ
259 set START_DAY = ${TODAYJ}
260 set STOP_DAY = `date -ud "${END_DATE}" +%Y%j` #> Convert YYYY-MM-DD to YYYYJJJ
261 set NDAYS = 0
262 set STKCASEE = 10 # 12US1_cmaq_cb6_2016ff_16j # Stack Emission Version Label
263
264 while ($TODAYJ <= $STOP_DAY ) #>Compare dates in terms of YYYYJJJ
265
266 set NDAYS = `echo "${NDAYS} + 1" | bc -l`
267
268 #> Retrieve Calendar day Information
269 set YYYYMMDD = `date -ud "${TODAYG}" +%Y%m%d` #> Convert YYYY-MM-DD to YYYYMMDD
270 set YYYYMM = `date -ud "${TODAYG}" +%Y%m` #> Convert YYYY-MM-DD to YYYYMM
271 set YYMMDD = `date -ud "${TODAYG}" +%y%m%d` #> Convert YYYY-MM-DD to YYMMDD
272 set YYYYJJJ = $TODAYJ
273
274 #> Calculate Yesterday's Date
275 set YESTERDAY = `date -ud "${TODAYG}-1days" +%Y%m%d` #> Convert YYYY-MM-DD to YYYYJJJ
276
277 # =====================================================================
278 #> Set Output String and Propagate Model Configuration Documentation
279 # =====================================================================
280 echo ""
281 echo "Set up input and output files for Day ${TODAYG}."
282
283 #> set output file name extensions
284 setenv CTM_APPL ${RUNID}_${YYYYMMDD}_${GRID_NAME}_${STKCASEE}
285
286 #> Copy Model Configuration To Output Folder
287 if ( ! -d "$OUTDIR" ) mkdir -p $OUTDIR
288 cp $BLD/CCTM_${VRSN}.cfg $OUTDIR/CCTM_${CTM_APPL}.cfg
289
逐日之IC/BC/臭氧/光學/氣象檔案執行前逐一準備好檔案。目錄、檔名「不建議」修改,以避免錯誤。 設定內容 包括:初始檔案:RUN5會讀ICON檔案(須事先處理,見CMAQ初始及邊界條件設定 ),其他批次則會讀前一天的CCTM_CGRID檔案,或前一批次最後天的CCTM_CGRID檔案邊界檔案:BCON檔案,見CMAQ初始及邊界條件設定氣象檔案:並不逐日變化,但mcip與wrfout有1小時之差異,須以add_firstHr補滿第1小時。(見筆記全球空品模擬結果作為CMAQ之初始及邊界條件 附錄)範例290 # =====================================================================
291 #> Input Files (Some are Day-Dependent)
292 # =====================================================================
293
294 #> Initial conditions
295 if ($NEW_START == true || $NEW_START == TRUE ) then
296 setenv ICFILE ICON_${VRSN}_${APPL}_${GRID_NAME}_${ICTYPE}_${YYYYMMDD}
297 if ( $MO != '01' ) then
298 set PYM = `date -ud "20${APYM}01 -1 month" +%y%m`
299 set PRUNID = ${VRSN}_${compilerString}_${PYM}_run12
300 set POUTDIR = ${OUTDIR}/../output_CCTM_${PRUNID}
301 ln -sf ${POUTDIR}/CCTM_CGRID_${PRUNID}_${YESTERDAY}_${GRID_NAME}_${STKCASEE}.nc \
302 $ICpath/$ICFILE
303 endif
304 setenv INIT_MEDC_1 notused
305 setenv INITIAL_RUN Y #related to restart soil information file
306 else
307 set ICpath = $OUTDIR
308 setenv ICFILE CCTM_CGRID_${RUNID}_${YESTERDAY}_${GRID_NAME}_${STKCASEE}.nc
309 setenv INIT_MEDC_1 $ICpath/CCTM_MEDIA_CONC_${RUNID}_${YESTERDAY}_${GRID_NAME}_${STKCASEE}.nc
310 setenv INITIAL_RUN N
311 if ( $TODAYJ == $START_DAY ) then
312 @ PRUN = $RUN - 1
313 set PRUNID = ${VRSN}_${compilerString}_${APYM}_run${PRUN}
314 set POUTDIR = ${OUTDIR}/../output_CCTM_${PRUNID}
315 ln -sf ${POUTDIR}/CCTM_CGRID_${PRUNID}_${YESTERDAY}_${GRID_NAME}_${STKCASEE}.nc \
316 $OUTDIR/CCTM_CGRID_${RUNID}_${YESTERDAY}_${GRID_NAME}_${STKCASEE}.nc
317 ln -sf ${POUTDIR}/CCTM_MEDIA_CONC_${PRUNID}_${YESTERDAY}_${GRID_NAME}_${STKCASEE}.nc \
318 $OUTDIR/CCTM_MEDIA_CONC_${RUNID}_${YESTERDAY}_${GRID_NAME}_${STKCASEE}.nc
319 endif
320 endif
321
322 #> Boundary conditions
323 set BCFILE = BCON_${VRSN}_${APPL}_${ICTYPE}_${STD_YMD}_${GRID_NAME}
324
325 #> Off-line photolysis rates
326 #set JVALfile = JTABLE_${YYYYJJJ}
327
328 #> Ozone column data
329 set OMIfile = OMI_1979_to_2017.dat
330
331 #> Optics file
332 set OPTfile = PHOT_OPTICS.dat
333
334 #> MCIP meteorology files
335 setenv GRID_BDY_2D $METpath/GRIDBDY2D_${APPL}.nc # GRID files are static, not day-specific
336 setenv GRID_CRO_2D $METpath/GRIDCRO2D_${APPL}.nc
337 setenv GRID_CRO_3D $METpath/GRIDCRO3D_${APPL}.nc
338 setenv GRID_DOT_2D $METpath/GRIDDOT2D_${APPL}.nc
339 setenv MET_CRO_2D $METpath/METCRO2D_${APPL}.nc
340 setenv MET_CRO_3D $METpath/METCRO3D_${APPL}.nc
341 setenv MET_DOT_3D $METpath/METDOT3D_${APPL}.nc
342 setenv MET_BDY_3D $METpath/METBDY3D_${APPL}.nc
343 setenv LUFRAC_CRO $METpath/LUFRAC_CRO_${APPL}.nc
344
逐日排放相關檔案之設定為模式修改的重點位置。應注意版本的差異 設定內容 包括:地面及高空排放源都可以分群組提供,以解析各群組的貢獻,進一步解析污染源敏感性。如不提供,只須在等號右邊保持空白即可。面源須每網格提供,點源檔案分為常數部分(const)及變數部分(timvar),且可以每日不同煙囪個數。建議單月仍然保持相同煙囪個數,以降低複雜度。OBS系列變數:為測站測值,進行濃度調整之用。其檔案格式與點源相同。如不執行修正,可免提供。非官方程式碼,使用請自行注意。雖然逐日提供排放量檔案可以詳細在CTM檔案中進行檢查,但是否一定須逐日提供,是否可以像mcip檔案批次提供,此點還有待實證。EMISSCTRL_NML之設定此檔案會對排放源檔案做對照、加權不必再重做SMOKE,直接修改此檔,即可得特定污染源的到增、減量效果逐日排放相關檔案之設定範例345 #> Emissions Control File
346 setenv EMISSCTRL_NML ${BLD}/EmissCtrl_${MECH}.nml
347
348 #> Spatial Masks For Emissions Scaling
349 setenv CMAQ_MASKS # $SZpath/12US1_surf_bench.nc #> horizontal grid-dependent surf zone file
350
351 #> Gridded Emissions Files
352 setenv N_EMIS_GR 1
353 set EMISfile = area_${GRID_NAME}.${YYYYMMDD}.nc
354 setenv GR_EMIS_001 ${EMISpath}/${APYM}/${EMISfile}
355 setenv GR_EMIS_LAB_001 GRIDDED_EMIS
356
357 set EMISfile =
358 setenv GR_EMIS_002 ${EMISpath2}/${EMISfile}
359 setenv GR_EMIS_LAB_002 GRIDDED_RWC
360
361 #> In-line point emissions configuration
362 setenv N_EMIS_PT 1 #> Number of elevated source groups
363 setenv N_OBSS_PT 1 #> Number of elevated source groups
364 setenv STN_GRPS_001 $CMAQ_DATA/sites/${APYM}/const.nc
365 setenv STN_OBSS_001 $CMAQ_DATA/sites/${APYM}/EPA57.${YYYYMMDD}.timvar.nc_12
366 setenv STN_OBSS_LAB_001 TWN_EPA_57
367 setenv STN_OBSS_DIAG_001
368
369 set STKCASEG = teds # 12US1_2016ff_16j # Stack Group Version Label
370 # set STKCASEE = 10 # 12US1_cmaq_cb6_2016ff_16j # Stack Emission Version Label
371
372 # Time-Independent Stack Parameters for Inline Point Sources
373 setenv STK_GRPS_001 $IN_PTpath/${STKCASEG}${STKCASEE}.${APYM}.const.nc
374 setenv LAYP_STTIME $STTIME
375 setenv LAYP_NSTEPS $NSTEPS #umber of time steps for calculating elevated-point-source emissions.
376
377 # Emission Rates for Inline Point Sources
378 setenv STK_EMIS_001 $IN_PTpath/${APYM}/${STKCASEG}${STKCASEE}.${YYYYMMDD}.timvar.nc_12
379 setenv LAYP_STDATE $YYYYJJJ
380
381 # Label Each Emissions Stream
382 setenv STK_EMIS_LAB_001 POINT_NONEGU
383 setenv STK_EMIS_LAB_002 # POINT_EGU
384 setenv STK_EMIS_LAB_003 # POINT_OTHER
385 setenv STK_EMIS_LAB_004 # POINT_AGFIRES
386 setenv STK_EMIS_LAB_005 # POINT_FIRES
387 setenv STK_EMIS_LAB_006 # POINT_OTHFIRES
387 setenv STK_EMIS_LAB_006 # POINT_OTHFIRES
388 setenv STK_EMIS_LAB_007 # POINT_OILGAS
389 setenv STK_EMIS_LAB_008 # POINT_CMV
390
391 # Stack emissions diagnostic files
392 setenv STK_EMIS_DIAG_001 # 2DSUM
393 #setenv STK_EMIS_DIAG_002 2DSUM
394 #setenv STK_EMIS_DIAG_003 2DSUM
395 #setenv STK_EMIS_DIAG_004 2DSUM
396 #setenv STK_EMIS_DIAG_005 2DSUM
397
398
在線排放之設定呼應前述科學設定,此處提供每日數據。在線排放會逐時讀取氣象數據(雲雨、風速、氣溫等)來計算當時的排放量,此部分尚未測試。設定內容 包括:雷電NOx生物源風吹沙海洋飛沫雙向氨氣吸收及排放(似一定要啟動)範例399 #> Lightning NOx configuration
400 if ( $CTM_LTNG_NO == 'Y' ) then
401 setenv LTNGNO "InLine" #> set LTNGNO to "Inline" to activate in-line calculation
402
403 #> In-line lightning NOx options
404 setenv USE_NLDN N #> use hourly NLDN strike file [ default: Y ]
405 if ( $USE_NLDN == Y ) then
406 setenv NLDN_STRIKES ${IN_LTpath}/NLDN.12US1.${YYYYMMDD}_bench.nc
407 endif
408 setenv LTNGPARMS_FILE ${IN_LTpath}/LTNG_AllParms_12US1_bench.nc #> lightning parameter file
409 endif
410
411 #> In-line biogenic emissions configuration
412 setenv BIOSW_YN N
413 if ( $CTM_BIOGEMIS == 'Y' ) then
414 set IN_BEISpath = ${INPDIR}/land
415 setenv GSPRO $BLD/gspro_biogenics.txt
416 # setenv B3GRD $IN_BEISpath/b3grd_bench.nc
417 setenv B3GRD $IN_BEISpath/b3gts_l.20190920-28.1.sChina_27k.teds10_16.ncf
418 setenv BIOSW_YN Y #> use frost date switch [ default: Y ]
419 setenv BIOSEASON $IN_BEISpath/bioseason.cmaq.201909_sChina.ncf
420 #> ignore season switch file if BIOSW_YN = N
421 setenv SUMMER_YN Y #> Use summer normalized emissions? [ default: Y ]
422 setenv PX_VERSION Y #> MCIP is PX version? [ default: N ]
423 setenv SOILINP $OUTDIR/CCTM_SOILOUT_${RUNID}_${YESTERDAY}.nc
424 #> Biogenic NO soil input file; ignore if INITIAL_RUN = Y
425 endif
426
427 #> Windblown dust emissions configuration
428 if ( $CTM_WB_DUST == 'Y' ) then
429 # Input variables for BELD3 Landuse option
430 setenv DUST_LU_1 $LUpath/beld3_12US1_459X299_output_a_bench.nc
431 setenv DUST_LU_2 $LUpath/beld4_12US1_459X299_output_tot_bench.nc
432 endif
433
434 #> In-line sea spray emissions configuration
435 setenv OCEAN_1 $SZpath/12US1_surf_bench.nc #> horizontal grid-dependent surf zone file
436
437 #> Bidirectional ammonia configuration
438 if ( $CTM_ABFLUX == 'Y' ) then
439 setenv E2C_SOIL ${LUpath}/${APYM}/2016_${GRID_NAME}_soil_bench${APYM}.nc
440 setenv E2C_CHEM ${LUpath}/${APYM}/2016_${GRID_NAME}_time${YYYYMMDD}_bench.nc
441 setenv E2C_CHEM_YEST ${LUpath}/${APYM}/2016_${GRID_NAME}_time${YESTERDAY}_bench.nc
442 setenv E2C_LU ${LUpath}/beld4.${GRID_NAME}.ncf
443 endif
444
程序分析之設定程序分析是空氣品質管理很重要的模擬工作,包括模式程序分析、來源分配解析、以及硫份追蹤3個部分,目前都還未開啟。 設定內容 包括:CTM_PROCAN 程序分析CTM_ISAM 來源分配解析STM_SO4TRACK硫份追蹤範例445 #> Inline Process Analysis
446 setenv CTM_PROCAN N #> use process analysis [ default: N]
447 if ( $?CTM_PROCAN ) then # $CTM_PROCAN is defined
448 if ( $CTM_PROCAN == 'Y' || $CTM_PROCAN == 'T' ) then
449 #> process analysis global column, row and layer ranges
450 # setenv PA_BCOL_ECOL "10 90" # default: all columns
451 # setenv PA_BROW_EROW "10 80" # default: all rows
452 # setenv PA_BLEV_ELEV "1 4" # default: all levels
453 setenv PACM_INFILE ${NMLpath}/pa_${MECH}.ctl
454 setenv PACM_REPORT $OUTDIR/"PA_REPORT".${YYYYMMDD}
455 endif
456 endif
457
458 #> Integrated Source Apportionment Method (ISAM) Options
459 setenv CTM_ISAM N
460 if ( $?CTM_ISAM ) then
461 if ( $CTM_ISAM == 'Y' || $CTM_ISAM == 'T' ) then
462 setenv SA_IOLIST ${WORKDIR}/isam_control.txt
463 setenv ISAM_BLEV_ELEV " 1 1"
464 setenv AISAM_BLEV_ELEV " 1 1"
465
466 #> Set Up ISAM Initial Condition Flags
467 if ($NEW_START == true || $NEW_START == TRUE ) then
468 setenv ISAM_NEW_START Y
469 setenv ISAM_PREVDAY
470 else
471 setenv ISAM_NEW_START N
472 setenv ISAM_PREVDAY "$OUTDIR/CCTM_SA_CGRID_${RUNID}_${YESTERDAY}.nc"
473 endif
474
475 #> Set Up ISAM Output Filenames
476 setenv SA_ACONC_1 "$OUTDIR/CCTM_SA_ACONC_${CTM_APPL}.nc -v"
477 setenv SA_CONC_1 "$OUTDIR/CCTM_SA_CONC_${CTM_APPL}.nc -v"
478 setenv SA_DD_1 "$OUTDIR/CCTM_SA_DRYDEP_${CTM_APPL}.nc -v"
479 setenv SA_WD_1 "$OUTDIR/CCTM_SA_WETDEP_${CTM_APPL}.nc -v"
480 setenv SA_CGRID_1 "$OUTDIR/CCTM_SA_CGRID_${CTM_APPL}.nc -v"
481
482 #> Set optional ISAM regions files
483 # setenv ISAM_REGIONS /work/MOD3EVAL/nsu/isam_v53/CCTM/scripts/input/RGN_ISAM.nc
484
485 endif
486 endif
487
488
489 #> Sulfur Tracking Model (STM)
490 setenv STM_SO4TRACK N #> sulfur tracking [ default: N ]
491 if ( $?STM_SO4TRACK ) then
492 if ( $STM_SO4TRACK == 'Y' || $STM_SO4TRACK == 'T' ) then
493
494 #> option to normalize sulfate tracers [ default: Y ]
495 setenv STM_ADJSO4 Y
496
497 endif
498 endif
499
輸出檔案名稱、層數檔名因與後處理有關,「不建議」任何修改高空層數如用做下一網格子系統的邊界值,「必須」有3D的輸出如果沒有下一網格子系統(IC/BC會需要完整層數),可以只輸出第一層,以減省磁碟機空間。無設定(整句加註):將會輸出所有層數1 1:只輸出第1層、地面層影響所及之設定VariableSettingdefaultLine#CONC_BLEV_ELEV” 1 1”all layers151ACONC_BLEV_ELEV” 1 1”all layers156APMDIAG_BLEV_ELEV(according to ACONC)NLAYS226PA_BLEV_ELEV“1 4”all levels452ISAM_BLEV_ELEV” 1 1”all levels463AISAM_BLEV_ELEV” 1 1”all levels464範例500 # =====================================================================
501 #> Output Files
502 # =====================================================================
503
504 #> set output file names
505 setenv S_CGRID "$OUTDIR/CCTM_CGRID_${CTM_APPL}.nc" #> 3D Inst. Concentrations
506 setenv CTM_CONC_1 "$OUTDIR/CCTM_CONC_${CTM_APPL}.nc -v" #> On-Hour Concentrations
507 setenv A_CONC_1 "$OUTDIR/CCTM_ACONC_${CTM_APPL}.nc -v" #> Hourly Avg. Concentrations
508 setenv MEDIA_CONC "$OUTDIR/CCTM_MEDIA_CONC_${CTM_APPL}.nc -v" #> NH3 Conc. in Media
509 setenv CTM_DRY_DEP_1 "$OUTDIR/CCTM_DRYDEP_${CTM_APPL}.nc -v" #> Hourly Dry Deposition
510 setenv CTM_DEPV_DIAG "$OUTDIR/CCTM_DEPV_${CTM_APPL}.nc -v" #> Dry Deposition Velocities
511 setenv B3GTS_S "$OUTDIR/CCTM_B3GTS_S_${CTM_APPL}.nc -v" #> Biogenic Emissions
512 setenv SOILOUT "$OUTDIR/CCTM_SOILOUT_${CTM_APPL}.nc" #> Soil Emissions
513 setenv CTM_WET_DEP_1 "$OUTDIR/CCTM_WETDEP1_${CTM_APPL}.nc -v" #> Wet Dep From All Clouds
514 setenv CTM_WET_DEP_2 "$OUTDIR/CCTM_WETDEP2_${CTM_APPL}.nc -v" #> Wet Dep From SubGrid Clouds
515 setenv CTM_PMDIAG_1 "$OUTDIR/CCTM_PMDIAG_${CTM_APPL}.nc -v" #> On-Hour Particle Diagnostics
516 setenv CTM_APMDIAG_1 "$OUTDIR/CCTM_APMDIAG_${CTM_APPL}.nc -v" #> Hourly Avg. Particle Diagnostics
517 setenv CTM_RJ_1 "$OUTDIR/CCTM_PHOTDIAG1_${CTM_APPL}.nc -v" #> 2D Surface Summary from Inline Photolysis
518 setenv CTM_RJ_2 "$OUTDIR/CCTM_PHOTDIAG2_${CTM_APPL}.nc -v" #> 3D Photolysis Rates
519 setenv CTM_RJ_3 "$OUTDIR/CCTM_PHOTDIAG3_${CTM_APPL}.nc -v" #> 3D Optical and Radiative Results from Photolysis
520 setenv CTM_SSEMIS_1 "$OUTDIR/CCTM_SSEMIS_${CTM_APPL}.nc -v" #> Sea Spray Emissions
521 setenv CTM_DUST_EMIS_1 "$OUTDIR/CCTM_DUSTEMIS_${CTM_APPL}.nc -v" #> Dust Emissions
522 setenv CTM_IPR_1 "$OUTDIR/CCTM_PA_1_${CTM_APPL}.nc -v" #> Process Analysis
523 setenv CTM_IPR_2 "$OUTDIR/CCTM_PA_2_${CTM_APPL}.nc -v" #> Process Analysis
524 setenv CTM_IPR_3 "$OUTDIR/CCTM_PA_3_${CTM_APPL}.nc -v" #> Process Analysis
525 setenv CTM_IRR_1 "$OUTDIR/CCTM_IRR_1_${CTM_APPL}.nc -v" #> Chem Process Analysis
526 setenv CTM_IRR_2 "$OUTDIR/CCTM_IRR_2_${CTM_APPL}.nc -v" #> Chem Process Analysis
527 setenv CTM_IRR_3 "$OUTDIR/CCTM_IRR_3_${CTM_APPL}.nc -v" #> Chem Process Analysis
528 setenv CTM_DRY_DEP_MOS "$OUTDIR/CCTM_DDMOS_${CTM_APPL}.nc -v" #> Dry Dep
529 setenv CTM_DRY_DEP_FST "$OUTDIR/CCTM_DDFST_${CTM_APPL}.nc -v" #> Dry Dep
530 setenv CTM_DEPV_MOS "$OUTDIR/CCTM_DEPVMOS_${CTM_APPL}.nc -v" #> Dry Dep Velocity
531 setenv CTM_DEPV_FST "$OUTDIR/CCTM_DEPVFST_${CTM_APPL}.nc -v" #> Dry Dep Velocity
532 setenv CTM_VDIFF_DIAG "$OUTDIR/CCTM_VDIFF_DIAG_${CTM_APPL}.nc -v" #> Vertical Dispersion Diagnostic
533 setenv CTM_VSED_DIAG "$OUTDIR/CCTM_VSED_DIAG_${CTM_APPL}.nc -v" #> Particle Grav. Settling Velocity
534 setenv CTM_LTNGDIAG_1 "$OUTDIR/CCTM_LTNGHRLY_${CTM_APPL}.nc -v" #> Hourly Avg Lightning NO
535 setenv CTM_LTNGDIAG_2 "$OUTDIR/CCTM_LTNGCOL_${CTM_APPL}.nc -v" #> Column Total Lightning NO
536 setenv CTM_VEXT_1 "$OUTDIR/CCTM_VEXT_${CTM_APPL}.nc -v" #> On-Hour 3D Concs at select sites
537
538 #> set floor file (neg concs)
539 setenv FLOOR_FILE ${OUTDIR}/FLOOR_${CTM_APPL}.txt
540
541 #> look for existing log files and output files
542 ( ls CTM_LOG_???.${CTM_APPL} > buff.txt ) >& /dev/null
543 ( ls ${LOGDIR}/CTM_LOG_???.${CTM_APPL} >> buff.txt ) >& /dev/null
544 set log_test = `cat buff.txt`; rm -f buff.txt
545
546 set OUT_FILES = (${FLOOR_FILE} ${S_CGRID} ${CTM_CONC_1} ${A_CONC_1} ${MEDIA_CONC} \
547 ${CTM_DRY_DEP_1} $CTM_DEPV_DIAG $B3GTS_S $SOILOUT $CTM_WET_DEP_1\
548 $CTM_WET_DEP_2 $CTM_PMDIAG_1 $CTM_APMDIAG_1 \
549 $CTM_RJ_1 $CTM_RJ_2 $CTM_RJ_3 $CTM_SSEMIS_1 $CTM_DUST_EMIS_1 $CTM_IPR_1 $CTM_IPR_2 \
550 $CTM_IPR_3 $CTM_IRR_1 $CTM_IRR_2 $CTM_IRR_3 $CTM_DRY_DEP_MOS \
551 $CTM_DRY_DEP_FST $CTM_DEPV_MOS $CTM_DEPV_FST $CTM_VDIFF_DIAG $CTM_VSED_DIAG \
552 $CTM_LTNGDIAG_1 $CTM_LTNGDIAG_2 $CTM_VEXT_1 )
553 if ( $?CTM_ISAM ) then
554 if ( $CTM_ISAM == 'Y' || $CTM_ISAM == 'T' ) then
555 set OUT_FILES = (${OUT_FILES} ${SA_ACONC_1} ${SA_CONC_1} ${SA_DD_1} ${SA_WD_1} \
556 ${SA_CGRID_1} )
557 endif
558 endif
559 set OUT_FILES = `echo $OUT_FILES | sed "s; -v;;g" | sed "s;MPI:;;g" `
560 ( ls $OUT_FILES > buff.txt ) >& /dev/null
561 set out_test = `cat buff.txt`; rm -f buff.txt
562
563 #> delete previous output if requested
564 if ( $CLOBBER_DATA == true || $CLOBBER_DATA == TRUE ) then
565 echo
566 echo "Existing Logs and Output Files for Day ${TODAYG} Will Be Deleted"
567
568 #> remove previous log files
569 foreach file ( ${log_test} )
570 #echo "Deleting log file: $file"
571 /bin/rm -f $file
572 end
573
574 #> remove previous output files
575 foreach file ( ${out_test} )
576 #echo "Deleting output file: $file"
577 /bin/rm -f $file
578 end
579 /bin/rm -f ${OUTDIR}/CCTM_EMDIAG*${RUNID}_${YYYYMMDD}.nc
580
581 else
582 #> error if previous log files exist
583 if ( "$log_test" != "" ) then
584 echo "*** Logs exist - run ABORTED ***"
585 echo "*** To overide, set CLOBBER_DATA = TRUE in run_cctm.csh ***"
586 echo "*** and these files will be automatically deleted. ***"
587 exit 1
588 endif
589
590 #> error if previous output files exist
591 if ( "$out_test" != "" ) then
592 echo "*** Output Files Exist - run will be ABORTED ***"
593 foreach file ( $out_test )
594 echo " cannot delete $file"
595 end
596 echo "*** To overide, set CLOBBER_DATA = TRUE in run_cctm.csh ***"
597 echo "*** and these files will be automatically deleted. ***"
598 exit 1
599 endif
600 endif
601
執行有關之環境設定 控制程式之環境變數設定起始日時、時間長度與步階、IC/BC/OMI及光解數據污染物名稱定義光化學反應與係數此段應無修改的必要,只要確定連結得到正確的檔案範例602 #> for the run control ...
603 setenv CTM_STDATE $YYYYJJJ
604 setenv CTM_STTIME $STTIME
605 setenv CTM_RUNLEN $NSTEPS
606 setenv CTM_TSTEP $TSTEP
607 setenv INIT_CONC_1 $ICpath/$ICFILE
608 setenv BNDY_CONC_1 $BCpath/$BCFILE
609 setenv OMI $OMIpath/$OMIfile
610 setenv OPTICS_DATA $OMIpath/$OPTfile
611 #setenv XJ_DATA $JVALpath/$JVALfile
612 set TR_DVpath = $METpath
613 set TR_DVfile = $MET_CRO_2D
614
615 #> species defn & photolysis
616 setenv gc_matrix_nml ${NMLpath}/GC_$MECH.nml
617 setenv ae_matrix_nml ${NMLpath}/AE_$MECH.nml
618 setenv nr_matrix_nml ${NMLpath}/NR_$MECH.nml
619 setenv tr_matrix_nml ${NMLpath}/Species_Table_TR_0.nml
620
621 #> check for photolysis input data
622 setenv CSQY_DATA ${NMLpath}/CSQY_DATA_$MECH
623
624 if (! (-e $CSQY_DATA ) ) then
625 echo " $CSQY_DATA not found "
626 exit 1
627 endif
628 if (! (-e $OPTICS_DATA ) ) then
629 echo " $OPTICS_DATA not found "
630 exit 1
631 endif
632
記憶體與路徑設定、執行程式執行程式所需電腦的環境資源,一般用unlimit控制,須視機器的情況來設定。此外mpi的路徑也不同,須考量EXE編譯程式版本配套使用。設定內容 包括:unlimit/limitmpi 路徑執行程式(line 656)完成當日工作,準備次日日期範例633 # ===================================================================
634 #> Execution Portion
635 # ===================================================================
636
637 #> Print attributes of the executable
638 if ( $CTM_DIAG_LVL != 0 ) then
639 ls -l $BLD/$EXEC
640 size $BLD/$EXEC
641 # unlimit
642 # limit
643 endif
644
645 #> Print Startup Dialogue Information to Standard Out
646 echo
647 echo "CMAQ Processing of Day $YYYYMMDD Began at `date`"
648 echo
649
650 #> Executable call for single PE, uncomment to invoke
651 #( /usr/bin/time -p $BLD/$EXEC ) |& tee buff_${EXECUTION_ID}.txt
652
653 #> Executable call for multi PE, configure for your system
654 set MPI = /usr/local/bin
655 set MPIRUN = $MPI/mpirun
656 ( /usr/bin/time -p mpirun --oversubscribe -np $NPROCS $BLD/$EXEC ) |& tee buff_${EXECUTION_ID}.txt
657
658 #> Harvest Timing Output so that it may be reported below
659 set rtarray = "${rtarray} `tail -3 buff_${EXECUTION_ID}.txt | grep -Eo '[+-]?[0-9]+([.][0-9]+)?' | head -1` "
660 rm -rf buff_${EXECUTION_ID}.txt
661
662 #> Abort script if abnormal termination
663 if ( ! -e $OUTDIR/CCTM_CGRID_${CTM_APPL}.nc ) then
664 echo ""
665 echo "**************************************************************"
666 echo "** Runscript Detected an Error: CGRID file was not written. **"
667 echo "** This indicates that CMAQ was interrupted or an issue **"
668 echo "** exists with writing output. The runscript will now **"
669 echo "** abort rather than proceeding to subsequent days. **"
670 echo "**************************************************************"
671 break
672 endif
673
674 #> Print Concluding Text
675 echo
676 echo "CMAQ Processing of Day $YYYYMMDD Finished at `date`"
677 echo
678 echo "\\\\\=====\\\\\=====\\\\\=====\\\\\=====/////=====/////=====/////=====/////"
679 echo
680
681 # ===================================================================
682 #> Finalize Run for This Day and Loop to Next Day
683 # ===================================================================
684
685 #> Save Log Files and Move on to Next Simulation Day
686 mv CTM_LOG_???.${CTM_APPL} $LOGDIR
687 if ( $CTM_DIAG_LVL != 0 ) then
688 mv CTM_DIAG_???.${CTM_APPL} $LOGDIR
689 endif
690
691 #> check for first run and first day and d02
692 #>duplicate the CCTM_ACON result for preceding day for preparing BCON_d04
693 if ( $TODAYJ == $START_DAY && $RUN == 5 && $DM == 'd02' ) then
694 cd $OUTDIR
695 python ../icon/ACONC31-30.py CCTM_ACONC_${CTM_APPL}.nc
696 cd $CMAQ_HOME
697 endif
698
699 #> The next simulation day will, by definition, be a restart
700 setenv NEW_START false
701
702 #> Increment both Gregorian and Julian Days
703 set TODAYG = `date -ud "${TODAYG}+1days" +%Y-%m-%d` #> Add a day for tomorrow
704 set TODAYJ = `date -ud "${TODAYG}" +%Y%j` #> Convert YYYY-MM-DD to YYYYJJJ
705
706 end #Loop to the next Simulation Day
批次時間報表707
708 # ===================================================================
709 #> Generate Timing Report
710 # ===================================================================
711 set RTMTOT = 0
712 foreach it ( `seq ${NDAYS}` )
713 set rt = `echo ${rtarray} | cut -d' ' -f${it}`
714 set RTMTOT = `echo "${RTMTOT} + ${rt}" | bc -l`
715 end
716
717 set RTMAVG = `echo "scale=2; ${RTMTOT} / ${NDAYS}" | bc -l`
718 set RTMTOT = `echo "scale=2; ${RTMTOT} / 1" | bc -l`
719
720 echo
721 echo "=================================="
722 echo " ***** CMAQ TIMING REPORT *****"
723 echo "=================================="
724 echo "Start Day: ${START_DATE}"
725 echo "End Day: ${END_DATE}"
726 echo "Number of Simulation Days: ${NDAYS}"
727 echo "Domain Name: ${GRID_NAME}"
728 echo "Number of Grid Cells: ${NCELLS} (ROW x COL x LAY)"
729 echo "Number of Layers: ${NZ}"
730 echo "Number of Processes: ${NPROCS}"
731 echo " All times are in seconds."
732 echo
733 echo "Num Day Wall Time"
734 set d = 0
735 set day = ${START_DATE}
736 foreach it ( `seq ${NDAYS}` )
737 # Set the right day and format it
738 set d = `echo "${d} + 1" | bc -l`
739 set n = `printf "%02d" ${d}`
740
741 # Choose the correct time variables
742 set rt = `echo ${rtarray} | cut -d' ' -f${it}`
743
744 # Write out row of timing data
745 echo "${n} ${day} ${rt}"
746
747 # Increment day for next loop
748 set day = `date -ud "${day}+1days" +%Y-%m-%d`
749 end
750 echo " Total Time = ${RTMTOT}"
751 echo " Avg. Time = ${RTMAVG}"
752 #source run_cctmMM_RR_DMX.csh 01 ${RUN} d04
753 exit
Referencehttps://forum.cmascenter.org/t/cctm-stops-running-on-long-simulation/292/12Example https://github.com/sinotec2/cmaq_relatives/blob/master/run_cctmMM_RR_DM.cshNotesCMAQ compilationsCMAQ初始及邊界條件設定Here: CCTM之設定Back to topCopyright © 2021- Dr. Yunchuan Kuang. Distributed by an MIT license. views site_visits visitors Page last modified: Apr 20 2022 at 08:28 PM. Edit this page at GitHub
-1.5%����
1 0 obj
<>>>
endobj
2 0 obj
<>
endobj
3 0 obj
<>/ProcSet[/PDF/Text/ImageB/ImageC/ImageI] >>/MediaBox[ 0 0 595.32 841.92] /Contents 4 0 R/Group<>/Tabs/S/StructParents 0>>
endobj
4 0 obj
<>
stream
x��YKo�6����>�]��C Hj)�bX�i��1�L�O����SE��"%Qm<��cU��:4�}����[u<�����|���th���}��������r����ח�o�_O��=wU�|w�
����wX�Ê U�E��?���*��p�����u�3�P��x����C5[Y�?uUu�Y�����\A�Ъ�&(Ц�T�#d
���#8����T��he�Na���?��z�{E"��͑vħ�9Si���3�W���7#�Q�ޞ��6��-ȯ֦�džl/��Y��z��s��k�*G�)=�t_�V9�o\=]��[�A��W��ڸ"�B}.B�U����Di4�, ���BRkKa��M���U���a�ܖ��̴A�rޫP_�8nAiÂm��5zci<��] -x �H6���.�;P��A����~y�Z=|y}��spv ���N�/�y�ȵB蘡���W���vs�#b�hoaWF�����U����������J(�����[z���
![�6��B qE����N)�eA�sJv�& jM�� ����g�U��`�g��&��[�w�� .4`��з�p���86qp���$�z��`ϼ�����7m���j��\'?yy�ڀ�g��ֈ�0t!��z�_����j��a ca�B��h��7 �C ێ��3�}��a�T�����xq>{yY��(8�,<�MlJ$�r�rc�
˾f
*�#�p�GtY��ly���-�R3$$,B�-��ir��R��r��D�` �v2 nR F��g=��C��Ѽ=N����l�f����:�3���Y4��kwT���#hf�b�`;���d+Y˻M��D<�� bB[��p涘�&�~��z�a<�����P��%��09T)N��o,���&�<<����}|"D���[�����h� ��l��"�YAဪ�MQ������X�]C�hNz5fz��㖰n�])�mt1��'G�(��_ق(��t6T*��k`�
Jί)�5 ָ��(��0Í[�ռԇ$�Qq ���C�:��{�z,�YTXXMY�5W1\p��BM�7�"ؚ�L1�8=; �,�!�oW��a4���((����'��6�qﯵN|��!ڮ
e0co
�%�5w_�9�2[.�V�}��L|V���X���s�� ǔ���z��$�PHPGF�v%4Zu3vC����jN��P<Z�5�'��,�@�nо4UsP�m�v/8�gHmcQ��[L2d��s�Ӹ�H�� do��n�kY�n�W�c^��u[s �Y0���<����XJ��{�C��(�\ݮ��h�?֖�e��Zl��� 2Y�� r�v�AN��4� ��9�x%�����D�4o�J �&K=�� :E~�ߒ��_1�D�f�#m���'tc~�6���B:Cy��x���$)���a|�f���d!��=���NrI��I6�8_�ilFF8����'�e��q0�
���W��\����h������F��326�8��zJ�8~�tz�T��a-�&j���S����.�@=7��� ���`��&R����٨m��~�pI˽��9���In��3]���=̆�C脁i�-���#%�- :[>)�E,7�!�$_��с�7�6��5����#r�6�kü���:�U�+[�5.�R���
yC��P�ʼ�ɛ���Sì7�-��S�V�
S+J�8���#F��gT�y��0t���N8�[Ӧ��Л�f�v'm&k�#i3��PڜS�ϕ��DZ��is����o(\Ǵ�l������.Nz�q�M��T.����%[��/!�'�9[;��P�^Dec�~ݼ]������ٯ�ɠ ��H}�ǤQj�Jϥ���8"�ƀ?'X����:{/�rC}]s:��|� dm�
endstream
endobj
5 0 obj
<>
endobj
6 0 obj
<>
endobj
7 0 obj
<>
endobj
8 0 obj
[ 9 0 R]
endobj
9 0 obj
<>
endobj
10 0 obj
<>
endobj
11 0 obj
<>
endobj
12 0 obj
<>
endobj
13 0 obj
<>
endobj
14 0 obj
<>
endobj
15 0 obj
[ 16 0 R]
endobj
16 0 obj
<>
endobj
17 0 obj
<>
endobj
18 0 obj
<>
endobj
19 0 obj
<>
endobj
20 0 obj
<>
endobj
21 0 obj
<>
endobj
22 0 obj
[ 23 0 R]
endobj
23 0 obj
<>
endobj
24 0 obj
<>
endobj
25 0 obj
<>
endobj
26 0 obj
<>/ProcSet[/PDF/Text/ImageB/ImageC/ImageI] >>/MediaBox[ 0 0 595.32 841.92] /Contents 27 0 R/Group<>/Tabs/S/StructParents 1>>
endobj
27 0 obj
<>
stream
x��YKoG���sp��w��~.6H ���A�:^"�v��C��g�{�ݚ �aK��E���c���^����饹��k_^�>��wݷ��o��|���>�}x�������燇���{�����I��/�7� ��� ���7����@�������w������ooF�d�LP�2�5��f��o��l�?����01�������{� ��Ո:=�{H>B�=�S|x�4�^�[����Ү#õݒt�5I�l*���"ȾnC���S-Ҽ��y�D��ZGO�U^85�8���i�X�~���O�_�w������w����
�m�+�7�$��:�6��A~y֮�r:H&-%��ד�^꼇~���� �N�������s�tߎ%�6���ݷ2*l| �*�5�, ��bt�tfzx�!�#=.��Le�O{j4heS���E��0�*����_��v�T�QH�M��+@oP��S7"`?�{�v�H��
��-��㸀 [+)�-��X6��;jH!�y|��t�@�rg0<�H�8��4������;�
��g��FJKw\̭.Fsҟ�7ӌ0�~TR`Tj�w4Dt݉���*�%y��T�o+d���Vs7���T���- ���̓%��Vz����6T�B�憋�����%�+��q��v���T܄@�,A��3j�y�T���w=L��������S]�ۻ��oF!M]qqwA!�@ؠY�7�����g���K1�R/��ؗ���@,��"MeW=�]zyY���Z�4��g�>��_�}�F�ϙ�����6�\8P�).7��M�z\<�d��n�N*������As��4�0Թ�;��FË�*T[�