490 likes | 636 Vues
DieCAST. V 1.0 March 3 rd 2009. outline. What is DieCAST Features Roadmap Questions. Features. F77, serial Multiple grids coupling (example in NORPAC2) Open/closed boundary options
E N D
DieCAST V 1.0 March 3rd 2009
outline • What is DieCAST • Features • Roadmap • Questions
Features • F77, serial • Multiple grids coupling (example in NORPAC2) • Open/closed boundary options • Turbulence parameterization options (in realwinds.f)Non-hydrostatic and hydrostatic dynamics (see mon.f) • Bio-geochemistry dynamics (in realwinds_bgc.f under • Immersed boundary methods for topography (see mon.f) • Lagrangian tracers releases capacity (example in NORPAC2) • hindcast wind forcing capability (in REALWINDS.f) • NetCDF and customized I/O interface • Ensemble simulation capability
Roadmap F90,MPI option configuration customized I/O Ensemble simulation F77, serial Multiple grids coupling Open/closed boundary options Turbulence parameterization Non-hydrostatic/hydrostatic Bio-geochemistry dynamics Immersed boundary methods Lagrangian tracers releases hindcast wind forcing
architecture • Physical • Mathematical • Ghost zone • Technical
Numerical structure Ghost zone: I=1 I=I0 J=1 J=J0 Model zone boundary shelf Cell value dynamics control
Code structure • Metgen: create x,y,z level coordinate • Inmet: • interpolate bathymetry data onto grid points • Indata • Interpolate S,T climatological data onto grid points • Prep • Manipulate S,T,Depth yielding to specified control condition, such as boundary condition. • DieCAST • Calculate internal U,V,W,S,T fluxes, • Output graphic data • Postprocessor • Create graph according to model product
Z-coordinate • Configuration • K0 • ZTOP • TLZ • Coordinate • Z(K0+K1) • K1=K0-1
Z-coordinate • Coordinate Types: • Linear exponential: • Z=a+b*exp(c*zz)+d*zz, 0 < zz < tlz • zz=0:Z=0; zz=tlz:Z=tlz • Control factor: c,d • Stretched: • DZ(1)=ZTOP • DZ=GX*DZ • User defined:
XY coordinate Configuration • I0 • J0 • DYDX • DXMNUT • Y0DEG
MetGen-Product domain size I0,J0,K0 domain location WESTDEG,Y0DEG cell size DXMNUT,DX,DY,DZ(z(k+1)-z(k)) cell location YDEG, YVDEG,Z(k) Storage YZGRID
MetGen-Product The same on the same latitude Cell is empty
Inmet-Bathymetry data • Bathy data format (current version) • Etopo5 • Resolution: 5 min, METERS(4320,2160) • Format: binary, netcdf • Interpolated grid data • RMETS(I0,J0) • Available bathymetric data • Etopo2,Etopo1 • http://www.ngdc.noaa.gov/mgg/global/global.html
Inmet-product Unique location In the world Cell is empty
Indata-S,T,Wind • Levitus • 360*180*33 • TLEV, SLEV • Seasonal data • TLEV(1,1) : • (0.5 deg E, -89.5 deg) • TLEV(360,180,) • (359.5 deg E, 89.5 deg ) • Grid “anualevitus” • S1,T1
Winds-Hellerman coordinate • I=0.5, 359.5,1 • J=-89.5,89.5,1 • Grid winds (“winds”) • TAUX(360,180) • TAUY(360,180)
Grid S T Winds Just for refference!!
Prep-Procedure Configuration Manipulating DEPTH CALL INITFS CALL BOUNDS
Prep_name list • IRUN, DAODT, DM0, DE0, • DMZ0, RZMX, TAU, • TAUN, DRAG, FLTW • KTRM • Flags • LOPEN, LMOVI, MVFREQ, LEVP • LRSTRT
Prep_depth • Read “depth” • At least 2 layer deep limit • Artificial shelf (NPB) • Avoid spurious vortex stretching at lat boundaries • Coupling option: depth=0 at n&s bound of overlap • Generate logical depth: KB(I,J) • Land/sea mask • IN=0,1,2 • Specify definite computational zone
Prep_initfs • Coupling option: (TAI) • Find easternmost and westernmost wet points • Set initial vertical diffusivities • EV,HV • Optional (NPB): HBK,VBK (WINDS) • Preprocessor for elliptic solver • Used for incompressibility
PREP-bounds • Read: “anualevitus” • S,T: • check for unstable points (Dan Wright scheme) • Create “anualcyc” • Seasonal boundary conditions ( “boundaries”) • TAVG,SAVG,SSURF,TSURF,TXYCLI,SXYTLI • “Anualevitus” “ anualcyc “ boundaries”
Prep-2 Uniquely exist In the world Cell is filled with Initialized S,T Static domain
Prep-summary • 1 separate set of preprocessing codes (metgen,inmet,indat,prep) for each domain • Input • bathymetry climatological data • Output of preprocessors • Domain structure • Initialized computational cell • Domain boundary
Prep-summary (continued) • Input • Controlling condition • Bathymetric data • Climatological S,T,Wind • Output • “YZGRID”,RUNDATA,BOUNDARIES,EVP,ANUALCYC
DieCAST-INITIALIZATION • Preconfiguration • LRSTRT,MPASS • READ • “RUNDATA”,ZKB,WINDS,EVP,BOUNDARIES,ANUALCYC • WINDMIX, • KBI1J
DieCAST-INITFSNPB/TAI • READ • “RUNDATA”,ZKB,WINDS,EVP,BOUNDARIES,ANUALCYC • WINDMIX, • KBI1J • PROCESS INPUT DATA • EVP ELLIPTIC SOLVER
DieCAST-FSNPB/TAI • CALCULATE LATERAL DRAG • HYDROSTATIC EUATIONS • TOP HALF LAYER APPROXIMATION • CALCULATE INTERNAL FLUXES
DieCAST-FS- INTERNAL FLUXES • 4TH-ORDER PRESSURE • VERTICAL FLUXES • LONGITUDINAL FLUXES • LATITUDINAL FLUXES • CONSERVATION EUATIONS • SLINITY • SERFACE EFFECT • OPEN BOUNDARY • SURFACE RESTORING • SURFACE FLUXES • BOTTON STRESS • A-C GRID INTERPOLATION • CREATE NON-DIVERGENT U,V,W • CHECK IN COMPRESSIBILITY • UPDATE USING FLTW METHOD
CHECK AND DIAGNOSTICS • CHECK FOR INSTABILITY • MEAN DIAGNOSTICS • HEAT • SURFACE PRESSURE • INTEGRATE VELOCITY • INTERPOLATE STREAMFUNCTION • TIME MEAN STREAMFUNCTION • MEAN U,V • MEAN S,T, LATITUDINAL VOLUME
DieCAST-COUPLING_NPB • NPB COUPLING TO TAI (U,V,S,T)_TAI(I0,J,K) =(U2,V2,S2,T2)_NPB(3,(J-2)/2+127,K)
DieCAST-FSTAI • INITIALIZING I0 ZONE: COUPLING • CALL FSTAI • CHECK AND DIAGNOSTICS • PASS INTERFACE DATA TO NPB • OUTPUT
DieCAST-FSTAI • INITIALIZING I0 ZONE: COUPLING • The same process as in NPB • PASS INTERFACE DATA TO NPB
Postprecessor • Software • Product • I/O format
Namelist example • &NML1 • DAODT=96., • MXIT=552960, • / • &FLAG • LRSTRT=1, • LMOVI=1, • LEVP=1, • / • &TIME_STEP • IRUN=50, • DAODT=12 • /
case • Global • Norpac1 • norpac2
run • Edit source file • Compile • run