Data pre-processing device and method for domestic ocean mode earth coupling system

Abstract

The invention discloses a data pre-processing device and method for a domestic ocean mode earth coupling system, and aims at data pre-processing between an LICOM3 ocean mode and a CESM earth coupling system. The data preprocessing device comprises an ocean mode input file automatic generation module, a grid coupling infrastructure file automatic generation module and an automatic detection module, and the ocean mode input file automatic generation module generates an ocean mode input file based on an original data file after receiving a first user instruction. The grid coupling infrastructure file automatic generation module receives a second user instruction and then generates a grid coupling infrastructure file based on the original data file, the automatic detection module detects whether the original data file, the ocean mode input file and the grid coupling infrastructure file are abnormal or not, and if yes, a termination instruction is output.

Description

Technical Field

[0001] The present invention belongs to the technical field of numerical simulation of the Earth system and data preprocessing. Specifically, it particularly relates to a data preprocessing device and method for the domestic ocean numerical model Earth coupling system CESM+LICOM3. Background Art

[0002] The Earth coupling system solution CESM+LICOM3 uses the LASG / IAP climate system ocean model version 3 (LICOM3) to replace the original ocean component POP2 in the Community Earth System Model (CESM), so as to enhance the ability of LICOM3 for high-resolution simulation and model tuning.

[0003] LICOM is an ocean circulation model independently developed by the ocean model team of the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics (LASG), Institute of Atmospheric Physics (IAP), Chinese Academy of Sciences (CAS), and is widely used in ocean dynamics research and ocean environment prediction. Since the LICOM3 model adopts a three-level grid structure, including horizontal grids and vertical stratifications, and can simulate the evolution of physical variables such as temperature, salinity, and flow fields in the ocean, its input data such as grids, terrain fields, initial fields, and forcing fields all need to be remade specifically for the LICOM3 model itself.

[0004] As the ocean model in the CESM coupling system, in order to perform data exchange with models such as the atmosphere and runoff, it is also necessary to produce mapping files (storing the mapping weights from the source model grid to the target model grid) and domain files (describing the geometric attributes and boundary conditions of the grids) between grids of different model resolutions, so as to achieve a full-coupling experiment of atmosphere-land-runoff-ocean-sea ice, and simulate the whole process from atmospheric precipitation to land runoff, then to river entry into the sea, and the energy and material exchange between the ocean and sea ice, for studying complex problems such as global climate change, water resource management, and ecosystem changes.

[0005] Currently, in the simulation application of the Earth coupling system solution CESM+LICOM3, data preprocessing operations generally require the following steps: (1) Generate the LICOM3 model horizontal grid file: Researchers use a custom Fortran script program to generate the original files of the horizontal grid and vertical grid of the orthogonal curvilinear grid (i.e., the three-level grid) introduced for LICOM3, and convert the horizontal grid file into a remap file in the SCRIP format required for making mapping.

[0006] (2)Generate the initial temperature and salinity fields and LICOM model forcing field files: According to the requirements of different resolution grids of the LICOM model, researchers manually interpolate through the CDO tool to generate the initial temperature and salinity fields of the LICOM three-level grid, including horizontal and vertical stratifications. In addition, the forcing field files required for the LICOM model operation (mainly including internal wave energy, chlorophyll concentration distribution, and reference sea surface salinity) need to be manually produced through the NCO and CDO tools respectively.

[0007] (3)Produce the mapping and domain files: When the correspondence between the atmospheric grid and the LICOM ocean grid is clear, researchers need to manually write SCRIP program scripts to generate the mapping and domain files required for the CESM system based on the atmospheric grid and the LICOM three-level grids with different resolutions.

[0008] (4)Detect the correctness of the files: Currently, researchers do not have an effective method to directly detect whether the generated mapping file is valid. Generally, as long as the mapping file can be normally generated through the SCRIP script, it is considered directly usable. Often, it is only when an abnormality occurs during the actual simulation operation of the CESM system that it is determined that there may be a problem with the generated mapping file or the original grid file. At this time, the original grid needs to be remade and the mapping file needs to be generated according to the abnormality prompt. The same is true for the initial temperature and salinity field files. After production, effective detection cannot be carried out. It is only when abnormalities such as NaN values appear in the integration results during the actual simulation operation of the CESM system that it is realized that there may be an error in the initial field file and it needs to be remade.

[0009] Existing data preprocessing tools usually lack special optimization for specific domestic ocean numerical models (such as LICOM3). When generating the initial temperature and salinity fields and forcing field files, a large amount of manual intervention is often required, which is not only inefficient but also prone to introducing human errors. Currently, researchers related to the LICOM3 model mainly use self-defined fortran script programs to generate the LICOM3 model grid and terrain field files, use the CDO and NCO tools to generate the initial fields and the forcing field files required for the LICOM3 model, and use the SCRIP tool to generate the mapping and domain files between the LICOM3 ocean model grid and the atmospheric model grid required for the CESM+LICOM3 coupling system. However, the combined use of these tools is rather cumbersome and complex, resulting in low efficiency and error-proneness in building the CESM+LICOM3 coupling system. The main disadvantages are summarized as follows: (1)Low efficiency and poor repeatability: Manually writing scripts to generate files is time-consuming and laborious, especially for complex three-level grid structures, which requires a large amount of manual intervention; and it is prone to introducing human errors, resulting in poor repeatability of the results and affecting the reliability of scientific research.

[0010] (2) Lack of automatic detection: There is a lack of an automatic detection function for the generated files in the prior art, and it is impossible to detect abnormalities in the files in a timely manner, which easily leads to the failure of the coupled system operation.

[0011] (3) Poor user experience: Users need to have a high programming ability and an in-depth understanding of the mode to complete data preprocessing. This increases the usage threshold and limits the wide application of the LICOM and CESM coupled systems.

[0012] Therefore, it is urgent to develop a data preprocessing device and method for a domestic ocean model earth coupling system that overcomes the above defects. Summary of the Invention

[0013] In view of the above problems, the present invention provides a data preprocessing device for a domestic ocean model earth coupling system. Among them, for the data preprocessing between the LICOM3 ocean model and the CESM earth coupling system, the data preprocessing device includes: An automatic generation module for ocean model input files, which generates ocean model input files based on the original data files after receiving the first user instruction; An automatic generation module for grid coupling infrastructure files, which generates grid coupling infrastructure files based on the original data files after receiving the second user instruction; An automatic detection module, which detects whether there are abnormalities in the original data files, the ocean model input files, and the grid coupling infrastructure files. If there are abnormalities, it outputs a termination instruction.

