Scribd
Upload
630 views7 pages

Gradsusr Moist Con Transp

1. The user is asking how to calculate vertically integrated moisture transport over a specific longitude and latitude region for selected years using GrADS. 2. A sample GrADS script is provided that calculates vertically integrated moisture transport using pressure, zonal wind, meridional wind, and specific humidity data. 3. However, the user has some questions about specifics in the script, including whether msl refers to pressure data, if shum refers to specific humidity, and if their time specification may cause issues.

Uploaded by

akhmadfadholi
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOC, PDF, TXT or read online on Scribd
630 views7 pages

Gradsusr Moist Con Transp

1. The user is asking how to calculate vertically integrated moisture transport over a specific longitude and latitude region for selected years using GrADS. 2. A sample GrADS script is provided that calculates vertically integrated moisture transport using pressure, zonal wind, meridional wind, and specific humidity data. 3. However, the user has some questions about specifics in the script, including whether msl refers to pressure data, if shum refers to specific humidity, and if their time specification may cause issues.

Uploaded by

akhmadfadholi
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOC, PDF, TXT or read online on Scribd
You are on page 1/ 7

[gradsusr] how to use vint function in

GrADS
sushant puranik sushantpuranik at gmail.com
Fri Jul 23 01:34:01 EDT 2010
 Previous message: [gradsusr] how to use vint function in GrADS
 Next message: [gradsusr] how to use vint function in GrADS

 Messages sorted by: [ date ] [ thread ] [ subject ] [ author ]

Thank you for suggestion. Actually i want to calculate following integral.

1000

Q = (1/g) ∫ (V q)dp.

400

I used this kind of command in grads

'define ps=1000'
'define a=vint(ps,vwnd*shum,400)'
'd a/9.81'
with this command i am able to make a plot.

whether it is correct or not ?

Thanking you

Sushant

On Fri, Jul 23, 2010 at 10:38 AM, MY <funnytyphoon at gmail.com> wrote:

> Hi,
>
> please try
>
> 'define ps=1000'
> 'd vint(ps,vwnd*shum,400)'
>
> best wishes,
> ++++++++++++++++++++++++++++++++++++++++++++++++++++
> YING, Ming
> Shanghai Typhoon Institute of China Meteorological Administration
> 166 Puxi Road, Xuhui District, Shanghai, 200030
> PR China
> ++++++++++++++++++++++++++++++++++++++++++++++++++++
>
>
>
>
>
> 2010/7/23 sushant puranik <sushantpuranik at gmail.com>:
> > Hi
> > I am trying to find vertical integral in GrADS using "vint" command.
For
> > that purpose i opened vwnd and specific humidity file as first and
second
> > respectively. The q file command for both the variable are as
> >
> > ga-> q file
> > File 1 : mean daily NMC reanalysis (2008)
> > Descriptor: /home/Sushant/Desktop/10may08v.nc
> > Binary: /home/Sushant/Desktop/10may08v.nc
> > Type = Gridded
> > Xsize = 29 Ysize = 25 Zsize = 8 Tsize = 1
> > Number of Variables = 1
> > vwnd 8 -999 mean Daily V wind
> >
> > ga-> q file 2
> > File 2 : mean daily NMC reanalysis (2008)
> > Descriptor: /home/Sushant/Desktop/10may08q.nc
> > Binary: /home/Sushant/Desktop/10may08q.nc
> > Type = Gridded
> > Xsize = 29 Ysize = 25 Zsize = 8 Tsize = 1
> > Number of Variables = 1
> > shum 8 -999 mean Daily specific humidity
> >
> > after this i used
> >
> > ga-> d vint(1000,vwnd*shum.2,400)
> > gxstart=0, gxstop=29, gystart=0, gystop=25.
> > gzstart=0, gzstop=1, gtstart=0, gtstop=1.
> > gxstart=0, gxstop=29, gystart=0, gystop=25.
> > gzstart=0, gzstop=1, gtstart=0, gtstop=1.
> > Error from VINT: Incompatible grids.
> > Contouring: -0 to 2.4e+36 interval 3e+35
> > Segmentation fault
> >
> > after showing this message, the grads window gets closed.
> >
> > Can anyone has solution for this perticular problem.
> >
> > Thanks in advance
> >
> >
> > --
> > Sushant
> >
[gradsusr] Vertical Integrated Moisture
Transport (VIMT)
Eric Altshuler ela at cola.iges.org
Mon May 12 17:47:29 EDT 2014
 Previous message: [gradsusr] Vertical Integrated Moisture Transport
(VIMT)
 Next message: [gradsusr] I need shapefile of Congo basin

 Messages sorted by: [ date ] [ thread ] [ subject ] [ author ]

Hi Muhammad,

I think there are two issues here, but the problem might not be caused by
either of them. First, you close file 1 after defining a variable. I don't
think this is a good idea, since it leaves grads without a default file and
the subsequent behavior of grads might be strange or incorrect (I ran into
something like this in the past). Similarly, if there are several files
open, the current default file should not be closed, for the same reason
(e.g. if you have 'set dfile 2', don't close file 2 until you've set the
default file to something other than 2). In your script, you should leave
file 1 open and use 2,3,4 to refer to the files you open later.

