资源描述:
《add a physics scheme into wrf model》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、AddaPhysicsSchemeintoWRFModelShu-huaChenUCDavis/AFWAPhysicsimplementationfeaturesAddingaphysicsschemeDomainsize:ids,ide,jds,jde,kds,kdeMemorysize:ims,ime,jms,jme,kms,kmeTilesize:its,ite,jts,jte,kts,kteThreeSetsofDimensionsRulesforWRFPhysicsNamingrulesxxx=individualschemeex,module_cu_kf.FNa
2、mingRulesmodule_yy_xxx.F(module)yy=raisforradiationblisforPBLcuisforcumulusmpisformicrophysics.YY=raisforradiationblisforPBLcuisforcumulusNamingRulesRXXYYTEN(tendencies)XX=variable(th,u,v,qv,qc,…)ex,RTHBLTENCodingrules(later)RulesforWRFPhysicsNamingrulesOneschemeonemoduleWRFPhysicsFeatures
3、REAL,PARAMETER::r_d=287.REAL,PARAMETER::r_v=461.6REAL,PARAMETER::cp=7.*r_d/2.REAL,PARAMETER::cv=cp-r_d..•Unifiedglobalconstatnts(module_model_constants.F)•Unifiedglobalconstatnts(module_model_constants.F)•Verticalindex(kmsisatthebottom)•Unifiedcommoncalculations(saturationmixingratio)WRFPh
4、ysicsFeaturesWRFLanguage•4DMoisturefield,moist(i,k,j,?)?=P_QV(watervapor)P_QC(cloudwater)P_QI(cloudice)P_QR(rain)P_QS(snow)P_QG(graupel)(module_state_description.F)WRFLanguage•4DMoisturefield,moist(i,k,j,?)•PARAM_FIRST_SCALARIF(P_QI.ge.PARAM_FIRST_SCALAR)(thememoryofcloudiceisallocated)...
5、PrepareyourcodeCreateanewmoduleDeclarenewvariablesandanewpackageinRegistryModifysolve_em.F(solve_eh.F)ImplementanewphysicsschemeDoinitializationModifynamelistModifyphy_prep(module_em.F)Modifyphysics_drive.intModifycumulus_driver.FModifyMakefileCompileandtestModifycalculate_phy_ten(module_e
6、m.F)ImplementanewphysicsschemeModifyphy_cu_ten(module_physics_addtendc.F)Replacecontinuationcharactersinthe6thcolumnwithf90continuation`&‘atendofpreviousline1.F90PrepareyourcodeSubroutinekessler(QV,T,&its,ite,jts,jte,kts,kte,&ims,ime,jms,jme,kms,kme,&ids,ide,jds,jde,kds,kde)F77Subroutineke
7、ssler(QV,T,...&its,ite,jts,jte,kts,kte,&ims,ime,jms,jme,kms,kme,&ids,ide,jds,jde,kds,kde)F90Replacecontinuationcharactersinthe6thcolumnwithf90continuation`&‘atendofpreviousline1.F90b)Replacethe1stcolumn`C`forcommentwith`!`PrepareyourcodecThisisatestF77!Thisisa