[0014] The above-mentioned data preprocessing device, wherein the automatic generation module for ocean model input files includes: A first reading and parsing unit, which reads and parses the input parameters of the LICOM3 horizontal grid file set by the user, and then fills the land points of the LICOM3 horizontal grid file to form a new horizontal grid file; An interpolation unit, which performs interpolation processing on the original data file to obtain a plurality of first temporary files, and respectively performs interpolation processing on the temperature temporary file and the salinity temporary file among the plurality of first temporary files to obtain two second temporary files; A merging unit, which merges the two second temporary files to obtain a third temporary file; A transposing unit, which performs north-south pole transposition on the third temporary file and the remaining first temporary files to obtain the ocean model input files, and the ocean model input files include a temperature-salinity initial field file and a forcing field file; A deletion unit, which deletes all the remaining temporary files and then outputs the temperature-salinity initial field file and the forcing field file.

[0015] The above data preprocessing device, wherein the interpolation unit performs horizontal bilinear interpolation on the original data file and remaps it onto the horizontal grid specified by the new horizontal grid file to obtain a plurality of first temporary files; The interpolation unit performs vertical interpolation on the temperature temporary file and the salinity temporary file respectively to convert them to the specified vertical levels, and correspondingly obtains the second temporary files.

[0016] The above data preprocessing device, wherein the grid coupling infrastructure file includes a weight mapping file and a domain file, and the grid coupling infrastructure file automatic generation module includes: A second reading and parsing unit that reads and parses the CESM horizontal grid file input parameters set by the user to obtain the atmospheric horizontal grid file, the ocean horizontal grid file, the atmospheric horizontal grid name, and the ocean horizontal grid name; A weight mapping file creation unit that constructs the weight mapping file based on the atmospheric horizontal grid file, the ocean horizontal grid file, the atmospheric horizontal grid name, and the ocean horizontal grid name through the ESMF_RegridWeightGen tool; A domain file creation unit that generates the domain file with the corresponding horizontal resolution according to the weight mapping file through the gen_domain tool.

[0017] The above data preprocessing device, wherein the weight mapping file creation unit includes: a) Define the longitude and latitude coordinates, cell area, and mask in the input source grid and target grid; b) Calculate the spatial intersection of the source grid and the target grid according to the input interpolation method, traverse all the source grids overlapping with each target, and calculate the weight matrix; c) Store the calculated weight matrix in a sparse format; d) Generate the weight mapping file containing the coordinates, area, mask of the source / target grid, the index of the sparse matrix, and the weight value based on the weight matrix stored in the sparse format.

[0018] The above data preprocessing device, wherein the domain file creation unit includes: a) Parse the weight mapping file and the grid information, and extract the weight matrix and the coordinates, cell boundaries, masks, and cell areas of the ocean and atmospheric grids; b) Copy the ocean grid mask, read or calculate the ocean cell area from the weight mapping file, and record the longitude and latitude coordinates of the ocean cell boundary to obtain an independent domain file of the pure ocean grid; c) Map the ocean mask from the ocean grid to the atmospheric grid through conserved interpolation weights, and mark the areas in the atmospheric grid that need to exchange data with the ocean model, so as to obtain the mask and regional information file of the ocean model on the atmospheric grid for the coupler in the CESM3 system to determine the ocean-atmosphere exchange area; d) Read the non-exchange area from the weight mapping file and adjust the cell area of the non-exchange area to obtain the mask and regional information file of the land model on the atmospheric grid.

[0019] The above data preprocessing device, wherein the automatic detection module includes: The first detection unit detects the correctness of the ocean and atmospheric horizontal grid files input to the ocean model input file automatic generation module and the grid coupling infrastructure file automatic generation module, including confirming whether the coordinate variables exist or are out of bounds and have reasonable dimensions, and verifying whether they are in a legal SCRIP or ESMF format; The second detection unit performs anomaly detection on each step of the ocean model input file automatic generation module and performs data integrity detection on the ocean model input file, including confirming whether the data file variables, units, and time ranges meet expectations, whether they are damaged or metadata is missing, and whether the variable dimensions are consistent; The third detection unit performs anomaly detection on each step of the grid coupling infrastructure file automatic generation module and performs inspection and verification on the grid coupling infrastructure file, including weight matrix inspection, grid consistency analysis and consistency verification. When an anomaly is detected, an anomaly prompt and the output termination instruction are output, and the automated execution program is terminated.

[0020] The above data preprocessing device, wherein the second detection unit is built into the ocean model input file automatic generation module; and / or the third detection unit is built into the grid coupling infrastructure file automatic generation module.

[0021] The above data preprocessing device, wherein the weight mapping file includes: The mapping file for the atmospheric-to-ocean conserved interpolation method, the mapping file for the atmospheric-to-ocean bilinear interpolation method, the mapping file for the atmospheric-to-ocean polar correction interpolation method, the mapping file for the ocean-to-atmosphere conserved interpolation method, the mapping file for the ocean-to-atmosphere bilinear interpolation method; Wherein, the domain file creation unit constructs the independent domain file of the pure ocean grid, the mask and regional information file of the ocean model on the atmospheric grid, and the mask and regional information file of the land model on the atmospheric grid based on the mapping file for the ocean-to-atmosphere conserved interpolation method.

[0022] The present invention also provides a data pre - processing method for a domestic ocean - atmosphere coupled system. For the data pre - processing between the LICOM3 ocean model and the CESM earth - atmosphere coupled system, the data pre - processing method includes: Automatic generation step of ocean model input files: After receiving a first user instruction, generate ocean model input files based on the original data files; Automatic generation step of grid coupling infrastructure files: After receiving a second user instruction, generate grid coupling infrastructure files based on the original data files; Automatic detection step: Detect whether there are abnormalities in the original data files, the ocean model input files, and the grid coupling infrastructure files. If there are abnormalities, output a termination instruction.

