NVAP-Next Generation (NVAP-NG) Readme 1.0 Introduction This readme file contains information about the NASA WATER VAPOR PROJECT NEXT GENERATION (NVAP-NG) data sets. The NVAP-NG data are available through the Langley Atmospheric Sciences Data Center (ASDC) as the following 13 data sets: NVAP_NG_CLW_12HOUR NVAP_NG_CLW_ANNUAL NVAP_NG_CLW_DAILY NVAP_NG_CLW_MONTHLY NVAP_NG_LWP_12HOUR NVAP_NG_LWP_ANNUAL NVAP_NG_LWP_DAILY NVAP_NG_LWP_MONTHLY NVAP_NG_PWC_12HOUR NVAP_NG_PWC_ANNUAL NVAP_NG_PWC_DAILY NVAP_NG_PWC_MONTHLY NVAP_NG_PWC_STATS This readme file includes the following six other sections: Section 2.0 - Data Set Description Section 3.0 - Format and Packaging Section 4.0 - Science Parameters Information Section 5.0 - Description and Implementation of Sample Read Software Section 6.0 - Sample Output Section 7.0 - Additional Information If users have questions while using the NVAP-NG sample read software, please contact the NASA Langley ASDC Science, Users and Data Services Office at: Atmospheric Sciences Data Center Science, Users and Data Services Office Mail Stop 157D 2 South Wright Street NASA Langley Research Center Hampton, Virginia 23681-2199 U.S.A. E-mail: larc@eos.nasa.gov Phone: (757)864-8656 FAX: (757)864-8807 URL: http://eosweb.larc.nasa.gov 2.0 Data Set Description 2.1 Data Set Information The NVAP Next Generation (NVAP-NG) data sets are the latest phase of NVAP, covering years 2000 and 2001. Built on the successful foundation of NVAP, the NVAP-NG products are a result of increased spatial and temporal resolution. The NVAP-NG data sets are global products at 1/2 degree and twice daily resolution with 5 vertical layers. These data sets are made up of packed integer*2 (16 bit) records after an ASCII header of 144 bytes (8 bits). Each file can hold many global, 0.5 degree x 0.5 degree grids [720x360]. Coordinate [1,1] is at the North Pole and 0 degrees longitude. For information about the NVAP & NVAP-NG Products see: http://www.cira.colostate.edu/Climate/NVAP/climate_wvsci.html Scientific information, results and discussions are found in the "Science and Technology Corp. Technical Report 3333" at: http://eosweb.larc.nasa.gov/PRODOCS/nvap/sci_tech_report_3333.pdf 2.2 File Naming Convention The NVAP-NG has been catagorized into 13 data sets (see Section 1.0 for complete names). There are 3 types of data recorded during the 4 time periods (12 HOUR, DAILY, MONTHLY, and ANNUAL), and a Statistics file. CLW Cloud Liquid Water LWP Liquid Water Path PWC Precipitable Water Column PWC_STATS PWC Statistics 2.3 Science Point of Contact Prof. Thomas Vonder Haar Princinal Investigator CIRA Colorado State University Ft. Collins, CO 80523 E-mail: vonderhaar@cira.colostate.edu Dr. David L. Randel CIRA Colorado State University Ft. Collins, CO 80523 E-mail: randel@cira.colostate.edu 3.0 Format and Packaging The following files are defined for this Sample Read Software distribution: Filename Format ----------------------- ---------------- nvap_all_2ascii.f Fortran Source Code nvap_all_glbave.f Fortran Source Code nvap_all_idl.pro IDL Source Code nvap_all_listheaders.f Fortran Source Code nvap_all_sample.f Fortran Source Code nvap_all_subroutines.f Fortran Source Code nvap_ng_read_ds.f Fortran Source Code nvap_ng_removefill.f Fortran Source Code nvap_ng_resize.f Fortran Source Code nvap_ng_readme This readme file 3.1 Zip Information To read the files, first unzip them with Info-Zip software. Information about downloading and using this freeware is available at the URL: http://www.info-zip.org/pub/infozip/Info-ZIP.html 3.2 File Naming Conversion The names of the NVAP-NG data files referenced in the "Science and Technology Corp. (STS) Technical Report 3333", have been renamed for conformity and ease of use at the ASDC. http://eosweb.larc.nasa.gov/PRODOCS/nvap/sci_tech_report_3333.pdf 3.2.1 The naming conversion used for these files are: nvap_ng_ggg_rrr_yyyymm.nat where: ggg gridded product clw Cloud Liquid Water lwp Liquid Water Path rrr recording period 12hr 12 hour daily daily mon monthly ann annual yyyy year (4-digit) mm month (2-digit) for example: nvap_ng_clw_daily_200001.nat This example is a NVAP-NG Cloud Liquid Water file with a daily recording for January 2000 in native format. 3.2.2 Annual file names do not contain the month field as part of the date. for example: nvap_ng_lwp_ann_2000.nat This example is a NVAP-NG Liquid Water Path file with an annual recording for the year 2000 in native format. 3.2.3 Precipitable Water Column (PWC) file names contain a file type specifier before the date: nvap_ng_ggg_rrr_tt_yyyymm.nat where: ggg gridded product pwc Precipitable Water Column pwc_stats PWC Statistics rrr recording period 12hr 12 hour daily daily mon monthly ann annual tt file type specifier dslwv data source for lwv dsq data source for level q dstc data source for tc lwv layer water vapor tc total column precipitable water q level q statistic yyyy year (4-digit) mm month (2-digit) for example: nvap_ng_pwc_daily_dstc_200001.nat This example is a NVAP-NG Precipitable Water Column file with a daily data source total column recordings for January 2000 in native format. 3.3 Browse Images Currently, there are a total of 168 Browse images in JPEG format displaying monthly averages of: Liquid Water Path ( 24 images) Total Column Water Vapor ( 24 images) Layered Water Vapor (120 images) L1: 850-1000mb L2: 700-850mb L3: 500-700mb L4: 300-500mb L5: 200-300mb * five layers per data month * 4.0 Science Parameters Information Parameters used for each Data Type in the NVAP-NG ASDC archive __________________________________________________________________________________________ Data Types: WV (Water Vapor), LW (Liquid Water), CW (Cloud Water), SH (Specific Humidity) Param Data Number Type Data Set Data File Ext Layer Range Units ------------------------------------------------------------------------------------------- 1 WV NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_tc_YYYYMM.nat TOTAL 0-11Z mm 2 WV NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_tc_YYYYMM.nat TOTAL 12-23Z mm 3 WV NVAP_NG_PWC_DAILY nvap_ng_pwc_daily_tc_YYYYMM.nat TOTAL 0-23Z mm 4 LW NVAP_NG_LWP_12HOUR nvap_ng_lwp_12hr_YYYYMM.nat TOTAL 0-11Z g/kg 5 LW NVAP_NG_LWP_12HOUR nvap_ng_lwp_12hr_YYYYMM.nat TOTAL 12-23Z g/kg 6 LW NVAP_NG_LWP_DAILY nvap_ng_lwp_daily_YYYYMM.nat TOTAL 0-23Z g/kg 7 CW NVAP_NG_CLW_12HOUR nvap_ng_clw_12hr_YYYYMM.nat TOTAL 0-11Z g/kg 8 CW NVAP_NG_CLW_12HOUR nvap_ng_clw_12hr_YYYYMM.nat TOTAL 12-23Z g/kg 9 CW NVAP_NG_CLW_DAILY nvap_ng_clw_daily_YYYYMM.nat TOTAL 0-23Z g/kg 12 SH NVAP_NG_PWC_DAILY nvap_ng_pwc_daily_q_YYYYMM.nat 300 hPa 0-23Z g/kg 13 SH NVAP_NG_PWC_DAILY nvap_ng_pwc_daily_q_YYYYMM.nat 500 hPa 0-23Z g/kg 14 SH NVAP_NG_PWC_DAILY nvap_ng_pwc_daily_q_YYYYMM.nat 700 hPa 0-23Z g/kg 15 SH NVAP_NG_PWC_DAILY nvap_ng_pwc_daily_q_YYYYMM.nat 850 hPa 0-23Z g/kg 16 SH NVAP_NG_PWC_DAILY nvap_ng_pwc_daily_q_YYYYMM.nat 1000 hPa 0-23Z g/kg 18 SH NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_q_YYYYMM.nat 300 hPa 0-11Z g/kg 19 SH NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_q_YYYYMM.nat 500 hPa 0-11Z g/kg 20 SH NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_q_YYYYMM.nat 700 hPa 0-11Z g/kg 21 SH NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_q_YYYYMM.nat 850 hPa 0-11Z g/kg 22 SH NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_q_YYYYMM.nat 1000 hPa 0-11Z g/kg 24 SH NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_q_YYYYMM.nat 300 hPa 12-23Z g/kg 25 SH NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_q_YYYYMM.nat 500 hPa 12-23Z g/kg 26 SH NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_q_YYYYMM.nat 700 hPa 12-23Z g/kg 27 SH NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_q_YYYYMM.nat 850 hPa 12-23Z g/kg 28 SH NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_q_YYYYMM.nat 1000 hPa 12-23Z g/kg 41 WV NVAP_NG_PWC_DAILY nvap_ng_pwc_daily_lwv_YYYYMM.nat 200-300 hPa 0-23Z mm 42 WV NVAP_NG_PWC_DAILY nvap_ng_pwc_daily_lwv_YYYYMM.nat 300-500 hPa 0-23Z mm 43 WV NVAP_NG_PWC_DAILY nvap_ng_pwc_daily_lwv_YYYYMM.nat 500-700 hPa 0-23Z mm 44 WV NVAP_NG_PWC_DAILY nvap_ng_pwc_daily_lwv_YYYYMM.nat 700-850 hPa 0-23Z mm 45 WV NVAP_NG_PWC_DAILY nvap_ng_pwc_daily_lwv_YYYYMM.nat 850-1000 hPa 0-23Z mm 46 WV NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_lwv_YYYYMM.nat 200-300 hPa 0-11Z mm 47 WV NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_lwv_YYYYMM.nat 300-500 hPa 0-11Z mm 48 WV NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_lwv_YYYYMM.nat 500-700 hPa 0-11Z mm 49 WV NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_lwv_YYYYMM.nat 700-850 hPa 0-11Z mm 50 WV NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_lwv_YYYYMM.nat 850-1000 hPa 0-11Z mm 51 WV NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_lwv_YYYYMM.nat 200-300 hPa 12-23Z mm 52 WV NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_lwv_YYYYMM.nat 300-500 hPa 12-23Z mm 53 WV NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_lwv_YYYYMM.nat 500-700 hPa 12-23Z mm 54 WV NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_lwv_YYYYMM.nat 700-850 hPa 12-23Z mm 55 WV NVAP_NG_PWC_12HOUR nvap_ng_pwc_12hr_lwv_YYYYMM.nat 850-1000 hPa 12-23Z mm ------------------------------------------------------------------------------------------- Data source Bit : add 100 to parameter number (12hour and daily grids) Data Population : add 200 to parameter number (tc(1-9) & 500hPa grids only) Data Range : add 300 to parameter number (12hour grids only) Standard Deviation: add 400 to parameter number (12hour grids only) Time test GRID : parm=991, 992 ___________________________________________________________________________________________ For more information about the NVAP-NG science parameters see: http://www.cira.colostate.edu/Climate/NVAP/climate_wvsci.html or view the "Science and Technology Corp. Technical Report 3333" at: http://eosweb.larc.nasa.gov/PRODOCS/nvap/sci_tech_report_3333.pdf 5.0 Description and Implementation of Sample Read Software 5.1 Interactive Data Language (IDL) source file IDL applications can only be run from within the IDL environment, so anyone who wishes to run IDL applications must have access to an IDL license. This NVAP-NG software package contains one IDL read software program used to visualize CIRA Climate Data Archive (CCDA) formatted grid fields. It is: nvap_all_idl.pro 5.1.1 Compile and execute using the interactive IDLDE Invoke IDL by typing "idlde" at the prompt. Select "Open" from the File menu. Select "nvap_all_idl.pro" from the Files window, click "OK". Select "Compile nvap_all_idl.pro" from the Run menu. Select "Run nvap_all_idl" from the Run menu. A pop-up window will appear. Select "OPEN and select or enter a data file (*nat). The Filter window will need to be modified to the path of the data files you would like to plot. Select a parameter from the listing and IDL will display the image. 5.1.2 Compile and execute at the command line NOTE The > character in front of a command means it's a unix command line prompt. Invoke IDL by typing "idl" at the prompt. > idl [then enter return] IDL Version 5.3 (IRIX mipseb). (c) 1999, Research Systems, Inc. Installation number: 17659. Licensed for use by: NASA - Langley Research Center IDL> nvap_all_idl [ enter this name ] % Compiled module: NVAP_NG_IDL. % Compiled module: CONGRID. % Compiled module: XMANAGER. A window will appear at the top left corner of the screen. Click the OPEN FILE button to select a filename. Select a file from the Files window. Click the ok button. Select the data field you wish to plot. Click the Return button to return back to the field window. From this window select Quit to return to the Quit window. Click the Quit button to exit the program. 5.2 FORTRAN77 source files All of these routines were developed in the FORTRAN77 language. These programs have been tested on the Sun4 and the SGI UNIX based workstations. All of these programs require the nvap_all_subroutines.f program to be compiled along with it. The nvap_all_subroutines.f program contains the read routines and global average routines needed by the other seven programs. 5.2.1 Creating FORTRAN77 Source Executable File To build the executable file for all the NVAP-NG source programs the user can compile the FORTRAN77 version of the software by invoking the FORTRAN77 compiler f77 program_name at the command line prompt. 5.2.2 Compile for Sun (nvap_all_glbave.f used as an example for all programs) To get these programs to run on a sun you must edit the nvap_all_subroutines.f program. Locate the line with "subroutine nvap_read(flname,igrd,istat)" in it. It should be around line 103. Then scan down to line 161 which reads: irecl = 720 / 4 ! bytes / 4 (this may need to be just bytes) change it to read: irecl = 720 Make the above changes to line 166 also. Locate at or near line 232 : Subrountine nvap_byteswap(bytebuff, isize) and uncomment this line. Compile and execute the program. f77 -o read_nvap nvap_all_glbave.f nvap_all_subroutines.f read_nvap (where read_nvap can be any name you select) 5.2.3 Compile for SGI f77 -o read_nvap nvap_all_glbave.f nvap_all_subroutines.f read_nvap (where read_nvap can be any name you select) 5.2.4 Invoking the Executable File This section illustrates the invocation of the read_nvap executable file. Section 4.1 describes how to run the read_nvap program for the NVAP-NG data sets. To execute the program, enter executable file name at the prompt. You will then be prompted to enter a data file name: %read_nvap %Input filename: nvap_ng_clw_12hr_200001.nat fld# source param strt date end date n.hemis s.hemis global %GOOD 1 900 7 2000/001 00 2000/001 11 0.87 0.87 0.87 100.0 2 900 8 2000/001 12 2000/001 23 0.87 0.87 0.87 100.0 3 900 7 2000/002 00 2000/002 11 0.88 0.88 0.88 100.0 4 900 8 2000/002 12 2000/002 23 0.86 0.86 0.86 100.0 5 900 7 2000/003 00 2000/003 11 0.88 0.88 0.88 100.0 etc... 6.0 Sample Output This NVAP-NG software package contains eight separate FORTRAN77 programs. They are: nvap_all_2ascii.f nvap_all_glbave.f nvap_all_listheaders.f nvap_ng_read_ds.f nvap_ng_removefill.f nvap_ng_resize.f nvap_all_sample.f nvap_all_subroutines.f 6.1 Program nvap_all_2ascii.f: ----------------------------- THIS ROUTINE PRINTS OUT A CCDA FORMATTED FIELD TO AN ASCII TEXT FILE. f77 -w -o read_nvap nvap_all_2ascii.f nvap_all_subroutines.f read_nvap - executable from being compiled To run it type read_nvap at the prompt: %read_nvap program prompts for the input file: Input ccda filename: nvap_ng_clw_12hr_200001.nat Output ascii filename: out do all fields in file? (y/n) [y] : y 900 7 2000 1 0 2000 1 11 720 360 900 8 2000 1 12 2000 1 23 720 360 900 7 2000 2 0 2000 2 11 720 360 900 8 2000 2 12 2000 2 23 720 360 900 7 2000 3 0 2000 3 11 720 360 900 8 2000 3 12 2000 3 23 720 360 etc... 6.2 Program nvap_all_glbave.f: ----------------------------- THIS ROUTINE READS THE CCDA FORMATTED GRID FILES AND WRITES OUT A SIMPLE GLOBAL AVERAGE OF ALL THE POINTS. f77 -w -o read_nvap nvap_all_glbave.f nvap_all_subroutines.f read_nvap - executable from being compiled To run it type read_nvap at the prompt: %read_nvap program prompts for input file: Input filename: nvap_ng_clw_12hr_200001.nat fld# source param strt date end date n.hemis s.hemis global %GOOD 1 900 7 2000/001 00 2000/001 11 0.87 0.87 0.87 100.0 2 900 8 2000/001 12 2000/001 23 0.87 0.87 0.87 100.0 3 900 7 2000/002 00 2000/002 11 0.88 0.88 0.88 100.0 4 900 8 2000/002 12 2000/002 23 0.86 0.86 0.86 100.0 5 900 7 2000/003 00 2000/003 11 0.88 0.88 0.88 100.0 etc... 6.3 Program nvap_all_listheaders.f: ---------------------------------- THIS ROUTINE LIST ALL OF THE HEADERS FOUND ON A ccda FORMATTED GRID FILE. f77 -w -o read_nvap nvap_all_listheaders.f nvap_all_subroutines.f read_nvap - executable from being compiled To run it type read_nvap at the prompt: %read_nvap program prompts for the input file: Input filename: nvap_ng_clw_12hr_200001.nat FLD SOURCE PARAM START TIME END TIME GRID SIZE LABEL 1 900 7 2000 1 0 2000 1 11 720 360 NVAP-NG CLOUD LIQUID WATER 0-11Z 2 900 8 2000 1 12 2000 1 23 720 360 NVAP-NG CLOUD LIQUID WATER 12-23Z 3 900 7 2000 2 0 2000 2 11 720 360 NVAP-NG CLOUD LIQUID WATER 0-11Z 4 900 8 2000 2 12 2000 2 23 720 360 NVAP-NG CLOUD LIQUID WATER 12-23Z 5 900 7 2000 3 0 2000 3 11 720 360 NVAP-NG CLOUD LIQUID WATER 0-11Z etc... 6.4 Program nvap_ng_read_ds.f: ------------------------------ THIS ROUTINE CONVERTS DATA SOURCE BITMAP FILES TO CCDA FORMAT. f77 -w -o read_nvap nvap_ng_read_ds.f nvap_all_subroutines.f read_nvap - executable from being compiled To run it type read_nvap at the prompt: %read_nvap program prompts for the input file: Input Data Source Bitmap filename? : nvap_ng_clw_12hr_200001.nat Output CCDA formatted filename? : out finished writing 16 grids for field : 1 - " finished writing 16 grids for field : 2 - output is written to the screen finished writing 16 grids for field : 3 - " The contents of the output file out: ccda900 73 2020 0 1 0 0 111 720 360 0.50 -0.50 89.750 0.2500.15000E+030.15259E-02-.32768E+05 DSB SSM/I Leftovers binary data ccda900 73 2020 0 1 0 0 111 720 360 0.50 -0.50 89.750 0.2500.15000E+030.15259E-02-.32768E+05 DSB SSM/T2 Sat 2 binary data etc... 6.5 Program nvap_ng_removefill.f: --------------------------------- THIS ROUTINE TAKES THE DATA SOURCE BITMAP FILE INFORMATION AND REMOVES SPATIAL AND TEMPORAL DATA FROM THE CCDA GRID. f77 -w -o read_nvap nvap_ng_removefill.f nvap_all_subroutines.f read_nvap - executable from being compiled To run it type read_nvap at the prompt: %read_nvap program prompts for input files: Input CCDA filename: nvap_ng_pwc_daily_tc_200001.nat Input DSB filename: nvap_ng_pwc_daily_dstc_200001.nat Output Filename: out USE ALL FIELDS IN FILE? (Y/N) [Y] : y 2000 1 NVAP-NG WATER VAPOR 0-23Z 2000 2 NVAP-NG WATER VAPOR 0-23Z 2000 3 NVAP-NG WATER VAPOR 0-23Z 2000 4 NVAP-NG WATER VAPOR 0-23Z 2000 5 NVAP-NG WATER VAPOR 0-23Z etc... 6.6 Program nvap_ng_resize.f: ----------------------------- THIS ROUTINE CONVERTS 1/2 DEGREE FILE INTO A 1 DEGREE FILE. f77 -w -o read_nvap nvap_ng_read_resize.f nvap_all_subroutines.f read_nvap - executable from being compiled: To run it type read_nvap at the prompt: %read_nvap program prompts for the input file: Input 1/2 degree filename? : nvap_ng_clw_12hr_200001.nat Output 1 degree filename? : out Resize all fields in file? (y/n) [y] : y 2000 1 resized field written to: out 2000 1 resized field written to: out 2000 2 resized field written to: out 2000 2 resized field written to: out etc... The contents of the output file out: ccda900 73 2020 0 1 0 0 111 360 180 1.00 -1.00 89.500 0.5000.85123E+000.22815E-04-.32768E+05 NVAP-NG CLOUD LIQUID WATER 0-11Z binary data 0 73 2020 0 1 0 0 111 360 180 1.00 -1.00 89.500 0.5000.85123E+000.22815E-04-.32768E+05 NVAP-NG CLOUD LIQUID WATER 0-11Z binary data etc ... 6.7 Program nvap_all_sample.f: ----------------------------- THIS ROUTINE READS A CCDA FORMATTED GRID FILE AND PERFORMS A SIMPLE (NO WEIGHTING) GRID AVERAGE. f77 -w -o read_nvap nvap_all_sample.f nvap_all_subroutines.f read_nvap - executable from being compiled: To run it type read_nvap at the prompt: %read_nvap program prompts for the input file: Input filename: nvap_ng_clw_12hr_200001.nat grid size = 720 360 simple (no weighting) grid ave = 1 0.8705708 simple (no weighting) grid ave = 2 0.8701813 simple (no weighting) grid ave = 3 0.8733956 simple (no weighting) grid ave = 4 0.8598542 simple (no weighting) grid ave = 5 0.8744755 etc... 6.8 Program nvap_all_subroutines.f: ---------------------------------- THIS PROGRAM CONTAINS SUBROUTINES CALLED BY THE OTHER PROGRAMS. IT MUST BE COMPILED ALONG WITH ALL THE OTHER FORTRAN ROUTINES. see section 5.2.1 for details * this file must be edited * 7.0 Additional Information 7.1 Contributions from: * Thomas H. Vonder Haar * Principal Investigator John M. Forsythe Darren McKague David L. Randel Benjamin C. Ruston Shannon Woo Science and Technology Corporation Fort Collins, Colorado STC-METSAT A Division of Science and Technology Corp. 226 RemingtonStreet, Fort Collins, CO 80524 Tel: (970) 221-5420; Fax: (970) 493-3410 7.2 Last Updated June 1, 2004