[an error occurred while processing this directive]

themes

 
 
 
 
 

Physical Processes Group

From Daniele Rosa
Date June 21, 2012

Hi

reading the SPCAM code it seems that local and non-local diffusion is computed by the standard CAM 3.5 PBL parameterization and then the CRM is allowed to do more mixing on constituents in the PBL layers.

So SPCAM does not replace the conventional PBL mixing, instead it adds some.

Do I get it right?

Thanks
Daniele


From Mike Pritchard
Date June 21, 2012

Hi Daniele,

Awesome question.

My understanding from closing the T,Q,U & V budgets in the SPCAM3.0 code was that the PBL routine is bypassed for Q & T diffusion but is used to do U & V diffusion (at least in the standard configuration of a 2D CRM; if the CRM is 3D the PBL scheme is entirely bypassed and the resolved momentum transport from the 3D CRM is used instead). I can't fully recall if this was also the case for SPCAM3.5, but I thought so.

Actually, I think this might be an important lingering issue for some of the many SPCAM runs that have been run on greater than 26 vertical levels - Rich Neale has mentioned in the past that the CAM3.5.X PBL scheme has bugs that make it behave unphysically on anything other than the L26 vertical grid. So it's quite possible that SPCAM3.5 simulations that have used grids other than L26 may have suffered from unphysical diffusive momentum transport.

Cheers,

-Mike.

From Marat Khairoutdinov
Date June 21, 2012

Daniel,

The vertical countergradient PBL transport is switched off in src/physics/cam/diffusion_solver.F90 when CRM is on. Some CAM grid diffusion is left to maintain CAM's advection stability and for physical reasons (as for stable stratification, SAM is not that good in mixing and diffusing things around).

Marat

From Daniele Rosa
Date June 22, 2012

Marat:

How is it done in the code? I can't find any "#ifdef CRM" or anything else.

Is there an equivalent PBL diffusion in SAM when it is run by itself?

Daniele

From Marat Khairoutdinov
Date June 22, 2012

try to look for #ifndef CRM, or just CRM. -Marat



[an error occurred while processing this directive]