[0023] In summary, the advantages of the present invention compared with the prior art are as follows: Aiming at the data pre - processing requirements between the LICOM3 (LASG IAP Ocean Model) ocean model and the CESM (Community Earth System Model) earth - atmosphere coupled system, the present invention provides an automated and intelligent data pre - processing solution. It can automatically generate the temperature - salinity initial field and forcing field files required by the LICOM3 ocean model, as well as the weight mapping file and domain file required by the CESM system, and has an automatic detection function, which can perform quality checks on the generated files to ensure the integrity and consistency of the data, aiming to improve the automation degree of data pre - processing and enhance the operation efficiency of the coupled system and the user experience.

[0024] Other features and advantages of the present invention will be described in the subsequent description, and some of them will become obvious from the description or be understood by implementing the present invention. The objectives and other advantages of the present invention can be achieved and obtained through the structures pointed out in the description, claims, and drawings. Brief Description of the Drawings

[0025] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or the prior art. Obviously, the following drawings are some embodiments of the present invention. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0026] Figure 1 It is a schematic structural diagram of the data pre - processing device of the present invention; Figure 2 It is Figure 1 The working flowchart of the automatic generation module of the ocean model input files in Figure 3 is Figure 1 The working flowchart of the automatic generation module for the medium grid coupling infrastructure file; Figure 4 is Figure 1 The functional block diagram of the automatic detection module in; Figure 5 The flowchart of the data preprocessing method of the present invention; Figure 6 is the code structure diagram of the automatic generation module for the ocean model input file. Specific embodiments

[0027] To make the objectives, technical solutions and advantages of the embodiments of the present invention clearer, the following will clearly and completely describe the technical solutions in the embodiments of the present invention with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are some, but not all, of the embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0028] The schematic embodiments of the present invention and their descriptions are used to explain the present invention, but not to limit the present invention. Additionally, elements / components using the same or similar reference numerals in the drawings and embodiments are used to represent the same or similar parts.

[0029] Regarding the "first", "second", "S1", "S2",... used herein, they do not particularly refer to the order or sequence, nor are they used to limit the present invention. They are only used to distinguish elements or operations described with the same technical terms.

[0030] Regarding the directional terms used herein, such as: up, down, left, right, front or back, etc., they are only references to the directions in the drawings. Therefore, the directional terms used are for explanation and not for limiting this creation.

[0031] Regarding the "including", "comprising", "having", "containing", etc. used herein, they are all open-ended terms, that is, they are meant to include but not be limited to.

[0032] Regarding the "and / or" used herein, it includes any one or all combinations of the described things.

[0033] Regarding "multiple" herein, it includes "two" and "more than two"; regarding "multiple groups" herein, it includes "two groups" and "more than two groups".

[0034] Regarding terms such as "substantially" and "about" used herein, they are used to modify any quantity or error that can vary slightly, but these slight variations or errors do not change their essence. Generally speaking, the range of such slight variations or errors modified by such terms can be 20% in some embodiments, 10% in some embodiments, 5% in some embodiments, or other values. Those skilled in the art should understand that the aforementioned values can be adjusted according to actual needs and are not limited thereto.

[0035] Certain terms used to describe this application will be discussed below or elsewhere in this specification to provide additional guidance to those skilled in the art regarding the description of this application.

[0036] Please refer to Figure 1 , Figure 1 which is a schematic structural diagram of the data preprocessing device of the present invention. As Figure 1 shown, a data preprocessing device for a domestic ocean model earth coupling system of the present invention, wherein, for data preprocessing between the LICOM3 ocean model and the CESM earth coupling system, the data preprocessing device includes: an ocean model input file automatic generation module 11, a grid coupling infrastructure file automatic generation module 12, and an automatic detection module 13; the ocean model input file automatic generation module 11 generates an ocean model input file based on the original data file after receiving a first user instruction; the grid coupling infrastructure file automatic generation module 12 generates a grid coupling infrastructure file based on the original data file after receiving a second user instruction, and the grid coupling infrastructure file includes a weight mapping file and a domain file; the automatic detection module 13 detects whether there are abnormalities in the original data file, the ocean model input file, and the grid coupling infrastructure file, and outputs a termination instruction if there are abnormalities.

[0037] The data preprocessing device of the present invention mainly includes three modules: an automatic generation module 11 for ocean model input files, an automatic generation module 12 for grid coupling infrastructure files, and an automatic detection module 13. Among them, the automatic generation module 11 for ocean model input files generates an ocean model input file for the LICOM3 model with one key according to the user's requirements for grids of different resolutions by designing a set of automated scripts. The ocean model input file includes temperature and salinity initial field and forcing field files; the automatic generation module 12 for grid coupling infrastructure files generates the grid coupling infrastructure files required for the CESM coupling system with one key according to the horizontal grid files of different resolutions between different models specified by the user. The grid coupling infrastructure files include weight mapping files and domain files. In this embodiment, the automatic generation module 12 for grid coupling infrastructure files generates the weight mapping files (mapping files) corresponding to different interpolation methods (such as conserve|bilinear|patch) required for the CESM coupling system with one key, and then automatically generates three domain files (domain files) according to the weight mapping file (mapping file) in the ocn-to-atm direction generated just now; the function of the automatic detection module is actually integrated inside the previous two modules, including performing correctness detection on the input grid files, performing data integrity detection on the generated temperature and salinity initial field files, performing validity detection on the generated weight mapping files (mapping files), and performing exception detection and handling on each step in the automatic generation module process.

[0038] Please refer to Figure 2 and Figure 6 , Figure 2 is Figure 1 the workflow diagram of the automatic generation module for ocean model input files in Figure 6 the code structure diagram of the automatic generation module for ocean model input files. As Figure 2 and Figure 6 shown and please combine with Figure 1 , the automatic generation module 11 for ocean model input files includes: A first reading and parsing unit 111, which reads and parses the input parameters of the LICOM3 horizontal grid file set by the user, and fills the land points of the LICOM3 horizontal grid file to form a new horizontal grid file; An interpolation unit 112, which performs interpolation processing on the original data file to obtain multiple first temporary files, and performs interpolation processing on the temperature temporary file and the salinity temporary file in the multiple first temporary files respectively to obtain two second temporary files; among them, the interpolation unit 112 performs horizontal bilinear interpolation on the original data file and remaps it to the horizontal grid specified by the new horizontal grid file to obtain multiple first temporary files; The interpolation unit performs vertical interpolation on the temperature temporary file and the salinity temporary file respectively to convert them to specified vertical levels, and obtains the corresponding second temporary file. The merging unit 113 merges the two second temporary files to obtain a third temporary file; The transposing unit 114 performs north-south pole transposition on the third temporary file and the remaining first temporary files to obtain the ocean model input file, and the ocean model input file includes a temperature and salinity initial field file and a forcing field file; The deleting unit 115 deletes all the remaining temporary files and then outputs the temperature and salinity initial field file and the forcing field file.

