----------------------------- TO COMPILE MAKE SURE TO load module NetCDF/3.6.3 ----------------------------- All files are copied from /scratch1/portfolios/NCEPDEV/hwrf/save/Samuel.Trahan/waveonly-20140310-220344 created by Samuel Trahan On Mon, Mar 10, 2014 at 6:41 PM, Samuel Trahan - NOAA Affiliate wrote: Hello all, I finally have code and data for you. Sorry it took so long. It is in an HPSS archive due to the large size of the sample data (~2GB). /NCEPDEV/hpssuser/g01/wx20st/waveonly-20140310-220344.tar Most of the source code is a library that is part of the HWRF. The main file of interest to you is the "sample.f90" which is the main program. It is straightforward to run the program: htar -xvf /NCEPDEV/hpssuser/g01/wx20st/waveonly-20140310-220344.tar cd waveonly make cd data ./test.sh If you look in test.sh, you'll see that all it does is execute the program. If you want, you can also give it more OpenMP threads and it will run a bit faster. You'll get a few output files after you run it: native-file-3.dat - domain 3 native grid data in a binary format native-file-2.dat - domain 2 native grid data in a binary format native-file-1.dat - domain 1 native grid data in a binary format latlon.dat - lat-lon-interpolated data in a binary format hourly-track-3.atcf - hourly track and intensity hourly-track-2.atcf - duplicate of hourly-track-3.atcf in ATCF format (no analysis time MSLP though) You will probably want to modify the format of the various files to suit your model's pre-processor. The source code is in sample.f90: subroutine output_native_grid - outputs the native grid data. called once per forecast time (not for analysis time) subroutine output_latlon_grid - outputs the lat-lon grid data, called once per forecast time (not for analysis time) subroutine output_track - writes the hourly track file. Sincerely, Sam Trahan Yung, My cell phone number is 410-564-9592. I have a critical meeting with the hurricane center from 3:15-4:15 today, so please do not call during that time. Also, please do not give out the phone number. One other thing. In the program, for debugging, I added a function that can bring up a window with the domain data in it (in greyscale, up-side-down). To use it, uncomment the lines that call view_pgm. Also, to use it, you need to make sure you log in with X enabled (ssh -Y username@zeus.rdhpcs.noaa.gov) and have the $DISPLAY environment variable set. You can test that by running "xterm" on Zeus and see if an xterm window appears. Be warned that it is quite slow though, due to the network latencies between us and Zeus. Sincerely Sam Trahan On Wed, Mar 12, 2014 at 11:14 AM, Yung Chao - NOAA Federal wrote: Sam, Could you give me your phone number, just in case I need to talk to you. I have run your program as the way it is. It runs OK. But I want to see what is the output looks like so I am being modify your program. My home phone is 301-236-0699301-236-0699, Cell 443-864-8504443-864-8504. I am working from home today too. More cmmunication with Sam through email and personal contact for deguging his program. The latest one is: YY, The new archive with the fix for the bad lat-lon dims is here: /NCEPDEV/hpssuser/g01/wx20st/waveonly-20140318-164100.tar Sincerely, Sam Trahan