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

Abstract

The invention discloses a construction method of an earth system coupling device and the earth system coupling device, and the construction method comprises the steps: a first component module construction step: constructing a plurality of first component modules of a global climate system based on full coupling, the plurality of first component modules comprise an atmospheric component module, a land component module, a sea ice component module, a ground ice component module, a river component module, a first ocean component module and a sea wave component module; a second component module construction step: constructing a second component module based on the ocean circulation simulation system, wherein the second component module is a second ocean component module; a replacement step: replacing the first ocean component module with the second ocean component module; and a coupling step: reconstructing a coupling application interface of the coupler of the global climate system based on full coupling, and carrying out data exchange and interaction among different component modules through the coupling application interface.

Description

Technical Field

[0001] The present invention relates to the field of earth system science and technology, and in particular 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 developed to LICOM3 (LASG / IAP Climate System Ocean Model Version 3 (LICOM3)). The model includes multiple sets of resolutions and parameterization schemes, with the highest resolution reaching 3-5km globally. LICOM, as an ocean sub-model, constructs an earth system model that can simulate past geological periods and current climate change, and estimate climate change in the next 10-30 years or even the next 100 years. The Community Earth System Model (CESM) is a fully coupled global climate model that provides users with the most advanced computer simulation of the earth's past, present and future climate states, and is an advanced technology for earth system simulation. The initial CESM1 version only includes 5 geophysical submodules, namely atmosphere (CAM), ocean (POP), land (CLM), sea ice (CICE) and ice sheet (CISM), and newly added river transport model (RTM) and 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 Intercomparison Project) with the Earth system coupling device of CESM1 (version 1.0.4) + LICOM3, and conducted global ocean-sea ice coupling experiments by providing a unified extra-atmospheric forcing data set and flux calculation scheme. 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: (1) The existing LICOM3 mode coupling scheme is no longer applicable to higher versions of the CESM framework; (2) There is a lack of automated parameter configuration methods for coupled experiments related to the LICOM3 ocean model. The steps are cumbersome and multiple configuration files and parameters need to be manually modified according to grids of different resolutions. This is prone to errors and has poor inheritance. (3) Previously, only ocean-sea ice coupling was applicable in the CESM1.0.4+LICOM3 coupling device, and it was impossible to realize the full coupling experiment of atmosphere-land-runoff-ocean-sea ice.

[0004] (4) The calculation results output by the previous LICOM3 mode are saved in binary form and need to be converted after post-processing before visualization analysis or post-processing. The process is complicated and consumes unnecessary resources.

[0005] Therefore, there is an urgent need to develop a construction method of an earth system coupling device and an earth system coupling device that overcomes the above-mentioned defects. Summary of the invention

[0006] In view of the above problems, the present invention provides a method for constructing an earth system coupling device, which comprises: First component module construction step: constructing a plurality of first component modules based on a fully coupled global climate system, wherein 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; Second component module construction step: constructing a second component module based on the ocean circulation simulation system, wherein the second component module is a second ocean component module; Substituting step: replacing the first ocean component module by the second ocean component module; Coupling step: reconstructing the coupling application interface of the coupler based on the fully coupled global climate system, and performing data exchange and interaction between different component modules through the coupling application interface.

[0007] In the above construction method, the second component module construction step further includes: Change 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 to correspond to the compilation control options corresponding to the compilers of different platforms; Added MACH as two machine configuration items for sunway and shuguang. You can choose to perform coupling experiments on different supercomputing platforms by setting different -mach options when executing create_newcase; An adaptive compilation option is added in the compilation configuration file, and the compilation option takes effect only when compiling the ocean component module.

[0008] The above-mentioned construction method, wherein, in the replacement step, the second ocean component module is preset into the global climate system as an equivalent mode of the first ocean component module as an optional second ocean component module configuration, which includes: integrating the mode function code of the second ocean component module into the directory, and adding instructions to replace the first ocean component module with the second ocean component module in the system configuration file as needed.

[0009] In the above construction method, 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.

[0010] In the above construction method, the step of reconstructing the component module application interface in the coupling step includes: Modify the name of the application interface file of the component module application interface corresponding to the second ocean component module; Reconstructing at least one function of the component module application interface corresponding to the second ocean component module; A source code subdirectory of the second ocean component module is added to the ocean model component source code directory corresponding to the second ocean component module in the global climate system.

