关于cesm自定义模拟文件的一些说明记录

记录时间2020/12/29

1. 关于 datm

aerosoldep_WACCM.ensmean_monthly_hist_1849-2018_0.9x1.25_CMIP6_c201223.nc

说明：
（1）NCAR 文件路径参考：CLMinputdata\atm\cam\chem\trop_mozart_aero\aero
（2）aerosoldep_WACCM.ensmean_monthly_hist_1849-2018_0.9x1.25_CMIP6_c201223.nc 文件是根据NCAR官方数据在时间上延展而得，
    原始时间为1849-2015年。后将2015年的数据复制到2016，2017，2018此3年中，得到1849-2018共170年数据；

2. 关于强迫场的 domain data

domain.lnd.360x720_ecmwf.c201217.nc

说明：
（1）不仅地表数据有 domain 文件，强迫场也有 domain 文件，且其分辨率需要和强迫场的分辨率一致，才能保证顺利插值。
（2）domain.lnd.360x720_ecmwf.c201217.nc 是根据ECMWF强迫场（0.5°x0.5°）有数据的值，同时用NCAR默认的hcru_hcru (0.5°x0.5°) 的地表的domain文件（domain.lnd.360x720_cruncep.100429.nc）进行补漏制作得到，目的是保证陆地表面上空都有强迫场数据。
（3）domain 中的变量 mask，1表示激活，0表示未激活。
（4）NCAR默认的强迫场一般都是全球范围内都有值，即所有网格的 mask 均为“1”； 但是ECMWF只是针对陆地，而海洋上空时没有强迫场的，因此海洋上空的 mask 为 “0”；
（5）由于NCAR的地表数据全球范围都有值，通过 mask 和 frac 来控制陆地和海洋的网格。mask中陆地为1，海洋为0。因此，如只需计算陆地，需要保证地表 mask 为1的地方上空强迫场必须有数值，即强迫场的domian文件中mask也为1。
（6）如此，只要保证强迫场 domain 为1的网格能完全覆盖 地表数据 mask为1的 网格即可（可以大于等于，不能小于）。

3. 关于地表的 domain data

说明：
（1）NCAR默认的（0.5°x0.5°）的domain文件（domain.lnd.360x720_cruncep.100429.nc）中，mask中陆地为1，海洋为0。
（2）domain.lnd.360x720_cruncep.100429.nc 中frac表示陆地面积占据网格面积的比例；
（3）由于地表数据surface data在每个网格都有数值，即不存在miss_value，因此需要通过地表的 domain 中的mask来确定是否激活该网格。

• 使用场景

情景1：使用自己的强迫场数据进行模拟时，需要制作 stream.txt* 文件，每个stream.txt* 文件中有对应的domaininInfo信息，如下。其中 domain 文件就是atm对应的domain

<?xml version="1.0"?>
<file id="stream" version="1.0">
<dataSource>
  GENERIC
</dataSource>
<domainInfo>
 <variableNames>
    time    time
       xc      lon
       yc      lat
       area    area
       mask    mask
 </variableNames>
 <filePath>
    /home/dayon/inputdata/share/domains/domain.clm
 </filePath>
 <fileNames>
    domain.lnd.360x720_ecmwf.c201217.nc
 </fileNames>
</domainInfo>
<fieldInfo>
  <variableNames>
    FLDS     lwdn
  </variableNames>
  <filePath>
    /home/dayon/inputdata/Ecmwf/LWdown
  </filePath>
  <fileNames>
   ECMWF_c2020_LWdown.2011-01.nc
ECMWF_c2020_LWdown.2011-02.nc
ECMWF_c2020_LWdown.2011-03.nc
ECMWF_c2020_LWdown.2011-04.nc
ECMWF_c2020_LWdown.2011-05.nc
ECMWF_c2020_LWdown.2011-06.nc
ECMWF_c2020_LWdown.2011-07.nc
ECMWF_c2020_LWdown.2011-08.nc
ECMWF_c2020_LWdown.2011-09.nc
ECMWF_c2020_LWdown.2011-10.nc
ECMWF_c2020_LWdown.2011-11.nc
ECMWF_c2020_LWdown.2011-12.nc
  </fileNames>
  <offset>
     0
  </offset>
</fieldInfo>
</file>

情景2：使用自己制作的surface data时，在 ./case.setup 之前，需要设置一下强迫场和地表的domain文件

./xmlchange ATM_DOMAIN_FILE = domain.lnd.360x720_ecmwf.c201217.nc     # 文件需要放到默认路径下
./xmlchange LND_DOMAIN_FILE = domain.lnd.360x720_cruncep.100429.nc     # 文件需要放到默认路径下

4. 关于地表的 surface data

说明：
（1）NCAR 默认的地表数据，每个变量是默认是不含miss_value，因此每个网格都有数值；
（2）surface data 文件中不包含mask，其输出结果的网格又地表对应的domain file 文件中的 mask 决定；

Variables:
    mxsoil_color           
           Size:       1x1
           Dimensions: 
           Datatype:   int32
           Attributes:
                       long_name = 'maximum numbers of soil colors'
                       units     = 'unitless'
    SOIL_COLOR             
           Size:       720x360
           Dimensions: lsmlon,lsmlat
           Datatype:   int32
           Attributes:
                       long_name = 'soil color'
                       units     = 'unitless'
    PCT_SAND               
           Size:       720x360x10
           Dimensions: lsmlon,lsmlat,nlevsoi
           Datatype:   double
           Attributes:
                       long_name = 'percent sand'
                       units     = 'unitless'
    PCT_CLAY               
           Size:       720x360x10
           Dimensions: lsmlon,lsmlat,nlevsoi
           Datatype:   double
           Attributes:
                       long_name = 'percent clay'
                       units     = 'unitless'
    ORGANIC                
           Size:       720x360x10
           Dimensions: lsmlon,lsmlat,nlevsoi
           Datatype:   double
           Attributes:
                       long_name = 'organic matter density at soil levels'
                       units     = 'kg/m3 (assumed carbon content 0.58 gC per gOM)'
    FMAX                   
           Size:       720x360
           Dimensions: lsmlon,lsmlat
           Datatype:   double
           Attributes:
                       long_name = 'maximum fractional saturated area'
                       units     = 'unitless'
    natpft                 
           Size:       15x1
           Dimensions: natpft
           Datatype:   int32
           Attributes:
                       long_name = 'indices of natural PFTs'
                       units     = 'index'
    cft                    
           Size:       2x1
           Dimensions: cft
           Datatype:   int32
           Attributes:
                       long_name = 'indices of CFTs'
                       units     = 'index'

5. 知识点

• 对于陆面模块（CLM），为保证其顺利运行，强迫场的domain中mask为“1”的范围一定要囊括住 surface 的domain中mask为“1”的范围；
  
• 模型计算输出的结果是以 surface 对应的 domain 中 mask=1的位置决定的；mask=1 的地方，网格才有结果；
• 当强迫场的空间分辨率和地表的空间分辨率一致时，根据前面的讲解，设置强迫场的domain文件时，也可使用surface 的domain文件替代。
• 当强迫场的空间分辨率和地表的空间分辨率不同时，需要强迫场在进行插值的时候，需要其插值的可能性，从而合理的制作domain文件及mask范围。