1 / 15

TIEGCM Code Structure Contents

TIEGCM Code Structure Contents. TIEGCM Main Init Model initialization Advance Advance model in time Dynamics Dynamics and Chemistry Aurora Aurora Parameterization Oplus O+ Minor Minor Species Composition DUV Neutral Winds DT Neutral Temperature

lolita
Download Presentation

TIEGCM Code Structure Contents

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. TIEGCM Code StructureContents • TIEGCM Main • Init Model initialization • AdvanceAdvance model in time • DynamicsDynamics and Chemistry • AuroraAurora Parameterization • Oplus O+ • MinorMinor Species Composition • DUV Neutral Winds • DT Neutral Temperature • Major Major Species Composition • Dynamo Electrodynamics • Outhist Write to output files • Filter Longitudinal filtering

  2. Contents TIEGCM Main init Model initialization input Read user input readsource Read source history nc_rdhist rdsource apxmall apxparm Call apex module apxmka advance Advance model in time

  3. Contents init is called by tgcm init Model initialization init_cons(cons.F) Set constants rd_soldata(soldata.F) Read SEE flux data consdyn (dynamo constants) setfft(util.F) Initialize fft init_4d (fields.F) Allocate 4d fields init_3d (fields.F) Allocate 3d fields allocdata allocate and initialize msc dynamic memory init_fsech(fields.F) Initialize secondary history fields

  4. advance is called by tgcm Contents Advance (advance.F) Advance the model in time sunloc (locate sun’s longitude) getgpi (read GPI data) getimf (read IMF data ) getgswm (read GSWM data) aurora_cons (update auroral constants) init_sflux (update solar flux data) ssflux efield (calculate electric field) Time-step loop addiag (calculate geopotential Z,ZG, mbar, etc.) lbc_gswm_addiag hdif1,hdif2 (horizontal diffusion) lsqdsq weimer05 (Weimer potential model) heelis (Heelis potential model) dynamics (dynamics and chemistry module) magpres_grav (magnetic pressure and gravity prep_dynamo (prepare for dynamo) dynamo (electrodynamics module) outhist (write to netCDF output files)

  5. Contents dynamics is called by advance dynamics Dynamics and chemistry driver swdot (omega) chemrates_tdep (chemical rates) bndcmp (composition lower boundary) ionvel (drift velocities) Minor species composition filter cpktkm (specific heat, molecular visc) chapman (line integrals) qinite (night-time ionization) qjnno, qjion (ion heating and O2 dissociation) filter qjoule_ti (ion joule heating) settei (ion and electron temperature) trsolv qrj (ionization and heating) newton, newto3p (implicit and explicit cooling) aurora (auroral parameterization) hdif3 (horizontal diffusion, stage 3) lsqdsq oplus (update O+) qjoule_tn (joule heating for TN) elden (electron density) trsolv dt (neutral temperature TN) trsolv lamdas (ion drag coefficients) filter duv (neutral winds U, V) filter vquart comp (major composition O2, O)

  6. Contents aurora is called by dynamics aurora Aurora parameterization aurora_cusp cusp aurora_heat heating aurora_ions ions all in aurora.F aion auroral electrons bion solar protons

  7. Contents oplus is called by dynamics oplus Update O+ oplus_flux O+ number flux divb divergence rrk Diffusion coefficients diffus (d/(h*dz)*tp+m*g/r)*en bdotdh (b(h)*del(h))*phi bdzdvb (bz*d/(h*dz)+divb)*phi trsolv Tridiagonal solver filter_op filter

  8. Contents dynamo is called by advance dynamo Calculate Electric Potential transf Field line integrals edges Polar boundary conditions rhspde Right-hand side of PDE divide Divide stencil calrhs_jrr J_rR contribution to dynamo mud Multi-grid PDE solver clearcee Init stencil arrays threed 3d electric potential stencmd 3x3 stencils mag2geo and geo2mag Grid transformation utilities stencil Modified stencil

  9. minor species composition routines called by dynamics. Contents comp_n2d (n2d) comp_n4s (n4s) comp_no (no) minor_n4s minor_no minor minor filter filter

  10. Contents minor Called by minor composition routines Minor (advance minor species) advec (horizontal advecion) trsolv Tridiagonal solver filter_minor filter

  11. Contents Major species composition (comp_o2o and comp are called by dynamics) comp_o2o (sources and sinks for O2, O) comp (advance O2, O) advecl (horizontal advection) smooth (Shapiro smoother) filter_o2o filter

  12. Contents duv is called by dynamics duv Neutral winds U, V advec (horizontal advection) glp (horizontal pressure forcing) dldp lat,lon derivatives smooth (Shapiro smoother) lbc_gswm_duv (GSWM data) blktri (block tridiagonal solver) filter_uv filter

  13. outhist is called by advance Contents output Write to netcdf output files output_hist Write to primary history file output_sechist Write to secondary history file copyhist Copy source history def_fsech Define secondary history define_hist Define history structure nc_define Define netCDF file nc_wrhist Write history to netCDF file savefile Save history file to MSS

  14. Contents • filter: Filter routines are called by: • comp Major composition • dt Neutral temperature • duv Neutral winds • minorMinor composition • oplus O+ • swdot Vertical motion filter (Called by dt, duv, swdot, oplus) filter2 (Called by comp,minor,oplus) fftrans (util.F) fft999 Fourier transform

  15. Contents dt is called by dynamics dt Neutral Temperature TN lbc_gswm_dt (GSWM data) advec (horizontal advection) advecv (vertical advection) smooth (Shapiro smoother) tsolv (tridiagonal solver) filter_tn filter

More Related