Python regrid_interp Examples

Example #1

File: hourly_regrid.py Project: DeneBowdalo/AtmosChem_Tools

gfdl_lon_diff =  360./len(gfdl_lon_c)

giss_lat_diff =  180./len(giss_lat_c)
giss_lon_diff =  360./len(giss_lon_c)

miroc_chem_lat_diff =  180./len(miroc_chem_lat_c)
miroc_chem_lon_diff =  360./len(miroc_chem_lon_c)


min_lat = np.min([cesm_cam_lat_diff,cmam_lat_diff,geos_ccm_lat_diff,geos_chem_lat_diff,gfdl_lat_diff,giss_lat_diff,miroc_chem_lat_diff])
min_lon = np.min([cesm_cam_lon_diff,cmam_lon_diff,geos_ccm_lon_diff,geos_chem_lon_diff,gfdl_lon_diff,giss_lon_diff,miroc_chem_lon_diff])

min_lat = 1.
min_lon = 1.

new_lat_e,new_lat_c,new_lon_e,new_lon_c,z_cesm_cam = modules.regrid_interp(cesm_cam_param,cesm_cam_lat_e,cesm_cam_lat_c,cesm_cam_lon_e,cesm_cam_lon_c,min_lat,min_lon,data_type,max_ph)
new_lat_e,new_lat_c,new_lon_e,new_lon_c,z_cmam = modules.regrid_interp(cmam_param,cmam_lat_e,cmam_lat_c,cmam_lon_e,cmam_lon_c,min_lat,min_lon,data_type,max_ph)
new_lat_e,new_lat_c,new_lon_e,new_lon_c,z_geos_ccm = modules.regrid_interp(geos_ccm_param,geos_ccm_lat_e,geos_ccm_lat_c,geos_ccm_lon_e,geos_ccm_lon_c,min_lat,min_lon,data_type,max_ph)
new_lat_e,new_lat_c,new_lon_e,new_lon_c,z_geos_chem = modules.regrid_interp(geos_chem_param,geos_chem_lat_e,geos_chem_lat_c,geos_chem_lon_e,geos_chem_lon_c,min_lat,min_lon,data_type,max_ph)
new_lat_e,new_lat_c,new_lon_e,new_lon_c,z_gfdl = modules.regrid_interp(gfdl_param,gfdl_lat_e,gfdl_lat_c,gfdl_lon_e,gfdl_lon_c,min_lat,min_lon,data_type,max_ph)
new_lat_e,new_lat_c,new_lon_e,new_lon_c,z_giss = modules.regrid_interp(giss_param,giss_lat_e,giss_lat_c,giss_lon_e,giss_lon_c,min_lat,min_lon,data_type,max_ph)
new_lat_e,new_lat_c,new_lon_e,new_lon_c,z_miroc_chem = modules.regrid_interp(miroc_chem_param,miroc_chem_lat_e,miroc_chem_lat_c,miroc_chem_lon_e,miroc_chem_lon_c,min_lat,min_lon,data_type,max_ph)

if (data_type == 'amp') or (data_type == 'ave'):
    z_std = np.std((z_cesm_cam,z_cmam,z_geos_ccm,z_geos_chem,z_gfdl,z_giss,z_miroc_chem),axis=0)
    z_ave = np.average((z_cesm_cam,z_cmam,z_geos_ccm,z_geos_chem,z_gfdl,z_giss,z_miroc_chem),axis=0)
    z_pc = (z_std/z_ave)*100.
if data_type == 'ph':
    data = np.array([z_cesm_cam,z_cmam,z_geos_ccm,z_geos_chem,z_gfdl,z_giss,z_miroc_chem])
    z_ave,z_std = modules.yamartino_ave_std(data,max_ph)
    z_pc = (z_std/z_ave)*100.

Example #2

miroc_chem_lon_diff = 360. / len(miroc_chem_lon_c)

min_lat = np.min([
    cesm_cam_lat_diff, cmam_lat_diff, geos_ccm_lat_diff, geos_chem_lat_diff,
    gfdl_lat_diff, giss_lat_diff, miroc_chem_lat_diff
])
min_lon = np.min([
    cesm_cam_lon_diff, cmam_lon_diff, geos_ccm_lon_diff, geos_chem_lon_diff,
    gfdl_lon_diff, giss_lon_diff, miroc_chem_lon_diff
])

min_lat = 1.
min_lon = 1.

new_lat_e, new_lat_c, new_lon_e, new_lon_c, z_cesm_cam = modules.regrid_interp(
    cesm_cam_param, cesm_cam_lat_e, cesm_cam_lat_c, cesm_cam_lon_e,
    cesm_cam_lon_c, min_lat, min_lon, data_type, max_ph)