[0011] The above construction method further includes a construction automatic configuration step: placing the namelist configuration file corresponding to the second ocean component module in the source code subdirectory.

[0012] In the above construction method, the step of building an automated configuration includes: Add the variables required by the source code of the second ocean component module to the namelist variable list compilation configuration file of the second ocean component module; In the namelist variable list definition configuration file of the second ocean component module, add the relevant attribute definitions required when the source code of the second ocean component module is running: In the variable list default value configuration file of the second ocean component module, default values ​​corresponding to the grid names of different resolutions of the second ocean component module are added.

[0013] In the above construction method, the second component module construction step further comprises: Modify the B1850_LICOM component set in the component set configuration file; Add res as the grid definition in the grid configuration file; Add pes default configuration for GRID in the parallel parameter configuration file; Added default configuration for namelist of grids inside the second ocean component module and the sea ice component module.

[0014] The above construction method further includes the steps of building data saving and reading: When constructing the PIO library in the global climate system, after enabling the PnetCDF method support of PIO through digital configuration, the ocean-related variable data obtained by the integral calculation of the second ocean component module is written into a preset file using the PnetCDF method.

[0015] The present invention further provides an earth system coupling device, wherein the earth system coupling device is obtained by any one of the construction methods described above, and comprises: A plurality of first component modules of a fully coupled global climate system, wherein 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 and an ocean wave component module; A second component module of the ocean circulation simulation system, wherein the second component module is a second ocean component module; A coupler of a fully coupled global climate system, through which data exchange and interaction between different component modules takes place.

[0016] In summary, the present invention has the following advantages over the prior art: (1) LICOM3 was coupled to the CESM1.3 high-resolution version system as a new ocean component module, and was compatible with multiple supercomputer platforms such as the domestic Shenwei supercomputer and X86 supercomputer; (2) Design a new LICOM3 coupling application interface with strong compatibility to facilitate the update and upgrade of subsequent CESM versions; (3) Design an automated experimental scheme for the ocean component module of LICOM3. The experimental steps are simple and clear, which greatly simplifies the operation of scientific experiments. (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; (5) Design a new parallel data standardized output scheme for the ocean component module of LICOM3, simplify the original post-processing scheme, and achieve high-speed and formatted output of output data.

[0017] Other features and advantages of the present invention will be described in the following description, and partly become apparent from the description, or understood by practicing the present invention. The purpose and other advantages of the present invention can be realized and obtained by the structures pointed out in the description, claims and drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0018] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0019] Figure 1 is a flow chart of the construction method of the present invention; Figure 2 It is a schematic diagram of the initial calculation and continued calculation process of the full coupling experiment case based on the coupling device of the present invention; Figure 3 This is the flowchart of the namelist-related automated processing of the CESM system; Figure 4 This is a schematic diagram of the namelist automation configuration structure for LICOM3 mode; Figure 5 This is a comparison diagram of the follow-up execution process before and after the output plan is modified; Figure 6 It is a schematic structural diagram of the earth system coupling device of the present invention. DETAILED DESCRIPTION

[0020] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

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

[0022] The terms “first”, “second”, “S1”, “S2”, etc. used in this document do not specifically refer to an 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.

[0023] The directional terms used in this document, such as up, down, left, right, front or back, etc., are only used to refer to the directions of the drawings. Therefore, the directional terms used are used to illustrate and not to limit the present invention.

[0024] The words “include,” “including,” “have,” “contain,” etc. used in this document are open-ended terms, meaning including but not limited to.

[0025] As used herein, "and / or" includes any or all combinations of the items described.

[0026] As used herein, “plurality” includes “two” and “more than two”; as used herein, “plurality groups” includes “two groups” and “more than two groups”.

[0027] The terms "substantially" and "approximately" used herein are used to modify any quantity or error that may vary slightly, but these slight changes or errors do not change their essence. Generally speaking, the range of slight changes or errors modified by such terms may 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.

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

