[gradsusr] Aggregating multiple data files (Sara Fenech)
Sestak, Michael L CIV 63134
michael.l.sestak at navy.mil
Mon Jun 2 10:16:16 EDT 2014
I think you need to look into the %ch template option and how to use CHSUB in your control file in order to use files with this sort of date ranges.
http://www.iges.org/grads/gadoc/templates.html
Michael Sestak
Fleet Numerical Meteorology and Oceanography Center
Monterey, CA
________________________________________
From: gradsusr-bounces at gradsusr.org [gradsusr-bounces at gradsusr.org] on behalf of gradsusr-request at gradsusr.org [gradsusr-request at gradsusr.org]
Sent: Monday, June 02, 2014 06:22
To: gradsusr at gradsusr.org
Subject: gradsusr Digest, Vol 52, Issue 3
Send gradsusr mailing list submissions to
gradsusr at gradsusr.org
To subscribe or unsubscribe via the World Wide Web, visit
http://gradsusr.org/mailman/listinfo/gradsusr
or, via email, send a message with subject or body 'help' to
gradsusr-request at gradsusr.org
You can reach the person managing the list at
gradsusr-owner at gradsusr.org
When replying, please edit your Subject line so it is more specific
than "Re: Contents of gradsusr digest..."
Today's Topics:
1. Aggregating multiple data files (Sara Fenech)
2. Scale changing (Felipe Correa)
----------------------------------------------------------------------
Message: 1
Date: Mon, 2 Jun 2014 13:05:11 +0200
From: Sara Fenech <sarafenech at gmail.com>
Subject: [gradsusr] Aggregating multiple data files
To: gradsusr at gradsusr.org
Message-ID:
<CAL2KrZ5jZ5+18MuziORivsR9bYY240Z9vq0AVXFzFmt504y4Gw at mail.gmail.com>
Content-Type: text/plain; charset="utf-8"
Dear all
I am trying to aggregate 4 .dat files using a data descriptor file using a
single DSET entry and setting option to template. I have the following 4
.dat files of hourly data for the corresponding months (2185 datasets in
all)
*dom2jun.dat* (from 00Z01Jun2011 to 23Z30Jun2011) *(720 files)*
*dom2jul.dat* (from 00Z01Jul2011 to 23Z30Jul2011) *(720 files)*
*dom2aug.dat * (from 00Z31Jul2011 to 23Z29Aug2011) *(720 files)*
*dom2sept.dat* (from 00Z30Aug2011 to 00Z31Aug2011) *(25 files)*
Using the following data descriptor file (*dom2JJA.ctl*) I need to open all
the .dat files one after the other so as to act as a single file and thus
can obtain different parameters for all the time range at once. The only
difference between the .dat files is that the time range varies.
The problem is that when I am opening this descriptor file in GrAds and I
try to access different parameters, the values are exponentials and do not
agree with the values I extract when opening the files individually.
Is there anyone who can spot what I'm doing wrong in my descriptor file
please? Thanks a lot and I appreciate all your help
Kind Regards
Sara
dset ^dom2%mc.dat
options byteswapped
undef 1.e30
title OUTPUT FROM WRF V3.2.1 MODEL
pdef 380 170 lcc 38.429 18.325 190.500 85.500 36.50000 36.50000
16.50000 9600.000 9600.000
xdef 1068 linear -4.69692 0.04324324
ydef 399 linear 28.84252 0.04324324
zdef 27 linear 1 1
tdef 2185 linear 00Z01JUN2011 60MN
VARS 84
LU_INDEX 1 0 LAND USE CATEGORY (-)
U 27 0 x-wind component (m s-1)
V 27 0 y-wind component (m s-1)
W 27 0 z-wind component (m s-1)
PH 27 0 perturbation geopotential (m2 s-2)
PHB 27 0 base-state geopotential (m2 s-2)
T 27 0 perturbation potential temperature (theta-t0) (K)
MU 1 0 perturbation dry air mass in column (Pa)
MUB 1 0 base state dry air mass in column (Pa)
NEST_POS 1 0 - (-)
P 27 0 perturbation pressure (Pa)
PB 27 0 BASE STATE PRESSURE (Pa)
P_HYD 27 0 hydrostatic pressure (Pa)
Q2 1 0 QV at 2 M (kg kg-1)
T2 1 0 TEMP at 2 M (K)
TH2 1 0 POT TEMP at 2 M (K)
PSFC 1 0 SFC PRESSURE (Pa)
U10 1 0 U at 10 M (m s-1)
V10 1 0 V at 10 M (m s-1)
QVAPOR 27 0 Water vapor mixing ratio (kg kg-1)
QCLOUD 27 0 Cloud water mixing ratio (kg kg-1)
QRAIN 27 0 Rain water mixing ratio (kg kg-1)
QICE 27 0 Ice mixing ratio (kg kg-1)
QSNOW 27 0 Snow mixing ratio (kg kg-1)
QGRAUP 27 0 Graupel mixing ratio (kg kg-1)
LANDMASK 1 0 LAND MASK (1 FOR LAND, 0 FOR WATER) (-)
TSLB 2 0 SOIL TEMPERATURE (K)
SMOIS 2 0 SOIL MOISTURE (m3 m-3)
SH2O 2 0 SOIL LIQUID WATER (m3 m-3)
SMCREL 2 0 RELATIVE SOIL MOISTURE (-)
SEAICE 1 0 SEA ICE FLAG (-)
XICEM 1 0 SEA ICE FLAG (PREVIOUS STEP) (-)
SFROFF 1 0 SURFACE RUNOFF (mm)
UDROFF 1 0 UNDERGROUND RUNOFF (mm)
IVGTYP 1 0 DOMINANT VEGETATION CATEGORY (-)
ISLTYP 1 0 DOMINANT SOIL CATEGORY (-)
VEGFRA 1 0 VEGETATION FRACTION (-)
GRDFLX 1 0 GROUND HEAT FLUX (W m-2)
ACGRDFLX 1 0 ACCUMULATED GROUND HEAT FLUX (J m-2)
SNOW 1 0 SNOW WATER EQUIVALENT (kg m-2)
SNOWH 1 0 PHYSICAL SNOW DEPTH (m)
CANWAT 1 0 CANOPY WATER (kg m-2)
SST 1 0 SEA SURFACE TEMPERATURE (K)
SSTSK 1 0 SKIN SEA SURFACE TEMPERATURE (K)
LAI 1 0 Leaf area index (area/area)
MAPFAC_M 1 0 Map scale factor on mass grid (-)
MAPFAC_MX 1 0 Map scale factor on mass grid, x direction (-)
MAPFAC_MY 1 0 Map scale factor on mass grid, y direction (-)
MF_VX_INV 1 0 Inverse map scale factor on v-grid, x direction (-)
F 1 0 Coriolis sine latitude term (s-1)
E 1 0 Coriolis cosine latitude term (s-1)
SINALPHA 1 0 Local sine of map rotation (-)
COSALPHA 1 0 Local cosine of map rotation (-)
HGT 1 0 Terrain Height (m)
TSK 1 0 SURFACE SKIN TEMPERATURE (K)
RAINC 1 0 ACCUMULATED TOTAL CUMULUS PRECIPITATION (mm)
RAINNC 1 0 ACCUMULATED TOTAL GRID SCALE PRECIPITATION (mm)
SNOWNC 1 0 ACCUMULATED TOTAL GRID SCALE SNOW AND ICE (mm)
GRAUPELNC 1 0 ACCUMULATED TOTAL GRID SCALE GRAUPEL (mm)
HAILNC 1 0 ACCUMULATED TOTAL GRID SCALE HAIL (mm)
CLDFRA 27 0 CLOUD FRACTION (-)
SWDOWN 1 0 DOWNWARD SHORT WAVE FLUX AT GROUND SURFACE (W m-2)
GLW 1 0 DOWNWARD LONG WAVE FLUX AT GROUND SURFACE (W m-2)
OLR 1 0 TOA OUTGOING LONG WAVE (W m-2)
XLAT 1 0 LATITUDE, SOUTH IS NEGATIVE (degree_north)
XLONG 1 0 LONGITUDE, WEST IS NEGATIVE (degree_east)
ALBEDO 1 0 ALBEDO (-)
ALBBCK 1 0 BACKGROUND ALBEDO (-)
EMISS 1 0 SURFACE EMISSIVITY (-)
NOAHRES 1 0 RESIDUAL OF THE NOAH SURFACE ENERGY BUDGET (W m{-2})
TMN 1 0 SOIL TEMPERATURE AT LOWER BOUNDARY (K)
XLAND 1 0 LAND MASK (1 FOR LAND, 2 FOR WATER) (-)
UST 1 0 U* IN SIMILARITY THEORY (m s-1)
PBLH 1 0 PBL HEIGHT (m)
HFX 1 0 UPWARD HEAT FLUX AT THE SURFACE (W m-2)
QFX 1 0 UPWARD MOISTURE FLUX AT THE SURFACE (kg m-2 s-1)
LH 1 0 LATENT HEAT FLUX AT THE SURFACE (W m-2)
ACHFX 1 0 ACCUMULATED UPWARD HEAT FLUX AT THE SURFACE (J m-2)
ACLHF 1 0 ACCUMULATED UPWARD LATENT HEAT FLUX AT THE SURFACE (J
m-2)
SNOWC 1 0 FLAG INDICATING SNOW COVERAGE (1 FOR SNOW COVER) (-)
SR 1 0 fraction of frozen precipitation (-)
POTEVP 1 0 accumulated potential evaporation (W m-2)
SNOPCX 1 0 snow phase change heat flux (W m-2)
SOILTB 1 0 bottom soil temperature (K)
ENDVARS
@ global String comment TITLE = OUTPUT FROM WRF V3.2.1 MODEL
@ global String comment START_DATE = 2011-06-01_00:00:00
@ global String comment SIMULATION_START_DATE = 2011-06-01_00:00:00
@ global String comment WEST-EAST_GRID_DIMENSION = 381
@ global String comment SOUTH-NORTH_GRID_DIMENSION = 171
@ global String comment BOTTOM-TOP_GRID_DIMENSION = 28
@ global String comment DX = 9600.00
@ global String comment DY = 9600.00
@ global String comment GRIDTYPE = C
@ global String comment DIFF_OPT = 1
@ global String comment KM_OPT = 4
@ global String comment DAMP_OPT = 0
@ global String comment DAMPCOEF = 0.20
@ global String comment KHDIF = 0.00
@ global String comment KVDIF = 0.00
@ global String comment MP_PHYSICS = 6
@ global String comment RA_LW_PHYSICS = 1
@ global String comment RA_SW_PHYSICS = 1
@ global String comment SF_SFCLAY_PHYSICS = 7
@ global String comment SF_SURFACE_PHYSICS = 7
@ global String comment BL_PBL_PHYSICS = 7
@ global String comment CU_PHYSICS = 1
@ global String comment SURFACE_INPUT_SOURCE = 1
@ global String comment SST_UPDATE = 0
@ global String comment GRID_FDDA = 0
@ global String comment GFDDA_INTERVAL_M = 0
@ global String comment GFDDA_END_H = 0
@ global String comment GRID_SFDDA = 0
@ global String comment SGFDDA_INTERVAL_M = 0
@ global String comment SGFDDA_END_H = 0
@ global String comment SF_URBAN_PHYSICS = 0
@ global String comment FEEDBACK = 0
@ global String comment SMOOTH_OPTION = 0
@ global String comment SWRAD_SCAT = 1.00
@ global String comment W_DAMPING = 0
@ global String comment MOIST_ADV_OPT = 1
@ global String comment SCALAR_ADV_OPT = 1
@ global String comment TKE_ADV_OPT = 1
@ global String comment DIFF_6TH_OPT = 0
@ global String comment DIFF_6TH_FACTOR = 0.12
@ global String comment OBS_NUDGE_OPT = 0
@ global String comment BUCKET_MM = -1.00
@ global String comment BUCKET_J = -1.00
@ global String comment PREC_ACC_DT = 0.00
@ global String comment OMLCALL = 0
@ global String comment ISFTCFLX = 0
@ global String comment WEST-EAST_PATCH_START_UNSTAG = 1
@ global String comment WEST-EAST_PATCH_END_UNSTAG = 380
@ global String comment WEST-EAST_PATCH_START_STAG = 1
@ global String comment WEST-EAST_PATCH_END_STAG = 381
@ global String comment SOUTH-NORTH_PATCH_START_UNSTAG = 1
@ global String comment SOUTH-NORTH_PATCH_END_UNSTAG = 170
@ global String comment SOUTH-NORTH_PATCH_START_STAG = 1
@ global String comment SOUTH-NORTH_PATCH_END_STAG = 171
@ global String comment BOTTOM-TOP_PATCH_START_UNSTAG = 1
@ global String comment BOTTOM-TOP_PATCH_END_UNSTAG = 27
@ global String comment BOTTOM-TOP_PATCH_START_STAG = 1
@ global String comment BOTTOM-TOP_PATCH_END_STAG = 28
@ global String comment GRID_ID = 2
@ global String comment PARENT_ID = 1
@ global String comment I_PARENT_START = 16
@ global String comment J_PARENT_START = 18
@ global String comment PARENT_GRID_RATIO = 5
@ global String comment DT = 48.00
@ global String comment CEN_LAT = 38.43
@ global String comment CEN_LON = 18.33
@ global String comment TRUELAT1 = 36.50
@ global String comment TRUELAT2 = 36.50
@ global String comment MOAD_CEN_LAT = 36.50
@ global String comment STAND_LON = 16.50
@ global String comment POLE_LAT = 90.00
@ global String comment POLE_LON = 0.00
@ global String comment GMT = 0.00
@ global String comment JULYR = 2011
@ global String comment JULDAY = 152
@ global String comment MAP_PROJ = 1
@ global String comment MMINLU = USGS
@ global String comment NUM_LAND_CAT = 24
@ global String comment ISWATER = 16
@ global String comment ISLAKE = -1
@ global String comment ISICE = 24
@ global String comment ISURBAN = 1
@ global String comment ISOILWATER = 14
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20140602/d0979433/attachment-0001.html
------------------------------
Message: 2
Date: Mon, 2 Jun 2014 10:22:22 -0300
From: Felipe Correa <felipe.nascimento.correa at gmail.com>
Subject: [gradsusr] Scale changing
To: gradsusr at gradsusr.org
Message-ID:
<CAMUN9_5ZyZ-mVX0c0sYOGDUbDgkCpw-MrzNb+Ai1VXBkL3MGUg at mail.gmail.com>
Content-Type: text/plain; charset="utf-8"
Hello guys! Cold morning for us here in the south hemisphere.
I need help changing a scale on grads. I have data using pressure levels
and I need to show it using altitude (meters or feet). How can I do that?
Thanks for all!
--
Felipe do Nascimento Corr?a
T?cnico em Meteorologia - ACA-m/ACA/IAE/
Departamento de Ci?ncia e Tecnologia Aeroespacial (DCTA)
For?a ?rea Brasileira - "Asas que protegem o pa?s"
Tel: (12) 3947-4548
Cel: (12) 9-9717-1058
E-mail1: felipefnc at iae.cta.br
E-mail2: felipe.nascimento.correa at gmail.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20140602/5277c547/attachment.html
------------------------------
_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org
http://gradsusr.org/mailman/listinfo/gradsusr
End of gradsusr Digest, Vol 52, Issue 3
***************************************
More information about the gradsusr
mailing list