Error running Spectral Element and Finite Volume dycores for 80 levels and PK Stratosphere

[ag4680@nyu_edu]

Hi,

Following the helpful guidance from Isla and Peter on running the Polvani Kushner Stratosphere and Topography respectively in NCAR dynamical cores, I was able to integrate the f19 and ne16np4 grids with 40 levels in the vertical.

Now I am aiming to increase the vertical resolution to 80 levels. However, when I do so, I get an error the goes like :

Error : GET_PBUF1D_FIELD_BY_INDEX : index (-1) out of range

This has something to do with the physics buffers and is generated in the src/physics/cam/physics_buffer.F90.in file. The 'index' variable associated with the buffer is -1 and is creating this error. My first guess is that this could be due to some issue related to undefined parameters for 80 levels in the namelist_defaults.xml file, but ofc I can be wrong.

I tried running the models with increased resolution in the horizontal but with same vertical levels (i.e ne30np4 and f09) and they run fine without generating this error. Hence this had got to do with the number of vertical levels I am using. Any help on the matter is deeply appreciated.


Regards,

Aman

 

[ag4680@nyu_edu]

Hi,I was suggested that this error was fixed in the physgrid branch which has now been committed to the cam/trunk and that I should use this instead.I checked out the cam/trunk and tried running the models for 80 vertical levels, after merging my modules. However, I am still getting the same error. Can you please suggest how to resovle this physics buffer issue?Any help on this issue is deeply appreciated.Regards,Aman

 

[eaton]

This turns out to be a bug in build-namelist which is turning on the gravity wave parameterization whenever the number of levels is >=60.  I will add a fix to the CAM trunk shortly, but in the meantime the simple fix is to change the conditional in components/cam/bld/build-namelist fromif ($waccm_phys or $cfg->get('nlev') >= 60) {
    # Spectral gravity waves are part of WACCM physics, and also drive the
    # QBO in the high vertical resolution configuration.
    add_default($nl, 'use_gw_front'     , 'val'=>'.true.');
    add_default($nl, 'use_gw_convect_dp', 'val'=>'.true.');
To: if ($waccm_phys or
    (!$simple_phys and $cfg->get('nlev') >= 60)) {
    # Spectral gravity waves are part of WACCM physics, and also drive the
    # QBO in the high vertical resolution configuration.
    add_default($nl, 'use_gw_front'     , 'val'=>'.true.');
    add_default($nl, 'use_gw_convect_dp', 'val'=>'.true.');

 

[ag4680@nyu_edu]

Hi Brian,Thank you very much for your response on this. This solved the problem.Regards,Aman