[0039] Specifically, the initial field is the initial description of the ocean state when the model starts, mainly including: the temperature and salinity field and the flow velocity field. And the initial flow velocity field is usually set to a stationary state (zero velocity field), or output based on the stable state after the model spin-up operation. Therefore, the present invention mainly deals with the temperature and salinity initial field for the LICOM3 model.

[0040] Currently, the initial temperature and salinity field file initialTS_{gridname}.nc used in the LICOM model is obtained by interpolating and converting the format of the original temperature and salinity datasets of PHC3 or CMEMS GLORYS for different horizontal resolution grids and vertical layer distributions of LICOM, and is used for the initialization of the temperature and salinity variables in the LICOM model. Among them, PHC3 (Polar Science Center Hydrographic Climatology 3.0) is a global ocean temperature and salinity climatology dataset developed by the Polar Science Center of the University of Washington in the United States. Its coverage is the global ocean, including polar and high-latitude regions, with a horizontal resolution of 1°×1°, 33 layers in the vertical direction (from the sea surface to the seabed), and including monthly and annual average climatology data of temperature and salinity. GLORYS is an important global ocean reanalysis and forecasting product under CMEMS (Copernicus Marine Environment Monitoring Service). Its coverage is also the global ocean area, providing a relatively high horizontal resolution (for example, the common version is 1 / 12°, about 8-9 km), 50 layers in the vertical direction, and including rich ocean physical variables such as three-dimensional temperature (Temperature), three-dimensional salinity (Salinity), three-dimensional zonal velocity (Zonal Velocity, U), three-dimensional meridional velocity (Meridional Velocity, V), sea surface height (Sea Surface Height, SSH), and mixed layer depth (Mixed Layer Depth, MLD). The time span is a long time series from the early 1990s (usually 1993) to the present, and will continue to be updated to provide near real-time data.

[0041] The forcing fields generally refer to the external input data that drive the operation of ocean models, mainly including wind stress, heat flux, freshwater flux, etc. These data generally come from atmospheric reanalysis datasets or the output of the atmospheric component of a coupled model. Since the present invention mainly focuses on the preprocessing device for the Earth Coupled System Project CESM+LICOM3, the conventional forcing field data of the LICOM model are output by the atmospheric component or other model components such as sea ice in the CESM system. The forcing field files referred to in the present invention mainly refer to the input data required for the internal parameterization scheme of the LICOM3 model, including the chlorophyll concentration distribution data file chl_monthly_{gridname}.nc required in the shortwave radiation penetration scheme dependent on chlorophyll concentration (Ohlmann, 2003), the tidal-related energy data file tidal_energy_{gridname}.nc required in the internal tide mixing parameterization scheme, and the reference sea surface salinity data file sss_phc3_monthly_{gridname}.nc required for sea surface salinity restoration forcing.

[0042] The automatic generation module 11 of the ocean model input files can generate four files, namely the initial temperature and salinity field initialTS_{gridname}.nc, the chlorophyll concentration distribution data file chl_monthly_{gridname}.nc, the tidal-related energy data file tidal_energy_{gridname}.nc, and the reference sea surface salinity data file sss_phc3_monthly_{gridname}.nc, with one click according to the instruction bash gen_cesm_maps.sh -scripgrid_file remap_licom_tripole_*.nc -res{gridname} input by the user through an automated script. In this automated script, the SCRIP-format LICOM3 horizontal grid remap_licom_tripole_{grid_size}.nc file is required as input, and the global ocean raw data that has been downloaded offline and converted in format needs to be prepared in advance. For example, the original temperature data file temp_jan_phc3.nc is obtained from the binary-format single-time-point PHC3 temperature dataset temp.ctl, the original salinity data file salt_jan_phc3.nc is obtained from the binary-format single-time-point PHC3 salinity dataset salt.ctl, the original sea surface salinity data file sss_month_phc3.nc is obtained from the binary-format 12-month PHC3 salinity dataset salt_month_phc3.ctl, the original chlorophyll concentration distribution data file smooth_chl_clim_monthly_licomgrid.nc is obtained by bilinear interpolation of the chlorophyll concentration distribution data at 12-month time points downloaded from the NASA GSFC (Goddard Space Flight Center) website, and the tidal dissipation data converted from the tidal_mixing_energy_{create_date}.nc estimated by the barotropic tidal model using the internal wave energy conversion parameterization method.

[0043] Code structure of the automated script for the automatic generation module 11 of the ocean model input files Figure 6As shown, the root directory includes the automation script file generate_licom_foring.sh, the raw data file sub-directory raw_data (including raw salinity data files, raw temperature data files, raw sea surface salinity data files, raw chlorophyll concentration distribution data files, raw tidal energy-related data files, etc.), and the SCRIP format horizontal grid file sub-directory scripgrids (including global horizontal grid files with resolutions of 100km, 10km, 5km, 3km, 2km, and 1km).