There is also an error in the calculation. The vertically integrated


moisture transport is NOT equal to the product of the vertical integrals of
wind and moisture. You need to calculate it this way:

'define zon=vint(surfpres,shum.2*uwnd.3,500)'
'define mer=vint(surfpres,shum.2*vwnd.4,500)'

Best regards,

Eric L. Altshuler
Research Scientist
Center for Ocean-Land-Atmosphere Studies
253 Research Hall, Mail Stop 6C5
George Mason University
4400 University Drive
Fairfax, VA 22030 USA

E-mail: ela at cola.iges.org


Phone: (703) 993-5725
Fax: (703) 993-5770

----- Original Message -----


From: "Muhammad Latif" <latif_space at yahoo.com>
To: gradsusr at gradsusr.org
Sent: Friday, May 9, 2014 4:06:19 PM
Subject: [gradsusr] Vertical Integrated Moisture Transport (VIMT)

Hi all,

I am trying to calculate Vertical Integrated Moisture Transport (VIMT) from


surface to 500 hPa in units of kg/m/sec using grads function
"vint(psexpr,expr,top)" only on the one time step.
I used the following script to get VIMT:

'sdfopen spr-daily.nc' ;* spr is the surface pressure in Pa.

'set time 20jun1970'


'define surfpres=pres/100' ;* /100 to convert in hPa
'close 1'

'sdfopen shum-daily.nc'
'sdfopen uwnd-daily.nc'
'sdfopen vwnd-daily.nc'

'set lat 0 45'


'set lon 50 100'

'set time 20jun1970'


'define vintq=vint(surfpres,shum.1,500)' ;*shum is the specific humidity in
kg/kg
'define vintu=vint(surfpres,uwnd.2,500)' ;*uwnd is the u component of wind
in m/sec
'define vintv=vint(surfpres,vwnd.3,500)' ;*vwnd is the v component of wind
in m/sec
'define zon=(vintq*vintu)'
'define mer=(vintq*vintv)'

'set gxout vector'


'd zon;mer'

Script runs successfully showing VIMT, but the magnitude of vector appears
(4e+06) in the vector scale bar.
I guess, this magnitude is not right. The value should lie around 20 or 40
(kg/m/sec).

I'll be very thankful if anyone could help me to sort out the issue
mentioned above.

Regards

Muhammad Latif
PhD Scholar
Department of Meteorology
CIIT, Islamabad

[gradsusr] Moisture flux


convergence/divergence
Jeffrey Duda jdduda at iastate.edu
Thu Sep 15 18:34:50 EDT 2011
 Previous message: [gradsusr] Moisture flux convergence/divergence
 Next message: [gradsusr] Accumulated precipitation script

 Messages sorted by: [ date ] [ thread ] [ subject ] [ author ]

Set the graphics output to vector ('set gxout vector'), then 'd zon;mer'
(note the semicolon instead of a comma).

Jeff Duda

2011/9/15 Hyacinth Nnamchi <hyacinth.1 at hotmail.com>

> Hello,
>
> This is not directly about GrADS but since most of address similar
> problems, I guess someone on this forum will be able to answer my
question.
>
> I understand that to calculate moisture transport, given zonal,
meridional
> winds and specific humidity, can simply do:
>
> 'define zon=uwnd.1*shum.3'
> 'define mer=vwnd.2*shum.3'
> 'd zon,mer'
>
> But, I want to plot vectors indicating field of net moisture flux (all
arrows converging/diverging) instead of just transport.
>
> Thanks for any ideas.
>
> Hyacinth
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>

--
Jeff Duda
Iowa State University
Meteorology Graduate Student
www.meteor.iastate.edu/~jdduda <http://www.meteor.iastate.edu/%7Ejdduda>
-------------- next part --------------
[gradsusr] moisture flux
ugo okoro uknac23 at yahoo.com
Thu Jan 3 10:10:22 EST 2013
 Previous message: [gradsusr] Built-in functions
 Next message: [gradsusr] Error in using skew.gs

 Messages sorted by: [ date ] [ thread ] [ subject ] [ author ]

Dear jennifer,

i want to compute the vertically integrated moisture transport over lon


(25W - 15E)
and lat (0 - 20N) at selected years (say 1994,1999,2003 and 2008 from NCEP
monthly)
using this script i saw online;

'sdfopen file-pressure.nc'
'sdfopen file-zonal-wind.nc'
'sdfopen file-meridional-wind.nc'
'define QiJAN=ave(vint(prmsl.1(z=1)/100,shum.1*uwnd.2,200),t=((t=181)+
(t=241)+(t=289)+(t=349))'
'define QjJAN=ave(vint(prmsl.1(z=1)/100,shum.1*vwnd.3,200),t=((t=181)+
(t=241)+(t=289)+(t=349))'
'set gxout vector'
'd QiJAN;QjJAN' but i dont get some details from the script-
1. is msl the pressure data?
2. is shum the specific humidity?
3. will my time be the problem with the script?

Thanks in anticipation.

ugo okoro
cmsr, iap-cas,
beijing

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy