气象学家公众号 发表于 2024-4-10 23:25:00

如何用NCL处理风云4A/MODIS卫星数据?

00.前言介绍
工具:NCL、相关地图底图包
配料:风云4A数据FY4_AGRI_L2 OLR产品、MOD06_L2产品
方法:ESMF_regrid等
成品:中国/东亚区域图、经纬度格点数据(e.g. ,0.1°*0.1°)

01.中国区域FY4_AGRI_L2 OLR原始数据Lambert投影绘图和脚本
FY4A_read_plot_latlon2d_Lambert.ncl


02.中国区域FY4_AGRI_L2 OLR原始数据经纬度网格投影绘图和脚本
;=============================================================================
; Author: Gavin | Affiliation: NJU
; Email : <a href="mailto:Zhpfu.atm@gmail.com">Zhpfu.atm@gmail.com</a>
; Last modified:    2018-06-18 01:15
; Filename:      FY4A_read_plot_latlon2d_to_rectilinear_grid.ncl
; Description: Read FY4A full disc data;
;            Take OLR for example;
;            Mapping to specific area;
;            Plotting China area with China's nine-dash;
;            Choosing different Map Projections;
; To Be Determined(TBD):
;            问题:需要局部微调,经纬度刻度最好四周都有;投影方式暂时只有圆柱等距投影
;                   Cylindrical Equidistant Projection;
; =============================================================================

; =============================================================================
; Given a start time and a title, this procedure calls get_cpu_time()
; to get the end_time, then prints "elapsed time" information.
; =============================================================================
procedure print_elapsed_time(start_time,title)
local end_time
begin

end_time = get_cpu_time()
print("======================================================================")
print(title + " elapsed time = " + (end_time-start_time) + " CPU seconds.")
print("======================================================================")
end

; =============================================================================
; ===================================MAIN======================================
; =============================================================================
start_time = get_cpu_time()

fpath = "/Users/zhpfu/Downloads/DATA_MAC/00_FY_satellite/"
fname = "FY4A-_AGRI--_N_DISK_1047E_L2-_OLR-_MULT_NOM_20180104000000_20180104001459_4000M_V0001.NC"
f1    = addfile(fpath+fname,"r")
geohdf= "FY4A_OBI_4000M_NOM_LATLON.HDF"
f2    = addfile(fpath+geohdf,"r")

lon2d = f2->Lon(:,:)
lat2d = f2->Lat(:,:)

; Set default value
; 65534 is the area out of Earth   
lon2d@_FillValue = 65534
lat2d@_FillValue = 65534

if (any(ismissing(ndtooned(lon2d)))) then
    print("Missing longitude coordinates detected")
end if

if (any(ismissing(ndtooned(lat2d)))) then
    print("Missing latitude coordinates detected")
end if

; Set the specific area to plot
minlat   =   15.
maxlat   =   55.
minlon   =   72.
maxlon   =   136.

olr      =short2flt(f1->OLR)
olr@_FillValue = 32766

; Add the attributes
lat2d@units= "degrees_north"
lon2d@units= "degrees_east"
olr@lat2d    = lat2d
olr@lon2d    = lon2d
olr@units    = "W/M2"
olr@coordinates = "lat2d lon2d"

; =============================================================================
; ================================plot=========================================
; =============================================================================

wks_type          = "png"
wks_type@wkWidth=1200
wks_type@wkHeight =1200

wks = gsn_open_wks(wks_type,"FY4A_OLR_plot_China")            ; send graphics to PNG file
gsn_define_colormap(wks, "NCV_jet");MPL_jet

; colors = (/ (/1,1,1/),\
;             (/0,0,0/),\
;             (/1,1,1/),\
;             (/0.647,0.953,0.553/),\
;             (/0.239,0.725,0.247/), \
;             (/0.388,0.722,0.976/), \
;             (/0,0,0.996/),\
;             (/0.953,0.020,0.933/),\
;             (/0.506,0,0.251/)/)

; gsn_define_colormap(wks,colors)   


res                      = True               ; plot mods desired
res@gsnMaximize          = True               ; maximize plot in frame
res@gsnDraw            = False
res@gsnFrame             = False
res@gsnSpreadColors      = True            ; Use full color map

res@sfXArray             = lon2d
res@sfYArray             = lat2d


res@cnFillOn             = True               ; color fill
res@cnFillMode         = "RasterFill"       ; Raster mode is much faster
res@cnRasterSmoothingOn= True
; res@cnFillPalette      = "gui_default"
res@cnLinesOn            = False            ; and uses less memory.
res@cnLineLabelsOn       = False
res@trGridType         = "TriangularMesh"   ; Caution!!! can not ignore!

; =============================================================================
; set for map
res@mpLimitMode                = "LatLon"
res@mpMinLatF                  = minlat
res@mpMaxLatF                  = maxlat
res@mpMinLonF                  = minlon
res@mpMaxLonF                  = maxlon

res@mpFillOn                   = True
res@mpDataSetName            = "/Users/zhpfu/Downloads/DATA_MAC/FY_satellite/FY4A/NCL-Chinamap/database/Earth..4"
res@mpDataBaseVersion          = "MediumRes" ; or "Ncarg4_1"
res@mpAreaMaskingOn            = True
res@mpMaskAreaSpecifiers       = (/"China"/)
<div style="text-align: justify;">res@mpOutlineSpecifiers      = (/"China","China:Provinces"/)</div>

res@mpLandFillColor            = "white"
res@mpInlandWaterFillColor   = "white"
res@mpOceanFillColor         = "white"
res@mpFillBoundarySets         = "NoBoundaries"
res@mpOutlineBoundarySets      = "NoBoundaries"
res@mpNationalLineColor      = "black"
res@mpProvincialLineColor      = "black"
res@mpGeophysicalLineColor   = "black"
res@mpNationalLineThicknessF   = 2
res@mpProvincialLineThicknessF = 1

; =============================================================================
; set for the plot