[0044] Among them, the processing flow of the automation script generate_licom_foring.sh of the ocean model input file automatic generation module 11 is as Figure 2 shown, and mainly includes the following steps: 1. The first reading and parsing unit 111 reads and parses the LICOM3 horizontal grid file input parameters set by the user, that is, the grid file path - scripgrid_file scripgrids / remap_licom_tripole_*.nc and the grid file name - res {gridname}, provided that this file actually exists in the corresponding path; 2. The first detection unit 131 of the automatic detection module 13 detects whether this global horizontal grid file is valid. If it is valid, land point filling is performed on this horizontal grid file to form a new fill_grid.nc file. Otherwise, an exception prompt A is given and the automation execution program is exited; 3. The interpolation unit 112 remaps all raw data files (including temperature, salinity, sea surface salinity, chlorophyll concentration, tidal energy) to the horizontal grid specified by fill_grid.nc using bilinear interpolation respectively, and saves them as new high-efficiency storage format (such as NetCDF4) temporary files temp1_temperature.nc, temp1_salinity.nc, temp1_sss.nc, temp1_chl.nc, temp1_tidal.nc in sequence; 4. The second detection unit 132 of the automatic detection module 13 judges each operation in the third step respectively. If the execution is successful, the interpolation unit 112 will perform vertical layer interpolation on the temperature and salinity temporary file data temp1_temperature.nc and temp1_salinity.nc after horizontal interpolation in the third step, convert them to the specified vertical layer (such as 80 layers), and save them as new temporary files temp2_temperature.nc and temp2_salinity.nc in the new efficient storage format in sequence. At the same time, delete temp1_temperature.nc and temp1_salinity.nc. If the execution fails, an exception prompt B will be given and the automated execution program will exit. 5. The automatic detection module 13 judges each operation in the fourth step respectively. If the execution is successful, the merging unit 113 will merge the temporary files temp2_temperature.nc and temp2_salinity.nc after horizontal and vertical interpolation in the fourth step into a temporary temperature-salinity initial field file temp_initialTS_{gridname}.nc. At the same time, delete temp2_temperature.nc and temp2_salinity.nc. If the execution fails, an exception prompt C will be given and the automated execution program will exit. 6. The second detection unit 132 of the automatic detection module 13 judges the operation in the fifth step. If the execution is successful, the transposing unit 114 will transpose the temporary files temp1_sss.nc, temp1_chl.nc, temp1_tidal.nc in the third step and the temporary file temp_initialTS_{gridname}.nc in the fifth step to the north and south poles in sequence. If the execution fails, an exception prompt D will be given and the automated execution program will exit. 7. The second detection unit 132 of the automatic detection module 13 judges each operation in the sixth step respectively. If the execution is successful, 4 files required by the LICOM model will be generated, namely the temperature-salinity initial field initialTS_{gridname}.nc, the chlorophyll concentration distribution data file chl_monthly_{gridname}.nc, the tidal-related energy data file tidal_energy_{gridname}.nc, and the reference sea surface salinity data file sss_phc3_monthly_{gridname}.nc. The deletion unit 115 deletes all the remaining temporary files temp*.nc. If the execution fails, an exception prompt E will be given and the automated execution program will exit.

[0045] Please refer to Figure 3 ,Figure 3 For Figure 1 the working flow chart of the intermediate grid coupling infrastructure file automatic generation module. As Figure 3 shown and please combine with Figure 1 , the intermediate grid coupling infrastructure file automatic generation module 12 includes: A second reading and parsing unit 121, which reads and parses the input parameters of the CESM horizontal grid file set by the user to obtain the atmospheric horizontal grid file, the ocean horizontal grid file, the atmospheric horizontal grid name, and the ocean horizontal grid name; A weight mapping file creation unit 122, which constructs the weight mapping file based on the atmospheric horizontal grid file, the ocean horizontal grid file, the atmospheric horizontal grid name, and the ocean horizontal grid name through the ESMF_RegridWeightGen tool; A domain file creation unit 123, which generates the domain file with the corresponding horizontal resolution according to the weight mapping file through the gen_domain tool.

[0046] Among them, the weight mapping file creation unit 122 includes: a) Define the longitude and latitude coordinates, cell area, and mask in the input source grid and target grid; b) Calculate the spatial intersection of the source grid and the target grid according to the input interpolation method, traverse all the source grids overlapping with each target, and calculate the weight matrix; c) Store the calculated weight matrix in a sparse format; d) Generate the weight mapping file including the coordinates, area, mask of the source / target grid, the index of the sparse matrix, and the weight value based on the weight matrix stored in the sparse format.

[0047] Among them, the domain file creation unit 123 includes: a) Parse the weight mapping file and grid information, and extract the weight matrix and the coordinates, cell boundaries, masks, and cell areas of the ocean and atmospheric grids; b) Copy the ocean grid mask, read or calculate the ocean cell area from the weight mapping file, and record the longitude and latitude coordinates of the ocean cell boundary to obtain the independent domain file of the pure ocean grid; c) Map the ocean mask from the ocean grid to the atmospheric grid through the conserved interpolation weight, and mark the areas in the atmospheric grid that need to exchange data with the ocean model to obtain the mask and area information file of the ocean model on the atmospheric grid for the coupler in the CESM3 system to determine the ocean-atmosphere exchange area; d) Read the non-exchange area from the weight mapping file and adjust the non-exchange area cell area to obtain the mask and area information file of the land model on the atmospheric grid.

[0048] Specifically, in an Earth system coupling model (such as CESM), the mapping weight mapping file and the domain file are important files used to define grid attributes and data interpolation relationships, and they play a crucial role when coupling different subsystems such as the atmosphere and the ocean. The mapping file is used to define the data mapping relationship between different grids, providing the interpolation weights from one grid to another to ensure the accuracy of data during the conversion process. In the CESM coupling system, the atmospheric and ocean models usually use different grid resolutions and structures, so it is necessary to interpolate the output data of one model onto the grid of another model through the mapping file to achieve seamless data transfer. The domain file is used to describe the geometric attributes and boundary conditions of the grid, providing detailed information about the grid for the coupling system, including the horizontal and vertical resolutions of the grid, that is, the size and hierarchical structure of the grid, specifying the boundary range of the grid, such as the land boundary of the ocean grid, and the grid weights used to calculate area averages or other statistics.

[0049] In CESM, different types of physical quantities (such as state variables and fluxes) have different physical properties, so different interpolation methods need to be adopted to ensure the accuracy of the simulation and physical conservation. State variables (such as temperature, humidity, etc.) usually represent the system state, and spatial continuity and accuracy are emphasized during interpolation. Common interpolation methods include bilinear interpolation or nearest neighbor interpolation to minimize errors.

[0050] Fluxes (such as heat flux, momentum flux, water flux, etc.) represent the flow of energy or matter, and conservation needs to be emphasized during interpolation to avoid introducing artificial sources / sinks. Therefore, conservative interpolation is generally used to ensure that the total flux does not change during the grid conversion process. Specific types of vector variables (such as wind speed, flow velocity, etc.) usually represent physical fields with direction and magnitude, such as the wind field in the atmosphere or the water flow in the ocean. Direction and magnitude accuracy need to be emphasized during interpolation to avoid introducing artificial errors or direction distortions. Therefore, the commonly used interpolation method is patch interpolation, which takes into account the topological structure of the grid to ensure that the physical properties of the vector (such as direction conservation and magnitude continuity) are correctly preserved during the grid conversion process.

