AMSR IDL (read_amsr_day_v5.pro): .pro. : PRO example (means this script starts here and it s name is example.pro)

Transcription

1 IDL 2007/1/17~18, C202 - REMSS (free, ) http//wwwremsscom AMSR IDL (read_amsr_day_v5pro) HomePage -> AMSR -> Download Data (via ftpssmicom/amsre) -> FTP server -> support/idl/read_amsr_day_v5pro AMSR HomePage -> AMSR -> Download Data (via ftpssmicom/amsre) -> FTP server -> bmaps_v05 -> -> ( 3 ) IDL PRO filename (filename pro ) PRO example (means this script starts here and it s name is examplepro),, END ***!? IDL IDL Command Input window pro Run -> Compile filenamepro (Ctrl+F5) Run -> Run filename (F5), Run -> Stop (Ctrl+R) Output Log window,,, IDL 1

2 print, data (1), openw, 2, D\testtxt (2) printf, 2, data (3) close, 2 (1) openr, 2, D\testtxt (2) readf, 2, b FOR FOR parameter = start, end DO BEGIN IF ( ) THEN BEGIN, ENDIF ELSE BEGIN ENDELSE (1) array name = INTARR( ) a = INTARR(2,2) (2) array name = FLTARR( ) array name[flag] a[0]=1 (1), print, a[0] (2) print, a ***! IDL 2

3 AMSR HomePage -> AMSR -> See the data description for details The center of the first cell of the 1440 column and 720 row map is at 0125 E longitude and latitude The center of the second cell is 0375 E longitude, latitude *** AMSR ASCII IDL (read_amsr_day_v5pro) ( ) inclass_1_testpro (Compile All), (txt or asc) ,, AMSR 251 = missing sst or wind speed due to sun glint; missing sst or wind speed near sea ice; missing sst or wind speed due to rain; missing water vapor due to heavy rain; missing sst due to high wind speed (> 20 m/s) 252 = sea ice 253 = observations exist, but are bad (not used in composite maps) 254 = no observations 255 = land mass (1) GMT ps (2) HomePage -> AMSR -> Daily Passes (Dynamic Data Imaging) -> -> Update Display -> GMT GSview \\ \gmt\ 4 GMT (1) GMT netcdf C\ (2), netcdf path, ( path_gmt344txt) etcdf = c\netcdf path = c\gmt\bin;c\netcdf\bin Ghostscript GSview ( gs_view ) IDL 3

4 IDL ( GMT )! 2, 250, SST data inclass_2_testpro GMT GMT totalbat ( ) runbat ( ) runbat GMT ps ps GSview QSCAT HomePage -> QSCAT -> Archive of Past Tropical Storms -> HomePage -> QSCAT -> Download QuikScat Data (via ftp) -> bmaps_v03a/y2004/m06/ -> gz QSCAT HomePage -> QSCAT -> Download QuikScat Data (via ftp) -> scatterometer_bmap_support/idl/ -> get_scat_daily_v03pro inclass_3_testpro ( ) IDL, GMT ps IDL 4

5 AMSR IDL - ASCII, PRO amsre_ a ;*** determine file and file path *** path = 'D\IDL\' name = 'amsre_ v5gz' ;*** call subroutine and return geophysical data *** file_name = path + name read_amsr_day_v5, file_name, time, sst, wind, vapor, cloud, rain ;*** give the output file *** data_as = path + 'amsre_ _ascendingtxt' data_ds = path + 'amsre_ _descendingtxt' OPENW, 5, data_as OPENW, 6, data_ds ;*** output data *** ;0 ascending (Daytime) ;1 descending (Nighttime) PRINTF, 5, ' Longitude Latitude PRINTF, 6, ' Longitude Latitude SST UTC time UTC time(fractional hours)' SST UTC time UTC time(fractional hours)' FOR ilat=420,500 DO BEGIN ;15125N ~ 35125N FOR ilon = 440,560 DO BEGIN ;110125E ~ E lon = * ilon lat = * ilat mt1 = FIX(time[ilon,ilat,*]) ;float to integer mt2 = FIX((time[ilon,ilat,*] - mt1) * 60 ) ; ex 148 hour -> 1448 PRINTF, 5, lon, lat, sst[ilon,ilat,0], mt1[0],'', mt2[0], time[ilon,ilat,0], FORMAT='(3F123,I92,A1,I02,F155)' PRINTF, 6, lon, lat, sst[ilon,ilat,1], mt1[1],'', mt2[1], time[ilon,ilat,1], FORMAT='(3F123,I92,A1,I02,F155)' CLOSE, /ALL END IDL 5