[0029] Please refer to Figure 1 Flow chart of the construction method of the present invention. Figure 1 As shown, a method for constructing an earth system coupling device of the present invention comprises: First component module construction step S1: constructing a plurality of first component modules based on a fully coupled global climate system, wherein 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; Second component module construction step S2: constructing a second component module based on the ocean circulation simulation system, wherein the second component module is a second ocean component module; Substitution step S3: replacing the first ocean component module with the second ocean component module; Coupling step S4: reconstructing a coupling application interface of a coupler based on a fully coupled global climate system, through which data exchange and interaction between different component modules are performed.

[0030] Specifically, based on the modular framework of the CESM system, after constructing a main body consisting 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 constructed coupler (CPL7) is used to manage the data information exchange and mode operation between the modules. 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 RTM (River Transport Model) and the wave component module WW3 (WaveWatch III), wherein the first ocean component module is replaced by the second component module based on the ocean circulation simulation system; Therefore, 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 a CESM1.3+LICOM3-based earth system coupling device, so that users can use the simplest case creation, compilation and execution steps to conduct coupling experiments related to the LICOM3 ocean model.

[0031] Wherein, in the replacement step, the second ocean component module is preset into the global climate system as an equivalent mode of 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 the instruction of replacing 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 the COMPSET replacing the default ocean component of POP2 as LICOM3 is added as needed in the system configuration file scripts / ccsm_utils / Case.template / config_compsets.xml, so as to directly create a coupled experiment with the ocean model as LICOM3.

[0032] Furthermore, the second component module construction step also includes: Change 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 to correspond to the compilation control options corresponding to the compilers of different platforms; Added MACH as two machine configuration items for sunway and shuguang. You can choose to perform coupling experiments on different supercomputing platforms by setting different -mach options when executing create_newcase; An adaptive compilation option is added in the compilation configuration file, and the compilation option takes effect only when compiling the ocean component module.

[0033] Specifically, in response to some error messages encountered during the compilation and operation of the second ocean component module on the new version of CESM1.3, some codes of the second ocean component module were modified. For example, all the places where use esmf was declared in the second ocean component module LICOM3 were deleted, and only the places where it was actually called in the ocn_comp_mct.F90 file were retained. However, because only one variable ESMF_Clock was used, use emsf was changed to use ESMF_ClockMod, only :ESMF_Clock. According to the design principles of the CESM system, different machine definitions can be added to the machine environment configuration file (scripts / ccsm_utils / Machines / config_machines.xml) to define the compilation control options corresponding to different platform compilers, so that the CESM system can be adapted to run on different platforms. In this solution, two machine configuration items with MACH as sunway and shuguang are added for the domestic Shenwei supercomputer and X86 supercomputer respectively. When executing create_newcase, different -mach options can be set to select different supercomputer platforms for coupling experiments, and adaptive compilation options are added for the second ocean component module in the compilation configuration file (scripts / ccsm_utils / Machines / config_compilers.xml).<ADD_FFLAGS MODEL="licom3"> -fdefault-double-8-fdefault-real-8< / ADD_FFLAGS> This option only works 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.

[0034] In this embodiment, the second component module construction step further includes: Modify the B1850_LICOM component set in the component set configuration file (scripts / ccsm_utils / Case.template / config_compsets.xml); Add res as grid definition in the grid configuration file (scripts / ccsm_utils / Case.template / config_grid.xml); Add the pes default configuration for GRID in the parallel parameter configuration file (scripts / ccsm_utils / Machines / config_pes.xml); Added default configuration for namelist of grids inside the second ocean component module and the sea ice component module.

