Construction method of earth system coupling device and earth system coupling device

Abstract

The present invention discloses a construction method of an Earth system coupling device and an Earth system coupling device. The construction method includes: a first component module construction step: constructing a plurality of first component modules of a fully coupled global climate system, the plurality of first component modules including an atmospheric component module, a land component module, a sea ice component module, a land ice component module, a river component module, a first ocean component module, and a wave component module; a second component module construction step: constructing a second component module based on an ocean circulation simulation system, the second component module being a second ocean component module; a replacement step: replacing the first ocean component module with the second ocean component module; a coupling step: reconstructing a coupling application interface of a coupler of the fully coupled global climate system, and performing data exchange and interaction between different component modules through the coupling application interface.

Description

Technical Field

[0001] The present invention relates to the technical field of earth system science and technology. Specifically, it particularly relates to a construction method of an earth system coupling device and an earth system coupling device. Background Art

[0002] LICOM is an ocean circulation model independently developed by the Chinese Academy of Sciences (CAS). After more than 30 years of continuous research and development, it has now evolved to the LICOM3 (LASG / IAP Climate System Ocean Model Version 3 (LICOM3)) version. The model includes multiple sets of resolutions and parameterization schemes, with the highest resolution reaching 3 - 5 km globally. As an ocean sub - model for constructing an earth system model, LICOM can simulate past geological periods and present climate changes, and predict climate changes in the next 10 - 30 years or even the next century. The public earth system model CESM (Community Earth System Model) is a fully coupled global climate model that provides users with the most advanced computer simulations of the earth's past, present, and future climate states, and is an advanced technology for earth system simulation. The initial CESM1 version only included 5 geophysical sub - modules, namely the atmosphere (CAM), ocean (POP), land (CLM), sea ice (CICE), and ice sheet (CISM), and newly added a river transport model (RTM) and a wave numerical model WAVE WATCH III (WW3).

[0003] In addition to participating in CMIP6 as the ocean component of the climate system model, LICOM3 also participated in OMIP (Ocean Model Inter - comparison Project) with the earth system coupling device of CESM1 (version 1.0.4) + LICOM3. By providing a unified atmospheric external forcing dataset and flux calculation scheme, global ocean - sea ice coupling experiments were carried out. However, in actual use, it was found that the existing earth system coupling device of CESM1 (version 1.0.4) + LICOM3 has the following defects:

[0004] (1) The existing LICOM3 model coupling scheme is no longer applicable to the high - version CESM framework;

[0005] (2) There is a lack of an automated parameter configuration method for conducting coupling experiments related to the LICOM3 ocean model. The steps are cumbersome. According to grids of different resolutions, multiple configuration files and parameters need to be manually modified, which is error - prone and has poor inheritance;

[0006] (3) Previously, only ocean - sea ice coupling was applicable in the CESM1.0.4 + LICOM3 coupling device, and full - coupling experiments of atmosphere - land - runoff - ocean - sea ice could not be achieved.

[0007] (4)Previously, the operation results output by the LICOM3 mode itself were saved in binary form and needed to be post-processed and converted for visual analysis or post-processing, which was a complex process and consumed extra resources.

[0008] Therefore, there is an urgent need to develop a construction method for an Earth system coupling device and an Earth system coupling device that overcome the above defects. Summary of the Invention

[0009] In view of the above problems, the present invention provides a construction method for an Earth system coupling device, which includes:

[0010] Step of constructing the first component modules: Construct a plurality of first component modules based on the fully coupled global climate system, and the plurality of first component modules include an atmospheric component module, a land component module, a sea ice component module, a land ice component module, a river component module, a first ocean component module, and a sea wave component module;

[0011] Step of constructing the second component module: Construct a second component module based on the ocean circulation simulation system, and the second component module is a second ocean component module;

[0012] Replacement step: Replace the first ocean component module with the second ocean component module;

[0013] Coupling step: Reconstruct the coupling application interface of the coupler based on the fully coupled global climate system, and perform data exchange and interaction between different component modules through the coupling application interface.

[0014] In the above construction method, the step of constructing the second component module further includes:

[0015] Modify use emsf to use ESMF_ClockMod, only : ESMF_Clock to introduce the single variable ESMF_Clock in the form of ESMF_ClockMod, only : ESMF_Clock;

[0016] Add different machine definitions in the machine environment configuration file to correspond to the compilation control options for different platform compilers;

[0017] Add two machine configuration items with MACH being sunway and shuguang, and different -mach options can be set when executing create_newcase to select coupling experiments on different supercomputer platforms;

[0018] Add adaptive compilation options in the compilation configuration file, and the compilation options only take effect when compiling the ocean component module.

[0019] The above construction method, wherein, in the replacement step, the second ocean component module is preset into the global climate system in the same mode as the first ocean component module as an optional second ocean component module configuration, including: integrating the mode function code of the second ocean component module into the directory, and adding instructions for replacing the first ocean component module with the second ocean component module in the system configuration file as needed.

[0020] The above construction method, wherein the coupling application interface includes a top-level control interface, a middle-level function interface, and a component module application interface corresponding to each component module, and the component module application interface is reconstructed in the coupling step.

