（P2）PALM软件的基本使用流程

上一篇介绍了PALM在Ubuntu系统上的下载安装及安装后的目录文件结构，这一篇以自带的案例介绍软件的基本操作流程。（请转至example_cbl项目）

输入

声明：官网均有各个变量的说明，以下代码来自软件自带案例
palm的最小输入文件是p3d文件，以下是example_cbl_p3d的详细说明。（文件语言格式Fortran，请自行先学习下语言的基本语法。笔者也是新学者，有非常多的问题，很多参数也不能准确的说出是什么意思，还请不吝赐教，谢谢！！！）

！几点说明
	！1、该运行文件基本由三大块组成：1）initialization_parameters，2）、runtime_parameters，3）调用的模块（自带案例没有调用模块，所以此处并没有涉及，我会在后面自定义项目中介绍）
	！2、一个块由“&”开始，以“/”结束
	！3、！是注释符号
!-------------------------------------------------------------------------------
!-- INITIALIZATION PARAMETER NAMELIST
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/inipar
!-------------------------------------------------------------------------------
&initialization_parameters   
!
!-- grid parameters   ！模型大小设置
!-------------------------------------------------------------------------------
    nx                         = 39, ! Number of gridboxes in x-direction (nx+1)
    ny                         = 39, ! Number of gridboxes in y-direction (ny+1)
    nz                         = 40, ! Number of gridboxes in z-direction (nz)     ！模型的xyz网格数量，注意x和y需要设置为实际网格数量减去1（说明：网格数设置需要是计算机核心数的倍数）

    dx                         = 50.0, ! Size of single gridbox in x-direction
    dy                         = 50.0, ! Size of single gridbox in y-direction
    dz                         = 50.0, ! Size of single gridbox in z-direction       ！各个方向上的分辨率

    dz_stretch_level           = 1225.0, ! Height (in m) where stretching starts
    dz_stretch_factor          = 1.08,   ! dz(k+1) = dz(k) * dz_stretch_factor     ！z方向上的网格拉伸
!
!-- initialization     ！模型初始化设置
!-------------------------------------------------------------------------------
    initializing_actions       = 'set_constant_profiles', ! initial conditions     ！初始化条件，这里使用的是恒定条件

    ug_surface                 = 0.0, ! u-comp of geostrophic wind at surface      
    vg_surface                 = 0.0, ! v-comp of geostrophic wind at surface    !地转风的分量（有没有谁能告诉我地转风怎么计算?）

    pt_surface                 = 300.0, ! initial surface potential temp   !初始化的潜在温度

    pt_vertical_gradient       =    0.0,
                                    1.0, ! piecewise temp gradients
    pt_vertical_gradient_level =    0.0,
                                  800.0, ! height level of temp gradients    ！温度梯度


!
!-- boundary conditions    ！边界条件
!-------------------------------------------------------------------------------
    surface_heatflux           = 0.1, ! sensible heat flux at the bottom surface   ！底边的显热通量
    bc_pt_b                    = 'neumann', ! required with surface_heatflux    ！
!
!-- numerics
!-------------------------------------------------------------------------------
    fft_method                 = 'temperton-algorithm',   ! build-in fft method

/ ! end of initialization parameter namelist

!-------------------------------------------------------------------------------
!-- RUNTIME PARAMETER NAMELIST
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/d3par
!-------------------------------------------------------------------------------
&runtime_parameters      !运行时参数
!
!-- run steering   
!-------------------------------------------------------------------------------
    end_time                   = 1800.0, ! simulation time of the 3D model    ！运行时间

    create_disturbances        = .TRUE., ! randomly perturbate horiz. velocity  ！扰乱
    dt_disturb                 = 150.0,  ! interval for random perturbations   
    disturbance_energy_limit   = 0.01,   ! upper limit for perturbation energy

    data_output_2d_on_each_pe  = .FALSE., ! don't do 2D output on each MPI rank

!
!-- data output    ！数据输出控制
！说明
！数据输出按时间分类：瞬时输出，时间平均输出
！数据输出按纬度分类：3d输出，截面输出，一维序列输出
!-------------------------------------------------------------------------------
    netcdf_data_format         = 2, ! use NetCDF3  

    dt_run_control             = 0.0,    ! output interval for run control  ！运行控制文件输出间隔
    dt_data_output             = 900.0,  ! output interval for general data  ！3d数据输出时间间隔
    dt_data_output_av          = 1800.0, ! output interval for averaged data ！时间平均输出数据的输出间隔
    dt_dopr                    = 900.0,  ! output interval for profile data

    data_output                = 'w_xy',
                                 'w_xz',
                                 'w_xz_av',
                                 'theta_xy',
                                 'theta_xz',      ！输出的数据，_xy的意思是输出xy截面，_av的意思是输出时间平均数据

    data_output_pr             = 'wtheta', 'w"theta"', 'w*theta*',
                                 'w*2',
                                 '#theta', 'theta*2',         

    section_xy                 = 2,
                                 10, ! grid index for 2D XY cross sections      ！xy截面输出的高度选择
    section_xz                 = 20, ! grid index for 2D XZ cross sections

    averaging_interval         = 900.0, ! averaging interval general data 
    dt_averaging_input         = 10.0,  ! averaging general data sampling rate

    averaging_interval_pr      = 600.0, ! averaging interval profile data
    dt_averaging_input_pr      = 10.0,  ! averaging profile data sampling rate

/ ! end of runtime parameter namelist

运行

该运行文件不用做修改，主要输出的是温度信息。运行的方法如下：
在current_version下打开终端，输入运行命令：palmrun -r example_cbl -c “default” -a “d3#” -X “4”
其中r后是项目名称，c是配置文件，x核心数量
输入y确定，开始运行，输出

输出文件介绍

输出文件保存在项目文件夹下的OUTPUT文件夹下，我的结果如下（我的后缀是002，是因为之前我已经运行过两次了）：

对应p3d文件的设置会有以上输出，p3d设置参考上述data_output部分

ts文件一直会输出
data_output                = 'w_xy', 'w_xz',  'w_xz_av', 'theta_xy', 'theta_xz',      
该设定会输出xy截面文件，xz截面文件，av_xz文件
data_output_pr             = 'wtheta', 'w"theta"', 'w*theta*','w*2', '#theta', 'theta*2',  
该设定会输出pr文件       

以example_cbl_xy.002.nc为例
1）我们打开它应该会看到两个变量w_xy和theta_xy
2）设定的输出时间间隔是900s，那应该输出900,1800两个时间点的数据
3）设定的输出高度是2和10（此处数字意思是纵向网格的索引，我们设置的网格分辨率是50，所以2代表第二个网格，也就是50-100m，所以其输出的截面是75m处的数据，同理10代表475m处）
用ncview工具快速浏览example_cbl_xy.002.nc，如下（背景是海贼王剧场版打金狮子的那一部黑衣出场的一幕）

至此，基本的使用流程就介绍完毕，之后会介绍如何创建自己的项目，如何引用更多的模块，引入用户自定义的地形地表，建筑物等等