[0051] The grid coupling infrastructure file automatic generation module 12 in the present invention can, according to the user input The generate_mapping_domain.sh -fatm {atm_gridname}*.nc -focn {ocn_gridname}*.nc -nocn {ocn_gridname} -natm {atm_gridname} command can generate 5 mapping files with one click by extracting the horizontal grid files and names of the atmosphere and ocean from the input parameters. These files include the mapping file for the conservative interpolation method from the atmosphere to the ocean (map_{atm_gridname}_TO_{ocn_gridname}_aave.nc), the bilinear interpolation mapping file from the atmosphere to the ocean (map_{atm_gridname}_TO_{ocn_gridname}_blin.nc), the polar correction interpolation mapping file from the atmosphere to the ocean (map_{atm_gridname}_TO_{ocn_gridname}_patc.nc), the mapping file for the conservative interpolation method from the ocean to the atmosphere (map_{ocn_gridname}_TO_{atm_gridname}_aave.nc), the bilinear interpolation mapping file from the ocean to the atmosphere (map_{ocn_gridname}_TO_{atm_gridname}_blin.nc), and 3 domain files, including the mask and regional information file of the land model on the atmospheric grid domain.lnd.{atm_gridname}_{ocn_gridname}.nc, the mask and regional information file of the ocean model on the atmospheric grid domain.ocn.{atm_gridname}_{ocn_gridname}.nc, and the independent domain file of the pure ocean grid domain.ocn.{ocn_gridname}.nc.

[0052] After compiling and adapting the underlying ESMF library (mainly including the ESMF_RegridWeightGen tool) and the internal gen_domain program through the tools / mapping program in this module, an automated script generate_mapping_domain.sh is added on the outermost layer for one-click processing. The specific processing flow is as Figure 3 shown, mainly including the following steps: 1. The second reading and parsing unit 121 reads and parses the parameters input by the user, namely the atmospheric horizontal grid file - fatm {atm_gridname}*.nc, the ocean horizontal grid file - focn {ocn_gridname}*.nc, the atmospheric horizontal grid name - nocn {ocn_gridname}, and the ocean horizontal grid name - natm {atm_gridname}, provided that the two grid files actually exist in the corresponding paths; 2. The first detection unit 131 of the automatic detection module 13 detects the correctness of the atmospheric horizontal grid file {atm_gridname}*.nc and the ocean horizontal grid file {ocn_gridname}*.nc, including confirming whether their coordinate variables (such as grid_corner_lon, grid_corner_lat, grid_dims) exist, are out of bounds, and have reasonable dimensions, and verifying whether they are in a legal SCRIP or ESMF format. If it is correct, continue to execute downward; otherwise, give an exception prompt A and exit the automated execution program; 3. The weight mapping file creation unit 122 uses the ESMF_RegridWeightGen tool to create a set of mapping weight files for the CESM system, namely 5 mapping files including map_{atm_gridname}_TO_{ocn_gridname}_aave.nc, map_{atm_gridname}_TO_{ocn_gridname}_blin.nc, map_{atm_gridname}_TO_{ocn_gridname}_patc.nc, map_{ocn_gridname}_TO_{atm_gridname}_aave.nc, and map_{ocn_gridname}_TO_{atm_gridname}_blin.nc. The specific implementation steps of the ESMF_RegridWeightGen tool include: a) defining the longitude and latitude coordinates, cell area, and mask in the input source grid and target grid; b) calculating the spatial intersection of the source grid cells and the target grid cells according to the input interpolation method. For each target cell, traverse all the source cells overlapping with it and calculate the weight value (such as the area ratio of conservative interpolation); c) storing the calculated weight matrix in a sparse format, recording the row (target cell), column (source cell), and weight value; d) generating a mapping weight mapping file in NetCDF format containing the coordinates, area, mask of the source / target grid, the index and weight value of the sparse matrix, and other metadata (interpolation method, version, etc.); 4. The third detection unit 133 of the automatic detection module 13 checks and validates the above-generated 5 mapping files respectively, including weight matrix check, i.e., verifying the weight value range, sparsity, and non-zero element distribution; grid consistency analysis, i.e., checking the mask mask, coverage range, and area statistics of the source grid and the target grid; conservation and consistency verification, i.e., ensuring that the total error of the conserved interpolation (such as area weighting) meets the expectation; and metadata integrity check, i.e., confirming whether the auxiliary information such as grid coordinates and area is complete. If the detection is normal, continue to execute downward; otherwise, give an exception prompt B and exit the automatic execution program. 5. The third detection unit 133 of the domain file creation unit 123 generates 3 domain files corresponding to the horizontal resolution through the gen_domain program according to the ocean-to-atmosphere conserved interpolation mapping file, i.e., map_{ocn_gridname}_TO_{atm_gridname}_aave.nc, which are the mask and regional information file domain.lnd.{atm_gridname}_{ocn_gridname}.nc of the land model on the atmospheric grid, the mask and regional information file domain.ocn.{atm_gridname}_{ocn_gridname}.nc of the ocean model on the atmospheric grid, and the independent domain file domain.ocn.{ocn_gridname}.nc of the pure ocean grid. The specific implementation steps of the gen_domain program include: a) Parsing the input mapping file and grid information, extracting the weight matrix and the coordinates, cell boundaries, mask (such as ocean / land marking), and cell area of the ocean and atmospheric grids; b) Copying the ocean grid mask, reading or calculating (such as spherical polygon area) the cell area from the ocean grid file and recording the longitude and latitude coordinates of the cell boundaries, and saving it as the ocean grid domain file domain.ocn.{ocn_gridname}.nc; c) Mapping the ocean mask from the ocean grid to the atmospheric grid through the conserved interpolation weight, marking which areas in the atmospheric grid need to exchange data with the ocean model for the coupler in the CESM system to determine the ocean-atmosphere exchange area, and saving the result as the file domain.ocn.{atm_gridname}_{ocn_gridname}.nc; d) Reading the original land mask from the atmospheric grid file, marking non-exchange areas such as glaciers and lakes, and adjusting the cell area to ensure that the land mask does not overlap with the ocean area in domain.ocn.{atm_gridname}_{ocn_gridname}.nc, and saving it as the file domain.lnd.{atm_gridname}_{ocn_gridname}.nc for the land model to use. 6. The third detection unit 133 of the automatic detection module 13 judges the operation in the fifth step. If the execution is successful, it will end normally; if the execution fails, it will give an exception prompt C and exit the automated execution program.