[0021] The above construction method, wherein the steps of reconstructing the component module application interface in the coupling step include:

[0022] Modifying the name of the application interface file of the component module application interface corresponding to the second ocean component module;

[0023] Reconstructing at least one function function of the component module application interface corresponding to the second ocean component module;

[0024] Adding a source code sub-directory of the second ocean component module to the ocean model component source code directory corresponding to the second ocean component module in the global climate system.

[0025] The above construction method, which further includes a construction automation configuration step: placing the namelist configuration file corresponding to the second ocean component module in the source code sub-directory.

[0026] The above construction automation configuration step includes:

[0027] Adding variables required for the operation of the source code of the second ocean component module to the namelist variable list compilation configuration file of the second ocean component module;

[0028] Adding relevant attribute definitions required for the operation of the source code of the second ocean component module to the namelist variable list definition configuration file of the second ocean component module:

[0029] Adding default values corresponding to the grid names with different resolutions of the second ocean component module to the variable list default value configuration file of the second ocean component module.

[0030] The above construction method, wherein the second component module construction step further includes:

[0031] Modify the B1850_LICOM component set in the component set configuration file;

[0032] Add a definition of res for the grid in the grid configuration file;

[0033] Add a default pes configuration for GRID in the parallel parameter configuration file;

[0034] Add a default namelist configuration for the grid inside the second ocean component module and the sea ice component module.

[0035] The above construction method, further comprising steps of constructing data saving and reading:

[0036] After enabling the PnetCDF mode support of PIO by numerical configuration when constructing the PIO library in the global climate system, write the ocean-related variable data obtained by integrating the second ocean component module using the PnetCDF mode into a preset file.

[0037] The present invention further provides an earth system coupling device, which is obtained by the construction method described in any one of the above, and the earth system coupling device includes:

[0038] Multiple first component modules of a fully coupled global climate system, and the multiple first component modules include an atmospheric component module, a land component module, a sea ice component module, a land ice component module, a river component module, and a wave component module;

[0039] A second component module of an ocean circulation simulation system, and the second component module is a second ocean component module;

[0040] A coupler of the fully coupled global climate system, and data exchange and interaction between different component modules are carried out through the coupler.

[0041] In summary, the effects of the present invention compared with the prior art are as follows:

[0042] (1) Coupled LICOM3 as a new ocean component module into the CESM1.3 high-resolution version system, and is compatible with multiple supercomputer platforms such as domestic Sunway supercomputers and X86 supercomputers;

[0043] (2) Design a brand-new coupling application interface for LICOM3, with strong compatibility and being convenient for subsequent updates and upgrades of the CESM version;

[0044] (3) Design an automated experimental scheme for the ocean component module of LICOM3, and the execution process of the experimental steps is simple and clear, greatly simplifying the operation of scientific experiments;

[0045] (4) Develop different - resolution grids for the ocean component module of LICOM3, and add a fully - coupled experimental configuration of atmosphere - land - runoff - ocean - sea ice based on the LICOM3 model.

[0046] (5) Design a new parallel data standard output scheme for the ocean component module of LICOM3, simplify the original post - processing scheme, and achieve high - speed and formatted output of the output data.

[0047] Other features and advantages of the present invention will be described in the following specification, and, in part, will be obvious from the specification, or will be understood by implementing the present invention. The objectives and other advantages of the present invention can be realized and obtained by the structures pointed out in the specification, claims, and drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0048] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or in 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 drawings in the following description are some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can be obtained based on these drawings.

[0049] Figure 1 is a flowchart of the construction method of the present invention;

[0050] Figure 2 is a schematic diagram of the initial calculation and continued calculation processes of the fully - coupled experiment case of the coupling device based on the present invention;

[0051] Figure 3 is a flowchart of the automated processing related to the CESM system namelist;

[0052] Figure 4 is a schematic diagram of the automated configuration structure of the LICOM3 model namelist;

[0053] Figure 5 is a schematic diagram of the comparison of the continued calculation execution processes before and after modifying the output scheme;

[0054] Figure 6 is a schematic diagram of the structure of the earth system coupling device of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0055] 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.

[0056] The schematic embodiments of the present invention and the descriptions thereof are used to explain the present invention, but do not limit the present invention. In addition, the same or similar reference numerals of elements / components used in the drawings and embodiments are used to represent the same or similar parts.

[0057] Regarding the "first", "second", "S1", "S2",... used in this article, 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.

[0058] Regarding the directional terms used in this article, 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 illustration and not for limiting this creation.

[0059] Regarding the "comprising", "including", "having", "containing", etc. used in this article, they are all open-ended terms, that is, they mean including but not limited to.

[0060] Regarding the "and / or" used in this article, it includes any one or all combinations of the described things.

[0061] Regarding the "plurality" in this article, it includes "two" and "more than two"; regarding the "multiple groups" in this article, it includes "two groups" and "more than two groups".

[0062] Regarding the terms "substantially", "about", etc. used in this article, they are used to modify any quantity or error that can vary slightly, but these slight variations or errors will not change their essence. Generally, 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.

[0063] 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.