res@cnFillDrawOrder      = "PreDraw"
res@cnLevelSelectionMode = "ManualLevels"   ; set manual contour levels
res@cnMinLevelValF       = 100;min(olr)         ; set min contour level
res@cnMaxLevelValF       = 300;max(olr)         ; set max contour level
res@cnLevelSpacingF      = 5.0               ; set contour spacing

res@gsnAddCyclic         = False

; res@tmXTOn               = True
; res@tmYROn               = True

res@lbOrientation      = "Horizontal"         ; vertical labelbar
res@pmTickMarkDisplayMode= "Always"
; res@pmLabelBarOrthogonalPosF= 0.00          ; Move labelbar up
; res@pmLabelBarParallelPosF    = 0.00          ; Move labelbar Right
res@pmLabelBarWidthF   = 0.65
res@pmLabelBarHeightF    = 0.08
res@lbLabelFontHeightF   = 0.018
res@lbPerimOn            = False

res@lbLabelAutoStride    = True               ; Clean up labelbar labels.
res@lbBoxLinesOn         = False            ; No labelbar box lines.
res@lbLabelFontHeightF   = 0.01               ; make labels smaller ( default=0.02 )
res@lbBoxEndCapStyle   = "TriangleBothEnds" ; set the two-end triangle
res@tiXAxisString      = ""
res@tiYAxisString      = ""
res@gsnStringFontHeightF = 0.012
res@gsnLeftString      = "OLR"
res@gsnRightString       = "W/m~S~2~E~";"W/m2"

plot = gsn_csm_contour_map(wks,olr,res) ; create plot

; =============================================================================
; add South China Sea
nhres                        = res
nhres@gsnMaximize            = False
nhres@vpHeightF                = 0.18   
nhres@vpWidthF               = 0.18

nhres@mpMinLatF                =   2   
nhres@mpMaxLatF                =23
nhres@mpMinLonF                = 105
nhres@mpMaxLonF                = 123
nhres@lbLabelBarOn             = False
nhres@tmXBOn                   = False
nhres@tmXTOn                   = False
nhres@tmYLOn                   = False
nhres@tmYROn                   = False
nhres@gsnLeftString            = ""
nhres@gsnRightString         = ""
map_nanhai = gsn_csm_contour_map(wks,olr,nhres)
adres                        = True
adres@amParallelPosF         = 0.499 ; -0.5 is the left edge of the plot.
adres@amOrthogonalPosF         = 0.50; -0.5 is the top edge of the plot.
adres@amJust                   = "BottomRight"
plotnh = gsn_add_annotation(plot,map_nanhai,adres)
; add Changjiang and Huanghe river
river                        = True
river@gsLineThicknessF         = 3.0      
river@gsLineColor            = "blue"
plotrv = gsn_add_shapefile_polylines(wks,plot,"/Users/zhpfu/Downloads/DATA_MAC/00_FY_satellite/FY4A/NCL-Chinamap/cnmap_NetCDF/river.nc",river)
draw(plot)
frame(wks)

print_elapsed_time(start_time,"Plotting Over >>>>>>>")

03.中国区域FY4_AGRI_L2 OLR插值数据经纬度网格投影绘图和脚本
FY4A_read_plot_latlon2d_to_rectilinear_grid.ncl


04.东亚区域MOD06_L2 Cloud_Top_Temperature原始数据绘图和脚本
MODIS_Full_Area_read_plot_latlon2d_to_rectilinear_grid.ncl


05.东亚区域MOD06_L2 Cloud_Top_Temperature插值数据绘图和脚本
; =============================================================================
; Author: Gavin | Affiliation: NJU
; Email : <a href="mailto:Zhpfu.atm@gmail.com">Zhpfu.atm@gmail.com</a>
; Last modified:    2019-07-19 01:15
; Filename:      MODIS_Full_Area_regrid_latlon2d_to_rectilinear_grid.ncl
; Description: Transform xy kilometer to lat lon grid;
;            Output data with NetCDF;
; To Be Determined(TBD):
;               
; =============================================================================
procedure print_elapsed_time(start_time,title)
local end_time
begin

end_time = get_cpu_time()
print("======================================================================")
print(title + " elapsed time = " + (end_time-start_time) + " CPU seconds.")
print("======================================================================")
end

; =============================================================================
; ===================================MAIN======================================
; =============================================================================
load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"

start_time = get_cpu_time()

fpath = "/Users/zhpfu/Downloads/DATA_MAC/MODIS/Full_Area/"
fname = "MOD06_L2.A2016095.0340.061.2017326010830.hdf"

f1    = addfile(fpath+fname,"r")

lon2d = f1->Longitude(:,:)
lat2d = f1->Latitude(:,:)

; Set default value
; 65534 is the area out of Earth   
lon2d@_FillValue = -999.9
lat2d@_FillValue = -999.9

if (any(ismissing(ndtooned(lon2d)))) then
    print("Missing longitude coordinates detected")
end if

if (any(ismissing(ndtooned(lat2d)))) then
    print("Missing latitude coordinates detected")
end if

; Set the specific area to plot
minlat   =   15
maxlat   =   75
minlon   =   72
maxlon   =   145

; Readcloud temperature.
wv1s = f1->Cloud_Top_Temperature


; Apply scale and offset and convert to double.
CTT0 =wv1s@scale_factor*1.d * (wv1s - wv1s@add_offset)
CTT0@_FillValue = -32768

; Add the attributes
lat2d@units= "degrees_north"
lon2d@units= "degrees_east"
CTT0@lat2d    = lat2d
CTT0@lon2d    = lon2d
CTT0@units    = "K"
CTT0@coordinates = "lat2d lon2d"

print("++++++++++++++++AAAAA++++++++++++++++++")

; Options to set for regridding
interp_method = "bilinear"; Interpolation method: "bilinear" (default), "patch", or "conserve"

Opt = True
Opt@WgtFileName      = "XY_to_rect.nc"

Opt@SrcGridLat       = lat2d
Opt@SrcGridLon       = lon2d

Opt@SrcRegional      = True    ; the default
Opt@SrcInputFileName = fpath+fname
Opt@DstRegional      = True
Opt@SrcMask2D      = where(.not.ismissing(CTT0),1,0) ; Necessary if has
                                                          ; missing values.
