[ncl-talk] YYYYMMDDHH conversion

Paolina Bongioannini Cerlini paolina.cerlini at unipg.it
Thu May 11 12:06:23 MDT 2017


Hi NCL people,
I read some data and since I have days and hours I build a
 time as YYYYMMDDHH then I  convert it in yrfrac:

; ***********************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
;************************************************
begin
;************************************************
; read in data
;************************************************
;
   diri = "./"
   fili = systemfunc("cd "+diri+" ; ls rce-bret-2005-rad.*.nc")
   nfili= dimsizes(fili)
   print("nfili="+nfili)
;
   f    = addfiles(diri+fili,"r") 
   ListSetType (f, "join") 
;
;    vNam = getfilevarnames (f)  ; all variables on file   
;    print(vNam)
;
    vNam1 = (/"x","y","z","p","raing","prcate1","tke","radfrc","radflx","ptsflx","qvsflx"/)  ; manually specify
    print(vNam1)
;
;	  
     x   = f[0]->x
     y   = f[0]->y
     z   = f[0]->z 
     p  = f[:]->p
     PREC= f[:]->raing
     RADQR= f[:]->rnflx
     LHF = f[:]->ptsflx
     SHF = f[:]->qvsflx
;
    printVarSummary(x)
    printVarSummary(y)
    printVarSummary(z)           ; [74]
    printVarSummary(p)      
;
; create time 
      TIME3D = yyyymmddhh_time(2017, 2017, 1, "integer")
      printVarSummary(TIME3D)
;      time3D = TIME3D({2017040100:2017051000})    ; coordinate subscripting
      time3D = TIME3D({2017040100:2017040203})
      time3D!0="time"
      time3D&time=time3D
      printVarSummary(time3D)  
;
;
;      print(time3D) 
;
      yrfrac = yyyymmddhh_to_yyyyfrac(time3D, 0)
      print(yrfrac)
; 
 but it gives the following 
error:

[paolina at cloud outputcheck]$ ncl energy_ARPS-PROVA_0.ncl>output.txt
[paolina at cloud outputcheck]$ more output.txt
 Copyright (C) 1995-2014 - All Rights Reserved
 University Corporation for Atmospheric Research
 NCAR Command Language Version 6.2.0
 The use of this software is governed by a License Agreement.
 See http://www.ncl.ucar.edu/ for more details.
(0)	nfili=28


Variable: vNam1
Type: string
Total Size: 88 bytes
            11 values
Number of Dimensions: 1
Dimensions and sizes:	[11]
Coordinates: 
(0)	x
(1)	y
(2)	z
(3)	p
(4)	raing
(5)	prcate1
(6)	tke
(7)	radfrc
(8)	radflx
(9)	ptsflx
(10)	qvsflx

Variable: x
Type: float
Total Size: 780 bytes
            195 values
Number of Dimensions: 1
Dimensions and sizes:	[fakeDim0 | 195]
Coordinates: 
Number Of Attributes: 3
  hdf_name :	x
  units :	m
  comment :	x coordinate

Variable: y
Type: float
Total Size: 780 bytes
            195 values
Number of Dimensions: 1
Dimensions and sizes:	[fakeDim1 | 195]
Coordinates: 
Number Of Attributes: 3
  hdf_name :	y
  units :	m
  comment :	y coordinate

Variable: z
Type: float
Total Size: 296 bytes
            74 values
Number of Dimensions: 1
Dimensions and sizes:	[fakeDim2 | 74]
Coordinates: 
Number Of Attributes: 3
  hdf_name :	z
  units :	m
  comment :	z coordinate

Variable: p
Type: float
Total Size: 315151200 bytes
            78787800 values
Number of Dimensions: 4
Dimensions and sizes:	[ncl_join | 28] x [fakeDim39 | 74] x [fakeDim40 | 195] x [fakeDim41 | 195]
Coordinates: 
Number Of Attributes: 6
  hdf_name :	p
  stag_dim :	0
  units :	Pascal
  comment :	Pressure
  hdf_comp_prm :	6
  hdf_comp_code :	4

Variable: TIME3D
Type: integer
Total Size: 35040 bytes
            8760 values
Number of Dimensions: 1
Dimensions and sizes:	[time | 8760]
Coordinates: 
            time: [2017010100..2017123123]
Number Of Attributes: 2
  long_name :	time
  units :	YYYYMMDDHH

Variable: time3D
Type: integer
Total Size: 112 bytes
            28 values
Number of Dimensions: 1
Dimensions and sizes:	[time | 28]
Coordinates: 
            time: [2017040100..2017040203]
Number Of Attributes: 2
  units :	YYYYMMDDHH
  long_name :	time


Variable: time3D
Type: integer
Total Size: 112 bytes
            28 values
Number of Dimensions: 1
Dimensions and sizes:	[time | 28]
Coordinates: 
            time: [2017040100..2017040203]
Number Of Attributes: 2
  units :	YYYYMMDDHH
  long_name :	time
(0)	2017040100
(1)	2017040101
(2)	2017040102
(3)	2017040103
(4)	2017040104
(5)	2017040105
(6)	2017040106
(7)	2017040107
(8)	2017040108
(9)	2017040109
(10)	2017040110
(11)	2017040111
(12)	2017040112
(13)	2017040113
(14)	2017040114
(15)	2017040115
(16)	2017040116
(17)	2017040117
(18)	2017040118
(19)	2017040119
(20)	2017040120
(21)	2017040121
(22)	2017040122
(23)	2017040123
(24)	2017040200
(25)	2017040201
(26)	2017040202
(27)	2017040203
fatal:Variable (nDay) is undefined
fatal:["Execute.c":8567]:Execute: Error occurred at or near line 7879 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl

fatal:["Execute.c":8567]:Execute: Error occurred at or near line 7937 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl

fatal:["Execute.c":8567]:Execute: Error occurred at or near line 63 in file energy_ARPS-PROVA_0.ncl
-------------------------------------------------------------------------------------------------------------------------------------------
In attachment the files.
Thank you for your help.

Paolina



-- 
Paolina Bongioannini Cerlini Ph.D. 
University of Perugia 
CIRIAF -room 8 
06125 - Perugia 

Tel.+39 075 5853576 
e-mail: paolina.cerlini at unipg.it 
skype: pcerlini 
--
-------------- next part --------------
A non-text attachment was scrubbed...
Name: energy_ARPS-PROVA_0.ncl
Type: application/octet-stream
Size: 6475 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170511/5fd8f29c/attachment.obj 
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: output.txt
Url: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170511/5fd8f29c/attachment.txt 


More information about the ncl-talk mailing list