[0035] 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 < / 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: (1) Add the following line to the component set configuration file (scripts / ccsm_utils / Case.template / config_compsets.xml) for compset: B1850_LICOM adds the following definition: <compset sname="B_1850_LICOM" alias="B1850_LICOM"> 1850_CAM4_CLM40%SP_CICE_LICOM3_RTM_SGLC_SWAV< / compset> (2) Add res to the ne30np4_lt1v1 grid related definition in the grid configuration file (scripts / ccsm_utils / Case.template / config_grid.xml): <grid sname="ne30np4_lt1v1" alias="ne30_lt1v1"> a%ne30np4_l%ne30np4_oi%lt1v1_r%r05_m%lt1v1_g%null_w%null< / grid> <gridhorz name="lt1v1"> <nx> 360< / nx> <ny> 218< / ny> <desc> lt1v1 is tripole 1-deg grid for LICOM low resolution:< / desc> < / gridhorz> <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> <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> < / griddom> <gridmap atm_grid="ne30np4" ocn_grid="lt1v1"> <ATM2OCN_FMAPNAME>cpl / gridmaps / ne30np4 / map_ne30np4_to_lt1v1_aave.20240809.nc< / ATM2OCN_FMAPNAME> <ATM2OCN_SMAPNAME>cpl / gridmaps / ne30np4 / map_ne30np4_to_lt1v1_blin.20240809.nc< / ATM2OCN_SMAPNAME> <ATM2OCN_VMAPNAME>cpl / gridmaps / ne30np4 / map_ne30np4_to_lt1v1_patc.20240809.nc< / ATM2OCN_VMAPNAME> <OCN2ATM_FMAPNAME>cpl / gridmaps / lt1v1 / map_lt1v1_to_ne30np4_aave.20240809.nc< / OCN2ATM_FMAPNAME> <OCN2ATM_SMAPNAME>cpl / gridmaps / lt1v1 / map_lt1v1_to_ne30np4_aave.20240809.nc< / OCN2ATM_SMAPNAME> < / gridmap> <gridmap rof_grid="r05" ocn_grid="lt1v1"><ROF2OCN_FMAPNAME>cpl / gridmaps / r05 / map_r05_TO_lt1v1_aave.240809.nc< / ROF2OCN_FMAPNAME> < / gridmap> (3) Add the following line for GRID in the parallel parameter configuration file (scripts / ccsm_utils / Machines / config_pes.xml): The default configuration of pes for "a%ne30np4_l%ne30np4_oi%lt1v1": <pes grid="a%ne30np4_l%ne30np4_oi%lt1v1" mach="sunway"> <NTASKS_ATM>80< / NTASKS_ATM><NTHRDS_ATM>1< / NTHRDS_ATM><ROOTPE_ATM>0< / ROOTPE_ATM><NINST_ATM>1< / NINST_ATM> <NTASKS_LND>$NTASKS_ATM< / NTASKS_LND><NTHRDS_LND>$NTHRDS_ATM< / NTHRDS_LND><ROOTPE_LND>$ROOTPE_ATM< / ROOTPE_LND><NINST_LND>1< / NINST_LND> <NTASKS_ICE>$NTASKS_ATM< / NTASKS_ICE><NTHRDS_ICE>$NTHRDS_ATM< / NTHRDS_ICE><ROOTPE_ICE>$ROOTPE_ATM< / ROOTPE_ICE><NINST_ICE>1< / NINST_ICE> <NTASKS_OCN>$NTASKS_ATM< / NTASKS_OCN><NTHRDS_OCN>$NTHRDS_ATM< / NTHRDS_OCN><ROOTPE_OCN>$ROOTPE_ATM< / ROOTPE_OCN><NINST_OCN>1< / NINST_OCN> <NTASKS_CPL>$NTASKS_ATM< / NTASKS_CPL><NTHRDS_CPL>$NTHRDS_ATM< / NTHRDS_CPL><ROOTPE_CPL>$ROOTPE_ATM< / ROOTPE_CPL> <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> <PES_LEVEL>0< / PES_LEVEL> < / pes> .

[0036] (4) At the same time, the namelist default configuration for the grid "lt1v1" needs to be added in 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 properties and files corresponding to ocn_grid="lt1v1" in models / ocn / licom3 / bld / namelist_files / namelist_defaults_licom3.xml; Please refer to Figure 2 , Figure 2 This is a schematic diagram of the initial calculation and continued calculation process of the fully coupled experimental case based on the coupling device of the present invention. This is a schematic diagram of the initial calculation and continued calculation process of the fully coupled experimental case with -compset as B1850_LICOM and -res as ne30_lt1v1 after adding different resolution grids of the second ocean component module and the fully coupled experimental configuration of atmosphere-land-runoff-ocean-sea ice based on the LICOM3 mode in the present invention. The process is concise and clear, and the experimental process can be executed normally without any special manual configuration, which is convenient for researchers to conduct various scientific experiments based on this version.

[0037] Furthermore, 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.

