Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

Cumulus scheme support

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

neel14

Member
Hi,
I have some confusion regarding the convective transport by the schemes. Do Grell-Freitas scheme support convective transport of chem species in WRF-Chem?
There is a pull request for implementation of GF scheme: https://github.com/wrf-model/WRF/pull/1121

However if you see this part of code from chem_driver.F

if( config_flags%cu_physics>0 .and. config_flags%chem_conv_tr>0 &
.and. config_flags%cu_physics/=kfcupscheme ) then
call wrf_debug(15,'calling conv transport for chemical species')
if(config_flags%chem_opt >0 )then

DO nv=PARAM_FIRST_SCALAR,num_chem_ct
chem_old(its:ite,kts:kte,jts:jte,nv) = chem(its:ite,kts:kte,jts:jte,chem_ct_indices(nv))
ENDDO
call grelldrvct(grid%DT,ktau,grid%DX, &
rho,grid%RAINCV_B,chem,
- - - - - - - - - - - - - - - - - --

chem_conv_tr=1 supported by Grell-freitas scheme, calls the grelldrvct from module_ctrans_grell.F which is for convective transport, wet deposition,etc.

Hence there is a conflict in my understanding whether GF scheme is implemented for convective transport,aqueous chem or not, as the pull request for its implementation and the code seem to be contradictory. It would be helpful if someone could clarify a bit. Appreciate your help & time.

Best
Neeldip
 
Hi,
It would be really helpful if someone from the admins/moderators of the forum could clarify on the issue.

Thanks.
 
Hi Neeldip,

This paper explains the changes: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6999733/

This improved scheme will be in the next release, but of course you can grab the code to use yourself.

Jordan
 
Hi Neeldip,

It has not yet been merged into the official release, but you should be able to grab the branch/files directly from the PR.

Jordan
 
Top