GSI Fundamentals (3): Diagnostics
|
|
- Adam Jayson Briggs
- 5 years ago
- Views:
Transcription
1 24 GSI Community Tutorial July 4-6, 24. Boulder, CO GSI Fundamentals (3): Diagnostics Kathryn Newman Ming Hu & Chunhua Zhou Developmental Testbed Center (DTC)
2 Outline GSI fundamentals (): Setup and Compilation GSI fundamentals (2): Run and Namelist GSI fundamentals (3): Diagnostics Standard output Observation innovation statistics and diagnostic files Convergence information Conventional Observation Errors Background Error Covariance (BE) Analysis increments GSI fundamentals (4): Applications ü This talk is tailored to Chapter 4 of the GSI User s Guide for Community Release v3.3
3 Standard Output (stdout) Details in User s Guide Section 4. ü Critical information about the GSI analysis can be obtained ü Users can check:. Did GSI successfully complete? 2. Do optimal iterations look correct? 3. Are the background and analysis fields reasonable? ü Helpful in understanding where and why GSI may have failed
4 stdout: structure The structure of stdout is as follows:. Read in data and prepare analysis: Read in configuration (namelist) Read in background Read in observations Partition domain and data for parallel analysis Read in constant fields (fixed files) 2. Optimal iteration (analysis) 3. Save analysis result *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. PROGRAM GSI_ANL HAS BEGUN. COMPILED ORG: NP23 STARTING DATE-TIME JUN 4,24 6:58: WEN Indicates the start of the GSI analysis
5 Stdout: read in configuration state_vectors*init_anasv: 2D-STATE VARIABLES ps sst state_vectors*init_anasv: 3D-STATE VARIABLES u v tv tsen q oz cw p3d control_vectors*init_anacv: 2D-CONTROL VARIABLES ARE ps sst control_vectors*init_anacv: 3D-CONTROL VARIABLES ARE sf vp t q oz cw Read in: General set-up for GSI 3D-var analysis ü anavinfo ü State and control variables ü 4D-var setup (for 4dvar interface) ü namelist &SETUP GENCODE = 78., FACTQMIN =.E+, FACTQMAX =.E+, FACTV =., DELTIM = 2., DTPHYS = 36., BIASCOR = -., BCOPTION =, DIURNALBC =.E+, NDAT = 87, NITER =, 2*, 48*, NITER_NO_QC = 5*, MITER = 2, QOPTION = 2, NHR_ASSIMILATION = 3, MIN_OFFSET = 8, IOUT_ITER = 22, / &GRIDOPTS &BKGERR &ANBKGERR &JCOPTS &STRONGOPTS &OBSQC &SUPEROB_RADAR &LAG_DATA &HYBRID_ENSEMBLE &RAPIDREFRESH_CLDSURF &CHEM Details covered in GSI fundamentals (2)
6 stdout: check background input dh = 3 iy,m,d,h,m,s= 2 nlon,lat,sig_regional= rmse_var=t ordering=xyz WrfType,WRF_REAL= 4 4 ndim= 3 Check that bk fields have been staggering= N/A read in successfully start_index= end_index= k,max,min,mid T= k,max,min,mid T= k,max,min,mid T= k,max,min,mid T= k,max,min,mid T= k,max,min,mid T= K Maximum Minimum Repeats for all fields at each vertical level: P_TOP, ation m r o f n i e Tim field sion of bk Dimen Variable name in netcdf file Central grid ZNU, ZNW, RDX, RDY, MAPFAC_M, XLAT, XLONG, MUB, MU, PHB, QVAPOR, U, V, LANDMASK, SEAICE, SST,IVGTYP, ISLTYP, VEGFRA, SNOW, U, V, SMOIS, TSLB, TSK
7 stdout: Check domain and convinfo Domain partition information: Task ID general_deter_subdomain: general_deter_subdomain: general_deter_subdomain: general_deter_subdomain: task,istart,jstart,ilat,jlon= task,istart,jstart,ilat,jlon= task,istart,jstart,ilat,jlon= task,istart,jstart,ilat,jlon= 2 3 Ob type Is Example using 4 processors (j-num) Js Read in convinfo: Js j-num i-num Is (i-num) usage READ_CONVINFO: tcp READ_CONVINFO: ps E-3.. READ_CONVINFO: ps E-3.. READ_CONVINFO: t 2 3..E-5.. READ_CONVINFO: gps assimilate monitor
8 stdout: observation ingest GSI resets file status depending on observation time & checks for consistency with usage in satinfo files read_obs_check: read_obs_check: read_obs_check: read_obs_check: read_obs_check: read_obs_check: bufr bufr bufr bufr bufr bufr file file file file file file date rw date date date is is is is data type hirs2_n4 data type hirs3_n prepbufr ps Not available at ob time not available satwnd not available radarbufr prepbufr t Ob time match prepbufr q hirs3bufr hirs3 n7 Not used in satinfo file not used in info file -- do not read file hirs2bufr not used in info file -- do not read file hirs3bufr List of observations to be read in and processors used to read them: READ_OBS: READ_OBS: READ_OBS: READ_OBS: read read read read hirs3 mhs ps t hirs3_n7 mhs_n8 ps t using using using using ntasks= ntasks= ntasks= ntasks= Loop over all data files to be read in observations, show GPS obs outside time window: READ_GPS: time outside window skip this report READ_PREPBUFR: file=prepbufr type= ithin= rmesh=2. isfcalc= ndata= READ_BUFRTOVS: file=amsuabufr type=amsua ithin= 2 rmesh= 6. isfcalc= ndata= sis= ntask= sis=amsua_n ntask= nread= 993 nread= Data file to read Number of observations in file GPS observation outside time window Grid for data thinning Data in analysis domain and time-window
9 stdout: Check Observations Input Observation distribution in an analysis using 4 processors OBS_PARA: OBS_PARA: OBS_PARA: OBS_PARA: OBS_PARA: OBS_PARA: OBS_PARA: OBS_PARA: OBS_PARA: OBS_PARA: OBS_PARA: OBS_PARA: OBS_PARA: OBS_PARA: OBS_PARA: OBS_PARA: OBS_PARA: OBS_PARA: ps t q pw sst gps_ref hirs3 hirs4 amsua amsua amsua amsub mhs mhs hirs4 amsua mhs n7 metop-a n5 n8 metop-a n7 n8 metop-a n9 n9 n Conventional obs fairly evenly distributed Task 2 3
10 stdout: observation innovation Set up RHS of analysis equation GLBSOI: jiter,jiterstart,jiterlast,jiterend= 2 First outer analysis loop Calculate observation innovation for each data type in first outer loop: SETUPALL:,obstype,isis,nreal,nchanl=ps ps 22 Following SETUPALL:,obstype,isis,nreal,nchanl=t t 24 output will repeat multiple SETUPALL:,obstype,isis,nreal,nchanl=gps_ref gps 6 times in SETUPALL:,obstype,isis,nreal,nchanl=amsua amsua_n5 33 stdout 5 INIT_CRTM: crtm_init() on path "./" ACCoeff ReadFile (Binary)(INFORMATION) : FILE:./amsua_n5.SpcCoeff.bin; ^M ACCoeff RELEASE.VERSION:.4 N_FOVS=3 N_CHANNELS=5 Read_ODPS_Binary(INFORMATION) : FILE:./amsua_n5.TauCoeff.bin; ^M ODPS RELEASE.VERSION: 2. N_LAYERS= N_COMPONENTS=2 N_ABSORBERS= N_CHANNELS=5 N_COEFFS=26 IRwaterCoeff_ReadFile(INFORMATION) : FILE:./Nalli.IRwater.EmisCoeff.bin; ^M IRwaterCoeff RELEASE.VERSION: 3.2 N_ANGLES= 76 N_FREQUENCIES= 2223 N_WIND_SPEEDS= SETUPRAD: write header record for amsua_n to file pe.amsua_n5_ radiance obs innovation computation followed by CRTM coefficients info (only for st outer loop)
11 stdout: Check outer loop and inner iteration Print Jo components at beginning of the inner loop: Begin Jo table outer loop Observation Type surface pressure temperature wind moisture gps radiance Jo Global End Jo table outer loop Nobs Nobs Jo E E E E E E+5 Jo E+5 Jo/n Jo/n.999 Print cost function values for each inner iteration: (st) outer loop Inner iteration GLBSOI: START pcgsoi jiter= Initial cost function = E+5 Initial gradient norm = E+3 cost,grad,step,b,step? = E E-4.E+ good cost,grad,step,b,step? = E E E- good cost,grad,step,b,step? = E E E- good cost,grad,step,b,step? = E E E- good E E E E+3
12 stdout: Check outer loop and inner iteration Print diagnostics about guess field after adding analysis increment: ================================================================================ Status Var Mean Min Max analysis U E E E+ analysis V E E E+ analysis TV E E E+2 analysis Q E-3.E E-2 analysis TSEN E E E+2 analysis OZ.E+.E+.E+ analysis DUMY.E+.E+.E+ analysis DIV.E+.E+.E+ analysis VOR.E+.E+.E+ analysis PRSL E E E+2 analysis PS E E E+2 analysis SST E E E+2 analysis radb E E E+ analysis pcpb.e+.e+.e+ analysis aftb.e+.e+.e+ ================================================================================ wwww u E E E+ wwww v E E E+ wwww tv E E E+ wwww tsen E E E+ wwww q E E E-3 wwww oz.e+.e+.e+ wwww cw.e+.e+.e+ wwww p3d E E E- wwww ps E E E- wwww sst E E E- Print diagnostics about analysis increment: GLBSOI: jiter,jiterstart,jiterlast,jiterend= 2 2 Repeat for 2nd outer loop
13 stdout: check Jo components Before st outer loop Observation Type surface pressure temperature wind moisture gps radiance Jo Global Nobs Nobs Jo Global Jo/n Jo/n.999 Jo After st outer loop surface pressure temperature wind moisture gps radiance Jo E E E E E E+5 Jo E Nobs E E E E E E+4 Jo E Jo/n Nobs E E E E E E+4 Jo E Jo/n.46 Shows fit of the analysis results to the data After 2nd outer loop surface pressure temperature wind moisture gps radiance Jo Global Outer loop set to 2: completion of outer loop is end of analysis
14 stdout: Check Analysis Result Output Save analysis result ordering=xy WrfType,WRF_REAL= 4 4 ndim= 2 Check that analysis fields look staggering= N/A reasonable start_index= end_index= k,max,min,mid T= k,max,min,mid T= k,max,min,mid T= k,max,min,mid T= rmse_var=t K Maximum Minimum Central grid QVAPOR, U, V, SEAICE, SST, TSK Variable name in netcdf file
15 Stdout: Final normal exit information ENDING DATE-TIME JUN 4,24 7:: WEN PROGRAM GSI_ANL HAS ENDED. *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. : *****************RESOURCE STATISTICS******************************* : : The total amount of wall time = : ( ) : *****************END OF RESOURCE STATISTICS*************************
16 Observation Fitting Statistics and Diagnostic Files Details in User s Guide Section 4.5 and A.2
17 Why need to check fitting statistics Data Analysis: adjust background fields based on observation data so that analysis fields fit the observation better. GSI has a series of text files to provide statistic information on how background or analysis fields fit to the observations in each outer loop (fort.2*) GSI also has a series of binary files to save diagnostic information for each observation (diag*)
18 Observation Innovation Statistics GSI User s Guide Table 4.4 Conventional observations File names are from fort.2 to fort.25 (or fit_variable.analysis_time) Each file is for one observation variable File name generated in run_gsi.ksh
19 Observation Innovation Statistics (Cont ) File name Variables in file fit of satellite radiance data, such as: fort.27 or fit_rad.analysis_time amsua_n5(, n6, n7, n8, metop-a, aqua, n9), amsub_n7, fort.28 fort.29 fort.2 fort.2 fort.22 fort.23 fort.24 fort.25 Ranges/units hirs3_n7, hirs4_n9 (, metop-a), mhs_n8(, metop-a, n9), Satellite atms_npp, cris_npp, iasi66_metop-a, airs28subset_aqua, radiance avhrr3_n8(, n9, metop-b), amsre_aqua, ssmi_f5, ssmis_f6(, f7,f8,f9,f2), sndrd?_g(, g2, g3, g4, g5), imgr_g(, g2, g3, g4, g5), seviri_m8(,m9 m) fit of pcp_ssmi, pcp_tmi fit of radar radial wind (rw) Radar radial wind fit of lidar wind (dw) fit of super wind data (srw) fit of GPS data fractional difference GPS fit of conventional sst data C Tropical cyclone central pressure Lagrangian data
20 Example: fit_p (fort.2) current fit of surface pressure data, ranges in mb pressure levels (hpa)=. 2. it obs type stype count bias ps ps ps ps all rms cpen qcpen rms cpen qcpen O-B current fit of surface pressure data, ranges in mb pressure levels (hpa)=. 2. it obs type stype count bias 3 ps ps ps ps all O-A Results from test case using 2 outer loops with inner iterations in each outer loop
21 Monitoring Rejected Data Used in Analysis Example: fit_w (fort.22) Vertical levels ptop it obs type styp pbot count bias rms cpen qcpen all count all bias All prepbufr types.93 all rms all cpen rej rej rej rej rej rej rej rej rej all all all all count bias rms cpen qcpen count bias rms cpen mon mon mon mon mon mon mon mon mon all all all all count bias rms cpen qcpen count bias rms cpen
22 Example: fit_w (fort.22) O-B O-A After st outer loop O-A After 2st outer loop ptop it obs type styp pbot count bias rms cpen ptop it obs type styp pbot count bias rms cpen ptop it obs type styp pbot count bias rms cpen
23 Example: fit_rad (fort.27) Channels in satinfo file RADINFO_READ: amsua_n5 2 amsua_n5 3 amsua_n5... Content of satinfo file jpch_rad= 268 chan= var= 3. varch_cld= 9. use= ermax= chan= 2 var= 2. varch_cld= 3.5 use= ermax= chan= 3 var= 2. varch_cld= 7. use= ermax= 4.5 b_rad= 4.5 b_rad= 4.5 b_rad=. pg_rad=. pg_rad=. pg_rad=. icld_det=. icld_det=. icld_det= RADINFO_READ: guess air mass bias correction amsua_n amsua_n it satellite instrument # read # keep # rad n8 amsua rad metop-a amsua rad aqua amsua coefficients below coefficients for all channels assim penalty qcpnlty penalty qcpnlty cpen qccpen O-B it satellite instrument 3 rad n8 amsua 3 rad metop-a amsua 3 rad aqua amsua O-A # read # keep # assim cpen Results from test case using 2 outer loops with inner iterations in each outer loop qccpen
24 Example: fit_rad (fort.27) Summary as a function of observation type sat n5... type amsua penalty qcpenalty nobs 6523 qc 2447 iland isnoice icoast ireduce qc2 qc3 qc4 qc ivarl nlgross qc6 qc7 6 6 radiance ob may include mult. channels, not all channels are used in analysis Summary as a function of channel: amsua_n5 2 2 amsua_n5 3 3 amsua_n5 4 4 amsua_n5 channel nobs nobs tossed... it satellite instrument rad n4 hirs2 rad n6 hirs3 rad n7 hirs3 rad n8 hirs4 rad metop-a hirs variance Bias (before BC) Bias (after BC) penalty O-B (w/bc) Standard dev Final summary for each observation type: # read # keep # assim 478 penalty qcpnlty cpen qccpen
25 Diagnostic files (User s Guide A.2) Files include observation departure for each obs: diag_amsua_metop-a_anl diag_amsua_metop-a_ges diag_amsub_n7_anl diag_amsub_n7_ges diag_amsua_n5_anl diag_amsua_n5_ges diag_conv_anl diag_conv_ges To get these files, turn write_diag on: write_diag()=.true.,write_diag(2)=.false.,write_diag(3)=.true., To read this binary information: Code to read these files (util/analysis_utilities/read_diag) read_diag_conv.f9 (diag_conv*) read_diag_rad.f9 (diag_amsub_n7* ) Compile:./make two executables: read_diag_conv.exe read_diag_rad.exe
26 Observation departure for each ob Run read_diag_conv.exe: (namelist.conv needed) &iosetup infilename='./diag_conv_anl.23222', outfilename='./results_conv_anl', / GSI diagnostic file Text file to save the content of diagnostic file Example content of results_conv_anl station ID 2997 : PADK : obs type obs time obs latitude obs longitude obs pressure usage obs value O-B Run read_diag_rad.exe: (namelist.rad needed) &iosetup infilename='./diag_amsua_n8_ges.23222', outfilename='./results_amsua_n5_anl', / Same as conventional
27 Convergence information Details in User s Guide Section 4.6 and A.3
28 Find convergence information From stdout (see example in stdout of this talk) From fort.22 Includes many details Cost function and gradient Contribution from background and each data type DTC provides a ksh script to filter this file util/analysis_utilities/plot_cost_grad/filter_fort22.ksh grep 'cost,grad,step,b' fort.22 sed -e 's/cost,grad,step,b,step? = //g' sed -e 's/good//g' > cost_gradient.txt outer loop e E E-3.E+ Inner iteration
29 Convergence information cost_gradient.txt outer loop Inner iteration e e e-2.e e e e e E E E E E E E E E E E E E E E E E E E E+ Penalty (cost function) Norm of gradient
30 Plot convergence information DTC provides a NCL script to make plot util/analysis_utilities/plot_cost_grad/gsi_cost_gradient.ncl load "$NCARG_ROOT/lib/ncarg/nclex/gsun/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" begin # of iterations in st outer loop nloop=5 # of iterations in 2nd outer loop nloop2=5 step = stringtofloat(systemfunc("cut c7-9./cost_gradient.txt")) cost = stringtofloat(systemfunc("cut c-35./cost_gradient.txt")) gradient = stringtofloat(systemfunc("cut c36-6./cost_gradient.txt")) titles = new(4,string) titles()="cost outer " titles()="gradient outer " titles(2)="cost outer 2" titles(3)="gradient outer 2" plot = new(4,graphic) figure file name and format xwks = gsn_open_wks("pdf","gsi_cost_gradient )
31 Check convergence: example figure Cost function Norm of gradient st outer loop with 5 inter iteration 2nd outer loop with 5 inter iteration Iteration step Iteration step
32 Conventional observation errors Details in User s Guide Section 4.7
33 Conventional Observation errors ü Each observation has its own observation errors Ø Global: oberrflg=true: obs errors generated from external obs error table oberrflg=false: read from PrepBUFR file Ø Regional: forced use of external obs error table regardless of oberrflg Observation error files in release (~/comgsi_v3.3/fix ): nam_errtable.r3dv (default) prepobs_errtable.global rtma/new_rtma_nam_errtable.r3dv Prepbufr Ob 2 OBSERVATION TYPE.E+4.267E+.5E+4.332E+.E+4.47E+.95E E+.9E E+.85E E+ P T.563E+.6326E E+.8635E E+.2E+ q.e+.e+.e+.e+.e+.e+.685e+.685e+.685e+.737e e e+.e+.e+.e+.e+.e+.e+ UV Ps Pw Observation error file easy to edit to tune error values
34 Gross check Gross check is a QC check to exclude questionable obs that degrade the analysis Users can adjust thresholds for data types in convinfo: otype ps ps t t type sub iuse twindow numgrp ngroup nmiter gross ermax ermin var_b var_pg ithin Gross check controlled by these columns ratio = (Observation-Background)/max(ermin,min(ermax,obserror)) If ratio > gross: observation rejected Rejected observations will show up in fit files under rej
35 Background Error Covariance Details in User s Guide 4.8 and Advanced User s Guide
36 Background Error Lecture tomorrow morning on BE May also use domain-specific BEs rather than pre-computed Pre-computed BEs available in v3.3 release (~/comgsi_v3.3/fix/*_endian): : regional background error statistics, computed using forecasts from the NCEP s NAM model covering North America. nam_glb_berror.f77.gcv : contains the global background errors based on the NCEP s GFS model, a global forecast model. nam_nmmstat_na.gcv ü Also included in this release package is the background error matrix for RTMA GSI: new_rtma_regional_nmm_berror.f77.gcv Be sure to use appropriate BE for domain Regional BE is latitude (not longitude) dependent can be used on any NH domain
37 anavinfo The background error covariance matrix used in the GSI analysis is the fixed berror multiplied by several factors set by the namelist and anavinfo vs scale factor for vertical correlation lengths for background error hzscl(3) scale factor for three scales specified for horizontal smoothing hswgt(3) weights to apply to each horizontal scales ü as/tsfc_sdv in the control_vector section of anavinfo are factors for tuning the variance of each analysis control variable
38 Analysis Increment Details in User s Guide Section 4.9 and A.4
39 Analysis increment: plot for ARW case util/analysis_utilities/plots_ncl/analysis_increment.ncl util/analysis_utilities/plots_ncl/gsi_singleobs_arw.ncl load "$NCARG_ROOT/lib/ncarg/nclex/gsun/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" GSI analysis begin cdf_analysis = addfile(./gsiprd_23222_prepbufr/wrf_inout.cdf","r") cdf_bk = addfile(./dtc/na3km/bk/wrfinput_d_2-3-22_2::.cdf","r") Ta = cdf_analysis->t(,:,:,:) Tb = cdf_bk->t(,:,:,:) DT = Ta Tb Analysis increment = GSI analysis First guess First guess
40 Analysis increment: plot for NMM case util/analysis_utilities/plots_ncl/gsi_singleobs_nmm.ncl util/analysis_utilities/plots_ncl/fill_nmm_grid2.ncl load "$NCARG_ROOT/lib/ncarg/nclex/gsun/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" load "./fill_nmm_grid2.ncl" begin cdf_analysis = addfile(./nmm_2722/gsiprd/wrf_inout.cdf","r") cdf_bk = addfile(./nmm_netcdf/2722/wrfinput_d_nmm_netcdf.cdf","r") Ta = cdf_analysis->t(,:,:,:) Tb = cdf_bk->t(,:,:,:) DT = Ta Tb do k=, nz- fill_nmm_grid2(dt(k,:,:),nx,ny,f3d(k,:,:),) end do Convert from E grid to A grid In: E grid field Out: A grid field
41 Analysis increment: example ΔT ΔU ΔV Δq GSI analysis increment at the 5th level
42 usage GSI took 3 m.5 s to finish *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. PROGRAM GSI_ANL HAS BEGUN. COMPILED ORG: NP23 STARTING DATE-TIME JUN 4,24 6:58: WEN ENDING DATE-TIME JUN 4,24 7:: WEN PROGRAM GSI_ANL HAS ENDED. *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. *. Additional resource stats *****************RESOURCE STATISTICS******************************* The total amount of wall time = The total amount of time in user mode = The total amount of time in sys mode = The maximum resident set size (KB) = Number of page faults without I/O activity = Number of page faults with I/O activity = Number of times filesystem performed INPUT = Number of times filesystem performed OUTPUT = Number of Voluntary Context Switches = 7723 Number of InVoluntary Context Switches = 229 *****************END OF RESOURCE STATISTICS*************************
43 Summary
44 Many methods can be used: We can diagnose GSI analysis by checking: ü Standard output ü Observation fitting statistics and diagnosis files ü Convergence information ü Analysis increment There are more: Single observation test Plot observation innovation overlay analysis increment Forecast
45 Questions?
GSI Fundamentals (5): Review and Applications
27 GSI/EnKF Community Tutorial July -4, 27, College Park, MD GSI Fundamentals (5): Review and Applications Jeff Beck, Ming Hu, Hui Shao2, Chunhua Zhou2, and Kathryn Newman2 Developmental Testbed Center
More informationGSI Fundamentals (2) Run and Namelist
2 GSI Community Tutorial June 29-July, 2, Boulder, CO GSI Fundamentals (2) Run and Namelist Ming Hu and Hui Shao Developmental Testbed Center Outlines GSI fundamentals (): Setup and Compilation GSI fundamentals
More informationGSI fundamentals (4): Background Error Covariance and Observation Error
27 GSI Community Tutorial July 2, 27, College Park, MD GSI fundamentals (4): Background Error Covariance and Observation Error Ming Hu Developmental Testbed Center Outline GSI fundamentals (): Setup and
More informationEnKF Fundamentals (2b): Applications
2015 EnKF Community Tutorial August 13-14, 2015. Boulder, CO EnKF Fundamentals (2b): Applications Kathryn Newman Ming Hu, and Chunhua Zhou Developmental Testbed Center (DTC) Outline EnKF fundamentals (1):
More informationEnKF Fundamentals (1): Configuration and Run
2017 Joint DTC-EMC-JCSDA GSI-EnKF Tutorial July 11-14, 2017, NCWCP, College Park, MD EnKF Fundamentals (1): Configuration and Run Chunhua Zhou * Kathryn Newman * and Ming Hu ** Developmental Testbed Center
More informationGSI Fundamentals (2) Run and Namelist
24 GSI Community Tutorial July 4-6, 24, NCAR, Boulder GSI Fundamentals (2) Run and Namelist Ming Hu and Hui Shao Developmental Testbed Center Outlines GSI fundamentals (): Setup and Compilation GSI fundamentals
More informationGSI Setup, Run and Namelist
GSI Setup, Run and Namelist Hui Shao GSI Community Tutorial, June 28-30, 2010, Boulder, CO Outline Installing GSI Running GSI Basic runtime options This talk is tailored based on GSI Community Release
More informationGSI Setup, Run and Namelist
GSI Setup, Run and Namelist Hui Shao GSI Community Tutorial, June 28-30, 2010, Boulder, CO Observation Error Observation Observation processing and assimilation PrepBUFR and BUFR processing: 06/29, Tue
More informationGSI Software Design. Donald Stark National Center for Atmospheric Research (NCAR) The Developmental Testbed Center (DTC) 24 July, 2010
GSI Software Design Donald Stark National Center for Atmospheric Research (NCAR) The Developmental Testbed Center (DTC) 24 July, 2010 1 Outline Tour of the directory structure The build system The rest
More informationWRF Post-Processing and Visualization
NCAR Earth System Laboratory National Center for Atmospheric Research NCAR is Sponsored by NSF and this work is partially supported by the Willis Research Network and the Research Partnership to Secure
More informationSyed RH Rizvi.
Community Tools: gen_be Syed RH Rizvi National Center For Atmospheric Research NCAR/ESSL/MMM, Boulder, CO-80307, USA rizvi@ucar.edu 0 Talk overview What is gen_be? How it works? Some technical details
More informationEnsemble Kalman Filter
NOAA Ensemble Kalman Filter Beta Release v1.0 compatible with GSI community release v3.3 User s Guide January 2015 Developmental Testbed Center National Center for Atmospheric Research National Centers
More informationThe inclusion of cloudy radiances in the NCEP GSI analysis system
he inclusion of cloudy radiances in the NCEP GSI analysis system 1,2,3 Min-Jeong Kim 1 Fuzhong Weng and 2 John Derber 1.NOAA/NESDIS/SAR 2. Joint Center for Satellite Data Assimilation (JCSDA) 3. CIRA,
More informationInstituting an observation database (ODB) capability in the GSI
Instituting an observation database (ODB) capability in the GSI Jeff Whitaker, Scott Gregory, and Tom Hamill NOAA / ESRL Physical Sciences Division Presentation to Blueprints for Next-Generation Data Assimilation
More informationInstituting an observation database
Instituting an observation database capability in the NCEP GSI Tom Hamill, Jeff Whitaker, Scott Gregory NOAA / ESRL Physical Sciences Division Presentation to DAOS, Exeter England, April 2016 Our intent:
More informationGeneration of WRF-ARW Background Errors (BE) For GSI
Generation of WRF-ARW Background Errors (BE) For GSI 1. Introduction: Syed RH Rizvi, Zhiquan Liu and Xiang-Yu Huang NCAR/ESSL/MMM, Boulder, CO Currently NCEP global background errors are used for running
More information2011 BUFR Tutorial, 13 November 2011, Boulder, CO. GSI BUFR Interface. Ming Hu** and Ruifang Li* *NCAR/MMM **Developmental Testbed Center
2011 BUFR Tutorial, 13 November 2011, Boulder, CO GSI BUFR Interface Ming Hu** and Ruifang Li* *NCAR/MMM **Developmental Testbed Center 1 Data processing procedure in GSI Observation ingestion and processing
More informationUser s Guide Version 1.2
User s Guide Version 1.2 - Compatible with GSI community release v3.6 September 2017 Hui Liu National Center for Atmospheric Research (NCAR) Ming Hu National Oceanic and Atmospheric Administration (NOAA)/Earth
More informationWRF-NMM Standard Initialization (SI) Matthew Pyle 8 August 2006
WRF-NMM Standard Initialization (SI) Matthew Pyle 8 August 2006 1 Outline Overview of the WRF-NMM Standard Initialization (SI) package. More detailed look at individual SI program components. SI software
More informationCOMMUNITY VERSION 3.3. User s Guide. June Developmental Testbed Center
COMMUNITY VERSION 3.3 User s Guide June 2014 Developmental Testbed Center National Center for Atmospheric Research National Centers for Environmental Prediction, NOAA Global Systems Division, Earth System
More informationUser s Guide Version 1.3
User s Guide Version 1.3 - Compatible with GSI community release v3.7 November 2018 Hui Liu National Center for Atmospheric Research (NCAR) Ming Hu, Guoqing Ge National Oceanic and Atmospheric Administration
More informationGridpoint Statistical Interpolation (GSI) Version 3.0 User s Guide
Gridpoint Statistical Interpolation (GSI) Version 3.0 User s Guide Developmental Testbed Center National Center for Atmospheric Research National Centers for Environmental Prediction, NOAA Global Systems
More informationCOMMUNITY VERSION 3.4. User s Guide. July Developmental Testbed Center
COMMUNITY VERSION 3.4 User s Guide July 2015 Developmental Testbed Center National Center for Atmospheric Research National Centers for Environmental Prediction, NOAA Global Systems Division, Earth System
More informationHARMONIE ATOVS data assimilation and coordinated impact study
HARMONIE ATOVS data assimilation and coordinated impact study ALADIN / HIRLAM 23rd Workshop / All-Staff Meeting Reykjavik, Iceland, 15-19 April, 2013 Magnus Lindskog, Mats Dahlbom, Sigurdur Thorsteinsson,
More informationBias correction of satellite data at ECMWF
Bias correction of satellite data at ECMWF Thomas Auligne Dick Dee, Graeme Kelly, Tony McNally ITSC XV - October 2006 Slide 1 Motivation for an adaptive system Simplify the bias correction process of manual
More informationGSI Fundamentals (1): Setup and Compilation
2012 GSI Summer Tutorial, Boulder, CO GSI Fundamentals (1): Setup and Compilation Donald Stark Na-onal Center for Atmospheric Research (NCAR) The Developmental Testbed Center (DTC) Wednesday 21 August,
More informationAbout the SPEEDY model (from Miyoshi PhD Thesis):
SPEEDY EXPERIMENTS. About the SPEEDY model (from Miyoshi PhD Thesis): The SPEEDY model (Molteni 2003) is a recently developed atmospheric general circulation model (AGCM) with a spectral primitive-equation
More informationUser s Guide Version 3.7
User s Guide Version 3.7 November 2018 Ming Hu, Guoqing Ge National Oceanic and Atmospheric Administration (NOAA)/Earth System Research Laboratory Cooperative Institute for Research in Environmental Sciences
More informationARPS-3dvar Program. 1. Introduction
1. Introduction The ARPS variational data assimilation program 3dvar is designed to assimilate the observation data to ARPS model. The program takes an ARPS forecast as background, and surface data, multi-level
More informationObservation feedback archiving in MARS. Acknowledgement:
Observation feedback archiving in MARS Piotr Kuchta Data and Services Section, ECMWF P.Kuchta@ecmwf.int Acknowledgement: Slide 1 Manuel Fuentes, Baudouin Rault, Erik Andersson, Anne Fouilloux, Lars Isaksen,
More informationConverging Remote Sensing and Data Assimilation through Data Fusion
Converging Remote Sensing and Data Assimilation through Data Fusion Benjamin T. Johnson, Ph.D. AER @ NOAA/STAR/JCSDA benjamin.t.johnson@noaa.gov The Team: Sid Boukabara Kevin Garrett Eric Maddy Ling Liu
More informationWRF-Var System. WRF-Var Tutorial July 21-22, Xin Zhang and Syed RH Rizvi. Michael Duda
System Tutorial July 21-22, 2008 Xin Zhang and Syed RH Rizvi Michael Duda Outline in the WRF Modeling System Software Code Overview 1 in the WRF Modeling System 2 in the WRF Modeling System Blue --> Supported
More informationCollabora've Development
Collabora've Development Ricardo Todling NASA GMAO 2013 Joint DTC- EMC- JCSDA GSI Workshop 1 In the Beginning NCEP: SSI (late 1980 s) Spectral formula'on of background error cov Direct assimila'on of radiances
More information3-D Tomographic Reconstruction
Mitglied der Helmholtz-Gemeinschaft 3-D Tomographic Reconstruction of Atmospheric Trace Gas Concentrations for Infrared Limb-Imagers February 21, 2011, Nanocenter, USC Jörn Ungermann Tomographic trace
More informationNWPSAF 1D-Var User Manual
NWPSAF 1D-Var User Manual Software Version 1.1.1 Fiona Smith, Met Office, Exeter, UK NWPSAF-MO-UD-032 Revision History Document revision Date Author Notes Version 1.1 06/12/16 F. Smith Draft Version valid
More informationProgress in the assimilation of GPS-RO at ECMWF
Progress in the assimilation of GPS-RO at ECMWF Sean Healy Mike Rennie, Milan Dragosavac, Paul Poli, Carla Cardinali, Peter Bauer. Impact of GPS-RO The impact of GPS-RO NWP and climate reanalysis applications
More informationFile Formats and Pre-Processing
File Formats and Pre-Processing File Formats Pre-processing Tools Useful Links Presenter: Tara Jensen Supported File Formats Forecasts GRIB GRIdded Binary file format (*version 1) NetCDF MET specific network
More informationBUFR/PrepBUFR File Processing
2011 BUFR Tutorial, 13 December 2011, Boulder, CO BUFR/PrepBUFR File Processing Ruifang Li* and Ming Hu** *NCAR/MMM **Developmental Testbed Center 1 Topics Covered Basic actions and steps of BUFR file
More informationNEMO data assimilation and PDAF cooperation. Lars Axell & Ye Liu, SMHI
NEMO data assimilation and PDAF cooperation Lars Axell & Ye Liu, SMHI Outline NEMO data assimilation An introduction NEMO DA 3D/4D Var 3D/4D EnVar PDAF Different ensemble filters NEMO-Nordic implementation
More informationObserva(on Processing. Nancy Collins or
Observa(on Processing Nancy Collins nancy@ucar.edu or dart@ucar.edu Roadmap What s in an Observa(on Provided tools and capabili(es Observa(on sources Observa(on and representa(veness error Types, Kinds,
More informationLevenberg-Marquardt minimisation in ROPP
Ref: SAF/GRAS/METO/REP/GSR/006 Web: www.grassaf.org Date: 4 February 2008 GRAS SAF Report 06 Levenberg-Marquardt minimisation in ROPP Huw Lewis Met Office, UK Lewis:Levenberg-Marquardt in ROPP GRAS SAF
More informationWRFDA-3DVar Setup, Run, and Diagnostics
WRFDA-3DVar Setup, Run, and Diagnostics Craig Schwartz Hui-Chuan Lin, Yong-Run Guo, Syed Rizvi, Xin Zhang, and Hans Huang WRFDA Tutorial, July 2013 Outline Installing WRFDA-3DVar Setting-up a WRFDA-3DVar
More informationStat Analysis Tool. Filtering Summarizing Aggregating. of Grid-Stat, Point-Stat, & Wavelet-Stat output
Stat Analysis Tool Filtering Summarizing Aggregating of Grid-Stat, Point-Stat, & Wavelet-Stat output Presenter: Tara Jensen What can Stat Analysis do for you? Can I get Q: Overall statistics for all gridded
More informationUser s Guide Version 3.5
User s Guide Version 3.5 August 2016 Ming Hu National Oceanic and Atmospheric Administration (NOAA)/Earth System Research Laboratory Cooperative Institute for Research in Environmental Sciences (CIRES)
More informationParalleliza(on Challenges for Ensemble Data Assimila(on
Paralleliza(on Challenges for Ensemble Data Assimila(on Helen Kershaw Institute for Mathematics Applied to Geophysics, National Center for Atmospheric Research Email: hkershaw@ucar.edu What am I going
More informationSURFEX LDAS March 2012
SURFEX LDAS March 2012 Alina Barbu Combined assimilation of satellite-derived soil moisture and LAI 2 Motivation of our work GEOLAND 2 project Land Carbon Information Service (LCIS) on vegetation/land
More informationGridpoint Statistical Interpolation (GSI)
Gridpoint Statistical Interpolation (GSI) Version 3.2 User s Guide Developmental Testbed Center National Center for Atmospheric Research National Centers for Environmental Prediction, NOAA Global Systems
More informationAAPP status report and preparations for processing METOP data
AAPP status report and preparations for processing METOP data Nigel C Atkinson *, Pascal Brunel, Philippe Marguinaud and Tiphaine Labrot * Met Office, Exeter, UK Météo-France, Centre de Météorologie Spatiale,
More informationWRF Utilities. Cindy Bruyère. WRF Users' Tutorial Mesoscale & Microscale Meteorological Division / NCAR
WRF Utilities Cindy Bruyère 1 Utilities UG: Chapter 3 & 10 Graphics Designing a mdel dmain Data Input Intermediate Output grib 1&2 intermediate frmat netcdf netcdf netcdf tls Other utilities MET 2 Graphics
More informationGSI Fundamentals (1): Setup and Compilation
GSI Fundamentals (1): Setup and Compilation Donald Stark Na-onal Center for Atmospheric Research (NCAR) The Developmental Testbed Center (DTC) Tuesday 11 August, 2015 Outline GSI fundamentals (1): Setup
More informationWRF Data Assimilation System: Software and Compilation
WRF Data Assimilation System: Software and Compilation Michael Kavulich, Jr. July 24 July 26, 2017 National Center for Atmospheric Research Boulder, CO 1 WRFDA System Outline Introduction Compiling the
More informationLecture 24. Lidar Simulation
Lecture 24. Lidar Simulation q Introduction q Lidar Modeling via Lidar Simulation & Error Analysis q Functions of Lidar Simulation and Error Analysis q How to Build up Lidar Simulation? q Range-resolved
More informationFormat specification for the SMET Weather Station Meteorological Data Format version 1.1
Format specification for the SMET Weather Station Meteorological Data Format version 1.1 Mathias Bavay November 28, 2017 Abstract The goal of this data format is to ease the exchange of meteorological
More informationLecture 09. Lidar Simulation and Error Analysis Overview (1)
Lecture 09. Lidar Simulation and Error Analysis Overview (1) Introduction Lidar Modeling via Lidar Simulation & Error Analysis Functions of Lidar Simulation and Error Analysis How to Build up Lidar Simulation?
More informationMotion tracking and cloud height assignment methods for Himawari-8 AMV
Motion tracking and cloud height assignment methods for Himawari-8 AMV Kazuki Shimoji Japan Meteorological Agency / Meteorological Satellite Center 3-235, Nakakiyoto, Kiyose, Tokyo, Japan Abstract Japanese
More informationThe impact of background error on incomplete observations for 4-D Var data assimilation with the FSU GSM
The impact of background error on incomplete observations for 4-D Var data assimilation with the FSU GSM Zhuo Liu and I.M. Navon (Department of Mathematics and School of Computational Science and Information
More informationCS 490: Computer Vision Image Segmentation: Thresholding. Fall 2015 Dr. Michael J. Reale
CS 490: Computer Vision Image Segmentation: Thresholding Fall 205 Dr. Michael J. Reale FUNDAMENTALS Introduction Before we talked about edge-based segmentation Now, we will discuss a form of regionbased
More informationMetview BUFR Tutorial. Meteorological Visualisation Section Operations Department ECMWF
Meteorological Visualisation Section Operations Department ECMWF 05/03/2015 This tutorial was tested with Metview version 4.3.0 and some features might not work for previous versions. Copyright 2015 European
More informationToward assimilating radio occultation data into atmospheric models
Toward assimilating radio occultation data into atmospheric models Stig Syndergaard Institute of Atmospheric Physics The University of Arizona, Tucson, AZ Thanks to: Georg B. Larsen, Per Høeg, Martin B.
More informationUSERS GUIDE for the Community release of the GFDL Vortex Tracker
USERS GUIDE for the Community release of the GFDL Vortex Tracker November 2011 Version 3.3b The Developmental Testbed Center Shaowu Bao, NOAA/ESRL/GSD and CIRES/CU Donald Stark, NCAR/RAL/JNT Ligia Bernardet,
More informationThe impact of background error on incomplete observations for 4-D Var data assimilation with the FSU GSM
The impact of background error on incomplete observations for 4-D Var data assimilation with the FSU GSM Zhuo Liu School of Computational Science and Information Technology Florida State University, Tallahassee,
More informationCommunity Tools (1) - PrepBUFR/BUFR: Basic tools, NCEP data, and BUFR table
2012 GSI Summer Tutorial, 22 Aug. 2012, Boulder, CO Community Tools (1) - PrepBUFR/BUFR: Basic tools, NCEP data, and BUFR table Ruifang Li NCAR/MMM Topics covered NCEP observation data Observation data
More informationUnderstanding the difference between the UAH and RSS retrievals of satellite-based tropospheric temperature estimate. Carl A. Mears and Frank J.
Understanding the difference between the UAH and RSS retrievals of satellite-based tropospheric temperature estimate Workshop on Reconciling Vertical Temperature Trends NCDC, October 27-29, 2003 Carl A.
More informationOBAN Class Homework Assignment No. 4 Distributed on November 3, Due Thursday, December 1, 2016
OBAN Class Homework Assignment No. 4 Distributed on November 3, 2016 Due Thursday, December 1, 2016 Original document on ANALAB was written by Huang, Gustafsson and Robertson (2000), and can be found in
More informationTutorial 1. Linear Regression
Tutorial 1. Linear Regression January 11, 2017 1 Tutorial: Linear Regression Agenda: 1. Spyder interface 2. Linear regression running example: boston data 3. Vectorize cost function 4. Closed form solution
More informationA problem of strong scalability?
Scalability of BATOR : A problem of strong scalability? Ryad El Khatib (CNRM/GMAP) Aladin Workshop / Hirlam All Staff Meeting Norrköping, 05-08 April 2011 Plan INTRODUCTION Presentation of BATOR software
More informationCommunity Tools (1) - PrepBUFR/BUFR: Basic tools, NCEP data tank, and Obsproc
2011 GSI Summer Tutorial, 29 June 2011, Boulder, CO Community Tools (1) - PrepBUFR/BUFR: Basic tools, NCEP data tank, and Obsproc Ruifang Li*, Ming Hu**, and Yongrun Guo* *NCAR/MMM **Developmental Testbed
More informationBasic Calculator Functions
Algebra I Common Graphing Calculator Directions Name Date Throughout our course, we have used the graphing calculator to help us graph functions and perform a variety of calculations. The graphing calculator
More informationData Management in BMKG. Database Center, BMKG July 11, 2017
Data Management in BMKG Database Center, BMKG July 11, 2017 Let s think about 4 Everlasting used data The biggest data producer in Indonesia Many data resources Forecaster may false but Data can t! STATION
More informationPSL Power Standards Lab 2020 Challenger Drive Alameda, CA USA TEL FAX
Reading PQube Meters remotely using Modbus-over-TCP PSL Power Standards Lab 2020 Challenger Drive Alameda, CA 94501 USA TEL ++1-510-522-4400 FAX ++1-510-522-4455 www.powerstandards.com 1 October 2010 A
More informationIntroduction to ANSYS DesignXplorer
Lecture 4 14. 5 Release Introduction to ANSYS DesignXplorer 1 2013 ANSYS, Inc. September 27, 2013 s are functions of different nature where the output parameters are described in terms of the input parameters
More informationOcean Simulations using MPAS-Ocean
Ocean Simulations using MPAS-Ocean Mark Petersen and the MPAS-Ocean development team Los Alamos National Laboratory U N C L A S S I F I E D Slide 1 Progress on MPAS-Ocean in 2010 MPAS-Ocean is a functioning
More informationAssimilating GPS radio occultation measurements with two-dimensional bending angle observation operators
QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY Q. J. R. Meteorol. Soc. 133: 1213 1227 (2007) Published online 6 July 2007 in Wiley InterScience (www.interscience.wiley.com).63 Assimilating GPS radio
More informationBUFR/PrepBUFR User s Guide Version 1.0 Draft
BUFR/PrepBUFR User s Guide Version 1.0 Draft Developmental Testbed Center December, 2011 Foreword BUFR (Binary Universal Form for the Representation of meteorological data) is Table Driven Data Representation
More informationUnifying Verification through a Python-wrapped Suite of Tools
Support for MET is provided by NOAA, US Air Force, NSF and NCAR through the Developmental Testbed Center (DTC) Unifying Verification through a Python-wrapped Suite of Tools Tara Jensen, John Halley Gotway,
More informationNon-Newtonian Transitional Flow in an Eccentric Annulus
Tutorial 8. Non-Newtonian Transitional Flow in an Eccentric Annulus Introduction The purpose of this tutorial is to illustrate the setup and solution of a 3D, turbulent flow of a non-newtonian fluid. Turbulent
More informationIntroduction to LETKF Data Assimilation
Introduction to LETKF Data Assimilation David Kuhl and Eric Kostelich AOSC614 November 12, 2008 Introduction The goal of these exercises is to familiarize you with LETKF data assimilation. We will first
More informationRegCM-ROMS Tutorial: Introduction to ROMS Ocean Model
RegCM-ROMS Tutorial: Introduction to ROMS Ocean Model Ufuk Utku Turuncoglu ICTP (International Center for Theoretical Physics) Earth System Physics Section - Outline Outline Introduction Grid generation
More informationPrac%cal Session 3: Atmospheric Model Configura%on Op%ons. Andrew Ge>elman
Prac%cal Session 3: Atmospheric Model Configura%on Op%ons Andrew Ge>elman Overview Monday: Running the model Tuesday: namelist Control of the model Diagnosing the model Today: Different configura%on op%ons
More informationAkarsh Pokkunuru EECS Department Contractive Auto-Encoders: Explicit Invariance During Feature Extraction
Akarsh Pokkunuru EECS Department 03-16-2017 Contractive Auto-Encoders: Explicit Invariance During Feature Extraction 1 AGENDA Introduction to Auto-encoders Types of Auto-encoders Analysis of different
More informationPrinciples of Wireless Sensor Networks. Fast-Lipschitz Optimization
http://www.ee.kth.se/~carlofi/teaching/pwsn-2011/wsn_course.shtml Lecture 5 Stockholm, October 14, 2011 Fast-Lipschitz Optimization Royal Institute of Technology - KTH Stockholm, Sweden e-mail: carlofi@kth.se
More informationIntroduction to the ClimValDiagTool
Introduction to the ClimValDiagTool K. Gottschaldt & V. Eyring, 13. 2. 2013 1. General Info 2. Access miklip.dkrz.de 3. Get the code 4. Prepare data 5. Walk through an example 6. Modify the example 7.
More informationLecture 14. Lidar Data Inversion and Sensitivity Analysis
Lecture 14. Lidar Data Inversion and Sensitivity Analysis Data inversion process Nonlinearity of PMT and discriminator More considerations for Na Doppler lidar Definition of sensitivity Summary Office
More informationIntroduction to Programming in C Department of Computer Science and Engineering. Lecture No. #16 Loops: Matrix Using Nested for Loop
Introduction to Programming in C Department of Computer Science and Engineering Lecture No. #16 Loops: Matrix Using Nested for Loop In this section, we will use the, for loop to code of the matrix problem.
More informationJust on time to face new challenges with NEC super-computer at Meteo-France
Just on time to face new challenges with NEC super-computer at Meteo-France Agenda of the procurement Presentation of the first phase Installation phase (power supply, air cooling) Use of a porting machine
More informationIntroduction to MODE
Introduction to MODE Verifying with Objects Presenter: Tara Jensen copyright 2017, UCAR, all rights reserved MODE Object Identification Source: Davis 2006 (Intensity presented as vertical dimension) Smoothing
More informationMcIDAS-V Tutorial Displaying Gridded Data updated January 2016 (software version 1.5)
McIDAS-V Tutorial Displaying Gridded Data updated January 2016 (software version 1.5) McIDAS-V is a free, open source, visualization and data analysis software package that is the next generation in SSEC's
More informationGrADS for Beginners. Laura Mariotti
GrADS for Beginners Laura Mariotti mariotti@ictp.it Outline n What is GrADS and how do I get it? n GrADS essentials n Getting started n Gridded data sets n Displaying data n Script language n Saving your
More informationAn introduction to the GOCE error variance covariance products: A user s perspective. Rory Bingham Newcastle University
An introduction to the GOCE error variance covariance products: A user s perspective Rory Bingham Newcastle University rory.bingham@ncl.ac.uk The ESA GOCE Virtual Archive The variance-covariance matrices
More informationVersion 3 Updated: 10 March Distributed Oceanographic Match-up Service (DOMS) User Interface Design
Distributed Oceanographic Match-up Service (DOMS) User Interface Design Shawn R. Smith 1, Jocelyn Elya 1, Adam Stallard 1, Thomas Huang 2, Vardis Tsontos 2, Benjamin Holt 2, Steven Worley 3, Zaihua Ji
More informationLecture 22 The Generalized Lasso
Lecture 22 The Generalized Lasso 07 December 2015 Taylor B. Arnold Yale Statistics STAT 312/612 Class Notes Midterm II - Due today Problem Set 7 - Available now, please hand in by the 16th Motivation Today
More informationHow to Use GOCE Level 2 Products
How to Use GOCE Level 2 Products Thomas Gruber 1), Reiner Rummel 1), Radboud Koop 2) 1) Institute of Astronomical and Physical Geodesy, Technical University Munich 2) Netherlands Institute for Space Research
More informationIntroduction to Data Assimilation
Introduction to Data Assimilation Eric Kostelich and David Kuhl MSRI Climate Change Summer School July 21, 2008 Introduction The goal of these exercises is to familiarize you with LETKF data assimilation.
More informationCFOSAT and WindRad simulation. IOVWST Meeting Zhen Li, Ad Stoffelen
CFOSAT and WindRad simulation IOVWST Meeting Zhen Li, Ad Stoffelen 02-05-2017 1 Outline CFOSAT Rotating Fan-beam concept FY-3E WindRad dual frequency Ku/C-band L1B simulation result Beam Weighting method
More informationNEMOVAR, a data assimilation framework for NEMO
, a data assimilation framework for Arthur Vidard September, 19 th 2008 19/09/2008 1/15 Forecast is produced by integration of a model from an initial state combines in a coherent manner all the available
More informationLecture 13. Lidar Data Inversion
Lecture 13. Lidar Data Inversion Review Doppler lidar architecture Daytime capability of Na Doppler lidar Introduction Common raw data format Basic ideas (clues) for data inversion Preprocess Main process
More informationIce Cover and Sea and Lake Ice Concentration with GOES-R ABI
GOES-R AWG Cryosphere Team Ice Cover and Sea and Lake Ice Concentration with GOES-R ABI Presented by Yinghui Liu 1 Team Members: Yinghui Liu 1, Jeffrey R Key 2, and Xuanji Wang 1 1 UW-Madison CIMSS 2 NOAA/NESDIS/STAR
More informationTutorial 1. Introduction to Using FLUENT: Fluid Flow and Heat Transfer in a Mixing Elbow
Tutorial 1. Introduction to Using FLUENT: Fluid Flow and Heat Transfer in a Mixing Elbow Introduction This tutorial illustrates the setup and solution of the two-dimensional turbulent fluid flow and heat
More informationSNMP MIBs and Traps Supported
This section describes the MIBs available on your system. When you access your MIB data you will expose additional MIBs not listed in this section. The additional MIBs you expose through the process are
More informationImprovements to a GPS radio occultation ray-tracing model and their impacts on assimilation of bending angle
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 108, NO. D17, 4548, doi:10.1029/2002jd003160, 2003 Improvements to a GPS radio occultation ray-tracing model and their impacts on assimilation of bending angle H.
More information