SST-T500の水平分布

$ srcdump.sh sst-t500.up.gs sst-t500.gs circlat.gs circlon.gs tile.gs qdims.gsf plev.ctl sfc.ctl

aofd165.bio.mie-u.ac.jp

/work1/am/2016.PolarLow/ERA-I/SST-T500

Tue Sep 20 17:08:17 JST 2016

======================

sst-t500.up.gs

======================

'open sfc.ctl'

'open plev.ctl'

'q ctlinfo 2'

say result

header='ERA-I_Barents110116-110123'

headout='ERA-I_Barents_UP'

varname='SST-T500'

outdir='Fig'

'!mkdir -p 'outdir

#'set mproj latlon'

'set mproj nps'

'set lat 60 88'

'set lon -90 90'

i=1

while (i <= 32 )

'set t ' i

'cc'

'set mpdset mres'

'set frame off'

'set grid off' ;#on 3 15'

'q dims'

line=sublin(result,5)

datetime=subwrd(line,6)

'set strsiz 0.2 0.25'

'set string 1 c 5'

'draw string 4.5 8.5 'datetime

'mul 1 4 1 1 -yini 1'

var1=sst.1'-273'

'set gxout shade2'

'set grads off'

'color -1 18 2 -kind white->rainbow -xcbar 1 1.1 0.5 2.7 -fh 0.12 -fw 0.11 -edge triangle -line on'

'draw string 0.9 2.8 degC'

'draw string 1.8 0.6 SST'

'set clab off'

'd maskout(' var1 ',' var1 '+3)'

'set gxout contour'

'set cthick 3'

'set ccolor 0'

'set cmin -1'

'set cmax 18'

'set cint 2'

'd maskout(' var1 ',' var1 '+3)'

'circlat 10 90 0.1'

'circlon 10 0.1'

'tile ' ci.1 ' 0.2 1'

'set gxout contour'

levs1='960 965 970 975 980 985 990 995 1000 1005 1010 1015 1020 1025 1030'

'set clevs 'levs1

'set clab off'

'set cthick 10'

'set ccolor 0'

'd ' msl.1 '/100'

'color -gxout contour -kind purple->thistle -levs ' levs1 ;#' -xcbar 7.5 7.6 6.0 8.2 -fh 0.12 -fw 0.11 -edge triangle -line on'

# 'draw string 7.5 8.3 hPa'

'set cthick 5'

'd ' msl.1 '/100'

# 'set clevs 990 1000'

# 'set clab forced'

# 'set cstyle 0'

# 'set ccolor 0'

# 'set cthick 0'

# 'set clopts 1 3 0.08'

# 'color -gxout contour -kind purple->thistle -levs 990 1000'

# 'd ' msl.1 '/100'

plev=500

'mul 1 4 1 2 -yini 0.1'

var2='t.2(lev=' plev ')-273'

'set gxout shade2'

'set grads off'

'color -50 -30 5 -kind blue->white -xcbar 1 1.1 3.3 5.5 -fh 0.12 -fw 0.11 -edge triangle -line on'

'draw string 0.9 5.6 degC'

'draw string 1.8 3.4 T'plev

'set clab off'

'd 'var2

'set gxout contour'

'set cthick 3'

'set ccolor 0'

'set cmin -40'

'set cmax -20'

'set cint 5'

'set cstyle 1'

'd 'var2

'set gxout contour'

levs1='960 965 970 975 980 985 990 995 1000 1005 1010 1015 1020 1025 1030'

'set clevs 'levs1

'set clab off'

'set cthick 10'

'set ccolor 0'

'd ' msl.1 '/100'

'color -gxout contour -kind purple->thistle -levs ' levs1 ;#' -xcbar 7.5 7.6 6.0 8.2 -fh 0.12 -fw 0.11 -edge triangle -line on'

# 'draw string 7.5 8.3 hPa'

'set cthick 5'

'd ' msl.1 '/100'

# 'set clevs 990 1000'

# 'set clab forced'

# 'set cstyle 0'

# 'set ccolor 0'

# 'set cthick 0'

# 'set clopts 1 3 0.08'

# 'd ' msl.1 '/100'

'set gxout contour'

'set cthick 10'

'set ccolor 0'

'set clevs -39'

'set cstyle 1'

'circlat 20 90 0.1'

'circlon 10 0.1'

'mul 1 4 1 3 -yini 0.1'

dvar=var1 '+273-' t.2 '(lev=' plev ')'