[0038] The step of reconstructing the component module application interface in the coupling step includes: Modify the name of the application interface file of the component module application interface corresponding to the second ocean component module; Reconstructing at least one function of the component module application interface corresponding to the second ocean component module; A source code subdirectory of the second ocean component module is added to the ocean model component source code directory corresponding to the second ocean component module in the global climate system.

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

[0040] In the present invention, the second ocean component module replaces the first ocean component module as a new ocean component mode connected to the CPL7 coupler, and the codes of the top-level control interface and the middle-level functional interface do not need to be modified, because in these two parts, the initialization, operation and termination of the ocean mode itself are irrelevant to the data exchange of 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 were placed in the same source directory as the dynamic framework, physical process simulation and other implementations within the LICOM mode. The structure was relatively chaotic and inconvenient for subsequent development. In addition, the function names for initialization, operation and termination of the coupling application interface defined were 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 since this name definition does not meet the standard, it is also necessary to manually modify the name of the function that calls the ocn component mode coupling application interface in ccsm_comp_mod.F90, from ocn_init_mct to licom_init_mct, ocn_run_mct to licom_run_mct, and ocn_final_mct to licom_final_mct, which in itself violates the design principles of the coupler interface. Therefore, in this solution, the licomcpl7.F90 file was modified to ocn_comp_mct.F90 and some of its internal functions were reconstructed. All application interfaces related to coupling were moved to the directory models / ocn / licom3 / drivers / cpl. At the same time, the code was adapted to the new CESM1.3 version to ensure that it can be compiled successfully.

[0041] Furthermore, the construction method further includes an automatic construction configuration step S5: placing the namelist configuration file corresponding to the second ocean component module in the source code subdirectory. The automatic construction configuration step includes: Add the variables required by the source code of the second ocean component module to the namelist variable list compilation configuration file of the second ocean component module; In the namelist variable list definition configuration file of the second ocean component module, add the relevant attribute definitions required when the source code of the second ocean component module is running: In the variable list default value configuration file of the second ocean component module, default values ​​corresponding to the grid names of different resolutions of the second ocean component module are added.

[0042] Specifically, namelist configuration in the CESM system is a key step used to modify model parameters without recompiling the model. Manually editing a large number of parameters is prone to errors. Automated configuration helps reduce input errors and can significantly reduce the time required to manually edit namelist files, ensuring parameter consistency and accuracy. In addition, automated configuration also helps to establish a standardized model operation process, facilitating team collaboration and model management.

[0043] The namelist generation process of CESM is highly encapsulated. Although you can directly modify the underlying namelist files in the RUNDIR / in directory, it is generally not recommended to do so. The previous CESM1.0.4+LICOM3 coupling solution adopted this method. Each time, you need to manually modify the configuration parameters of multiple files, which is a cumbersome and error-prone process. CESM recommends indirectly modifying these parameters through user_nl_* files and env_*.xml files. This process involves multiple steps, such as using cesm_setup or CASE.build to call the preview_namelists script before submitting the model. preview_namelists will automatically identify the *.buildnml.csh corresponding to the actual mode components (such as GIAF_LICOM mode components are DATM_SLND_CICE_LICOM3_DROF_SGLC_SWAV) according to the -compset set when create_newcase (such as DATM corresponds to datm.buildnml.csh, SLND corresponds to slnd.buildnml.csh, LICOM3 corresponds to licom3.buildnml.csh), and licom3.buildnml.csh will then automatically generate the final namelist, i.e., licom3_in file, according to the identified grid type (i.e., $OCN_GRID) by calling models / ocn / licom3 / bld / build-namelist. Figure 3 This is the flowchart of namelist-related automated processing in the CESM system.

[0044] In this solution, the namelist-related configurations for the second ocean component module LICOM3 are 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 the build-namelist dependency namelist_definition_licom3.xml and namelist_defaults_licom3.xml files. Figure 4 This is a schematic diagram of the LICOM3 mode namelist automatic configuration structure.

