Clavrx使用经验（非本人编写）

MissLH88

老师突然要求我朋友用这个模式处理静止卫星资料产出云分类产品，都是本科生，身边没有人问，我朋友研究了半个月才研究明白这个模式的使用方法。网上资料太少，我就让他写下来，一是万一我们自己以后也要使用方便查看，二是和更多可能要使用这个模式的人分享一些浅薄的经验。

Clavrx可用于识别云类，参考：

说明：http://cimss.ssec.wisc.edu/clavr/clavrx_docs/

Pdf文档：http://cimss.ssec.wisc.edu/clavr/documentation.html

参考文献：

https://journals.ametsoc.org/view/journals/apme/51/6/jamc-d-11-02.1.xml

https://journals.ametsoc.org/view/journals/apme/44/6/jam2236.1.xml

最近时间的计算结果可在此处获取（但此处无葵花结果）：

安装步骤

首先需要确保netcdf库与hdf5库正确安装，在编译hdf5库时需要在./configure后面加上：

--enable-fortran --enable-cxx --enable-fortran2003

在编译hdf4时需要加上

-DHAVE_NETCDF

否则hdf4可能与netcdf相互干扰。

下载源码：

svn checkout -q https://svn.ssec.wisc.edu/repos/ ... 02_10_2020_r3712/./

本研究中需使用葵花8卫星数据，因此还需要安装libHimawari,地址:

wget https://gitlab.ssec.wisc.edu/rayg/himawari/-/archive/master/himawari-master.tar.gz

安装libHimawari，解压文件：

tar -zxvf himawari-master.tar.gz

将其中src/makefile中所有编译器换为编译clavrx所用的编译器。

需修改第39，58，59，60，61行。

                              

更改保存后直接make即可。

编译clavrx：

此处以ifort编译器为例，如使用gfortran编译器即将参数中_ifort改为_gfortran，但需注意clavrx需与hdf4，hdf5，netcdf，libHimawari使用同一编译器。

./configure -hdf4_ifort=#填入hdf4安装路径，即bin文件夹所在位置

                -hdf5_ifort=#填入hdf5安装路径，即bin文件夹所在位置

             -netcdf_ifort=#填入netcdf安装路径，即bin文件夹所在位置

                -LIBHIM_IFORT=#填入libHimawari安装路径，即src文件夹所在位置

清除已安装的残余（其中ifort需根据所用编译器修改）：

make ARCH=ifort-debug distclean

安装：make ARCH=ifort-debugall_plus，注意第一次安装必须加上all_plus以编译全部文件。

资料获取：

下载静态需求文件地址：

https://bin.ssec.wisc.edu/pub/CSPP/hidden/CLAVRx/v3.0/CSPP_CLAVRX_V3.0_STATIC.tar.xz

下载luts资料地址：

ftp://ftp.ssec.wisc.edu/clavr/cld_ref_luts/

下载动态需求文件地址：

https://geodb.ssec.wisc.edu/ancillary/2021_07_13_194/

https://jpssdb.ssec.wisc.edu/cspp_v_2_0/ancillary/

卫星资料提前备好，按原下载格式（.nc，.bz2）即可

Goes卫星各观测像素点视角文件（葵花8该文件无已有结果，可根据该文件格式自行计算）：

ftp://ftp.ssec.wisc.edu/clavr/awalther/g16/goes_east_abi_conus_static_nav.nc

输入文件：

clavrx_options中expert mode设为0即可

file_list中第一行为原文件路径，第二行为结果存放路径，第三行开始为原数据文件名

level2_list中存放需要输出的结果。可参考手册中Level-2 File Contents部分。

注意问题：

默认Acha_Mode_User_Set为无效值，需要在user_options.f90文件第225行更改为有效值，参考pdf手册中acha_mode选项。

Ahi_hsd_reader中Him_Rad_Subsample函数Status变量未赋初值，需在该文件334行加上Status = Sym%SUCCESS

注意level2_mod.f90第3469及3508行if (Sds_Info(Var_Idx)%Scaling_Type /= 0) then需改为

                            if(trim(Sds_Info(Var_Idx)%Sds_Name)/="latitude".and.trim(Sds_Info(Var_Idx)%Sds_Name)/="longitude")then

2926行最后一个减号改为加号

一大碗年糕

顶