[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