Showing posts with label GIS. Show all posts
Showing posts with label GIS. Show all posts

Thursday, July 23, 2020

2m dew point, ERA5

In Kelvin

Temperature to which the air, at 2 metres above the surface of the Earth, would have to be cooled for saturation to occur.It is a measure of the humidity of the air. Combined with temperature and pressure, it can be used to calculate the relative humidity. 2m dew point temperature is calculated by interpolating between the lowest model level and the Earth's surface, taking account of the atmospheric conditions. Temperature measured in kelvin can be converted to degrees Celsius (°C) by subtracting 273.15.

https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-land?tab=overview


Low dew point == low humidity according to wikipedia, because less vapor is need for condensation
https://en.wikipedia.org/wiki/Dew_point

Dew pointRelative humidity at 32 °C (90 °F)
Over 26 °COver 80 °F73% and higher
24–26 °C75–80 °F62–72%
21–24 °C70–74 °F52–61%
18–21 °C65–69 °F44–51%
16–18 °C60–64 °F37–43%
13–16 °C55–59 °F31–36%
10–12 °C50–54 °F26–30%
Under 10 °CUnder 50 °F25% and lower
This seems to explain that the condensation on leaves is often seen in the early morning, but not during the daytime. So, if a mutant strain has a lower dew point, it tend be in drier locations. If a strain has a higher dew point, it tends to be more humid locations


Monday, June 15, 2020

CMIP5 daily data on single levels

CMIP5 daily data on single levels


https://cds.climate.copernicus.eu/cdsapp#!/dataset/projections-cmip5-daily-single-levels?tab=overview


This catalogue entry provides daily climate projections on single levels from a large number of experiments, models, members and time periods computed in the framework of the fifth phase of the Coupled Model Intercomparison Project (CMIP5). The term "single levels" is used to express that the variables are computed at one vertical level which can be surface (or a level close to the surface) or a dedicated pressure level in the atmosphere. Multiple vertical levels are excluded from this catalogue entry.
CMIP5 data are used extensively in the Intergovernmental Panel on Climate Change Assessment Reports (the latest one is IPCC AR5, which was published in 2014). The use of these data is mostly aimed at:
  • addressing outstanding scientific questions that arose as part of the IPCC reporting process;
  • improving the understanding of the climate system;
  • providing estimates of future climate change and related uncertainties;
  • providing input data for the adaptation to the climate change;
  • examining climate predictability and exploring the ability of models to predict climate on decadal time scales;
  • evaluating how realistic the different models are in simulating the recent past.
The term "experiments" refers to the three main categories of CMIP5 simulations:
  • Historical experiments which cover the period where modern climate observations exist. These experiments show how the GCMs performs for the past climate and can be used as a reference period for comparison with scenario runs for the future. The period covered is typically 1850-2005.;
  • Ensemble of experiments from the Atmospheric Model Intercomparison Project (AMIP), which prescribes the oceanic variables for all models and during all period of the experiment. This configuration removes the added complexity of ocean-atmosphere feedbacks in the climate system. The period covered is typically 1950-2005.
  • Ensemble of climate projection experiments following the Representative Concentration Pathways (RCP) 2.6, 4.5, 6.0 and 8.5. The RCP scenarios provide different pathways of the future climate forcing. The period covered is typically, 2006-2100 some extended RCP experimental data is available from 2100-2300.
In CMIP5, the same experiments were run using different GCMs. In addition, for each model, the same experiment was repeatedly done using slightly different conditions (like initial conditions or different physical parameterisations for instance) producing in that way an ensemble of experiments closely related. Note that CMIP5 GCM data can be also used as lateral boundary conditions for Regional Climate Models (RCMs). RCMs are also available in the CDS (see CORDEX datasets).
The data are produced by the participating institutes of the CMIP5 project. The latest CMIP GCM experiments will form the CMIP6 dataset, which will be published in the CDS in a later stage.
More details about the products are given in the Documentation section.
DATA DESCRIPTION
Data typeGridded
Horizontal coverageGlobal
Horizontal resolutionFrom 0.125°x0.125° to 5°x5° depending on the model
Vertical coverageVariables are provided in one single level, which may differ among variables
Temporal coverageFrom 1850 to 2300 (shorter for some experiments)
Temporal resolutionDay
File formatNetCDF
MAIN VARIABLES
NameUnitsDescription
10m wind speedm s-1Magnitude of the two-dimensional horizontal air velocity near the surface.
2m temperatureKTemperature of the air near the surface.
Daily near surface relative humidity%Amount of moisture in the air near the surface divided by the maximum amount of moisture that could exist in the air at a specific temperature and location.
Maximum 2m temperature in the last 24 hoursKDaily maximum near-surface air temperature.
Mean precipitation fluxkg m-2 s-1Amount of water per unit area and time.
Mean sea level pressurePaTime average of the air pressure at sea level.
Minimum 2m temperature in the last 24 hoursKDaily minimum near-surface air temperature.
Near surface specific humidityDimensionlessAmount of moisture in the air near the surface divided by amount of air plus moist at that location.
Snowfallkg m-2 s-1Mass of water in the form of snow precipitating per unit area.
Surface solar radiation downwardsW m-2Radiative shortwave flux of energy downward at the surface.
Record updated 2020-05-14 15:52:20 UTC