'set gxout shade2'

'set grads off'

'color 25 40 5 -kind white->yellow->red -xcbar 1 1.1 6.0 8.2 -fh 0.12 -fw 0.11 -edge triangle -line on'

'draw string 0.9 8.3 degC'

'draw string 1.8 6.1 SST-T'plev

'set clab off'

'd maskout(' dvar ',' var1 '+3)'

'set gxout contour'

'set cthick 3'

'set ccolor 0'

'set cmin 20'

'set cmax 40'

'set cint 5'

'd maskout(' dvar ',' var1 '+3)'

'tile ' ci.1 ' 0.2 1 -type 1'

'set gxout contour'

levs1='960 965 970 975 980 985 990 995 1000 1005 1010 1015 1020 1025 1030'

'set clevs 'levs1

'set clab off'

'set cthick 10'

'set ccolor 0'

'd ' msl.1 '/100'

'color -gxout contour -kind purple->thistle -levs ' levs1 ' -xcbar 7.5 7.6 2.3 6.5 -fh 0.12 -fw 0.11 -edge triangle -line on'

'draw string 7.5 6.6 hPa'

'set cthick 5'

'd ' msl.1 '/100'

# 'set clevs 990 1000'

# 'set clab forced'

# 'set cstyle 0'

# 'set ccolor 0'

# 'set cthick 0'

# 'set clopts 1 3 0.08'

# 'd ' msl.1 '/100'

'circlat 20 90 0.1'

'circlon 10 0.1'

num=''

if (i<10)

num='00'i

endif

if (i>=10 & i<100)

num='0'i

endif

if (i>100)

num=i

endif

figfile=outdir '/' headout '_' varname '_' num '_' datetime '.png'

say 'Output: ' figfile

'gxprint ' figfile

i=i+1

endwhile

'allclose'

#'quit'

----------------------

End of sst-t500.up.gs

----------------------

======================

sst-t500.gs

======================

'open sfc.ctl'

'open plev.ctl'

'q ctlinfo 2'

say result

header='ERA-I_Barents110116-110123'

headout='ERA-I_Barents'

varname='SST-T500'

outdir='Fig'

'!mkdir -p 'outdir

#'set mproj latlon'

'set mproj nps'

'set lat 35 88'

'set lon -90 90'

i=1

while (i <= 32 )

'set t ' i

'cc'

'set mpdset mres'

'set frame off'

'set grid off' ;#on 3 15'

'q dims'

line=sublin(result,5)

datetime=subwrd(line,6)

'set strsiz 0.2 0.25'

'set string 1 c 5'

'draw string 4.5 8.5 'datetime

'mul 1 4 1 1 -yini 1'

var1=sst.1'-273'

'set gxout shade2'

'set grads off'

'color -1 18 2 -kind white->rainbow -xcbar 1 1.1 0.5 2.7 -fh 0.12 -fw 0.11 -edge triangle -line on'

'draw string 0.9 2.8 degC'

'draw string 1.8 0.6 SST'

'set clab off'

'd maskout(' var1 ',' var1 '+3)'

'set gxout contour'

'set cthick 3'

'set ccolor 0'

'set cmin -1'

'set cmax 18'

'set cint 2'

'd maskout(' var1 ',' var1 '+3)'

'circlat 10 90 0.1'

'circlon 10 0.1'

'tile ' ci.1 ' 0.2 1'

'set gxout contour'

levs1='960 970 980 990 1000 1010 1020 1030'

'set clevs 'levs1

'set clab off'

'set cthick 10'

'set ccolor 0'

'd ' msl.1 '/100'

'color -gxout contour -kind purple->thistle -levs ' levs1 ;#' -xcbar 7.5 7.6 6.0 8.2 -fh 0.12 -fw 0.11 -edge triangle -line on'

# 'draw string 7.5 8.3 hPa'

'set cthick 6'

'd ' msl.1 '/100'

# 'set clevs 990 1000'

# 'set clab forced'

# 'set cstyle 0'

# 'set ccolor 0'

# 'set cthick 0'

# 'set clopts 1 3 0.08'

# 'color -gxout contour -kind purple->thistle -levs 990 1000'

# 'd ' msl.1 '/100'

plev=500

'mul 1 4 1 2 -yini 0.1'

var2='t.2(lev=' plev ')-273'

'set gxout shade2'

'set grads off'

'color -40 -20 5 -kind blue->white -xcbar 1 1.1 3.3 5.5 -fh 0.12 -fw 0.11 -edge triangle -line on'