[0064] Please refer to Figure 1 is the flowchart of the construction method of the present invention. As Figure 1As shown, a method for constructing an earth system coupling device of the present invention includes:

[0065] Step S1 of constructing the first component module: Construct a plurality of first component modules based on the fully coupled global climate system. The plurality of first component modules include an atmospheric component module, a land component module, a sea ice component module, a land ice component module, a river component module, a first ocean component module, and a wave component module;

[0066] Step S2 of constructing the second component module: Construct a second component module based on the ocean circulation simulation system. The second component module is the second ocean component module;

[0067] Replacement step S3: Replace the first ocean component module with the second ocean component module;

[0068] Coupling step S4: Reconstruct the coupling application interface of the coupler based on the fully coupled global climate system, and perform data exchange and interaction between different component modules through the coupling application interface.

[0069] Specifically, based on the modular framework of the CESM system, after constructing the main body composed of an atmospheric component module, a land component module, a sea ice component module, a land ice component module, a river component module, a first ocean component module, and a wave component module, the data information exchange and model operation between each module are managed through the constructed coupler (CPL7). Among them, in the present invention, the atmospheric component module adopts CAM (Community Atmosphere Model), the first ocean component module adopts POP2 (Parallel Ocean Program), the land component module adopts CLM (Community Land Model), the sea ice component module adopts CICE (Community Ice CodE), the land ice component module adopts CISM (CESM Ice Sheet Model), the river component module adopts RTM (River Transport Model), and the wave component module adopts WW3 (WaveWatch III). Among them, the first ocean component module is replaced by a second component module based on the ocean circulation simulation system;

[0070] Thus, the construction method of the present invention couples the domestic climate system ocean model version 3 (LICOM3) as a new ocean component module to the CESM1.3 high-resolution version system, and performs a series of automated adaptations to obtain an earth system coupling device based on CESM1.3 + LICOM3, so that users can perform coupling experiments related to the LICOM3 ocean model with the simplest case creation, compilation, and execution steps.

[0071] Among them, in the said alternative step, the second ocean component module is preset into the global climate system in the same mode as the first ocean component module as an optional second ocean component module configuration, including: integrating the mode function code of the second ocean component module into the directory, and at the same time adding an instruction to replace the first ocean component module with the second ocean component module in the system configuration file as needed. Specifically, in this step, the mode function code of the second ocean component module LICOM3 is integrated into the models / ocn / licom3 directory, and at the same time, a COMPSET that replaces the POP2 default ocean component with LICOM3 is added as needed in the system configuration file scripts / ccsm_utils / Case.template / config_compsets.xml, which is convenient for directly creating a coupled experiment with the ocean mode of LICOM3.

[0072] Furthermore, the second component module construction step also includes:

[0073] Modify use emsf to use ESMF_ClockMod, only : ESMF_Clock to introduce the single variable ESMF_Clock in the form of ESMF_ClockMod, only : ESMF_Clock;

[0074] Add different machine definitions in the machine environment configuration file to correspond to the compilation control options of different platform compilers;

[0075] Add two machine configuration items with MACH being sunway and shuguang, and different -mach options can be set when executing create_newcase to select to conduct coupled experiments on different supercomputer platforms;

[0076] Add adaptive compilation options in the compilation configuration file, and the compilation options only take effect when compiling the ocean component module.

[0077] Specifically, in response to some error messages encountered during the compilation and running of the second ocean component module on the new version of CESM1.3, some code of the second ocean component module was modified. For example, all places where use esmf is declared in the second ocean component module LICOM3 were deleted, and only the actual call places in the ocn_comp_mct.F90 file were retained. However, since only one variable ESMF_Clock was used, use emsf was changed to use ESMF_ClockMod, only:ESMF_Clock. According to the design principle of the CESM system, different compilation control options corresponding to compilers on different platforms can be added by defining different machines in the machine environment configuration file (scripts / ccsm_utils / Machines / config_machines.xml), enabling the CESM system to be adapted to run on different platforms. In this solution, two machine configuration items with MACH being sunway and shuguang were added for the domestic Sunway supercomputer and the X86 supercomputer respectively. During the execution of create_newcase, different -mach options can be set to select different supercomputer platforms for coupling experiments. Additionally, adaptive compilation options <ADD_FFLAGS MODEL="licom3">-fdefault-double-8 -fdefault-real-8< / ADD_FFLAGS> were added for the second ocean component module in the compilation configuration file (scripts / ccsm_utils / Machines / config_compilers.xml). This option only takes effect when compiling the component mode of the second ocean component module LICOM3, ensuring the normal operation of the second ocean component module LICOM3 without affecting the compilation of other modes.

[0078] In this embodiment, the step of constructing the second component module further includes:

[0079] Modify the B1850_LICOM component set in the component set configuration file (scripts / ccsm_utils / Case.template / config_compsets.xml);

[0080] Add a res grid definition in the grid configuration file (scripts / ccsm_utils / Case.template / config_grid.xml);

[0081] Add the pes default configuration for GRID in the parallel parameter configuration file (scripts / ccsm_utils / Machines / config_pes.xml);