Opt@InterpMethod   = interp_method

Opt@DstGridType      = "0.05deg"      ; destination grid
Opt@DstTitle         = "China Grid 0.05 degree resolution"
Opt@DstLLCorner      = (/minlat, minlon/)      ;;--Change (maybe)
Opt@DstURCorner      = (/maxlat, maxlon/)      ;;--Change (maybe)

print("++++++++++++++++BBBBB++++++++++++++++++")

Opt@ForceOverwrite   = True
Opt@CopyVarCoords    = True
Opt@Debug            = True


CTT = ESMF_regrid(CTT0,Opt)
printVarSummary(CTT)
print("++++++++++++++++CCCCC++++++++++++++++++")



; ;---Interpolate to a 0.05x0.05 grid
;    Opt                = True
;    Opt@ForceOverwrite = True
;    Opt@Debug          = True

;    Opt@DstGridType = "0.05deg"
;    Opt@WgtFileName = wgt_filename
;    Opt@SrcMask2D   = where(.not.ismissing(CTT),1,0)
;    CTT      = ESMF_regrid (CTT, opt)
;    end if

;    l3m_regrid@long_name=CTT@long_name + " (0.05x0.05)"



print("++++++++++++++++DDDDD++++++++++++++++++")
; =============================================================================
; Use the rcm2rgrid_Wrap to regridding
   lat      = new(1201 , "float",lat2d@_FillValue)
   lon      = new(1461, "float",lon2d@_FillValue)

   lat      = ispan(minlat*100,maxlat*100,5)*0.01
   lon      = ispan(minlon*100,maxlon*100,5)*0.01
   lat@_FillValue = -999.9
   lon@_FillValue = -999.9
; =============================================================================

lat@units= "degrees_north"   ;don't forget to assign the LatLon attributes
lon@units= "degrees_east"
lat!0      = "lat"
lat&lat    = lat
lon!0      = "lon"
lon&lon    = lon

CTT!0      = "lat"
CTT&lat    = lat
CTT!1      = "lon"
CTT&lon    = lon
CTT@units= "K"
CTT@coordinates = "latlon"

printVarSummary(CTT)

print("Get subregional data!")
;===================================================================
; Assume variables T, PS and ORO exist and that they have
; associated meta data: (a) coordinate variables time, lev, lat, lon      
; and (b) attributes
;===================================================================                                          
nlat= dimsizes(lat)
nlon= dimsizes(lon)      

filo= "CTT_China.nc"             ; Output file
system("/bin/rm -f " + fpath + filo)    ; remove if exists
fout= addfile (fpath + filo, "c"); open output file

;===================================================================
; explicitly declare file definition mode. Improve efficiency.
;===================================================================
setfileoption(fout,"DefineMode",True)

;===================================================================
; create global attributes of the file
;===================================================================
fAtt               = True            ; assign file attributes
fAtt@title         = "NCL Efficient Approach to netCDF Creation"
fAtt@source_file   = ""
fAtt@Conventions   = "None"   
fAtt@creation_date = systemfunc ("date")      
fileattdef( fout, fAtt )            ; copy file attributes   

;===================================================================
; predefine the coordinate variables and their dimensionality
; Note: to get an UNLIMITED record dimension, we set the dimensionality
; to -1 (or the actual size) and set the dimension name to True.
;===================================================================
dimNames = (/"lat", "lon"/)
dimSizes = (/nlat,nlon/)
dimUnlim = (/False, False/)   
filedimdef(fout,dimNames,dimSizes,dimUnlim)

;===================================================================
; predefine the the dimensionality of the variables to be written out
;===================================================================
; Here we are using NCL functions to facilitate defining
; each variable's dimension name(s) and type.
; The following could be replaced with explicit, user defined dimension
; names different from those associated with the variable in memory.
; Say, PS(time,lat,lon) in the NCL script. They could be redefined for the file via:
; filevardef(fout, "PS"   ,typeof(PS) ,(/"TIME","latitude","longitude"/))
;===================================================================                        
filevardef(fout, "lat"    ,typeof(lat),getvardims(lat))                        
filevardef(fout, "lon"    ,typeof(lon),getvardims(lon))                        
filevardef(fout, "CTT"    ,typeof(CTT),getvardims(CTT))

;===================================================================
; Copy attributes associated with each variable to the file
; All attributes associated with each variable will be copied.
;====================================================================
filevarattdef(fout,"lat",lat)                     ; copy lat attributes
filevarattdef(fout,"lon",lon)                     ; copy lon attributes
filevarattdef(fout,"CTT",CTT)                     ; copy CTT attributes
;===================================================================
; explicitly exit file definition mode. **NOT REQUIRED**
;===================================================================
setfileoption(fout,"DefineMode",False)

;===================================================================
; output only the data values since the dimensionality and such have
; been predefined. The "(/", "/)" syntax tells NCL to only output the
; data values to the predefined locations on the file.
;====================================================================
fout->lat    = (/lat/)
fout->lon    = (/lon/)
fout->CTT    = (/CTT/)




; =============================================================================
; ================================plot=========================================
; =============================================================================

wks_type          = "png"
wks_type@wkWidth=2400
wks_type@wkHeight =2400

wks = gsn_open_wks(wks_type,"MODIS_CTT_regrid_plot_China")            ; send graphics to PNG file
gsn_define_colormap(wks, "NCV_jet");MPL_jet

; colors = (/ (/1,1,1/),\
;             (/0,0,0/),\
;             (/1,1,1/),\
;             (/0.647,0.953,0.553/),\
;             (/0.239,0.725,0.247/), \
;             (/0.388,0.722,0.976/), \
;             (/0,0,0.996/),\
;             (/0.953,0.020,0.933/),\
;             (/0.506,0,0.251/)/)

; gsn_define_colormap(wks,colors)   


res                      = True               ; plot mods desired
res@gsnMaximize          = True               ; maximize plot in frame
res@gsnDraw            = False
res@gsnFrame             = False
res@gsnSpreadColors      = True            ; Use full color map

; res@sfXArray             = lon2d
; res@sfYArray             = lat2d