[0045] The following takes the automatic adaptation of the 10 km high-resolution grid (ocn_grid="lt0.1v1") in LICOM mode as an example. For the lt0.1v1 high-resolution grid, a new parameter "daily_accum" needs to be added to the namelist of the group "namctl". The following three files need to be modified: (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; (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: <entry id="daily_accum" type="logical" category="licom" group="namctl"> Default: .true. (3) Add the default values ​​corresponding to different ocn_grid in the variable list default value configuration file (models / ocn / licom3 / bld / namelist_files / namelist_defaults_licom3.xml) of the second ocean component module: <daily_accum> .false.< / daily_accum> <daily_accum ocn_grid="lt0.1v1">.true.< / daily_accum> For the parameters defined in namelist_definition_licom3.xml, if you want to set different values ​​according to different ocn_grids, you can directly make the following changes in namelist_defaults_licom3.xml: <horiz_grid_opt> file< / horiz_grid_opt> <horiz_grid_opt ocn_grid="lt1v1"> tripole< / horiz_grid_opt> <horiz_grid_opt ocn_grid="lt0.1v1"> tripole< / horiz_grid_opt> <horiz_grid_file ocn_grid="lt1v1"> ocn / licom / lt1v1 / grid / lt1v1_horiz_grid.da< / horiz_grid_file> <horiz_grid_file ocn_grid="lt0.1v1"> ocn / licom / lt0.1v1 / grid / lt0.1v1_horiz_grid.da< / horiz_grid_file> <vert_grid_file ocn_grid="lt1v1"> ocn / licom / lt1v1 / grid / lt1v1_vert_grid.da< / vert_grid_file> <vert_grid_file ocn_grid="lt0.1v1"> ocn / licom / lt0.1v1 / grid / lt0.1v1_vert_grid.da< / vert_grid_file> .

[0046] Furthermore, the construction method also includes constructing a data saving and reading step S6: after enabling the PnetCDF method support of PIO through digital configuration when constructing the PIO library in the global climate system, the ocean-related variable data obtained by the integral calculation of the second ocean component module is written into a preset file using the PnetCDF method.

[0047] Specifically, in the code of the second ocean component module LICOM3 used in the existing coupling device, the output mode operation result is saved in the form of binary I / O to the fort.22.*** file, which cannot be directly visualized, post-processed, or executed in a continued calculation experiment. A Fortran script conversion tool is required to read the variable data in fort.22.*** one by one and repackage them into a NetCDF file for visualization analysis and post-processing. At the same time, due to the problem of big-endian conversion in the storage of binary data, another Fortran script is required to convert the fort.22.*** file into big-endian before it can be used as the input for continued calculation to continue the continued calculation experiment. Therefore, the data storage method of the second ocean component module LICOM3 is too difficult to use and is time-consuming and labor-intensive.

[0048] In the present invention, the data saving and reading method of the second ocean component module is modified. Instead of using binary I / O, the method uses the parallel NetCDF (i.e., PnetCDF) method through the PIO library built into the CESM system. The PIO (Parallel I / O) library is developed to improve the input / output (I / O) performance of the CESM component model. It is automatically built as part of the CESM system build process. By default, CESM uses the serial NetCDF implementation of PIO, and PnetCDF is disabled by default in PIO. If you want to use PnetCDF, you must first ensure that the PnetCDF library has been successfully installed on the cluster machine used, and enable PnetCDF support for PIO through parameter configuration when building PIO.

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

[0050] The modification method is as follows: (1) Add the models / ocn / licom3 / source / io_types.F90 file. This module contains declarations of all required IO data types (such as io_dim, io_field_desc) and operation functions of 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 IO unit managers (such as get_unit, release_unit), as well as data file description structures (such as datafile) and their corresponding operation functions (such as construct_file, destroy_file, add_attrib_file, extract_attrib_file).

[0051] (2) Add the models / ocn / licom3 / source / io.F90 file. This module only provides 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. It will call different types of I / O interface implementations based on the file format ('bin' for binary, 'nc' for NetCDF).

[0052] (3) Add the models / ocn / licom3 / source / io_netcdf.F90 file. 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, and write_nstd_netcdf.

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

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

[0055] Figure 5 This is a comparison diagram of the subsequent calculation execution process before and after modifying the LICOM3 mode output scheme. It can be seen that the modified scheme no longer needs to use special script conversion tools to convert the binary format fort.22.*** file into big and small endian and rename it, and other tedious manual operations. The continued calculation experiment can be directly executed. The process is intuitive, clear, concise, efficient and easy to use.

[0056] Please refer to Figure 6 , Figure 6 FIG. 1 is a schematic diagram of the structure of the earth system coupling device of the present invention. Figure 6 As shown, an earth system coupling device of the present invention is obtained by any of the construction methods described above, and the earth system coupling device includes: A plurality of first component modules of a fully coupled global climate system, wherein the plurality of 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; A second component module of the ocean circulation simulation system, wherein the second component module is a second ocean component module LICOM3; The coupler CPL7 of the fully coupled global climate system is used to exchange data and interact with each other between the different component modules.

[0057] In summary, the present invention also has the following beneficial effects: (1) This plan uses the domestic ocean model LICOM3 to replace the foreign model, which will help promote the domestic ability to independently develop climate models and enhance the voice in international climate science research.

[0058] (2) This solution is compatible with different supercomputing platforms, which increases the applicability and flexibility of the model and facilitates research in different computing environments.

[0059] (3) This program can support a wider range of scientific research, including climate change, extreme weather events, ocean acidification, and other research areas.

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

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

[0062] Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent substitutions for some of the technical features therein; and these modifications or substitutions 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 method for constructing an earth system coupling device, characterized in that: include: First component module construction step: constructing a plurality of first component modules based on a fully coupled global climate system, wherein 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; Second component module construction step: constructing a second component module based on the ocean circulation simulation system, wherein the second component module is a second ocean component module; Substituting step: replacing the first ocean component module by the second ocean component module; Coupling step: reconstructing the coupling application interface of the coupler based on the fully coupled global climate system, and performing data exchange and interaction between different component modules through the coupling application interface.

2. The construction method according to claim 1, characterized in that: The second component module construction step also includes: Change 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 corresponding to the compilers of different platforms; Added MACH as two machine configuration items for sunway and shuguang. You can choose to perform coupling experiments on different supercomputing platforms by setting different -mach options when executing create_newcase; An adaptive compilation option is added in the compilation configuration file (scripts / ccsm_utils / Machines / config_compilers.xml). The compilation option takes effect only when compiling the ocean component module.

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

4. The construction method according to claim 1, characterized in that: 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. The component module application interface is reconstructed in the coupling step.

5. The construction method according to claim 4, characterized in that: The step of reconstructing the component module application interface in the coupling step includes: Modify the name of the application interface file of the component module application interface corresponding to the second ocean component module; Reconstructing at least one function of the component module application interface corresponding to the second ocean component module; A source code subdirectory of the second ocean component module is added to the ocean model component source code directory corresponding to the second ocean component module in the global climate system.

6. The construction method according to claim 5, characterized in that: The method also includes building an automated configuration step: placing a namelist configuration file corresponding to the second ocean component module in the source code subdirectory.

7. The construction method according to claim 6, characterized in that: The steps of building the automated configuration include: Add the variables required by the source code of the second ocean component module to the namelist variable list compilation configuration file of the second ocean component module; In the namelist variable list definition configuration file of the second ocean component module, add the relevant attribute definitions required when the source code of the second ocean component module is running: In the variable list default value configuration file of the second ocean component module, default values ​​corresponding to the grid names of different resolutions of the second ocean component module are added.

8. The construction method according to claim 2, characterized in that: The second component module construction step also includes: Modify the B1850_LICOM component set in the component set configuration file; Add res as the grid definition in the grid configuration file; Add pes default configuration for GRID in the parallel parameter configuration file; Added default configuration for namelist of grids inside the second ocean component module and the sea ice component module.

9. The construction method according to claim 1, characterized in that: It also includes building data saving and reading steps: When constructing the PIO library in the global climate system, after enabling the PnetCDF method support of PIO through digital configuration, the ocean-related variable data obtained by the integral calculation of the second ocean component module is written into a preset file using the PnetCDF method.

10. An earth system coupling device, characterized in that: Obtained by the construction method described in any one of claims 1 to 9 above, the earth system coupling device comprises: A plurality of first component modules of a fully coupled global climate system, wherein 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 and an ocean wave component module; A second component module of the ocean circulation simulation system, wherein the second component module is a second ocean component module; A coupler of a fully coupled global climate system, through which data exchange and interaction between different component modules takes place.

Images