资源描述
ENES and PRISM:A European approach to Earth System modellingSophie Valcke,CERFACSand the PRISM team across EuropeSophie Valcke,CERFACSOutlineENESThe PRISM project:partnersgoalsmodel componentsstandard physical interfacesarchitecture and GUIPRISM first coupler:Oasis3(08/2019):configuration communicationinterpolations/transformationsPRISM final coupler(12/2019)configuration,communication,interpolations/transformationsComputing in Atmospheric Sciences Workshop 20192OutlineENESComputing in AtmospENESClimate research in Europe:Societal/political needs in Europe high(IPCC,mitigation,)Recognised excellence;scientific diversity(models,approaches,)How to organise Earth System modelling in Europe?The one-big-centre-does-it-all not suitable:-expertise lies within national centres-diversity is key to research Need for shared infrastructures:-software(PRISM)-hardwareComputing in Atmospheric Sciences Workshop 20193ENESClimate research in EuropeENESENES:European Network for Earth System modelling Think tank to organize,plan and seek funding for efficient distributed Earth System modelling in Europe Follows a EuroClivar recommendation Open to any institute/industry(MoU)Coordinated by Guy Brasseur(MPI,Hamburg)50 members so far(enes.org)Computing in Atmospheric Sciences Workshop 20194ENESENES:European Network forA long term strategy:1.Jointly develop shared software infrastructure for Earth System modelling2.Provide European integrated service to access and use this infrastructure3.Provide and manage hyper-computing access by 2019PRISMENESComputing in Atmospheric Sciences Workshop 20195A long term strategy:Jointly dThe PRISM projectPRISM:PRogram for Integrated Earth System Modelling A European project,started December 2019,fundedfor 3 years by the European Commission(4.8 M)Coordinators:Guy Brasseur(MPI,Hamburg)Gerbrand Komen (KNMI,Amsterdam)PRISM Director:Reinhard Budich(MPI)Computing in Atmospheric Sciences Workshop 20196The PRISM projectPRISM:PRogra=22 partners:leading climate research institutes and computer vendorsMPG-IMET,GermanyKNMI,NetherlandsMPI-MAD,GermanyMet-Office,UKUREADMY,UKIPSL,FranceMto-France,FranceCERFACS,FranceDMI,DenmarkSHMI,SwedenNERSC,Norway ETH Zurich,Switzerland ING,ItalyMPI-BGC,GermanyPIK,GermanyECMWF,EuropeUCL-ASTR,BelgiumNEC EuropeFECIT/FujitsuSGI EuropeSUN EuropePRISM partnersComputing in Atmospheric Sciences Workshop 20197=22 partners:leading climatProvide software infrastructure to easily assemble Earth system coupled models based on existing state-of-art European component modelslaunch/monitor complex/ensemble Earth system simulationsPRISM goalsHelp climate modellers spend more time on science:Computing in Atmospheric Sciences Workshop 20198Provide software infrastructur Define and promote technical and scientific standards for ESM:Scientific standards:Physical interfaces between model componentsGlobal Earth System parameters Technical standards:Compiling,running,post-processing environmentArchitecture and Graphical User Interface Coupler and I/O softwareData and grid formatCoding and qualityInteraction with other groups(ESMF,ESG/NOMADS,CF,RPN?,.)PRISM goalsComputing in Atmospheric Sciences Workshop 20199 Define and promote technical Atmosphere:Mto-France(ARPEGE),MPG-IMET(ECHAM),IPSL(LMDZ),MetOffice(Unified Model),UREADMY,INGVAtmospheric Chemistry:MPG-IMET,UREADMY,IPSL,MetOffice,Mto-France,KNMILand Surface:IPSL(Orchide),MetOffice,MPG-IMET,UREADMY,Mto-France(ISBA)Sea Ice:NERSC,UCL-ASTR,MetOffice,IPSL,MPG-IMETOcean Biogeochemistry:MPI-BGC,IPSL,MPG-IMET,MetOfficeOcean:UREADMY,MetOffice(FOAM),MPI-M(HOPE),IPSL(OPA/ORCA)Regional Climate:SHMI,DMI,MetOfficeCoupler:CERFACS,NEC,CCRLE,FECIT,SGI,MPI-MAD PRISM model componentsComputing in Atmospheric Sciences Workshop 201910Atmosphere:Atmospheric ChemistOcean model1-Rainfall+int.energy2-Snowfall+int.energy3-Incoming solar radiat.4-Solar zenith angle5-Fraction of diffuse solar radiation6-Downward infrared radiation7-Sensitivity of atmos temp.&humidity to surf.fluxes1*-Sensible heat flux2*-Surf.emissivity 3*-Albedo,direct4*-Albedo,diffuse5*-Surf.radiative temp.6*-Evaporation+int.energy+Qlat7*-Wind stress8-Subgrid fractions1-Surface pressure2-4 Air temperature,humidity and wind5-Wind module6-Height of these 4 variables1*-rCd2*-rCe3*-rCh1x-Non solar heat flux2x-Solar radiation3x-Fresh water flux4x-Salt flux5x-Wind stress6x-U37x-Mass of snow and ice8-Subgrid fractions1*-Surf.Temp2*-Surf.Roughness3*-Displacement height4x-Surface velocity1-Continental runoff +internal Energy1-2 Temp./Salinity at sea-ice base3-Sea surface temperature4-Surf.radiative temp.5-Surface ocean current6-Sea surface salinity7-Surface height8-Absorbed solar radiation (in first oceanic layer)Iceberg parameters12345678 Ocean surface moduleSurface layer turbulenceSea ice modelwave model+3412Atmosphere model Land surface modelNote on subgrid fractiondependance:x-Sea Ice categories (incl.open ocean)*-Sea Ice or Land Surf.categoriesA proposal for PRISM standard O-A-SI physical interfaces:Computing in Atmospheric Sciences Workshop 201911Ocean model1-Rainfall+int.PRISM central server+PRISM local sitesGUI:adaptation of ECMWF prepIFS and SMS scheduler User/developerPrism centralserverPrism localsiteWebWebPRISM architecture and GUI:Computing in Atmospheric Sciences Workshop 201912PRISM central server+PRISM lBased on Oasis developed since 1991 in CERFACS to couple existing GCMs developed independently at the time:Models at relatively low resolution(10000-20000 pts)Small number of 2D coupling fields(10)Low coupling frequency(once/day)flexibility was very important,efficiency not so much!performs:synchronisation of the component models coupling fields exchange and interpolation I/O actionstested on VPP5000,NEC SX5,SGI Octane and O3000,Compaq Alpha cluster,Linux PC cluster(MPI-Globus)PRISM first coupler:Oasis3 Computing in Atmospheric Sciences Workshop 201913Based on Oasis developed sinceOasis regular users:CERFACSMETEO-FRANCE(France)IPSL-LODYC,LMD,LSCE(France)ECMWF(UK)UCL(Belgium)MPI-M&D(Germany)SMHI(Sweden)BMRC(Australia)IRI(USA)and:AWI(Germany)PIK(Germany)Met Office(UK)UGAMP(UK)KNMI(Netherlands)CSIRO(Australia)FSU/COAPS(USA)LASG(China)JAMSTEC(Japan)?PRISM project first coupler:Oasis3Computing in Atmospheric Sciences Workshop 201914Oasis regular users:and:PRISPRISM project first coupler:Oasis3Oasis3 configuration:In text file namcouple read by Oasis3 at the beginning of the run,e.g.total run timenumber and names of component modelsnumber and names of coupling fields;for each field:source and target symbolic namecoupling and/or I/O status,coupling or I/O periodtransformations/interpolations Component model grid(longitudes,latitudes,masks,mesh surfaces,mesh corner locations)must be available in binary or NetCDF files.Computing in Atmospheric Sciences Workshop 201915PRISM project first coupler:OPRISM project first coupler:Oasis3Oasis3 communication:New PRISM System model interface(PSMILe)based on MPI1 or MPI2 message passingParallel communication between parallel models and Oasis3 interpolation process A A A B B B A A Afile A A A O O O OOasis3Direct communication between models with same grid and partitioningI/O functionality(automatic switch between coupled and forced mode)Modularity:at each model time step,exchange is performed or not depending on users specifications in namcouple.Automatic time integration depending on users specificationComputing in Atmospheric Sciences Workshop 201916PRISM project first coupler:OOasis3 interpolations/transformationsPRISM project first coupler:Oasis3 A A A O O O OOasis3Oasis3=performed by separate sequential process=on 2D scalar fields only Interfacing with RPN Fast Scalar INTerpolator package nearest-neighbour,bilinear,bicubic for regular Lat-Lon grids Interfacing with SCRIP1.4 library(Los Alamos Software Release LACC 98-45):nearest-neighbour,1st and 2nd order conservative remapping for all grids bilinear and bicubic interpolation for logically-rectangular grids Bilinear and bicubic interpolation for reduced atmospheric grids Other spatial transformations:flux correction,merging,etc.General algebraic operations Computing in Atmospheric Sciences Workshop 201917PRISM project first coupler:OPRISM project final couplerHigher resolution,parallel and scalable modelsHigher coupling frequencies desirable Higher number of models and(3D)coupling fieldsNeed to optimise and parallelise the couplerThe final PRISM coupler will be composed of:a Drivera Transformera new PRISM System Model Interface LibraryComputing in Atmospheric Sciences Workshop 201918PRISM project final couplerHigPRISM project final couplerFinal coupler configuration(XML files):The user chooses the models through the GUI.Each component model comes with:an Application Description(AD)a Potential Model Input and Output Description(PMIOD).The user configures his particular coupled run through the GUI:total run time,etc.for each field described in the PMIOD:source or targetcoupling or I/O status,coupling or I/O periodtransformations/interpolations,etc.Based on the users choice,the GUI produces the XML configuration files.At run-time the Driver reads and distributes configuration information.The PSMILes and Transformer act accordingly to the users specifications.Computing in Atmospheric Sciences Workshop 201919PRISM project final couplerFinPRISM project final couplerFinal coupler communication:More elaborate PSMILe based on MPI1 or MPI2(grid definition transferred through the PSMILe API)Modularity as for Oasis3:at each model time step,exchange is performed or not depending on users specifications.As for Oasis3,automatic time integrationAs for Oasis3,I/O functionality(automatic switch between coupled and forced mode)Parallel communication:as for Oasis3+repartitioning.OB OB OBCCC O1 O1CCParallel calculation of interpolation weights and addresses in the source PSMILeExtraction of useful part of source field only.Computing in Atmospheric Sciences Workshop 201920PRISM project final couplerFinFinal coupler interpolations/transformations =as for Oasis3+Support of vector fieldsSupport of 3D fieldsMore flexibility for field combination/merging,etc.The PRISM project final couplerComputing in Atmospheric Sciences Workshop 201921The PRISM project final coupleConclusionsENES and PRISMPRISM first coupler:Oasis3,now available PRISM final coupler prototype due 11/2019PRISM final coupler due 12/2019 and after PRISM?Follow-on project re-submitted at next EU-call in 2019(CAPRI rejected)International interaction and collaboration essential in all cases!enes.org;http/cerfacs.fr/PRISM/prism.htmlvalckecerfacs.fr Computing in Atmospheric Sciences Workshop 201922ConclusionsComputing in AtmospPRISM project first coupler:Oasis3Oasis3 communication;PSMILe API:Initialization:call prism_init_comp()Retrieval of component model local communicatorcall prism_get_localcomm()Coupling or I/O field declarations(name,type,shape,local partition,)call prism_def_var(field_idx,)End of definitioncall prism_enddef()In model time stepping loop,coupling or I/O field exchangecall prism_put(field_id1,time,field_array1,ierror),call prism_get(field_id2,time,field_array2,ierror)=Automatic averaging/accumulation,coupling exchange,and/or I/O depending on time argument and users specifications in namcoupleTermination:call prism_terminate()Computing in Atmospheric Sciences Workshop 201923PRISM project first coupler:OPRISM project final couplerFinal coupler communication;PSMILe API:As for Oasis3 PSMILe+Definition of grid(1D,2D,3D)call prism_def_grid()call prism_set_corners()call prism_set_mask()Definition of grid for vector and bundle fieldscall prism_set_vector()Call prism_set_subgrid()Coupling or I/O field declarations support vector,bundles,1D,2D and 3D fieldsExtraction of SCC and SMIOC information:call prism_get_persist()Computing in Atmospheric Sciences Workshop 201924PRISM project final couplerFinATM SMIOC V1:from OCE,T1V2:to OCE,T2V3:to LANDuserLAND SMIOC V3:from ATMV4:from fileV4userOCE SMIOCV1:to ATM,T1V2:from ATM,T2userDriverTOCEATMLANDfileV4Definition PhaseOCE PMIODV1:out,metadata V1V2:in,metadata V2ATM PMIODV1:in,metadata V1V2:out,metadata V2V3:out,metadata V3LAND PMIODV3:in,metadata V3V4:in,metadata V4LAND ADATM ADOCE ADV2V1V1V3V4Deployment PhaseV2userComposition PhaseuserSCC ATM:.OCE:.LAND:.userComputing in Atmospheric Sciences Workshop 201925ATM SMIOC userLAND SMIOC userODeployment PhaseV3V6V2V2V1V1V4V5V5V7V7Mj SMIOC V1:cf SCCV4:cf SCCuserMk SMIOC V4:cf SCCV5:in,fileV5,TnlV5kuserMi SMIOCV1:cf SCCV2:cf SCCV3:in,fileV3,TliuserComposition PhaseSMIOC:Specific Model Input and Output Config.userSCC V1:Mi-Mj,Tli,TnlijV2:Mi-Mj,Tij(+V6)V4:Mj-MkuseruserSCC:Specific Coupling ConfigurationDriverTMiMjMkDefinition PhaseMi PMIODV1:out,metadata V1V2:out,metadata V2V3:in,metadata V3Mj PMIODV1:in,metadata V1V4:out,metadata V4Mk PMIODV4:in,metadata V4V5:in,metadata V5fileV6fileV3fileV5PMIOD:Potential Model Input and Output DescriptionMi:Model i T:TransformerComputing in Atmospheric Sciences Workshop 201926Deployment PhaseV3V6V2V2V1V1Coupling infrastructureSupporting softwareScientificcodeRunning environmentSoftware structure of an Earth System ModelShareComputing in Atmospheric Sciences Workshop 201927Coupling Supporting softwareScOn going PRISM/ESMF collaborationCoupling infrastructureSupporting softwareUser codeRunning environmentPRISMESMFEarth System ModelComputing in Atmospheric Sciences Workshop 201928On going PRISM/ESMF collabor谢谢你的阅读v知识就是财富v丰富你的人生谢谢你的阅读知识就是财富 Thank you拯畏怖汾关炉烹霉躲渠早膘岸缅兰辆坐蔬光膊列板哮瞥疹傻俘源拯割宜跟三叉神经痛-治疗三叉神经痛-治疗 拯畏怖汾关炉烹霉躲渠早膘岸缅兰辆坐蔬光膊列板哮瞥疹
展开阅读全文