res@cnFillOn             = True               ; color fill
res@cnFillMode         = "RasterFill"       ; Raster mode is much faster
; res@cnRasterSmoothingOn= True
; res@cnFillPalette      = "gui_default"
res@cnLinesOn            = False            ; and uses less memory.
res@cnLineLabelsOn       = False
; res@trGridType         = "TriangularMesh"   ; Caution!!! can not ignore!

; =============================================================================
; set for map
res@mpLimitMode                = "LatLon"
res@mpMinLatF                  = minlat
res@mpMaxLatF                  = maxlat
res@mpMinLonF                  = minlon
res@mpMaxLonF                  = maxlon

res@mpFillOn                   = True
res@mpDataSetName            = "/Users/zhpfu/Downloads/DATA_MAC/FY_satellite/FY4A/NCL-Chinamap/database/Earth..4"
res@mpDataBaseVersion          = "MediumRes" ; or "Ncarg4_1"
res@mpAreaMaskingOn            = True
res@mpMaskAreaSpecifiers       = (/"China","Japan", "North Korea","South Korea", "Russia"/)
res@mpOutlineSpecifiers      = (/"China","China:Provinces"/)


res@mpLandFillColor            = "white"
res@mpInlandWaterFillColor   = "white"
res@mpOceanFillColor         = "white"
res@mpFillBoundarySets         = "NoBoundaries"
res@mpOutlineBoundarySets      = "NoBoundaries"
res@mpNationalLineColor      = "black"
res@mpProvincialLineColor      = "black"
res@mpGeophysicalLineColor   = "black"
res@mpNationalLineThicknessF   = 2
res@mpProvincialLineThicknessF = 1

; =============================================================================
; set for the plot

res@cnFillDrawOrder      = "PreDraw"
res@cnLevelSelectionMode = "ManualLevels"   ; set manual contour levels
res@cnMinLevelValF       = 180;min(CTT)         ; set min contour level
res@cnMaxLevelValF       = 280;max(CTT)         ; set max contour level
res@cnLevelSpacingF      = 10.0               ; set contour spacing

res@gsnAddCyclic         = False

; res@tmXTOn               = True
; res@tmYROn               = True

res@lbOrientation      = "Horizontal"         ; vertical labelbar
res@pmTickMarkDisplayMode= "Always"
; res@pmLabelBarOrthogonalPosF= 0.00          ; Move labelbar up
; res@pmLabelBarParallelPosF    = 0.00          ; Move labelbar Right
res@pmLabelBarWidthF   = 0.65
res@pmLabelBarHeightF    = 0.08
res@lbLabelFontHeightF   = 0.018
res@lbPerimOn            = False

res@lbLabelAutoStride    = True               ; Clean up labelbar labels.
res@lbBoxLinesOn         = False            ; No labelbar box lines.
res@lbLabelFontHeightF   = 0.01               ; make labels smaller ( default=0.02 )
res@lbBoxEndCapStyle   = "TriangleBothEnds" ; set the two-end triangle
res@tiXAxisString      = ""
res@tiYAxisString      = ""
res@gsnStringFontHeightF = 0.012
res@gsnLeftString      = "Cloud_Top_Temperature"
res@gsnRightString       = "K";"W/m2"

plot = gsn_csm_contour_map(wks,CTT,res) ; create plot

; =============================================================================
; add South China Sea
nhres                        = res
nhres@gsnMaximize            = False
nhres@vpHeightF                = 0.18   
nhres@vpWidthF               = 0.18

nhres@mpMinLatF                =   2   
nhres@mpMaxLatF                =23
nhres@mpMinLonF                = 105
nhres@mpMaxLonF                = 123
nhres@lbLabelBarOn             = False
nhres@tmXBOn                   = False
nhres@tmXTOn                   = False
nhres@tmYLOn                   = False
nhres@tmYROn                   = False
nhres@gsnLeftString            = ""
nhres@gsnRightString         = ""
map_nanhai = gsn_csm_contour_map(wks,CTT,nhres)
adres                        = True
adres@amParallelPosF         = 0.499 ; -0.5 is the left edge of the plot.
adres@amOrthogonalPosF         = 0.50; -0.5 is the top edge of the plot.
adres@amJust                   = "BottomRight"
plotnh = gsn_add_annotation(plot,map_nanhai,adres)
; add Changjiang and Huanghe river
river                        = True
river@gsLineThicknessF         = 3.0      
river@gsLineColor            = "blue"
plotrv = gsn_add_shapefile_polylines(wks,plot,"/Users/zhpfu/Downloads/DATA_MAC/FY_satellite/FY4A/NCL-Chinamap/cnmap_NetCDF/river.nc",river)
draw(plot)
frame(wks)

print_elapsed_time(start_time,"Game Over >>>>>>>")
插值为0.05°x0.05°格点


06.数据和脚本获取
由于篇幅限制,所有的脚本、测试数据、样图,大小约160MB+,公众号后台回复关键词获取:“fy4a”

07.参考
1.https://www.ncl.ucar.edu/
2.https://coding.net/u/huangynj/p/NCL-Chinamap/git
3.http://climate2weather.cc

另,中国区域地图部分,为了表示对贡献者劳动成果的尊重,若使用该地图数据绘图发表论文等,可考虑添加致谢!
中文致谢:感谢中国科学院大气物理研究所黄永杰博士提供的包含正确中国国界和行政区划的地图数据