[0082] Added default configuration for namelist of grids inside the second ocean component module and the sea ice component module.

[0083] For example, B1850 is a specific climate model experiment, which is a climate simulation experiment conducted using CESM. It provides a climate benchmark, which is very important for understanding the natural variability of the climate system and anthropogenic climate change. The ocean model used by the CESM system's own compset B1850 is POP2, that is, <compset sname="B_1850" alias="B1850">1850_CAM4_CLM40%SP_CICE_POP2_RTM_SGLC_SWAV

[0084] < / compset> ,The following takes compset as B1850_LICOM and res as ne30np4_lt1v1 as an example to explain the relevant process of fully coupled experiment configuration. The name of the 100km low-resolution grid in LICOM3 mode is defined as "lt1v1", and the name of the 10km high-resolution grid is defined as "lt0.1v1". The following files are mainly modified:

[0085] (1) Add the following line to the component set configuration file (scripts / ccsm_utils / Case.template / config_compsets.xml) for compset:

[0086] B1850_LICOM adds the following definition:

[0087] <compset sname="B_1850_LICOM" alias="B1850_LICOM">1850_CAM4_CLM40%SP_CICE_LICOM3_RTM_SGLC_SWAV< / compset>

[0088] (2) Add res to the ne30np4_lt1v1 grid related definition in the grid configuration file (scripts / ccsm_utils / Case.template / config_grid.xml):

[0089] <grid sname="ne30np4_lt1v1" alias="ne30_lt1v1">a%ne30np4_l%ne30np4_oi%lt1v1_r%r05_m%lt1v1_g%null_w%null< / grid>

[0090] <gridhorz name="lt1v1">

[0091] <nx> 360< / nx> <ny> 218< / ny>

[0092] <desc>lt1v1 is tripole 1-deg grid for LICOM low resolution:< / desc>

[0093] < / gridhorz>

[0094] <griddom grid="lt1v1" mask="lt1v1"><ICE_DOMAIN_FILE>domain.ocn.lt1v1.240809.nc< / ICE_DOMAIN_FILE><OCN_DOMAIN_FILE>domain.ocn.lt1v1.240809.nc< / OCN_DOMAIN_FILE>< / griddom>

[0095] <griddom grid="ne30np4" mask="lt1v1"><ATM_DOMAIN_FILE>domain.lnd.ne30np4_lt1v1.240809.nc< / ATM_DOMAIN_FILE><LND_DOMAIN_FILE>domain.lnd.ne30np4_lt1v1.240809.nc< / LND_DOMAIN_FILE><ICE_DOMAIN_FILE>domain.ocn.ne30np4_lt1v1.240809.nc< / ICE_DOMAIN_FILE><OCN_DOMAIN_FILE>domain.ocn.ne30np4_lt1v1.240809.nc< / OCN_DOMAIN_FILE>

[0096] < / griddom>

[0097] <gridmap atm_grid="ne30np4" ocn_grid="lt1v1">

[0098] <ATM2OCN_FMAPNAME>cpl / gridmaps / ne30np4 / map_ne30np4_to_lt1v1_aave.20240809.nc< / ATM2OCN_FMAPNAME>

[0099] <ATM2OCN_SMAPNAME>cpl / gridmaps / ne30np4 / map_ne30np4_to_lt1v1_blin.20240809.nc< / ATM2OCN_SMAPNAME>

[0100] <ATM2OCN_VMAPNAME>cpl / gridmaps / ne30np4 / map_ne30np4_to_lt1v1_patc.20240809.nc< / ATM2OCN_VMAPNAME>

[0101] <OCN2ATM_FMAPNAME>cpl / gridmaps / lt1v1 / map_lt1v1_to_ne30np4_aave.20240809.nc< / OCN2ATM_FMAPNAME>

[0102] <OCN2ATM_SMAPNAME>cpl / gridmaps / lt1v1 / map_lt1v1_to_ne30np4_aave.20240809.nc< / OCN2ATM_SMAPNAME>

[0103] < / gridmap>

[0104] <gridmap rof_grid="r05" ocn_grid="lt1v1"><ROF2OCN_FMAPNAME>cpl / gridmaps / r05 / map_r05_TO_lt1v1_aave.240809.nc< / ROF2OCN_FMAPNAME>

[0105] < / gridmap>

[0106] (3) Add the default PES configuration for GRID as

[0107] "a%ne30np4_l%ne30np4_oi%lt1v1" to the parallel parameter configuration file (scripts / ccsm_utils / Machines / config_pes.xml):

[0108] <pes grid="a%ne30np4_l%ne30np4_oi%lt1v1" mach="sunway">

[0109] <NTASKS_ATM>80< / NTASKS_ATM><NTHRDS_ATM>1< / NTHRDS_ATM><ROOTPE_ATM>0< / ROOTPE_ATM><NINST_ATM>1< / NINST_ATM>

[0110] <NTASKS_LND>$NTASKS_ATM< / NTASKS_LND><NTHRDS_LND>$NTHRDS_ATM< / NTHRDS_LND><ROOTPE_LND>$ROOTPE_ATM< / ROOTPE_LND><NINST_LND>1< / NINST_LND>

