Example #1
0
def window_cloud_ceiling(f):
    dset = da.read_nc(f, ["cal333+cal05+cal20+cal80+csat", "cal333+cal05+cal20+cal80+csat_cprof", "nprof"])
    vcm, cprof, nprof = (
        dset["cal333+cal05+cal20+cal80+csat"],
        dset["cal333+cal05+cal20+cal80+csat_cprof"],
        dset["nprof"],
    )
    cf_lat, cf_lat2 = compute_cf(vcm, cprof, nprof)

    ceiling = tropic_width.cloud_cover_top(vcm.altitude, cf_lat)
    return ceiling, vcm.lat, cprof, nprof
Example #2
0
def window_cloud_ceiling(f):
    dset = da.read_nc(f, [
        'cal333+cal05+cal20+cal80+csat', 'cal333+cal05+cal20+cal80+csat_cprof',
        'nprof'
    ])
    vcm, cprof, nprof = dset['cal333+cal05+cal20+cal80+csat'], dset[
        'cal333+cal05+cal20+cal80+csat_cprof'], dset['nprof']
    cf_lat, cf_lat2 = compute_cf(vcm, cprof, nprof)

    ceiling = tropic_width.cloud_cover_top(vcm.altitude, cf_lat)
    return ceiling, vcm.lat, cprof, nprof
Example #3
0
def pcolor_zonal(x, y, vcmarray, title=None):

    z = vcmarray.copy()
    z = np.ma.masked_where(z==0, z)
    plt.contourf(x, y, z.T, np.r_[0:0.8:0.05])
    plt.colorbar()
    if title is not None:
        plt.title(title[10:-4])
    plt.clim(0,0.8)
    plt.xticks(np.r_[-90:90+30:30])
    ctop = tropic_width.cloud_cover_top(y, vcmarray, 0.05)
    latup, latdown, ceiling, height = tropic_width.tropic_width3(x, y, vcmarray)
    plt.axvline(x=latup, ls='--', color='grey')
    plt.axvline(x=latdown, ls='--', color='grey')
    plt.axhline(y=ceiling, ls='--', color='grey')
    plt.axhline(y=height, ls='--', color='grey')
    plt.plot(x, ctop, ls='--', color='w', alpha=0.5)
    plt.xlim(-82,82)
    plt.ylim(0, 20) 
    plt.ylabel('Altitude [km]')   
    plt.xlabel('Latitude')
Example #4
0
def pcolor_zonal(x, y, vcmarray, title=None):

    z = vcmarray.copy()
    z = np.ma.masked_where(z == 0, z)
    plt.contourf(x, y, z.T, np.r_[0:0.8:0.05])
    plt.colorbar()
    if title is not None:
        plt.title(title[10:-4])
    plt.clim(0, 0.8)
    plt.xticks(np.r_[-90:90 + 30:30])
    ctop = tropic_width.cloud_cover_top(y, vcmarray, 0.05)
    latup, latdown, ceiling, height = tropic_width.tropic_width3(
        x, y, vcmarray)
    plt.axvline(x=latup, ls='--', color='grey')
    plt.axvline(x=latdown, ls='--', color='grey')
    plt.axhline(y=ceiling, ls='--', color='grey')
    plt.axhline(y=height, ls='--', color='grey')
    plt.plot(x, ctop, ls='--', color='w', alpha=0.5)
    plt.xlim(-82, 82)
    plt.ylim(0, 20)
    plt.ylabel('Altitude [km]')
    plt.xlabel('Latitude')
Example #5
0
def pcolor_zonal(x, y, vcmarray, title=None):

    plt.figure(figsize=[15,5])
    plt.contourf(x, y, vcmarray.T, np.r_[0:0.8:0.05])
    plt.colorbar()
    if title is not None:
        plt.title(title)
    plt.clim(0,0.8)
    plt.xlim(-82,82)
    plt.xticks(np.r_[-90:90+30:30])
    
    cover_top = tropic_width.cloud_cover_top(y, vcmarray)
    plt.plot(x, cover_top, lw=2, color='w', alpha=0.5)
    
    latup, latdown, ceiling, height = tropic_width.tropic_width3(x, y, vcmarray)
    plt.figure(figsize=[15,5])
    plt.plot(x, cover_top)
    plt.axvline(x=latup)
    plt.axvline(x=latdown)
    plt.axhline(y=ceiling)
    plt.axhline(y=height)
        
    print 'lat range, second method : ', latup, latdown