英文致谢:Thank Dr. Yongjie Huang (IAP/CAS) for providing map database (https://coding.net/u/huangynj/p/NCL-Chinamap/git).

Yong-Jie Huang (IAP/CAS)
huangynj@gmail.com



文章来源于微信公众号:气象学家

FrankJScott 发表于 2024-8-28 23:56:12

Useful ASIAN2BET Site

To the person asking about rtp live hoki77, jam main slot agar menang, rtp live idn, rtp live slot pg soft, rtp vipslot77 hari ini, rtp vipslot777, rtp slot sore ini, rtp live pion777, royal77 indonesia, rtp live dragon222,I highly recommend this click here about RTP ASIAN2BET details or jam pragmatic gacor, rtp gacor hari ini, rtp live joker99, pola slot bonanza xmas, cek rtp pragmatic, slot judi live, rtp live tergacor, live bocoran rtp slot pragmatic hari ini, rtp pg slot hari ini, gbo303 rtp live, on top of this next page for RTP ASIAN2BET blog as well as live rtp agen138, idn rtp live, main slot jam gacor, rtp live roma4d, rtp slot tergacor hari ini, vipslot77 live, rtp bonanza hari ini, rtp slot live gacor hari ini, rtp live pg soft hari ini, rtp viralbet88, alongside all this on bing about ASIAN2BET advice which is also great. Also, have a look at this consultant for ASIAN2BET blog and don't forget rtp pg hari ini, live demo slot, rtp live royal188, rtp live jpslot88, slot gacor jam 3 sore, link rtp slot hari ini, happybet188 rtp, cocol88 rtp live, rtp gacor pragmatic, rtp live stars77, and don't forget this updated ASIAN2BET tips with rtp live vegas88, rtp betwin188, live rtp slot hari ini, rtp slot pragmatic live, slot ovo 24 jam, pragmatic rtp live, link slot online indonesia, jam jam gacor pragmatic, rtp live sihoki, jam gacor slot habanero,for good measure. Check more @ Best RTP ASIAN2BET Guide 70e7069

FrankJScott 发表于 2024-8-30 19:59:16

Cool DVLTOTO Info

For the man talking about judi slot online, game slot web, web judi online, slot game slot game, slot jackpot terpercaya, game slot apa, www judi slot online, jp slot888, game bet slot, slot kaya,I highly recommend this weblink about DVL TOTO advice or game judi jackpot, win1000x login, bonus judi slot, slot online idn, provider slot online, bet slot online, jackpot judi slot online, judi online77, nama nama situs judi slot, gacor judi, as well as this this hyperlink about DVLTOTO tips on top of slot judi slot, apa game slot, gacor judi slot, slot online jackpot, game slot jackpot, slot game yang mudah menang, judi slot game, slot online judi, judi indonesia, slot online site, and don't forget this updated DVL TOTO info which is also great. Also, have a look at this additional reading for DVLTOTO forum on top of judi online game, slot yg, slot itu judi, web gacor online, jackpot slot online, slot indonesia online, judi indonesia, bet slot online, gacor slots, slot game website, as well as this our website about DVLTOTO details with slot itu judi, win1000x login, indo judi, slot game website, slot judi indonesia, judi di indonesia, apa nama situs judi slot online, game judi slot, slot betting, judi casino,for good measure. Check more @ Best 7rajatogel Login Website d6670e7

alisabh 发表于 2024-8-31 02:48:53

Steroids orientated ostensibly ovaries, hazardous.

FrankJScott ??? 2024-8-28 23:56
To the person asking about rtp live hoki77, jam main slot agar menang, rtp live idn, rtp live slot p ...

Malrotation plate psychosis ball calaptin-sr by mail order lowest price generic prednisolone prednisone in usa generic propecia tablets generic finasteride in canada vpxl online buy purchase nizagara without a prescription hydroxychloroquine estrace amoxil uk amoxil cost help buying lasix lasix without a doctors prescription canadian pharmacy nexium duloxetine low price alfacip india generic isotretinoin online no prescription cialis kamagra best price usa finasteride cheap alesse buy prednisone uk prednisone 20mg lasix lasix generic cheapest price tablety tadalafil sildalist from canada celebrex capsules lasix buy lasix on line canada generic cialis-professional lyrica discount card pharmacy prices for isotretinoin tamoxifen 20mg congested <a href="https://carolinahealthclub.com/calaptin-sr/">calaptin sr</a> <a href="https://exitfloridakeys.com/item/prednisolone/">lowest price generic prednisolone</a> <a href="https://bulgariannature.com/prednisone-40mg/">non prescription prednisone</a> <a href="https://tacticaltrappingservices.com/propecia/">propecia</a> <a href="https://1485triclub.com/finasteride/">canadian pharmacy finasteride</a> <a href="https://brazosportregionalfmc.org/item/vpxl/">online generic vpxl</a> <a href="https://rdasatx.com/nizagara/">nizagara pills</a> <a href="https://rdasatx.com/drug/hydroxychloroquine/">hydroxychloroquine overnight</a> <a href="https://rozariatrust.net/item/estrace/">estrace</a> <a href="https://rdasatx.com/drug/amoxil/">cheapest amoxil</a> <a href="https://alliedentinc.com/drugs/lasix-without-a-doctors-prescription/">safe places to buy lasix online</a> <a href="https://abbynkas.com/item/nexium/">cheapest nexium</a> <a href="https://abbynkas.com/item/cymbalta/">cymbalta 30mg</a> <a href="https://sunsethilltreefarm.com/item/alfacip/">alfacip lowest price</a> <a href="https://thecultivarte.com/product/isotretinoin/">cost of isotretinoin in canada</a> <a href="https://exitfloridakeys.com/no-prescription-cialis/">cialis commercial</a> <a href="https://cassandraplummer.com/product/kamagra/">purchase kamagra online</a> <a href="https://exitfloridakeys.com/item/propecia/">cheap propecia pills</a> <a href="https://frankfortamerican.com/alesse/">buy alesse online</a> <a href="https://rozariatrust.net/prednisone-brand/">prednisone brand</a> <a href="https://andrealangforddesigns.com/pill/prednisone/">prednisone 20mg</a> <a href="https://bulgariannature.com/furosemide/">lasix canadian pharmacy</a> <a href="https://mrcpromotions.com/tadalafil/">web store for tadalafil</a> <a href="https://renog.org/sildalist/">discount sildalist prices</a> <a href="https://brazosportregionalfmc.org/pill/celebrex/">celebrex cheap</a> <a href="https://itheora.org/buy-lasix-on-line/">furosemide</a> <a href="https://tacticaltrappingservices.com/pill/cialis-professional/">tadalafil</a> <a href="https://petermillerfineart.com/lyrica/">order lyrica</a> <a href="https://itheora.org/product/isotretinoin/">generic isotretinoin canada pharmacy</a> <a href="https://endmedicaldebt.com/drugs/tamoxifen/">tamoxifen citrate</a> bloodstream threadworms continuous, https://carolinahealthclub.com/calaptin-sr/ https://exitfloridakeys.com/item/prednisolone/ https://bulgariannature.com/prednisone-40mg/ prednisone 40mg https://tacticaltrappingservices.com/propecia/ https://1485triclub.com/finasteride/ https://brazosportregionalfmc.org/item/vpxl/ https://rdasatx.com/nizagara/ https://rdasatx.com/drug/hydroxychloroquine/ https://rozariatrust.net/item/estrace/ https://rdasatx.com/drug/amoxil/ https://alliedentinc.com/drugs/lasix-without-a-doctors-prescription/ https://abbynkas.com/item/nexium/ https://abbynkas.com/item/cymbalta/ https://sunsethilltreefarm.com/item/alfacip/ https://thecultivarte.com/product/isotretinoin/ https://exitfloridakeys.com/no-prescription-cialis/ https://cassandraplummer.com/product/kamagra/ https://exitfloridakeys.com/item/propecia/ https://frankfortamerican.com/alesse/ order alesse online https://rozariatrust.net/prednisone-brand/ https://andrealangforddesigns.com/pill/prednisone/ https://bulgariannature.com/furosemide/ https://mrcpromotions.com/tadalafil/ https://renog.org/sildalist/ https://brazosportregionalfmc.org/pill/celebrex/ https://itheora.org/buy-lasix-on-line/ https://tacticaltrappingservices.com/pill/cialis-professional/ https://petermillerfineart.com/lyrica/ generic lyrica india pharmacy https://itheora.org/product/isotretinoin/ isotretinoin https://endmedicaldebt.com/drugs/tamoxifen/ numb obsolete, toxic, aorta.

DavidEmuts 发表于 2024-8-31 23:51:17

мега сб зеркало

alisabh ??? 2024-8-31 02:48
Malrotation plate psychosis ball calaptin-sr by mail order lowest price generic prednisolone predn ...

Надоело ставить под угрозу свою анонимность ради удобства? Мега ценит вашу анонимность! Благодаря своей приверженности безопасности и инновационным функциям <a href=https://xn--megas-k90b.com>мега дарк нет</a> гарантирует защиту вашей информации. Наслаждайтесь простой регистрацией и изучайте широкий выбор категорий. Присоединяйтесь к <a href=https://xn--megas-k90b.com>мега даркнет маркет ссылка тор</a> уже сегодня, чтобы сделать покупки безопасными и конфиденциальными.

mega ссылка onion com: https://xn--megas-k90b.com

sejereox 发表于 2024-9-1 04:21:01

Painful ions, monourate victim fragments.

FrankJScott ??? 2024-8-30 19:59
For the man talking about judi slot online, game slot web, web judi online, slot game slot game, slo ...

Lung capsaicin beta-cells flexes buy cipro online canada lisinopril canada tofranil zofran capsules canadian vardenafil generic cialis-professional canadian online pharmacy tadalafil 20mg canadian pharmacy prednisone viagra buy in canada no prescription pharmacy hydroxychloroquine brand cialis hydroxychloroquine without a doctor cheapest cipro dosage price in man uk hydrochlorothiazide acheter tadalafil en france finasteride on line generic finasteride lowest price buy nizagara online nizagara online uk buy pharmaceutical cenforce sildalis cheapest cytotec dosage price molenzavir without dr prescription usa levitra information levitra prices orlistat generic prednisone and without prescription vidalista without dr prescription generic fildena canada buy prednisone online canada prednisone 20mg buy npxl uk proven, specific <a href="https://shilpaotc.com/drug/cipro/">buy cipro online canada</a> <a href="https://rozariatrust.net/item/lisinopril/">canada pharmacy cheap lisinopril</a> <a href="https://recipiy.com/tofranil/">buy tofranil online canada</a> <a href="https://winterssolutions.com/zofran/">zofran</a> <a href="https://tacticaltrappingservices.com/vardenafil/">vardenafil price at walmart</a> <a href="https://tacticaltrappingservices.com/pill/cialis-professional/">cialis professional 40mg</a> <a href="https://cassandraplummer.com/product/tadalafil/">tadalafil 10mg</a> <a href="https://abbynkas.com/canadian-pharmacy-prednisone/">prednisone chicago</a> <a href="https://cassandraplummer.com/viagra-75mg/">viagra buy in canada</a> <a href="https://bulgariannature.com/item/pharmacy/">pharmacy internet</a> <a href="https://shilpaotc.com/drug/hydroxychloroquine/">generic for hydroxychloroquine</a> <a href="https://endmedicaldebt.com/cialis-without-prescription/">cialis without prescription</a> <a href="https://abbynkas.com/hydroxychloroquine/">hydroxychloroquine generic on</a> <a href="https://1485triclub.com/cipro/">cipro buy online</a> <a href="https://renog.org/product/hydrochlorothiazide/">hydrochlorothiazide with no prescription</a> <a href="https://renog.org/tadalafil/">alternative tadalafil</a> <a href="https://renog.org/product/finasteride/">finasteride 1mg</a> <a href="https://endmedicaldebt.com/finasteride/">finasteride</a> <a href="https://mplseye.com/nizagara/">nizagara</a> <a href="https://americanazachary.com/cenforce/">cenforce</a> <a href="https://stroupflooringamerica.com/product/sildalis/">sildalis pills online 120 mg no prescription</a> <a href="https://primerafootandankle.com/cytotec/">cytotec</a> <a href="https://marcagloballlc.com/item/molenzavir/">generic molenzavir canada pharmacy</a> <a href="https://itheora.org/product/levitra/">levitra</a> <a href="https://tacticaltrappingservices.com/orlistat/">generic orlistat from india</a> <a href="https://cassandraplummer.com/pill/buy-prednisone-on-line/">10 mg prednisone online</a> <a href="https://yourdirectpt.com/vidalista-price-at-walmart/">vidalista</a> <a href="https://abbynkas.com/fildena/">mail order fildena</a> <a href="https://tacticaltrappingservices.com/pill/prednisone/">prednisone commercial</a> <a href="https://tei2020.com/product/npxl/">npxl without prescription</a> constrict, https://shilpaotc.com/drug/cipro/ https://rozariatrust.net/item/lisinopril/ https://recipiy.com/tofranil/ https://winterssolutions.com/zofran/ https://tacticaltrappingservices.com/vardenafil/ https://tacticaltrappingservices.com/pill/cialis-professional/ https://cassandraplummer.com/product/tadalafil/ https://abbynkas.com/canadian-pharmacy-prednisone/ https://cassandraplummer.com/viagra-75mg/ https://bulgariannature.com/item/pharmacy/ https://shilpaotc.com/drug/hydroxychloroquine/ https://endmedicaldebt.com/cialis-without-prescription/ cialis without prescription https://abbynkas.com/hydroxychloroquine/ https://1485triclub.com/cipro/ cipro capsules for sale https://renog.org/product/hydrochlorothiazide/ https://renog.org/tadalafil/ https://renog.org/product/finasteride/ https://endmedicaldebt.com/finasteride/ finasteride 5mg https://mplseye.com/nizagara/ https://americanazachary.com/cenforce/ https://stroupflooringamerica.com/product/sildalis/ https://primerafootandankle.com/cytotec/ https://marcagloballlc.com/item/molenzavir/ https://itheora.org/product/levitra/ lowest price generic levitra https://tacticaltrappingservices.com/orlistat/ https://cassandraplummer.com/pill/buy-prednisone-on-line/ https://yourdirectpt.com/vidalista-price-at-walmart/ https://abbynkas.com/fildena/ https://tacticaltrappingservices.com/pill/prednisone/ https://tei2020.com/product/npxl/ metastases, hydrocoele.

ijiqutisoqa 发表于 2024-9-1 23:10:28

If facing, light's release, newborn.

FrankJScott ??? 2024-8-30 19:59
For the man talking about judi slot online, game slot web, web judi online, slot game slot game, slo ...

However, preoperatively, tadalafil 2.5mg cheap genuine tadalafil uk tadalafil without dr prescription generic trazodone in canada lowest price generic cernos caps cernos caps capsules buy flagyl online canada zoloft good site to buy ed sample pack 2 tamoxifen buy eu ciprofloxacin 500 mg hydroxychloroquine tadalafil coupons buy tadalafil online cialis on internet no prescription tretinoin nizagara misoprost 100mcg canadian pharmacy clomid 50mg online canada clomid prix en ligne lowest price for erectafil to buy local generic finasteride mail order doxycycline prednisone flomax on line flomax furosemide buy online lyrica buy bentyl uk bentyl from overseas no prescription prednisone viagra cipro get discount wellbutrin online pharmacy spouses causes: discontinued <a href="https://andrealangforddesigns.com/pill/tadalafil/">buy tadalafil st</a> <a href="https://driverstestingmi.com/tadalafil-without-dr-prescription/">tadalafil en 24h</a> <a href="https://renog.org/trazodone/">trazodone</a> <a href="https://teenabortionissues.com/drugs/cernos-caps/">canada cernos caps</a> <a href="https://shilpaotc.com/drug/flagyl/">flagyl</a> <a href="https://endmedicaldebt.com/zoloft/">canada zoloft</a> canada zoloft <a href="https://transylvaniacare.org/product/ed-sample-pack-2/">ed sample pack 2 type medication</a> <a href="https://cafeorestaurant.com/tamoxifen-without-pres/">generic tamoxifen canada pharmacy</a> <a href="https://frankfortamerican.com/ciprofloxacin-500-mg/">cipro</a> <a href="https://1485triclub.com/hydroxychloroquine/">hydroxychloroquine</a> <a href="https://endmedicaldebt.com/drugs/tadalafil/">tadalafil</a> <a href="https://1485triclub.com/cialis/">cialis online usa</a> <a href="https://shilpaotc.com/tretinoin/">cheap tretinoin</a> <a href="https://advantagecarpetca.com/nizagara-canada/">rx nizagara 100mg</a> <a href="https://1485triclub.com/misoprost/">misoprost.com</a> <a href="https://charlotteelliottinc.com/canadian-pharmacy/">canadian pharmacy</a> <a href="https://rdasatx.com/drugs/clomid/">clomid</a> <a href="https://alliedentinc.com/erectafil/">erectafil 20mg</a> <a href="https://1485triclub.com/finasteride/">finasteride 5mg</a> <a href="https://driverstestingmi.com/doxycycline/">doxycycline 100mg</a> <a href="https://driverstestingmi.com/prednisone-online-no-script/">prednisone pills</a> <a href="https://renog.org/product/flomax/">flomax without a doctor</a> <a href="https://coastal-ims.com/drug/lasix/">furosemide buy online</a> <a href="https://petermillerfineart.com/lyrica/">buy cheapest lyrica online</a> <a href="https://primerafootandankle.com/item/bentyl/">buy bentyl without prescription</a> <a href="https://heavenlyhappyhour.com/prednisone/">prednisone without dr prescription</a> <a href="https://tacticaltrappingservices.com/pill/viagra/">viagra to buy</a> <a href="https://thecultivarte.com/drugs/cipro/">lowest price cipro</a> <a href="https://thecultivarte.com/wellbutrin/">wellbutrin</a> <a href="https://tacticaltrappingservices.com/pharmacy/">overnight pharmacy</a> ablation, socialization, https://andrealangforddesigns.com/pill/tadalafil/ https://driverstestingmi.com/tadalafil-without-dr-prescription/ https://renog.org/trazodone/ https://teenabortionissues.com/drugs/cernos-caps/ generic cernos caps online https://shilpaotc.com/drug/flagyl/ https://endmedicaldebt.com/zoloft/ https://transylvaniacare.org/product/ed-sample-pack-2/ https://cafeorestaurant.com/tamoxifen-without-pres/ https://frankfortamerican.com/ciprofloxacin-500-mg/ https://1485triclub.com/hydroxychloroquine/ https://endmedicaldebt.com/drugs/tadalafil/ https://1485triclub.com/cialis/ https://shilpaotc.com/tretinoin/ https://advantagecarpetca.com/nizagara-canada/ https://1485triclub.com/misoprost/ https://charlotteelliottinc.com/canadian-pharmacy/ https://rdasatx.com/drugs/clomid/ https://alliedentinc.com/erectafil/ https://1485triclub.com/finasteride/ https://driverstestingmi.com/doxycycline/ https://driverstestingmi.com/prednisone-online-no-script/ https://renog.org/product/flomax/ https://coastal-ims.com/drug/lasix/ https://petermillerfineart.com/lyrica/ https://primerafootandankle.com/item/bentyl/ https://heavenlyhappyhour.com/prednisone/ https://tacticaltrappingservices.com/pill/viagra/ https://thecultivarte.com/drugs/cipro/ cipro 1000mg original https://thecultivarte.com/wellbutrin/ https://tacticaltrappingservices.com/pharmacy/ ventures decompressed dryer.

FrankJScott 发表于 2024-9-4 01:20:43

New Tajir4D Login Site

FrankJScott ??? 2024-8-28 23:56
To the person asking about rtp live hoki77, jam main slot agar menang, rtp live idn, rtp live slot p ...

In reply to the guy inquiring about judi 2d, t0gel online, cara main 4d, new toto 4d, slot judi toto, slot4d terbaru, online bandar, toto 4d toto, judi 4d, online 4d toto,I highly recommend this this contact form about TAJIR4D info or toto 4d online malaysia, bandar togel idn, slot4d singapore, idn live 4d, judi online 4d, bandar judi togel online, agen play, bandar togel bri 24 jam, bandar slot singapore, 4d toto online, on top of this full article for Tajir4D Login blog as well as slot4d singapore, all toto 4d, toto 4d online malaysia, play toto online, online toto 4d, toto judi slot, 4d 4d toto, toto game online, bandar judi togel, 4d toto online, as well as this top TAJIR4D tips which is also great. Also, have a look at this high rated Tajir4D Login info alongside all new toto 4d, dewatogel slot online login, 4d online, bandar toto, 4d 4d toto, bandar slot 4d, master wap, 4d org, toto judi, toto malaysia online, alongside all this going here about TAJIR4D blog with play toto online, agen macau, judi sgp, dewatogel mobile, judi toto slot, bandar 4d slot, 4d casino, agen slot 2d, idn live 4d, bandar togel idn,for good measure. Check more @ Cool 7RAJATOGEL Info 3_8c134

FrankJScott 发表于 2024-9-4 02:37:19

Top ASIAN2BET Login Blog

DavidEmuts ??? 2024-8-31 23:51
??????? ??????? ??? ?????? ???? ??????????? ???? ??? ...

To the person inquiring about bandar togel bonus deposit, judi slot4d, 4d toto, website 4d, keren 4d slot login, chat 4d, cara main 4d singapore, agen play, 4d online login, play 4d online,I highly suggest this killer deal for ASIAN2BET link or idn 4d slot, bandar judi terpercaya di indonesia, web judi togel, toto play online, toto slot judi, toto 2d, toto judi, toto judi slot, bandar togel bonus deposit, agen slot 2d, not to mention this what google did to me on ASIAN2BET site not to mention toto login, cara toto, new toto 4d, master wap, toto play online, bandar 4d slot, 11 toto, bandar judi hongkong, online toto 4d, www 4d net m, alongside all this her response about ASIAN2BET Login advice which is also great. Also, have a look at this great ASIAN2BET forum as well as dewatogel mobile, sdy slot 4d, toto judi 4d login, toto judi slot, bandar togel bri 24 jam, toto 4d online malaysia, judi sgp, whatsapp 4d, online toto 4d, pembayaran 4d, as well as this excellent ASIAN2BET Login tips with toto 4d 12, master wap, judi 4d, withdraw 4d, slot4d terbaru, judi hk online, bandar togel bonus deposit 100, situs bandar judi slot online, agen slot 2d, website 4d,for good measure. Check more @ Top Rated ASIAN2BET Blog 4cab878

FrankJScott 发表于 2024-9-6 00:56:14

Top Rated Daily Sport Predictions Guide

To the lady asking about winning baseball mlb picks strategy, daily predictions for table tennis, h2h stats picks each week, soccer mls predictions for tonight, basketball nba picks for today, tennis betting promo codes, best tennis prediction strategy, up to date nfl bets news, best nfl prediction model, nfl & nba predictions for tomorrow,I highly recommend this best excellent sport prediction tips or ice hockey nhl predictions today, daily picks for tennis h2h, tennis prediction news today, winning soccer mls picks strategy, winning nfl tips picks today, winning american football nfl predictions, best mls prediction strategy, mls predictions for tonight, american football nfl predictions each week, tenis pronosticos manana, alongside all this published here for cool sport prediction advice alongside all tennis predictions maker, best tennis prediction strategy, expert nfl tips predictions against the spread, best nfl prediction model, winning camel racing picks today, mls picks for today, table tennis news today, camel racing predictions today straight up, best table tennis predictions platform, best camel racing prediction strategy, as well as this she said on high rated sport prediction info which is also great. Also, have a look at this go here for excellent sport prediction site not to mention h2h stats predictions, best baseball mlb picks model, best ice hockey nhl betting picks site, mls predictions for this week, tennis news today, basketball nba picks straight up, free daily nfl bets predictions, mls picks today, baseball mlb predictions this weekend, best basketball nba picks model, not to mention this the original source on daily sport prediction advice with best ice hockey nhl picks & parlays, best mls picks & parlays, expert tennis predictions against the spread, free daily nfl bets predictions, basketball nba predictions for tomorrow, profitable american football nfl picks, table tennis picks this weekend, table tennis rumours & news, tennis computer predictions, latest h2h stats news,for good measure. Check more @ Top CUANSLOT Info e7069c2
页: [1] 2 3 4
查看完整版本: 如何用NCL处理风云4A/MODIS卫星数据?