Add a Physics Scheme into WRF Model.ppt

1、Add a Physics Scheme into WRF Model,Shu-hua Chen UC Davis/AFWA,Physics implementation features Adding a physics scheme,Domain size: ids, ide, jds, jde, kds, kde Memory size: ims, ime, jms, jme, kms, kme Tile size: its, ite, jts, jte, kts, kte,Three Sets of Dimensions,Rules for WRF Physics,Naming rul

2、es,xxx = individual scheme ex, module_cu_kf.F,Naming Rules,module_yy_xxx.F (module),yy = ra is for radiation bl is for PBL cu is for cumulus mp is for microphysics.,YY = ra is for radiation bl is for PBL cu is for cumulus,Naming Rules,RXXYYTEN (tendencies),XX = variable (th, u, v, qv, qc, ),ex, RTHB

3、LTEN,Coding rules (later),Rules for WRF Physics,Naming rules,One scheme one module,WRF Physics Features, Unified global constatnts (module_model_constants.F), Unified global constatnts (module_model_constants.F), Vertical index (kms is at the bottom), Unified common calculations (saturation mixing r

4、atio),WRF Physics Features,WRF Language, 4D Moisture field, moist(i,k,j,?),? = P_QV (water vapor) P_QC (cloud water) P_QI (cloud ice) P_QR (rain) P_QS (snow) P_QG (graupel) (module_state_description.F),WRF Language, 4D Moisture field, moist(i,k,j,?), PARAM_FIRST_SCALAR,IF ( P_QI .ge. PARAM_FIRST_SCA

5、LAR ) (the memory of cloud ice is allocated) . . .,Prepare your code,Create a new module,Declare new variables and a new package in Registry,Modify solve_em.F (solve_eh.F),Implement a new physics scheme,Do initialization,Modify namelist,Modify phy_prep (module_em.F),Modify physics_,Modify c

6、umulus_driver.F,Modify Makefile,Compile and test,Modify calculate_phy_ten (module_em.F),Implement a new physics scheme,Modify phy_cu_ten (module_physics_addtendc.F),Replace continuation characters in the 6th column with f90 continuation & at end of previous line,1.F90,Prepare your code,Replace conti

7、nuation characters in the 6th column with f90 continuation & at end of previous line,1.F90,b)Replace the 1st column C for comment with !,Prepare your code,1.F90,2.No common block,Prepare your code,1.F90,2.No common block,3.Use “ implicit none ”,4.Use “ intent ”,Prepare your code,Subroutine sub(T,q,p

8、, .) real,intent(out), & dimension(ims:ime,kms:kme,jms:jme): T real,intent( in), & dimension(ims:ime,kms:kme,jms:jme): q real,intent(inout), & dimension(ims:ime,kms:kme,jms:jme): p,1.F90,2.No common block,3.Use “ implicit none ”,4.Use “ intent ”,5.Variable dimensions,Prepare your code,Subroutine sub

9、(glo,.) real,intent(out), & dimension(ims:ime,kms:kme,jms:jme): glo real,dimension(its:ite,kts:kte,jts:jte): loc,1.F90,2.No common block,3.Use “ implicit none ”,4.Use “ intent ”,5.Variable dimensions,6.Do loops,Prepare your code,Implement a new physics scheme,Create a new module,ex, module_cu_chen.F

10、 (put all your codes in),Go Registry and declare a new package (and new variables) (WRFV1/Registry),package chenscheme cu_physics=3 - -,package kfscheme cu_physics=1 - - package bmjscheme cu_physics=2 - -,Implement a new physics scheme,Cloud microphysics package kesslerscheme mp_physics=1 - moist:qv

11、,qc,qr package linscheme mp_physics=2 - moist:qv,qc,qr,qi,qs,qg package ncepcloud3 mp_physics=3 - moist:qv,qc,qr package ncepcloud5 mp_physics=4 - moist:qv,qc,Create a new module,ex, module_cu_chen.F (put all your codes in),Go Registry and declare a new package (and new variables) (WRFV1/Registry),I

12、mplement a new physics scheme,Modify namelist.input and assign,cu_physics = 3,Create a new module,ex, module_cu_chen.F (put all your codes in),Go Registry and declare a new package (and new variables) (WRFV1/Registry),INIT,WRF,.,solve_em,phy_init,start_domain_em,cu_init,(dyn_em) (start_em.F),(phys)

13、(module_physics_init.F),(dyn_em),phys/module_physics_init.F,Pass new variables down to cu_init,INIT,WRF,.,solve_em,phy_init,start_domain_em,cu_init,(dyn_em) (start_em.F),(phys) (module_physics_init.F),(dyn_em),phys/module_physics_init.F,Go subroutine cu_init Include the new module and create a new S

14、ELECT case,Pass new variables down to cu_init,cps_select: SELECT CASE(config_flags%cu_physics) CASE (KFSCHEME) CALL kfinit(.) CASE (BMJSCHEME) CALL bmjinit(.) CASE DEFAULT END SELECT cps_select,Subroutine cu_init() . USE module_cu_kf USE module_cu_bmj .,CASE (CHENSCHEME) CALL cheninit(.),phys/module

15、_physics_init.F,USE module_cu_chen,WRF,solve_em,DYNAMICS,phy_init,INIT,.,.,phy_prep,moist_physics_prep,phy_prep/moist_physics_prep, Calculate required variables, Convert variables from C grid to A grid,WRF,solve_em,DYNAMICS,phy_init,INIT,.,.,microphysics_driver,radiation_driver,cumulus_driver,pbl_dr

16、iver,phy_prep,moist_physics_prep,chencps,Three-level structure,solve_em,Go physics driver (cumulus_driver.F) Include the new module and create a new SELECT CASE in driver Check available variables in drivers (variables are explained inside drivers),cumulus_driver.F,Subroutine cumulus_driver . USE mo

17、dule_cu_kf USE module_bmj_kf .,cumulus_driver.F,USE module_cu_chen,cps_select: SELECT CASE(config_flags%cu_physics) CASE (KFSCHEME) CALL KFCPS(.) CASE (BMJSCHEME) CALL BMJCPS(.) CASE DEFAULT END SELECT cps_select,CASE (CHENSCHEME) CALL CHENCPS(.),SUBROUTINE cumulus_driver(arg1, arg2, & newvar1, newv

18、ar2, & its,ite,jts,jte,kts,kte, & ims,ime,jms,jme,kms,kme, & ids,ide,jds,jde,kds,kde ) INTEGER, INTENT(IN) : its,ite,jts,jte,kts,kte, & ims,ime,jms,jme,kms,kme, & ids,ide,jds,jde,kds,kde REAL, INTENT(IN) : arg1, arg2 REAL, INTENT(OUT), DIMENSION(kms:kme) : & newvar1,newvar2,.,Physics_,solve_em,cumulus_driver,c