Thursday, June 11, 2020

USA states geo code coordinates


https://developers.google.com/public-data/docs/canonical/states_csv

statelatitudelongitudename
AK63.588753-154.493062Alaska
AL32.318231-86.902298Alabama
AR35.20105-91.831833Arkansas
AZ34.048928-111.093731Arizona
CA36.778261-119.417932California
CO39.550051-105.782067Colorado
CT41.603221-73.087749Connecticut
DC38.905985-77.033418District of Columbia
DE38.910832-75.52767Delaware
FL27.664827-81.515754Florida
GA32.157435-82.907123Georgia
HI19.898682-155.665857Hawaii
IA41.878003-93.097702Iowa
ID44.068202-114.742041Idaho
IL40.633125-89.398528Illinois
IN40.551217-85.602364Indiana
KS39.011902-98.484246Kansas
KY37.839333-84.270018Kentucky
LA31.244823-92.145024Louisiana
MA42.407211-71.382437Massachusetts
MD39.045755-76.641271Maryland
ME45.253783-69.445469Maine
MI44.314844-85.602364Michigan
MN46.729553-94.6859Minnesota
MO37.964253-91.831833Missouri
MS32.354668-89.398528Mississippi
MT46.879682-110.362566Montana
NC35.759573-79.0193North Carolina
ND47.551493-101.002012North Dakota
NE41.492537-99.901813Nebraska
NH43.193852-71.572395New Hampshire
NJ40.058324-74.405661New Jersey
NM34.97273-105.032363New Mexico
NV38.80261-116.419389Nevada
NY43.299428-74.217933New York
OH40.417287-82.907123Ohio
OK35.007752-97.092877Oklahoma
OR43.804133-120.554201Oregon
PA41.203322-77.194525Pennsylvania
PR18.220833-66.590149Puerto Rico
RI41.580095-71.477429Rhode Island
SC33.836081-81.163725South Carolina
SD43.969515-99.901813South Dakota
TN35.517491-86.580447Tennessee
TX31.968599-99.901813Texas
UT39.32098-111.093731Utah
VA37.431573-78.656894Virginia
VT44.558803-72.577841Vermont
WA47.751074-120.740139Washington
WI43.78444-88.787868Wisconsin
WV38.597626-80.454903West Virginia
WY43.075968-107.290284Wyoming

Tuesday, June 9, 2020

*** ERA5-Land hourly data from 1981 to present



https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-land?tab=overview

ERA5-Land hourly data from 1981 to present

Hong can only found Dec 1, 2020 data on Feb 19, 2021. 

2022-02-02. Hong found the title changed to "1950 to present", and resolution changed to 0.1 degree with a resolution about 9km. 

How to extract specific information from GRIB files?



https://gis.stackexchange.com/questions/173775/how-to-extract-specific-information-from-grib-files



https://github.com/jswhit/pygrib/issues/27

Coordinates, spatial grids
https://confluence.ecmwf.int/display/CKB/ERA5%3A+data+documentation#ERA5:datadocumentation-Spatialgrid

To find specific coordinates within this 720 x 360 array i just extrapolate the coordinates to the array. Ex: If I want to find data for New York New York coordinates ~40.5N 74 W I'll map those coordinates to the array to by doing 2x(90-40.5) and 2x (180-74) to find the following result:
grbs_old = pygrib.open('nomads.ncep.noaa.gov/cgi-bin/'+str(grib_filename))
data=grbs_old[6].values
data[99][212]
Grid box
https://confluence.ecmwf.int/display/CKB/Model+grid+box+and+time+step
Horizontal
In the horizontal the discrete points are arranged in a two dimensional grid and hence are referred to as grid points. The grid can be regular or irregular. An example of a regular latitude/longitude grid would be where the grid points are located every 1 degree of longitude in the east-west direction and every 1 degree of latitude in the north-south direction. Each grid point is associated with an area that either surrounds the grid point or lies between the grid points. This area is referred to as the "grid box".
Grid point values cannot properly represent variability on spatial scales smaller than the grid box.
(In addition to using grid points, the ECMWF Integrated Forecasting System (IFS) also uses an additional mathematical concept, spectral space, to represent horizontal space. This concept uses a set of wavy basis functions, spherical harmonics, to describe variations in the horizontal. The IFS switches between spectral space and grid point space, in order to perform specific computations.)


For Hong's Pygrib file,
grbs = pygrib.open('2019-2020June10.grib')
valuesTmp = grbs[1].values
valuesTmp.shape #return (721, 1440)

so, 720/180 = 4

1440/360 = 4

It turns out that GRIB use 0-360 longitude and many software use -180W, 180E longitude. 

Hong tried ( - 180W mod 360), and temperature results from DC matches AccuWeather results,





Tuesday, January 10, 2017

Python and GIS course

A student said that Python is used in GIS course. He was asked to review Python before taking the course.