'draw string 0.9 5.6 degC'

'draw string 1.8 3.4 T'plev

'set clab off'

'd 'var2

'set gxout contour'

'set cthick 3'

'set ccolor 0'

'set cmin -40'

'set cmax -20'

'set cint 5'

'set cstyle 1'

'd 'var2

'set gxout contour'

levs1='960 970 980 990 1000 1010 1020 1030'

'set clevs 'levs1

'set clab off'

'set cthick 10'

'set ccolor 0'

'd ' msl.1 '/100'

'color -gxout contour -kind purple->thistle -levs ' levs1 ;#' -xcbar 7.5 7.6 6.0 8.2 -fh 0.12 -fw 0.11 -edge triangle -line on'

# 'draw string 7.5 8.3 hPa'

'set cthick 6'

'd ' msl.1 '/100'

# 'set clevs 990 1000'

# 'set clab forced'

# 'set cstyle 0'

# 'set ccolor 0'

# 'set cthick 0'

# 'set clopts 1 3 0.08'

# 'd ' msl.1 '/100'

'circlat 20 90 0.1'

'circlon 10 0.1'

'mul 1 4 1 3 -yini 0.1'

dvar=sst.1 '-' t.2'(lev=' plev ')'

'set gxout shade2'

'set grads off'

'color 20 40 5 -kind white->yellow->red -xcbar 1 1.1 6.0 8.2 -fh 0.12 -fw 0.11 -edge triangle -line on'

'draw string 0.9 8.3 degC'

'draw string 1.8 6.1 SST-T'plev

'set clab off'

'd maskout(' dvar ',' var1 '+3)'

'set gxout contour'

'set cthick 3'

'set ccolor 0'

'set cmin 20'

'set cmax 40'

'set cint 5'

'd maskout(' dvar ',' var1 '+3)'

'tile ' ci.1 ' 0.2 1 -type 1'

'set gxout contour'

levs1='960 970 980 990 1000 1010 1020 1030'

'set clevs 'levs1

'set clab off'

'set cthick 10'

'set ccolor 0'

'd ' msl.1 '/100'

'color -gxout contour -kind purple->thistle -levs ' levs1 ' -xcbar 7.5 7.6 2.3 6.5 -fh 0.12 -fw 0.11 -edge triangle -line on'

'draw string 7.5 6.6 hPa'

'set cthick 6'

'd ' msl.1 '/100'

# 'set clevs 990 1000'

# 'set clab forced'

# 'set cstyle 0'

# 'set ccolor 0'

# 'set cthick 0'

# 'set clopts 1 3 0.08'

# 'd ' msl.1 '/100'

'circlat 20 90 0.1'

'circlon 10 0.1'

num=''

if (i<10)

num='00'i

endif

if (i>=10 & i<100)

num='0'i

endif

if (i>100)

num=i

endif

figfile=outdir '/' headout '_' varname '_' num '_' datetime '.png'

say 'Output: ' figfile

'gxprint ' figfile

i=i+1

endwhile

'allclose'

#'quit'

----------------------

End of sst-t500.gs

----------------------

======================

circlat.gs

======================

function circlat(args)

* usage: run circlat dlat (strsiz (font))

* mproj=nps もしくは sps において 緯度の数字を描く

* dlat は数字を描く経度の間隔で必須．

* lon は数字を描く経度で必須．

* strsiz は文字の大きさで必須．

* font はフォント番号でオプション

dlat=subwrd(args,1)

if (dlat=0); say '*circlat* dlat=0'; 'quit'; return; endif

lon0=subwrd(args,2)

strsiz=subwrd(args,3)

if (strsiz=""); say '*circlat* strsiz is not defined'; 'quit'; return; endif

if (strsiz=0); say '*circlat* strsiz=0'; 'quit'; return; endif

font=subwrd(args,4)

'q gxinfo'

cline=sublin(result,6)

mproj=subwrd(cline,3)

if (mproj!=3 & mproj!=4)

say ' *circlon* (mproj!=3:nps & mproj !=4:sps) is not supported'

say 'mproj='%mproj

'quit'

endif

if (mproj=3);* mproj=3 is nps

'query w2xy 0 90'

* say 'result='result

endif

if (mproj=4);* mproj=4 is sps

'query w2xy 0 -90'

endif

xpol=subwrd(result,3)

ypol=subwrd(result,6)

'q dims'

lonline=sublin(result,2)

