IC file problem | Polvani Kushner Strat on a Spectral-Element dycore

[ag4680@nyu_edu]

Hi,After successfully implementing a Polvani Kushner 2002 stratosphere(simple physics. HS94) with simple wave-2 topography for the spectral numerics(t42l40), I am trying to implement the same for spectral-element numerics(ne16p4). The repository has IC files for 30 levels but not for a required 40 hybrid levels and hence I tried using the 'interpic_new' tool provided in the CESM2Beta6 version in order to carry out vertical interpolation. However, when I run the executable ./interpic, I get a segmentation fault in the interp_driver() -> wrap_get_var_double() -> nf_get_var_double() function during the call wrap_get_var_double(ncidi, vari%y_vid, y_in)So, as an alternative, I commented the function calls in the y-direction following which, I get an error that reads : ENDRUN:LININTERP: Non-monotonic input coordinate array found.While the 1-D variables are copied successfully, this error is caused due only for 3-D variables. I think this is because the spectral element core stores the 2D horizotal grid as a single 'non monotonic vector' and the interpic script is not meant to handle such a vector.I now believe that the ./interpic tool might not be able to help me accomplish this task. Is there any workaround to get an initial condition with user defined vertical levels?I have no prior experience with spectral element cores and I hope I am not overlooking anything en route. Any help on this matter is much appreciated.Regards,Aman

 

[ag4680@nyu_edu]

Hi,Just a quick update. For now I have used a combination of default 30 level input files for ne16np4, some fortran coding and the ncks command to create a working initial condition file. The model seems to work just fine. Thanks.Regards,Aman

 

[islas]

Hi Aman,Glad it's working fine for you now.  In case it is of use to you or others, I have included the option of three spectral element resolutions (ne16, ne30 and ne120) in the NCL scripts used to generate initial condition files on the simpler models website here....http://www.cesm.ucar.edu/models/simpler-models/changeinitialconditions.htmlThis makes the isothermal, zero u and v initial conditions of Held and Suarez but could be modified for other initial conditions. Best wishes,Isla