[0053] It should be noted that in this embodiment, the weight mapping file includes: the atmosphere-to-ocean conservation interpolation method mapping file, the atmosphere-to-ocean bilinear interpolation method mapping file, the atmosphere-to-ocean polar correction interpolation method mapping file, the ocean-to-atmosphere conservation interpolation method mapping file, the ocean-to-atmosphere bilinear interpolation method mapping file; wherein, the domain file creation unit constructs the independent domain definition file of the pure ocean grid, the mask and regional information file of the ocean model on the atmospheric grid, and the mask and regional information file of the land model on the atmospheric grid based on the ocean-to-atmosphere conservation interpolation method mapping file.

[0054] Based on the automatic generation module, the present invention significantly improves the data preprocessing efficiency and reduces manual intervention; at the same time, users do not need to manually write scripts, which reduces the requirements for users' programming ability and improves the user experience; moreover, it reduces human errors, improves the repeatability of results, and contributes to the reliability of scientific research.

[0055] Please refer to Figure 4 , Figure 4 For Figure 1 the functional block diagram of the automatic detection module in Figure 4 shown, and in combination with Figure 1 , the automatic detection module 13 includes: The first detection unit 131 detects the correctness of the ocean and atmospheric horizontal grid files input into the ocean model input file automatic generation module and the grid coupling infrastructure file automatic generation module, including confirming whether the coordinate variables exist or are out of bounds and the dimensions are reasonable, and verifying whether it is in a legal SCRIP or ESMF format; The second detection unit 132 performs anomaly detection on each step of the ocean model input file automatic generation module and data integrity detection on the ocean model input file, including confirming whether the data file variables, units, and time ranges meet expectations, whether they are damaged or metadata is missing, and whether the variable dimensions are consistent; The third detection unit 133 performs anomaly detection on each step of the grid coupling infrastructure file automatic generation module and checks and verifies the grid coupling infrastructure file, including weight matrix check, grid consistency analysis and consistency verification. When an anomaly is detected, it outputs an anomaly prompt and the output termination instruction, and terminates the automated execution program.

[0056] Among them, in this embodiment, the second detection unit 132 is built into the ocean model input file automatic generation module; and / or; the third detection unit 133 is built into the grid coupling infrastructure file automatic generation module.

[0057] Specifically, the first detection unit 131 detects the correctness of the ocean and atmosphere horizontal grid files input to the ocean model input file automatic generation module 11 and the grid coupling infrastructure file automatic generation module 12, including confirming whether its coordinate variables (such as grid_corner_lon, grid_corner_lat, grid_dims) exist or are out of bounds and the dimensions are reasonable, and verifying whether it is in a legal SCRIP or ESMF format to ensure the effective execution of subsequent modules; The second detection unit 132 built into the ocean model input file automatic generation module 11 performs anomaly detection for each step to ensure that it can successfully generate the files required for the LICOM model to run. In addition, data integrity detection needs to be performed on the generated temperature and salinity initial field files, including confirming whether the data file variables, units, and time ranges meet expectations, whether they are damaged or metadata is missing, and whether the variable dimensions are consistent; Built into the mapping and domain file automatic generation module, it mainly checks and verifies the previously generated mapping files, including weight matrix check, that is, verifying the weight value range, sparsity, and non-zero element distribution, grid consistency analysis, that is, checking the mask masks, coverage, and area statistics of the source grid and the target grid, conservation and consistency verification, that is, ensuring that the total error of the conserved interpolation (such as area weighting) meets expectations, and metadata integrity check, that is, confirming whether the auxiliary information such as grid coordinates and areas is complete. For the detected abnormal mapping, an anomaly prompt is given, and the automated execution program is terminated to ensure the generation of correct mapping files before conducting scientific experiments in the Earth Coupled System CESM.

[0058] The present invention can timely detect file anomalies through the automatic detection module, reduce the error rate, and improve the operation success rate of the coupling system.

[0059] Please refer to Figure 5 , Figure 5 For the flowchart of the data preprocessing method of the present invention. As Figure 5 shown, a data preprocessing method for a domestic ocean model Earth coupling system of the present invention, wherein, for the data preprocessing between the LICOM3 ocean model and the CESM Earth coupling system, the data preprocessing method includes: Ocean model input file automatic generation step S11: Generate an ocean model input file based on the original data file after receiving the first user instruction; Grid coupling infrastructure file automatic generation step S12: After receiving the second user instruction, generate a grid coupling infrastructure file based on the original data file; Automatic detection step S13: Detect whether there are abnormalities in the original data file, the ocean model input file, and the grid coupling infrastructure file. If there are abnormalities, output a termination instruction.

[0060] In summary, the present invention aims to solve the problems of inefficiency, lack of automated processing, and poor user experience in the existing data preprocessing technology of the current Earth coupling system solution CESM+LICOM3. Its beneficial effects are as follows: (1) Through a carefully designed script, one-key automatic generation of the temperature-salinity initial field (including horizontal and vertical stratification) and forcing field files required for the three-level grid of the LICOM3 ocean model itself is achieved. This improvement simplifies the data preprocessing process of the LICOM ocean model, greatly reduces manual intervention, avoids human errors, and improves the efficiency and accuracy of data generation.

[0061] (2) On the premise that the user specifies the atmospheric grid and the LICOM3 ocean grid, relevant map and domain files can be automatically generated simultaneously through the script for the simulation operation of the CESM system. This makes the grid matching work efficient and accurate, meeting the personalized needs of different users.

[0062] (3) An automatic detection function is innovatively introduced. Each step of the above files can be automatically detected. Once an abnormality is found, a warning prompt is immediately given. This function not only saves the time and effort of manual error checking, but also improves the reliability and stability of the data, thereby enhancing the simulation operation efficiency of the CESM coupling system.

[0063] In short, through the above improvements, the present invention effectively solves the problems of low efficiency, poor accuracy, insufficient adaptability, and lack of automatic detection in the existing technology, providing a more efficient, reliable, and convenient data preprocessing tool for domestic ocean numerical simulation research.