lon1=subwrd(lonline,6)

lon2=subwrd(lonline,8)

latline=sublin(result,3)

lat1=subwrd(latline,6)

lat2=subwrd(latline,8)

lat=math_int(lat1/dlat)*dlat-dlat

'set line 0';* set color for recf

'set string 1 c'

'set strsiz '%strsiz

if (font!=""); ' set font '%font; endif

while(lat<lat2)

if (lat1<=lat)

'query w2xy '%lon0%' '%lat

x=subwrd(result,3)

y=subwrd(result,6)

if (lat<0); clat=-lat%'S'; endif

if (lat>0); clat=lat%'N'; endif

if (lat=0); clat='Eq'; endif

len=math_strlen(clat)

xwl=strsiz*len*1.1/2

xwr=strsiz*len*1.3/2

ywb=strsiz*1.4/2

ywt=strsiz*1.2/2

* 'draw recf '%x-xwl%' '%y-ywb%' '%x+xwr%' '%y+ywt

'draw string '%x%' '%y%' '%clat

endif

lat=lat+dlat

endwhile

'set string 1 bl 3 0'

return

----------------------

End of circlat.gs

----------------------

======================

circlon.gs

======================

function circlon(args)

* usage: run circlon dlon (strsiz (font))

* mproj=nps もしくは sps において longitudeの数字を描く

* dlon は数字を描く経度の間隔で必須．

* strsiz は文字の大きさでオプション

* font はフォント番号でオプション

dlon=subwrd(args,1)

strsiz=subwrd(args,2)

font=subwrd(args,3)

'q gxinfo'

cline=sublin(result,6)

mproj=subwrd(cline,3)

if (mproj!=3 & mproj!=4)

say ' *circlon* (mproj!=3:nps & mproj !=4:sps) is not supported'

say 'mproj='%mproj

'quit'

endif

if (mproj=3);* mproj=3 is nps

'query w2xy 0 90'

* say 'result='result

endif

if (mproj=4);* mproj=4 is sps

'query w2xy 0 -90'

endif

xpol=subwrd(result,3)

ypol=subwrd(result,6)

'q dims'

lonline=sublin(result,2)

latline=sublin(result,3)

if (mproj=3);* mproj=nps

lat0=subwrd(latline,6)

else;* mproj=sps

lat0=subwrd(latline,8)

endif

lon1=subwrd(lonline,6)

lon2=subwrd(lonline,8)

say 'lon1='lon1' lon2='lon2

'query w2xy '%lon1%' '%lat0

xd=subwrd(result,3)

yd=subwrd(result,6)

say 'xpol='xpol

say 'ypol='ypol

say 'xd='xd

say 'yd='yd

* 'run sqrt '%(xpol-xd)*(xpol-xd)+(ypol-yd)*(ypol-yd); rad=result;* for Grads 1.7 or earlier

rad=math_pow((xpol-xd)*(xpol-xd)+(ypol-yd)*(ypol-yd),0.5);

* say 'rad='rad

* say 'lon1='lon1' lon2='lon2

radd=rad+0.2

if (strsiz!="")

radd=rad+1.5*strsiz

'set strsiz '%strsiz

endif

if (font!=""); ' set font '%font; endif

lonm=(lon1+lon2)/2

*lon=0; while (lon<=360)

lon=lon1; while (lon<=lon2)

if (lon1<=lon & lon<lon2)

* say 'lon='lon

anglestr=lon-lonm

angleline=(anglestr-90)/180*3.14

* 'set string 1 c 3 '%anglestr

'set string 1 c'

cosangle=math_cos(angleline);

sinangle=math_sin(angleline);

* 'run cos '%angleline; cosangle=result;* for Grads 1.7 or earlier

* 'run sin '%angleline; sinangle=result;* for Grads 1.7 or earlier

xrim=xpol+rad*cosangle

yrim=ypol+rad*sinangle

* 'draw line '%xpol%' '%ypol%' '%xrim%' '%yrim

xchar=xpol+radd*cosangle

ychar=ypol+radd*sinangle

ew=''

if (-360<lon & lon<-180) ; ew='W'; clon=(360+lon)%ew; endif

if (-180<lon & lon<0) ; ew='W'; clon=-lon%ew ; endif

if ( 0<lon & lon<180) ; ew='E'; clon=lon%ew ; endif

if ( 180<lon & lon<360) ; ew='W'; clon=(360-lon)%ew; endif