new_lat_e, new_lat_c, new_lon_e, new_lon_c, z_cmam = modules.regrid_interp(
    cmam_param, cmam_lat_e, cmam_lat_c, cmam_lon_e, cmam_lon_c, min_lat,
    min_lon, data_type, max_ph)
new_lat_e, new_lat_c, new_lon_e, new_lon_c, z_geos_ccm = modules.regrid_interp(
    geos_ccm_param, geos_ccm_lat_e, geos_ccm_lat_c, geos_ccm_lon_e,
    geos_ccm_lon_c, min_lat, min_lon, data_type, max_ph)
new_lat_e, new_lat_c, new_lon_e, new_lon_c, z_geos_chem = modules.regrid_interp(
    geos_chem_param, geos_chem_lat_e, geos_chem_lat_c, geos_chem_lon_e,
    geos_chem_lon_c, min_lat, min_lon, data_type, max_ph)
new_lat_e, new_lat_c, new_lon_e, new_lon_c, z_gfdl = modules.regrid_interp(
    gfdl_param, gfdl_lat_e, gfdl_lat_c, gfdl_lon_e, gfdl_lon_c, min_lat,
    min_lon, data_type, max_ph)
new_lat_e, new_lat_c, new_lon_e, new_lon_c, z_giss = modules.regrid_interp(
    giss_param, giss_lat_e, giss_lat_c, giss_lon_e, giss_lon_c, min_lat,

Example #3

File: monthly_regrid.py Project: DeneBowdalo/AtmosChem_Tools

miroc_chem_lat_diff =  180./len(miroc_chem_lat_c)
miroc_chem_lon_diff =  360./len(miroc_chem_lon_c)

ncar_lat_diff =  180./len(ncar_lat_c)
ncar_lon_diff =  360./len(ncar_lon_c)

stochad_lat_diff =  180./len(stochad_lat_c)
stochad_lon_diff =  360./len(stochad_lon_c)

#min_lat = np.min([cesm_cam_lat_diff,cmam_lat_diff,emac_lat_diff,geos_ccm_lat_diff,geos_chem_lat_diff,gfdl_lat_diff,giss_lat_diff,hadgem2_lat_diff,miroc_chem_lat_diff,ncar_lat_diff,stochad_lat_diff])
#min_lon = np.min([cesm_cam_lon_diff,cmam_lon_diff,emac_lon_diff,geos_ccm_lon_diff,geos_chem_lon_diff,gfdl_lon_diff,giss_lon_diff,hadgem2_lon_diff,miroc_chem_lon_diff,ncar_lon_diff,stochad_lon_diff])

min_lat = 1.
min_lon = 1.

new_lat_e,new_lat_c,new_lon_e,new_lon_c,z_cesm_cam = modules.regrid_interp(cesm_cam_param,cesm_cam_lat_e,cesm_cam_lat_c,cesm_cam_lon_e,cesm_cam_lon_c,min_lat,min_lon,data_type,max_ph)
new_lat_e,new_lat_c,new_lon_e,new_lon_c,z_cmam = modules.regrid_interp(cmam_param,cmam_lat_e,cmam_lat_c,cmam_lon_e,cmam_lon_c,min_lat,min_lon,data_type,max_ph)
new_lat_e,new_lat_c,new_lon_e,new_lon_c,z_emac = modules.regrid_interp(emac_param,emac_lat_e,emac_lat_c,emac_lon_e,emac_lon_c,min_lat,min_lon,data_type,max_ph)
new_lat_e,new_lat_c,new_lon_e,new_lon_c,z_geos_ccm = modules.regrid_interp(geos_ccm_param,geos_ccm_lat_e,geos_ccm_lat_c,geos_ccm_lon_e,geos_ccm_lon_c,min_lat,min_lon,data_type,max_ph)
new_lat_e,new_lat_c,new_lon_e,new_lon_c,z_geos_chem = modules.regrid_interp(geos_chem_param,geos_chem_lat_e,geos_chem_lat_c,geos_chem_lon_e,geos_chem_lon_c,min_lat,min_lon,data_type,max_ph)
new_lat_e,new_lat_c,new_lon_e,new_lon_c,z_gfdl = modules.regrid_interp(gfdl_param,gfdl_lat_e,gfdl_lat_c,gfdl_lon_e,gfdl_lon_c,min_lat,min_lon,data_type,max_ph)
new_lat_e,new_lat_c,new_lon_e,new_lon_c,z_giss = modules.regrid_interp(giss_param,giss_lat_e,giss_lat_c,giss_lon_e,giss_lon_c,min_lat,min_lon,data_type,max_ph)
new_lat_e,new_lat_c,new_lon_e,new_lon_c,z_hadgem2 = modules.regrid_interp(hadgem2_param,hadgem2_lat_e,hadgem2_lat_c,hadgem2_lon_e,hadgem2_lon_c,min_lat,min_lon,data_type,max_ph)
new_lat_e,new_lat_c,new_lon_e,new_lon_c,z_miroc_chem = modules.regrid_interp(miroc_chem_param,miroc_chem_lat_e,miroc_chem_lat_c,miroc_chem_lon_e,miroc_chem_lon_c,min_lat,min_lon,data_type,max_ph)
new_lat_e,new_lat_c,new_lon_e,new_lon_c,z_ncar = modules.regrid_interp(ncar_param,ncar_lat_e,ncar_lat_c,ncar_lon_e,ncar_lon_c,min_lat,min_lon,data_type,max_ph)
new_lat_e,new_lat_c,new_lon_e,new_lon_c,z_stochad = modules.regrid_interp(stochad_param,stochad_lat_e,stochad_lat_c,stochad_lon_e,stochad_lon_c,min_lat,min_lon,data_type,max_ph)

if (data_type == 'amp') or (data_type == 'ave'):
    z_std = np.std((z_cesm_cam,z_cmam,z_emac,z_geos_ccm,z_geos_chem,z_gfdl,z_giss,z_hadgem2,z_miroc_chem,z_ncar,z_stochad),axis=0)
    z_ave = np.average((z_cesm_cam,z_cmam,z_emac,z_geos_ccm,z_geos_chem,z_gfdl,z_giss,z_hadgem2,z_miroc_chem,z_ncar,z_stochad),axis=0)
    z_pc = (z_std/z_ave)*100.

Example #4

miroc_chem_lon_diff = 360. / len(miroc_chem_lon_c)

ncar_lat_diff = 180. / len(ncar_lat_c)
ncar_lon_diff = 360. / len(ncar_lon_c)

stochad_lat_diff = 180. / len(stochad_lat_c)
stochad_lon_diff = 360. / len(stochad_lon_c)

#min_lat = np.min([cesm_cam_lat_diff,cmam_lat_diff,emac_lat_diff,geos_ccm_lat_diff,geos_chem_lat_diff,gfdl_lat_diff,giss_lat_diff,hadgem2_lat_diff,miroc_chem_lat_diff,ncar_lat_diff,stochad_lat_diff])
#min_lon = np.min([cesm_cam_lon_diff,cmam_lon_diff,emac_lon_diff,geos_ccm_lon_diff,geos_chem_lon_diff,gfdl_lon_diff,giss_lon_diff,hadgem2_lon_diff,miroc_chem_lon_diff,ncar_lon_diff,stochad_lon_diff])

min_lat = 1.
min_lon = 1.

new_lat_e, new_lat_c, new_lon_e, new_lon_c, z_cesm_cam = modules.regrid_interp(
    cesm_cam_param, cesm_cam_lat_e, cesm_cam_lat_c, cesm_cam_lon_e,
    cesm_cam_lon_c, min_lat, min_lon, data_type, max_ph)
new_lat_e, new_lat_c, new_lon_e, new_lon_c, z_cmam = modules.regrid_interp(
    cmam_param, cmam_lat_e, cmam_lat_c, cmam_lon_e, cmam_lon_c, min_lat,
    min_lon, data_type, max_ph)
new_lat_e, new_lat_c, new_lon_e, new_lon_c, z_emac = modules.regrid_interp(
    emac_param, emac_lat_e, emac_lat_c, emac_lon_e, emac_lon_c, min_lat,
    min_lon, data_type, max_ph)
new_lat_e, new_lat_c, new_lon_e, new_lon_c, z_geos_ccm = modules.regrid_interp(
    geos_ccm_param, geos_ccm_lat_e, geos_ccm_lat_c, geos_ccm_lon_e,
    geos_ccm_lon_c, min_lat, min_lon, data_type, max_ph)
new_lat_e, new_lat_c, new_lon_e, new_lon_c, z_geos_chem = modules.regrid_interp(
    geos_chem_param, geos_chem_lat_e, geos_chem_lat_c, geos_chem_lon_e,
    geos_chem_lon_c, min_lat, min_lon, data_type, max_ph)
new_lat_e, new_lat_c, new_lon_e, new_lon_c, z_gfdl = modules.regrid_interp(
    gfdl_param, gfdl_lat_e, gfdl_lat_c, gfdl_lon_e, gfdl_lon_c, min_lat,