[0064] Although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: they can still modify the technical solutions recorded in the foregoing embodiments, or perform equivalent replacements for some of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A data preprocessing device for a domestic ocean model earth coupling system, characterized in that For the data preprocessing between the LICOM3 ocean model and the CESM earth coupling system, the data preprocessing device includes: An ocean model input file automatic generation module, which generates an ocean model input file based on the original data file after receiving a first user instruction; A grid coupling infrastructure file automatic generation module, which generates a grid coupling infrastructure file based on the original data file after receiving a second user instruction; An automatic detection module, which detects whether there are abnormalities in the original data file, the ocean model input file, and the grid coupling infrastructure file. If there are abnormalities, it outputs a termination instruction.

2. The data preprocessing device according to claim 1, wherein The ocean model input file automatic generation module includes: A first reading and parsing unit, which reads and parses the LICOM3 horizontal grid file input parameters set by the user, and then fills the land points in the LICOM3 horizontal grid file to form a new horizontal grid file; An interpolation unit, which performs interpolation processing on the original data file to obtain multiple first temporary files, and respectively performs interpolation processing on the temperature temporary file and the salinity temporary file in the multiple first temporary files to obtain two second temporary files; A merging unit, which merges the two second temporary files to obtain a third temporary file; A transposing unit, which transposes the third temporary file and the remaining first temporary files between the north and south poles to obtain the ocean model input file. The ocean model input file includes a temperature-salinity initial field file and a forcing field file; A deletion unit, which deletes all the remaining temporary files and then outputs the temperature-salinity initial field file and the forcing field file.

3. The data preprocessing device according to claim 2, wherein The interpolation unit performs horizontal bilinear interpolation on the original data file and remaps it to the horizontal grid specified by the new horizontal grid file to obtain multiple first temporary files; The interpolation unit respectively performs vertical interpolation on the temperature temporary file and the salinity temporary file to convert them to the specified vertical levels to obtain the corresponding second temporary files.

4. The data preprocessing device according to claim 2, wherein The grid coupling infrastructure file includes a weight mapping file and a domain file. The grid coupling infrastructure file automatic generation module includes: A second reading and parsing unit, which reads and parses the CESM horizontal grid file input parameters set by the user to obtain the atmospheric horizontal grid file, the ocean horizontal grid file, the atmospheric horizontal grid name, and the ocean horizontal grid name; A weight mapping file creation unit, which constructs the weight mapping file based on the atmospheric horizontal grid file, the ocean horizontal grid file, the atmospheric horizontal grid name, and the ocean horizontal grid name through the ESMF_RegridWeightGen tool; A domain file creation unit, which generates the domain file with the corresponding horizontal resolution according to the weight mapping file through the gen_domain tool.

5. The data preprocessing device according to claim 4, wherein The weight mapping file creation unit includes: a) Define the longitude and latitude coordinates, cell area, and mask in the input source grid and target grid; b) Calculate the spatial intersection of the source grid and the target grid according to the input interpolation method, traverse all the source grids overlapping with each target, and calculate the weight matrix; c) Store the calculated weight matrix in a sparse format; d) Generating the weight mapping file that includes the coordinates, areas, masks of the source / destination grids, the indices of the sparse matrix, and the weight values based on the weight matrix stored in a sparse format.

6. The data preprocessing device according to claim 5, characterized in that, The domain file creation unit includes: a) Parsing the weight mapping file and the grid information, and extracting the weight matrix, and the coordinates, cell boundaries, masks, and cell areas of the ocean and atmosphere grids; b) Copying the ocean grid mask, reading or calculating the ocean cell areas from the weight mapping file, and recording the longitude and latitude coordinates of the ocean cell boundaries to obtain an independent domain file for the pure ocean grid; c) Mapping the ocean mask from the ocean grid to the atmosphere grid through conservative interpolation weights, and marking the areas in the atmosphere grid that need to exchange data with the ocean model, so as to obtain the mask and regional information file of the ocean model on the atmosphere grid for the coupler in the CESM3 system to determine the ocean-atmosphere exchange area; d) Reading the non-exchange area from the weight mapping file and adjusting the cell areas of the non-exchange area to obtain the mask and regional information file of the land model on the atmosphere grid.

7. The data preprocessing device according to claim 1, wherein The automatic detection module includes: The first detection unit, which detects the correctness of the ocean and atmosphere horizontal grid files input to the ocean model input file automatic generation module and the grid coupling infrastructure file automatic generation module, including confirming whether the coordinate variables exist or are out of bounds and have reasonable dimensions, and verifying whether they are in a legal SCRIP or ESMF format; The second detection unit, which performs anomaly detection on each step of the ocean model input file automatic generation module and performs data integrity detection on the ocean model input file, including confirming whether the data file variables, units, and time ranges meet the expectations, whether they are damaged or metadata is missing, and whether the variable dimensions are consistent; The third detection unit, which performs anomaly detection on each step of the grid coupling infrastructure file automatic generation module and performs inspection and verification on the grid coupling infrastructure file, including weight matrix inspection, grid consistency analysis and consistency verification, and outputs an anomaly prompt and the output termination instruction when detecting an anomaly, and terminates the automated execution program.

8. The data preprocessing device according to claim 7, wherein The second detection unit is built into the ocean model input file automatic generation module; and / or the third detection unit is built into the grid coupling infrastructure file automatic generation module.

9. The data preprocessing device according to claim 6, wherein The weight mapping file includes: The mapping file for the atmosphere-to-ocean conservative interpolation method, the mapping file for the atmosphere-to-ocean bilinear interpolation method, the mapping file for the atmosphere-to-ocean polar correction interpolation method, the mapping file for the ocean-to-atmosphere conservative interpolation method, the mapping file for the ocean-to-atmosphere bilinear interpolation method; Among them, the domain file creation unit constructs the independent domain file for the pure ocean grid, the mask and regional information file of the ocean model on the atmosphere grid, and the mask and regional information file of the land model on the atmosphere grid based on the mapping file for the ocean-to-atmosphere conservative interpolation method.

10. A data preprocessing method for a domestic ocean model earth coupling system, characterized in that, For the data preprocessing between the LICOM3 ocean model and the CESM earth coupling system, the data preprocessing method includes: Steps for automatically generating an ocean model input file: Generate an ocean model input file based on the original data file after receiving a first user instruction; Steps for automatically generating a grid coupling infrastructure file: Generate a grid coupling infrastructure file based on the original data file after receiving a second user instruction; Automatic detection step: Detect whether there are any abnormalities in the original data file, the ocean model input file, and the grid coupling infrastructure file, and output a termination instruction if there are any abnormalities.

Images