if ( 360<lon & lon<540) ; ew='E'; clon=(lon-360)%ew; endif

if ( 540<lon & lon<720) ; ew='W'; clon=(720-lon)%ew; endif

if ( -360=lon | 0=lon | 360=lon | 720=lon ) ; ew='' ; clon=0%ew; endif

if ( -180=lon | 180=lon | 540=lon ) ; ew=''; clon=180%ew; endif

'draw string '%xchar%' '%ychar%' '%clon

endif

lon=lon+dlon

endwhile

'set string 1 bl 3 0'

return

lat=lat1; while (lat<=lat2)

endwhile

----------------------

End of circlon.gs

----------------------

======================

tile.gs

======================

*

* Help is in the end of this script

*

function tile( args )

_version = '0.01b1'

rc = gsfallow( 'on' )

if( args = '' )

help()

return

endif

***** Default value *****

var = 'none'

vmin = 'none'

vmax = 'none'

int = 9

color = 1

type = 5

thickness = 3

varmin = 'none'

varmax = 'none'

***** Arguement *****

i = 1

while( 1 )

arg = subwrd( args, i )

i = i + 1;

if( arg = '' ); break; endif

while( 1 )

*** option

if( arg = '-min' ); vmin = subwrd(args,i); i=i+1; break; endif

if( arg = '-max' ); vmax = subwrd(args,i); i=i+1; break; endif

if( arg = '-int' ); int = subwrd(args,i); i=i+1; break; endif

if( arg = '-color' ); color = subwrd(args,i); i=i+1; break; endif

if( arg = '-type' ); type = subwrd(args,i); i=i+1; break; endif

if( arg = '-thickness' ); thickness = subwrd(args,i); i=i+1; break; endif

*** var, min, max

if( var != 'none' & vmin != 'none' & vmax = 'none' & valnum(arg) != 0 )

vmax = arg

break

endif

if( var != 'none' & vmin = 'none' & valnum(arg) != 0 )

vmin = arg

break

endif

if( var = 'none' & valnum(arg) = 0 )

var = arg

varmin = arg

break

endif

if( varmax = 'none' & valnum(arg) = 0 )

varmax = arg

break

endif

say 'syntax error : 'arg

say 'type "tile" for help'

return

endwhile

endwhile

if( vmin = 'none' ); vmin = -1e+30; endif

if( vmax = 'none' ); vmax = 1e+30; endif

'set tile 1 'type' 'int' 'int' 'thickness' 'color

'set rgb 16 tile 1'

* for 1-D

if( qdims( 'varying' ) = 1 )

'set gxout linefill'

'set lfcols -1 16'

'd 'varmin';'varmax

* for 2-D

else

'set gxout shaded'

'set clevs 'vmin' 'vmax

'set ccols -1 16 -1'

'd 'var

endif

return

*

* help

*

function help()

say ' Name:'

say ' tile '_version' - Draw tile (hatch).'

say ''

say ' Usage:'

say ' tile ( var ( min max | -min min | -max max ) | varmin varmax ) '

say ' [-type type] [-int int]'

say ' [-color color] [-thickness thickness]'

say ''

say ' var : Variable.'

say ' varmin,varmax : variable range (only for 1D chart)'

say ' min,max : value range to be drawn (default: [-1e+30:1e+30])'

say ' varmin,varmax : range of variable to be drawn (only for 1D chart)'

say ' int : integer tile interval (default: 9)'

say ' color : hatch color (default: 1)'

say ' type : hatch line type (default: 5)'

say ' thickness : hatch line thickness between 1-12 (default: 3)'

say ''

say ' Note:'

say ' [arg-name] : specify if needed'

say ' (arg1 | arg2) : arg1 or arg2 must be specified'

say ''

say ' Only for GrADS 2.1 or later.'

say ''

say ' Copyright (C) 2015-2015 Chihiro Kodama'

say ' Distributed under GNU GPL (http://www.gnu.org/licenses/gpl.html)'

say ''

return

----------------------

End of tile.gs

----------------------

======================

qdims.gsf

======================

* id = fnum

* varying/fixed : number of varying/fixed dimensions

*

* xtype

* lonmin, lonmax

* xmin, xmax

* ytype

* ...

* z

* t

* e

function qdims( id )

ret = ''

'q dims'

if( id = 'fnum' )

line = sublin( result, 1 )

ret = subwrd( line, 5 )

return ret

endif

if( id = 'varying' | id = 'fixed' )

ret = 0

i = 2