6 AMSR IDL -, PRO amsre_ b ;*** determine file and file path *** path = 'D\IDL\' name = 'amsre_ v5gz' ;*** call subroutine and return geophysical data *** file_name = path + name read_amsr_day_v5, file_name, time, sst, wind, vapor, cloud, rain ;*** give the output file *** data_as = path + 'amsre_ _ascendingasc' data_ds = path + 'amsre_ _descendingasc' data_as_miss = path + 'amsre_ _ascending_missasc' data_ds_miss = path + 'amsre_ _descending_missasc' OPENW, 5, data_as OPENW, 6, data_ds OPENW, 7, data_as_miss OPENW, 8, data_ds_miss ;*** output data *** ;0 ascending (Daytime) ;1 descending (Nighttime) miss = -999 FOR ilat=420,500 DO BEGIN ;15125N ~ 35125N FOR ilon = 440,560 DO BEGIN ;110125E ~ E lon = * ilon lat = * ilat IF (sst[ilon,ilat,0] LE 250) THEN BEGIN PRINTF, 5, lon, lat, sst[ilon,ilat,0] ENDIF ELSE BEGIN PRINTF, 7, lon, lat, miss ENDELSE IF (sst[ilon,ilat,1] LE 250) THEN BEGIN PRINTF, 6, lon, lat, sst[ilon,ilat,1] ENDIF ELSE BEGIN PRINTF, 8, lon, lat, miss ENDELSE CLOSE, /ALL END IDL 6

7 AMSR GMT <totalbat> makecpt -Crainbow -T2/32/1 -Z > acpt surface %1asc -G%1_sstgrd -I005 -R115/135/15/30 grdimage %1_sstgrd -R115/135/15/30 -Jq125/04 -Cacpt -X15 -Y18 -K > %1_sstps grdcontour %1_sstgrd -C1 -A2f9 -G3 -R115/135/15/30 -Jq125/04 -W3/255/255/255/255 -K-O >> %1_sstps psxy %1_missasc -R115/135/15/30 -Jq125/04 -G70/70/70 -K -O -Ss02 >> %1_sstps pscoast -R115/135/15/30 -Jq125/04 -Ba5f5g5/a5f5g5WeSn -Di -W2/0/0/0 -Ii -K -O >> %1_sstps psscale -Cacpt -D4/-055/10/014h -B"AMSR-SST data %1" -I+05 -L -K -O >> %1_sstps del %1_sstgrd 1 color bar makecpt ( cpt) C-, rainbow ; T- / / 2 ( ) surface ( asc grd) G- ; I- ( ); R- / / / ( ); 3 grd grdimage ( ps) J- (q) / C- ; X-X ( ); Y-Y ( ) 4 grdcontour ( ps) C- ; A- (f ); G- ; W- / / / 5 missing data psxy (, ps) G- / / ; S- (s ) 6 pscoast B- X /Y (a, f, q ) (,, WeSn ); D- (i ); W- / / ; Ii- 7 color bar psscale D-bar X / Y / / ; B- ; I-?; L-? 8 grd <runbat> call total amsre_ _ascending call total amsre_ _descending IDL 7

8 AMSR GMT PRO inclass_3_test ;*********************************************************************************************** u=fltarr(1440,720,2) v=fltarr(1440,720,2) ;*********************************************************************************************** path='d\idl\' filename=path+' gz' Get_Scat_Daily_v03,filename, mingmt,windspd,winddir,scatflag,radrain ;*********************************************************************************************** good=where(winddir NE -999) angle=!pi*(winddir[good]/180) ;convert degrees to!pi (!PI=314159) u[good]=windspd[good]*sin(angle) ;u and v has minus value where towards the West and South v[good]=windspd[good]*cos(angle) ;calculate wind vector ;see "http//wwwremsscom/qscat/qscat_descriptionhtml" for wind direction ;note wind speed are all positive values ;*********************************************************************************************** U_as=path+'U_ _astxt' U_ds=path+'U_ _dstxt' V_as=path+'V_ _astxt' V_ds=path+'V_ _dstxt' miss_as=path+'miss_ _astxt' miss_ds=path+'miss_ _dstxt' wpd_as=path+'wpd_ _astxt' wpd_ds=path+'wpd_ _dstxt' slg_as=path+'slg_ _astxt' slg_ds=path+'slg_ _dstxt' OPENW,2,U_as OPENW,3,U_ds OPENW,4,V_as OPENW,5,V_ds OPENW,6,miss_as OPENW,7,miss_ds OPENW,8,wpd_as OPENW,9,wpd_ds OPENW,10,slg_as OPENW,11,slg_ds miss=-999 FOR ilat=400,480 DO BEGIN ;10125N ~ 30125N FOR ilon=480,600 DO BEGIN ;120125E ~ E lon= *ilon lat= *ilat IDL 8

9 IF(windspd[ilon,ilat,0] NE -999) THEN BEGIN PRINTF,2, lon, lat, u[ilon,ilat,0], FORMAT='(3F153)' PRINTF,4, lon, lat, v[ilon,ilat,0], FORMAT='(3F153)' PRINTF,8, lon, lat, windspd[ilon,ilat,0], FORMAT='(3F153)' PRINTF,10, lon, lat, scatflag[ilon,ilat,0], FORMAT='(2F153,I12)' ENDIF ELSE BEGIN PRINTF,6, lon, lat, miss ENDELSE IF(windspd[ilon,ilat,1] NE -999) THEN BEGIN PRINTF,3, lon, lat, u[ilon,ilat,1], FORMAT='(3F153)' PRINTF,5, lon, lat, v[ilon,ilat,1], FORMAT='(3F153)' PRINTF,9, lon, lat, windspd[ilon,ilat,1], FORMAT='(3F153)' PRINTF,11, lon, lat, scatflag[ilon,ilat,1], FORMAT='(2F153,I12)' ENDIF ELSE BEGIN PRINTF,7, lon, lat, miss ENDELSE ;*********************************************************************** CLOSE, /ALL END IDL 9