[0111] <NTASKS_ICE>$NTASKS_ATM< / NTASKS_ICE><NTHRDS_ICE>$NTHRDS_ATM< / NTHRDS_ICE><ROOTPE_ICE>$ROOTPE_ATM< / ROOTPE_ICE><NINST_ICE>1< / NINST_ICE>

[0112] <NTASKS_OCN>$NTASKS_ATM< / NTASKS_OCN><NTHRDS_OCN>$NTHRDS_ATM< / NTHRDS_OCN><ROOTPE_OCN>$ROOTPE_ATM< / ROOTPE_OCN><NINST_OCN>1< / NINST_OCN>

[0113] <NTASKS_CPL>$NTASKS_ATM< / NTASKS_CPL><NTHRDS_CPL>$NTHRDS_ATM< / NTHRDS_CPL><ROOTPE_CPL>$ROOTPE_ATM< / ROOTPE_CPL>

[0114] <NTASKS_GLC>$NTASKS_ATM< / NTASKS_GLC><NTHRDS_GLC>$NTHRDS_ATM< / NTHRDS_GLC><ROOTPE_GLC>$ROOTPE_ATM< / ROOTPE_GLC><NINST_GLC>1< / NINST_GLC><PIO_NUMTASKS>-1< / PIO_NUMTASKS><PIO_STRIDE>-1< / PIO_STRIDE><PIO_TYPENAME>netcdf< / PIO_TYPENAME><PIO_ROOT>1< / PIO_ROOT>

[0115] <PES_LEVEL>0< / PES_LEVEL>

[0116] < / pes> .

[0117] (4) At the same time, add the default namelist configuration for the grid "lt1v1" inside the second ocean component module LICOM3 and the sea ice component module CICE. For example, add the decomp configuration corresponding to nproc="80" res="lt1v1" in models / ocn / licom3 / bld / licom_decomp.xml and models / ice / cice / bld / cice_decomp.xml; add the sea ice horizontal and vertical grid files corresponding to hgrid="lt1v1" in models / ice / cice / bld / namelist_files / namelist_defaults_cice.xml; add the ocean horizontal and vertical grid attributes and files corresponding to ocn_grid="lt1v1” in models / ocn / licom3 / bld / namelist_files / namelist_defaults_licom3.xml;

[0118] Please refer to Figure 2 , Figure 2 is a schematic diagram of the initial calculation and continuation calculation process of the fully coupled experiment case of the coupling device based on the present invention. After adding different resolution grids of the second ocean component module and the fully coupled experiment configuration of atmosphere-land-runoff-ocean-sea ice based on the LICOM3 model in the present invention, the schematic diagram of the initial calculation and continuation calculation process of the fully coupled experiment case with -compset as B1850_LICOM and -res as ne30_lt1v1 is shown. The process is simple and clear, and the experiment process can be normally executed without any special manual configuration, which is convenient for scientific researchers to conduct various scientific experiments based on this version.

[0119] Further, the coupling application interface includes a top-level control interface, a middle-level function interface, and a component module application interface corresponding to each component module, and the component module application interface is reconstructed in the coupling step.

[0120] Among them, the steps of reconstructing the component module application interface in the coupling step include:

[0121] Modify the name of the application interface file corresponding to the component module application interface of the second ocean component module;

[0122] Refactor at least one functional function corresponding to the component module application interface of the second ocean component module;

[0123] Add a source code subdirectory of the second ocean component module to the ocean model component source code directory corresponding to the second ocean component module in the global climate system.

[0124] Specifically, the CPL7 coupler is divided into three levels: the top-level control interface, the middle-level functional interface, and the component model application interface. Among them, the top-level control interface is responsible for controlling the coupling relationship and coupling process between component models. The middle-level functional interface is responsible for processing processes such as parallel grid mapping, rearrangement, communication, and data merging between component models. The component model application interface is responsible for providing coupling data for the component model, receiving and sending coupling data, and controlling the integration operation process of the model itself. Each component model application interface has initialization, running, and end functions. The top-level control interface also has these three types of functions. Each type of function calls the corresponding functions in the corresponding component model according to the specified coupling order.

[0125] In the present invention, the second ocean component module needs to replace the first ocean component module and be connected to the CPL7 coupler as the new ocean component mode. The code of the top-level control interface and the middle-level functional interface does not need to be modified because, in these two parts, the initialization, operation, and termination of the ocean mode are independent of data exchange with other component modes, and only the corresponding functional functions in the coupling application interface need to be implemented. In the previous LICOM3 version, all coupling application interfaces were implemented in the licomcpl7.F90 file and placed in the same source directory as the implementation of the dynamic framework and physical process simulation within the LICOM mode, resulting in a rather chaotic structure that is not convenient for subsequent development. Additionally, the function names for the initialization, operation, and termination of the defined coupling application interfaces are licom_init_mct, licom_run_mct, and licom_final_mct respectively. In theory, the top-level control interface code does not need to be modified, but due to this non-standard naming, the function names for calling the coupling application interface functions of the ocn component mode in ccsm_comp_mod.F90 need to be manually modified from ocn_init_mct to licom_init_mct, ocn_run_mct to licom_run_mct, and ocn_final_mct to licom_final_mct, which itself violates the design principles of the coupler interface. Therefore, in this solution, the licomcpl7.F90 file is modified to ocn_comp_mct.F90, and some of its internal functional functions are refactored. All coupling-related application interfaces are unified and moved to the directory models / ocn / licom3 / drivers / cpl, and code adaptation is performed for the new CESM1.3 version to ensure successful compilation.

