1 / 49

DieCAST

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

norman-baird
Download Presentation

DieCAST

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. DieCAST V 1.0 March 3rd 2009

  2. outline • What is DieCAST • Features • Roadmap • Questions

  3. 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

  4. 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

  5. architecture • Physical • Mathematical • Ghost zone • Technical

  6. Numerical structure Ghost zone: I=1 I=I0 J=1 J=J0 Model zone boundary shelf Cell value dynamics control

  7. DieCAST code structure

  8. 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

  9. Data & Directory structure

  10. Z-coordinate • Configuration • K0 • ZTOP • TLZ • Coordinate • Z(K0+K1) • K1=K0-1

  11. 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:

  12. Z coordinate-linear exponential

  13. Z coordinate-stretched

  14. Z coordinate-user defined

  15. XY coordinate Configuration • I0 • J0 • DYDX • DXMNUT • Y0DEG

  16. MetGen configuration-name list

  17. 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

  18. MetGen-Product The same on the same latitude Cell is empty

  19. 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

  20. Bathymetry data-Etopo5

  21. Interpolate climatological data to grid points

  22. Inmet-product Unique location In the world Cell is empty

  23. 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

  24. Winds-Hellerman coordinate • I=0.5, 359.5,1 • J=-89.5,89.5,1 • Grid winds (“winds”) • TAUX(360,180) • TAUY(360,180)

  25. Grid S T Winds Just for refference!!

  26. Prep-Procedure Configuration Manipulating DEPTH CALL INITFS CALL BOUNDS

  27. Prep_name list • IRUN, DAODT, DM0, DE0, • DMZ0, RZMX, TAU, • TAUN, DRAG, FLTW • KTRM • Flags • LOPEN, LMOVI, MVFREQ, LEVP • LRSTRT

  28. 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

  29. 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

  30. 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”

  31. Prep-2 Uniquely exist In the world Cell is filled with Initialized S,T Static domain

  32. 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

  33. Prep-summary (continued) • Input • Controlling condition • Bathymetric data • Climatological S,T,Wind • Output • “YZGRID”,RUNDATA,BOUNDARIES,EVP,ANUALCYC

  34. DieCAST

  35. DieCAST-INITIALIZATION • Preconfiguration • LRSTRT,MPASS • READ • “RUNDATA”,ZKB,WINDS,EVP,BOUNDARIES,ANUALCYC • WINDMIX, • KBI1J

  36. A & C grid

  37. DieCAST-INITFSNPB/TAI • READ • “RUNDATA”,ZKB,WINDS,EVP,BOUNDARIES,ANUALCYC • WINDMIX, • KBI1J • PROCESS INPUT DATA • EVP ELLIPTIC SOLVER

  38. DieCAST-FSNPB/TAI • CALCULATE LATERAL DRAG • HYDROSTATIC EUATIONS • TOP HALF LAYER APPROXIMATION • CALCULATE INTERNAL FLUXES

  39. 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

  40. 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

  41. DieCAST-3

  42. 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)

  43. DieCAST-FSTAI • INITIALIZING I0 ZONE: COUPLING • CALL FSTAI • CHECK AND DIAGNOSTICS • PASS INTERFACE DATA TO NPB • OUTPUT

  44. DieCAST-FSTAI • INITIALIZING I0 ZONE: COUPLING • The same process as in NPB • PASS INTERFACE DATA TO NPB

  45. I/O structure

  46. Postprecessor • Software • Product • I/O format

  47. Namelist example • &NML1 • DAODT=96., • MXIT=552960, • / • &FLAG • LRSTRT=1, • LMOVI=1, • LEVP=1, • / • &TIME_STEP • IRUN=50, • DAODT=12 • /

  48. case • Global • Norpac1 • norpac2

  49. run • Edit source file • Compile • run

More Related