[gradsusr] question to Jayakrishnan.P.R on opening MODIS hdf files
Jayakrishnan P.R
prjayakrishnan at gmail.com
Wed Oct 24 10:49:16 EDT 2012
Dear Sagar,
I herewith attach the programs as well as the method adopted
by me to plot swath data. Even though it is a very difficult task, you can
write all these commands into a single shell script. I had intended to
write a recipe in opengrads cookbooks about this thing...But I could not
do......If any more queries about this thing then kindly ask me.......
JKPR...
*Hello,
Here I am describing how to process MODIS HDF files using grads station
format. For this purpose I have attached some executables which will be
useful for the analysis
Firstly navigate to the folder in which the hdf files and other executables
are placed. Here comes the demonstration.... Go through it. *
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
jkpr at jkpr-desktop:~/Desktop/Modis data download/monsoon 2007/2007 10 th
july 20 45 IST$ ./ggg MAM07S0.A2007191.2040.002.2007194090750.hdf
Number of SDS arrays in file: 30
1) Latitude dimensions = 408 41
2) Longitude dimensions = 408 41
3) Scan_Start_Time dimensions = 408 41
4) Solar_Zenith dimensions = 408 41
5) Solar_Azimuth dimensions = 408 41
6) Sensor_Zenith dimensions = 408 41
7) Sensor_Azimuth dimensions = 408 41
8) Brightness_Temperature dimensions = 12 408 41
9) Cloud_Mask dimensions = 408 41
10) Surface_Temperature dimensions = 408 41
11) Surface_Pressure dimensions = 408 41
12) Surface_Elevation dimensions = 408 41
13) Processing_Flag dimensions = 408 41
14) Tropopause_Height dimensions = 408 41
15) Guess_Temperature_Profile dimensions = 20 408 41
16) Guess_Moisture_Profile dimensions = 20 408 41
17) Retrieved_Temperature_Profile dimensions = 20 408 41
18) Retrieved_Moisture_Profile dimensions = 20 408 41
19) Retrieved_Height_Profile dimensions = 20 408 41
20) Total_Ozone dimensions = 408 41
21) Total_Totals dimensions = 408 41
22) Lifted_Index dimensions = 408 41
23) K_Index dimensions = 408 41
24) Water_Vapor dimensions = 408 41
25) Water_Vapor_Direct dimensions = 408 41
26) Water_Vapor_Low dimensions = 408 41
27) Water_Vapor_High dimensions = 408 41
28) Quality_Assurance dimensions = 408 41 10
29) Quality_Assurance_Infrared dimensions = 408 41 5
30) Subset Starting Frame Indices 5km dimensions = 412
Enter total number of parameters to write out or 0 to exit program
3
Enter parameter numbers from above list separated by white space
1 2 22
The valid parameter numbers you Enter are: 1 2 22
****************************************
SDS name = Latitude
SDS type = 5
SDS rank = 2
SDS scale = 1.000000
SDS dims = 408 41
sds file MAM07S0.A2007191.2040.002.2007194090750.Latitude.
Binary file successfully created : MAM07S0.A2007191.2040.002.
2007194090750.Latitude
****************************************
SDS name = Longitude
SDS type = 5
SDS rank = 2
SDS scale = 1.000000
SDS dims = 408 41
sds file MAM07S0.A2007191.2040.002.2007194090750.Longitude.
Binary file successfully created : MAM07S0.A2007191.2040.002.
2007194090750.Longitude
****************************************
SDS name = Lifted_Index
SDS type = 22
SDS rank = 2
SDS scale = 0.010000
SDS dims = 408 41
sds file MAM07S0.A2007191.2040.002.2007194090750.Lifted_Index .
Binary file successfully created : MAM07S0.A2007191.2040.002.
2007194090750.Lifted_Index
*Now incorperate the output files
MAM07S0.A2007191.2040.002.2007194090750.Lifted_Index.
,MAM07S0.A2007191.2040.002.2007194090750.Latitude
,MAM07S0.A2007191.2040.002.2007194090750.Longitude to the fortran program
named test.f to get the output file result.txt, which will combine the
binary files to make an ascii file.
Now run the output of test.f which is result.txt to the input of the
fortran gradskutty.f, and we will get the output as passmosissep.txt
The output passmodis.txt can be used as the input of stnkutty.c, which is a
c program to generate the station binary file named passmodis.bin
After successful creation of the station binary file run the station map
utility using the command
stnmap -i stnexample.ctl, to create the station map
Then open the grads
open stnexample.ctl
set lat -10 45
set lon 60 90
set gxout stnmark
set cmark 11
d li
cbarn
This will give the required pass of the satellite in coloured form.
*
--
Sincerely
***********************************************
Jayakrishnan.P.R
CSIR-Senior Research Fellow
Department of Atmospheric Sciences
Cochin University of Science and Technology (CUSAT), Cochin-682 016
Kerala, India.
Mob: 09895417565
***********************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20121024/fcdddf2b/attachment-0003.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: fwdprocedurestoplothdfformattedmodisdatausingstati.zip
Type: application/zip
Size: 1149975 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20121024/fcdddf2b/attachment-0003.zip
More information about the gradsusr
mailing list