[0126] Furthermore, the construction method further includes an automated configuration step S5 of constructing: placing the namelist configuration file corresponding to the second ocean component module in the source code subdirectory. The automated configuration step of construction includes:

[0127] Adding variables required for the operation of the source code of the second ocean component module to the compilation configuration file of the namelist variable list of the second ocean component module;

[0128] Adding relevant property definitions required for the operation of the source code of the second ocean component module to the definition configuration file of the namelist variable list of the second ocean component module:

[0129] Adding default values corresponding to the grid names of different resolutions of the second ocean component module to the default value configuration file of the variable list of the second ocean component module.

[0130] Specifically, the namelist configuration in the CESM system is a crucial step for modifying model parameters without recompiling the model. Manually editing a large number of parameters is error-prone. Automated configuration helps reduce input errors and can significantly reduce the time required for manually editing the namelist file, ensuring parameter consistency and accuracy. Additionally, automated configuration also helps establish a standardized model operation process, facilitating team collaboration and model management.

[0131] The namelist generation process of CESM is highly encapsulated. Although the underlying namelist files in the RUNDIR / in directory can be directly modified, this is generally not recommended. In the previous CESM1.0.4+LICOM3 coupling scheme, this method was used, and each time it was necessary to manually modify the configuration parameters of multiple files, which was a cumbersome process and prone to errors. CESM recommends indirectly modifying these parameters through user_nl_* files and env_*.xml files. This process involves multiple steps. For example, before submitting the model, use cesm_setup or CASE.build to call the preview_namelists script. preview_namelists will automatically identify the corresponding *.buildnml.csh (such as datm.buildnml.csh for DATM, slnd.buildnml.csh for SLND, and licom3.buildnml.csh for LICOM3) of each model component actually used according to the -compset set when creating a new case (for example, the GIAF_LICOM model component is DATM_SLND_CICE_LICOM3_DROF_SGLC_SWAV). Then licom3.buildnml.csh will automatically generate the final namelist, that is, the licom3_in file, by calling models / ocn / licom3 / bld / build-namelist according to the identified grid type (i.e., $OCN_GRID). Figure 3 It is the flowchart for the automated processing related to the CESM system namelist.

[0132] In this solution, for the second ocean component module LICOM3, the namelist-related configuration is placed in the models / ocn / licom3 / bld / directory. The custom parameters used by the second ocean component module LICOM3 need to be defined and assigned through build-namelist depending on the namelist_definition_licom3.xml and namelist_defaults_licom3.xml files. Figure 4Schematic diagram of the automated configuration structure for the LICOM3 mode namelist.

[0133] The following takes the automated adaptation of the LICOM mode 10-kilometer high-resolution grid (ocn_grid = "lt0.1v1") as an example for illustration. For the lt0.1v1 high-resolution grid, new parameter "daily_accum" needs to be added to the namelist with group "namctl", and the following three files need to be modified:

[0134] (1) Add the variable add_default($nl, 'daily_accum') to the namelist variable list compilation configuration file (models / ocn / licom3 / bld / build-namelist) of the second ocean component module;

[0135] (2) Add the following definition to the namelist variable list definition configuration file (models / ocn / licom3 / bld / namelist_files / namelist_definition_licom3.xml) of the second ocean component module:

[0136] <entry

[0137] id = "daily_accum"

[0138] type = "logical"

[0139] category = "licom"

[0140] group = "namctl">

[0141] Default:.true.

[0142]

[0143] (3) Add the default values corresponding to different ocn_grid to the variable list default value configuration file (models / ocn / licom3 / bld / namelist_files / namelist_defaults_licom3.xml) of the second ocean component module:

[0144] <daily_accum>.false.< / daily_accum>

[0145] <daily_accum ocn_grid = "lt0.1v1">.true.< / daily_accum>

[0146] For the parameters already defined in namelist_definition_licom3.xml, if different values need to be set according to different ocn_grid, the following similar modifications can be directly made in namelist_defaults_licom3.xml:

[0147] <horiz_grid_opt>file< / horiz_grid_opt>

[0148] <horiz_grid_opt ocn_grid="lt1v1">tripole< / horiz_grid_opt>

[0149] <horiz_grid_opt ocn_grid="lt0.1v1">tripole< / horiz_grid_opt>

[0150] <horiz_grid_file ocn_grid="lt1v1">ocn / licom / lt1v1 / grid / lt1v1_horiz_grid.da< / horiz_grid_file>

[0151] <horiz_grid_file ocn_grid="lt0.1v1">ocn / licom / lt0.1v1 / grid / lt0.1v1_horiz_grid.da< / horiz_grid_file>