while( i <= 6 )

line = sublin( result, i )

stat = subwrd( line, 3 )

if( stat = id ) ; ret = ret + 1 ; endif

i = i + 1

endwhile

return ret

endif

if( id = 'xtype' | id = 'lonmin' | id = 'lonmax' | id = 'xmin' | id = 'xmax' )

line = sublin( result, 2 )

endif

if( id = 'ytype' | id = 'latmin' | id = 'latmax' | id = 'ymin' | id = 'ymax' )

line = sublin( result, 3 )

endif

if( id = 'ztype' | id = 'levmin' | id = 'levmax' | id = 'zmin' | id = 'zmax' )

line = sublin( result, 4 )

endif

if( id = 'ttype' | id = 'timemin' | id = 'timemax' | id = 'tmin' | id = 'tmax' )

line = sublin( result, 5 )

endif

if( id = 'etype' | id = 'ensmin' | id = 'ensmax' | id = 'emin' | id = 'emax' )

line = sublin( result, 6 )

endif

type = subwrd( line, 3 )

if( id = 'xtype' | id = 'ytype' | id = 'ztype' | id = 'ttype' | id = 'etype' )

ret = type

endif

if( id = 'lonmin' | id = 'latmin' | id = 'levmin' | id = 'timemin' | id = 'ensmin' )

ret = subwrd( line, 6 )

return ret

endif

if( id = 'lonmax' | id = 'latmax' | id = 'levmax' | id = 'timemax' | id = 'ensmax' )

if( type = 'fixed' ) ; ret = subwrd( line, 6 )

else ; ret = subwrd( line, 8 ) ; endif

return ret

endif

if( id = 'xmin' | id = 'ymin' | id = 'zmin' | id = 'tmin' | id = 'emin' )

if( type = 'fixed' ) ; ret = subwrd( line, 9 )

else ; ret = subwrd( line, 11 ) ; endif

return ret

endif

if( id = 'xmax' | id = 'ymax' | id = 'zmax' | id = 'tmax' | id = 'emax' )

if( type = 'fixed' ) ; ret = subwrd( line, 9 )

else ; ret = subwrd( line, 13 ) ; endif

return ret

endif

return ret

----------------------

End of qdims.gsf

----------------------

======================

plev.ctl

======================

dset /work4/data/ERA-Interim/ERA-I_Barents_110116-23_plev.nc

title

undef -9.99e+33

dtype netcdf

UNPACK scale_factor add_offset

OPTIONS yrev zrev

xdef 480 linear -180 0.75

ydef 81 linear 30 0.75

zdef 17 levels 1000 975 950 925 900 850 800 700

600 500 400 300 250 200 150 100 50

tdef 32 linear 00Z16JAN2011 360mn

vars 11

pv=>pv 17 t,z,y,x Potential vorticity

z=>z 17 t,z,y,x Geopotential

t=>t 17 t,z,y,x Temperature

q=>q 17 t,z,y,x Specific humidity

w=>w 17 t,z,y,x Vertical velocity

r=>r 17 t,z,y,x Relative humidity

clwc=>clwc 17 t,z,y,x Specific cloud liquid water content

ciwc=>ciwc 17 t,z,y,x Specific cloud ice water content

cc=>cc 17 t,z,y,x Fraction of cloud cover

u=>u 17 t,z,y,x U component of wind

v=>v 17 t,z,y,x V component of wind

endvars

----------------------

End of plev.ctl

----------------------

======================

sfc.ctl

======================

dset /work4/data/ERA-Interim/ERA-I_Barents_110116-23_sfc.nc

title

undef -9.99e+33

dtype netcdf

UNPACK scale_factor add_offset

OPTIONS yrev

xdef 480 linear -180 0.75

ydef 81 linear 30 0.75

zdef 1 linear 0 1

tdef 32 linear 00Z16JAN2011 360mn

vars 9

ci=>ci 0 t,y,x Sea-ice cover

sst=>sst 0 t,y,x Sea surface temperature

sp=>sp 0 t,y,x Surface pressure

msl=>msl 0 t,y,x Mean sea level pressure

u10=>u10 0 t,y,x 10 metre U wind component

v10=>v10 0 t,y,x 10 metre V wind component

t2m=>t2m 0 t,y,x 2 metre temperature

d2m=>d2m 0 t,y,x 2 metre dewpoint temperature

al=>al 0 t,y,x Albedo

endvars

----------------------

End of sfc.ctl

----------------------