[0152] <vert_grid_file ocn_grid="lt1v1">ocn / licom / lt1v1 / grid / lt1v1_vert_grid.da< / vert_grid_file>

[0153] <vert_grid_file ocn_grid="lt0.1v1">ocn / licom / lt0.1v1 / grid / lt0.1v1_vert_grid.da< / vert_grid_file>.

[0154] Furthermore, the construction method further includes constructing a data saving and reading step S6: After enabling the PnetCDF method support of PIO through several configurations when constructing the PIO library in the global climate system, use the PnetCDF method to write the ocean-related variable data obtained by integrating the second ocean component module into a preset file.

[0155] Specifically, in the code of the second ocean component module LICOM3 used in the existing coupling device, the output mode operation results are saved in the form of binary I / O to the fort.22.*** file. This file cannot be directly used for visualization analysis, post-processing, and continuation calculation experiments. A Fortran script conversion tool is needed to read the variable data in fort.22.*** one by one and repackage it into a NetCDF file for visualization analysis and post-processing. At the same time, due to the issue of endian conversion in the storage of binary data, another Fortran script is required to perform endian conversion on the fort.22.*** file before it can be used as the input for the continuation calculation to continue the continuation calculation experiment. Therefore, the data storage method of the second ocean component module LICOM3 is extremely inconvenient to use and time-consuming and laborious.

[0156] In the present invention, the data storage and reading method for the second ocean component module is modified. Instead of using binary I / O, the parallel NetCDF (i.e., PnetCDF) of the PIO library built into the CESM system is used. The PIO (Parallel I / O) library is developed to improve the input / output (I / O) performance of the CESM component models and is automatically built as part of the CESM system construction process. By default, CESM uses the serial NetCDF implementation of PIO, and PnetCDF is turned off by default in PIO. If PnetCDF is to be used, it is first necessary to ensure that the PnetCDF library has been successfully installed on the cluster machine in use, and PnetCDF support for PIO is enabled through parameter configuration when building PIO.

[0157] The present invention uses the PnetCDF method to write data such as ocean-related variables obtained from the integration calculation of the second ocean component module, such as sea surface height, meridional velocity, zonal velocity, sea surface temperature, sea surface salinity, and coupled flux temperature, salinity, velocity, etc. into the ***.nc file. At the same time, the continuation calculation processing logic is modified to read the variable data obtained from the last calculation from the file named ***.nc (such as GIAF_lower.licom.r.0001-01-06.nc) saved most recently by the last initial calculation run pointed to by rpointer.ocn, so as to continue the mode calculation, ensuring the consistency and stability of the long-time series integration operation.

[0158] The modification method is as follows:

[0159] (1) Add the file models / ocn / licom3 / source / io_types.F90. This module contains the declarations of all the necessary IO data types (such as io_dim, io_field_desc) and the operation functions for these data types (such as construct_io_field, destroy_io_field, construct_io_dim, add_attrib_io_field, extract_attrib_io_field, init_io). It also includes some global variables (such as stdin, stdout, stderr, max_units, in_use) used by various IO operations and the IO unit manager (such as get_unit, release_unit), as well as the data file description structure (such as datafile) and its corresponding operation functions (such as construct_file, destroy_file, add_attrib_file, extract_attrib_file).

[0160] (2) Add the file models / ocn / licom3 / source / io.F90. This module provides only a general parallel input / output interface for writing arrays, namely data_set. This routine is the main interface for array and file io functions, including reading, writing, opening, and closing, and it calls different types of I / O interfaces for implementation according to the judgment of the file format ('bin' for binary, 'nc' for NetCDF).

[0161] (3) Add the file models / ocn / licom3 / source / io_netcdf.F90. This module provides a general input / output interface for writing arrays in NetCDF format using the pio library, including public interfaces such as open_read_netcdf, open_netcd, close_netcdf, sync_netcdf, field_exists_netcdf, define_field_netcdf, read_field_netcd, write_field_netcdf, define_nstd_netcd, write_nstd_netcdf.

[0162] (4)Add the file models / ocn / licom3 / source / io_pio.F90. This module mainly defines the interfaces related to the initialization of PIO, namely io_pio_init and io_pio_initdecomp.

[0163] (5)Add the file models / ocn / licom3 / source / restart.F90. This module encapsulates interface functions such as the initialization of parallel NetCDF (init_restart), saving restart data files (write_restart), reading restart data files (read_restart), and generating restart file suffix names (create_restart_suffix) implemented through the PIO library.

[0164] Figure 5 For the comparison schematic diagram of the restart execution process before and after modifying the LICOM3 model output scheme, it can be seen that the modified scheme no longer requires cumbersome manual operations such as endian conversion and renaming of the binary format fort.22.*** files using special script conversion tools. Instead, the restart experiment can be directly executed, with an intuitive, clear, concise, and efficient process, which is convenient to use.

[0165] Please refer to Figure 6 , Figure 6 which is the structural schematic diagram of the Earth system coupling device of the present invention. As Figure 6 shown, an Earth system coupling device of the present invention is obtained by the construction method described in any one of the above. The Earth system coupling device includes:

[0166] Multiple first component modules of a fully coupled global climate system. The multiple first component modules include an atmospheric component module CAM, a land component module CLM, a sea ice component module CICE, a land ice component module CISM, a river component module RTM, and a wave component module WW3;

[0167] A second component module of an ocean circulation simulation system. The second component module is a second ocean component module LICOM3;

[0168] A coupler CPL7 of a fully coupled global climate system, through which data exchange and interaction between different component modules are carried out.

[0169] In summary, the present invention also has the following beneficial effects:

[0170] (1)This solution uses the domestic ocean model LICOM3 to replace foreign models, which helps to promote the independent R & D ability of domestic climate models and enhance the right to speak in international climate science research.

[0171] (2) This solution can be compatible with different supercomputer platforms simultaneously, increasing the applicability and flexibility of the model and providing convenience for research in different computing environments.

[0172] (3) This solution can support a wider range of scientific research, including multiple research fields such as climate warming, extreme weather events, and ocean acidification.

[0173] (4) This solution promotes the comparative study between different Earth system models, helps the scientific community better understand model differences and uncertainties, and promotes data sharing and model standardization.

[0174] In summary, the CESM1.3+LICOM3 coupling solution has significant beneficial effects in improving the overall performance of climate system simulation, promoting independent technological innovation, and enhancing the scope of model application.

[0175] 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 on some of the technical features; and these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A construction method of an earth system coupling device, characterized in that, Including: First component module construction step: Construct multiple first component modules based on a fully coupled global climate system. The multiple first component modules include an atmospheric component module, a land component module, a sea ice component module, a land ice component module, a river component module, a first ocean component module, and a wave component module; Second component module construction step: Construct a second component module based on an ocean circulation simulation system. The second component module is a second ocean component module; Replacement step: Replace the first ocean component module with the second ocean component module; Coupling step: Reconstruct the coupling application interface of the coupler based on the fully coupled global climate system, and perform data exchange and interaction between different component modules through the coupling application interface; Wherein, the second component module construction step further includes: Modify use emsf to use ESMF_ClockMod, only : ESMF_Clock to introduce the single variable ESMF_Clock in the form of ESMF_ClockMod,only : ESMF_Clock; Add different machine definitions in the machine environment configuration file (scripts / ccsm_utils / Machines / config_machines.xml) to correspond to the compilation control options of different platform compilers; Add two machine configuration items with MACH as sunway and shuguang, and different - mach options can be set when executing create_newcase to select different supercomputer platforms for coupling experiments; Add adaptive compilation options in the compilation configuration file (scripts / ccsm_utils / Machines / config_compilers.xml). The compilation options only take effect when compiling the ocean component module; Wherein, the coupling application interface includes a top - level control interface, a middle - layer function interface, and a component module application interface corresponding to each component module. The coupling step reconstructs the component module application interface; The steps of reconstructing the component module application interface in the coupling step include: Modify the name of the application interface file of the component module application interface corresponding to the second ocean component module; Reconstruct at least one function function of the component module application interface corresponding to the second ocean component module; Add a source code sub - directory of the second ocean component module in the ocean model component source code directory corresponding to the second ocean component module in the global climate system.

2. The construction method according to claim 1, characterized in that, In the replacement step, preset the second ocean component module into the global climate system as an optional second ocean component module configuration in the same mode as the first ocean component module. This includes integrating the mode function code of the second ocean component module into the directory, and at the same time adding instructions to replace the first ocean component module with the second ocean component module in the system configuration file as needed.

3. The construction method according to claim 1, characterized in that It further includes a construction automation configuration step: placing the namelist configuration file corresponding to the second ocean component module in the source code subdirectory.

4. The construction method according to claim 3, characterized in that The construction automation configuration step includes: Adding variables required during the runtime of the source code of the second ocean component module to the compilation configuration file of the namelist variable list of the second ocean component module; Adding relevant property definitions required during the runtime of the source code of the second ocean component module to the definition configuration file of the namelist variable list of the second ocean component module: Adding default values corresponding to the grid names with different resolutions of the second ocean component module to the default value configuration file of the variable list of the second ocean component module.

5. The construction method according to claim 1, wherein The second component module construction step further includes: Modifying the B1850_LICOM component set in the component set configuration file; Adding res as the grid definition in the grid configuration file; Adding the pes default configuration for GRID in the parallel parameter configuration file; Adding the namelist default configuration for the grid inside the second ocean component module and the sea ice component module.

6. The construction method according to claim 1, wherein It further includes a construction data saving and reading step: After enabling the PnetCDF mode support of PIO through numerical configuration when constructing the PIO library in the global climate system, using the PnetCDF mode to write the ocean-related variable data obtained by integrating the second ocean component module into a preset file.

7. An earth system coupling device, characterized in that, Obtained by the construction method according to any one of claims 1-6 above, the earth system coupling device includes: Multiple first component modules of a fully coupled global climate system, and the multiple first component modules include an atmospheric component module, a land component module, a sea ice component module, a land ice component module, a river component module, and a wave component module; A second component module of an ocean circulation simulation system, and the second component module is a second ocean component module; A coupler of a fully coupled global climate system, and data exchange and interaction between different component modules are carried out through the coupler.

Images