running GFDL vortex tracker with netcdf file

Submitted by hjshin on Wed, 05/26/2021 - 23:52

Hello :)

I'm trying to  run GFDL vortex tracker to get extratropical cyclones information using "tcgen" or "midlat" option with NetCDF file.

but I had found that there are no information about that in google

 

first, I tried "Tracker" mode with netcdf file.

I used GFDL vortex tracker v.3.9a.

and I succeeded run for tracking the TC (tracker mode) with the  test_data set downloaded from DTC.

 

but  when doing "tracker mode" with netcdf file,

It seems not to recognize netcdf file...

may be should've be manipulate netcdf file to have variables separated with level?   (example: v => v700hPa)

or some errors in the namelist file?

 

here are the namelist and  the log file I used...

please look through what the problem is.

 

sincerely,

hyun jin.

----<name list>-----------------------------------------------------

&datein
  inp%bcc = 20,
  inp%byy = 20,
  inp%bmm = 09,
  inp%bdd = 03,
  inp%bhh = 00,
  inp%model = 7,
  inp%modtyp = "global",
  inp%lt_units = "hours",
  inp%file_seq = "multi",
  inp%nesttyp = "fixed",
/

&atcfinfo
  atcfnum = 81,
  atcfname = "GKIM",
  atcfymdh = 2020090300,
  atcffreq = 300,
/

&trackerinfo
  trkrinfo%type = "tracker",
  trkrinfo%mslpthresh = 0.0015,
  trkrinfo%use_backup_mslp_grad_check = "y",
  trkrinfo%v850thresh = 1.5,
  trkrinfo%use_backup_850_vt_check = "y",
  trkrinfo%enable_timing = 1,
  trkrinfo%gridtype = "global",
  trkrinfo%contint = 100.0,
  trkrinfo%want_oci = T,
  trkrinfo%out_vit = "n",
  trkrinfo%use_land_mask = "n",
  trkrinfo%inp_data_type = "netcdf",
/

&phaseinfo
  phaseflag = "n",
  phasescheme = "n",
  wcore_depth = 1.0,
/

&structinfo
  structflag = "n",
  ikeflag = "n",
/

&fnameinfo
  gmodname = "GKIM",
  rundescr = "tracker",
  atcfdescr = "",
/

&waitinfo
  use_waitfor = "n",
  wait_min_age = 10,
  wait_min_size = 100,
  wait_max_wait = 3600,
  wait_sleeptime = 5,
  use_per_fcst_command = "y",
  per_fcst_command = "./deliver %[FHOUR] %[FMIN]",
/

&netcdflist
  netcdfinfo%lat_name = "lats",
  netcdfinfo%lmaskname = "slmask",
  netcdfinfo%lon_name = "lons",
  netcdfinfo%mslpname = "psl",
  netcdfinfo%netcdf_filename = "GKIM.tracker.2020090300.f %[FMIN]",
  netcdfinfo%num_netcdf_vars = 999,
  netcdfinfo%rv700name = "vorpot",
  netcdfinfo%rv850name = "vorpot",
  netcdfinfo%time_name = "time",
  netcdfinfo%time_units = "hours",
  netcdfinfo%tmean_300_500_name = "X",
  netcdfinfo%u500name = "u",
  netcdfinfo%u700name = "u",
  netcdfinfo%u850name = "u",
  netcdfinfo%usfcname = "u",
  netcdfinfo%v500name = "v",
  netcdfinfo%v700name = "v",
  netcdfinfo%v850name = "v",
  netcdfinfo%vsfcname = "v",
  netcdfinfo%z200name = "hgt",
  netcdfinfo%z300name = "hgt",
  netcdfinfo%z350name = "hgt",
  netcdfinfo%z400name = "hgt",
  netcdfinfo%z450name = "hgt",
  netcdfinfo%z500name = "hgt",
  netcdfinfo%z550name = "hgt",
  netcdfinfo%z600name = "hgt",
  netcdfinfo%z650name = "hgt",
  netcdfinfo%z700name = "hgt",
  netcdfinfo%z750name = "hgt",
  netcdfinfo%z800name = "hgt",
  netcdfinfo%z850name = "hgt",
  netcdfinfo%z900name = "hgt",
/

&parmpreflist
  user_wants_to_track_zeta700 = "y",
  user_wants_to_track_wcirc850 = "y",
  user_wants_to_track_wcirc700 = "y",
  user_wants_to_track_gph850 = "y",
  user_wants_to_track_gph700 = "y",
  user_wants_to_track_mslp = "y",
  user_wants_to_track_wcircsfc = "y",
  user_wants_to_track_zetasfc = "y",
  user_wants_to_track_thick500850 = "y",
  user_wants_to_track_thick200500 = "y",
  user_wants_to_track_thick200850 = "y",
  user_wants_to_track_zeta850 = "y",
/

&verbose
  verb = 3,
  verb_g2 = "X",
/

--------<log message>----------------------------

 TIMING: beginning ...  05:35:25

* . * . * . * . * . * . * . * . * . * . * . * . * . * . * . * . * . * . * . * . 
     PROGRAM GETTRK   HAS BEGUN. COMPILED 1999104.62     ORG: NP22   
     STARTING DATE-TIME  MAY 27,2021  05:35:25.224  147  THU   2459362

 just before trackerinfo read namelist
 just after trackerinfo read namelist
 in read_nlists, verb=            1
  
 After datein namelist in trak.f, namelist parms follow:
 Forecast initial year  = byy =           20
 Forecast initial month = bmm =            9
 Forecast initial day   = bdd =            3
 Forecast initial hour  = bhh =            0
 Forecast model identifier = model=            7
 Forecast model type = modtyp= global  
 Forecast model data lead time units= lt_units= hours  
 Forecast model data sequencing setup= file_seq= multi 
 Forecast model nest type = fixed   
  
 Values read in from atcfinfo namelist: 
ATCF ID = 81  ATCF Name = GKIM
ATCF date (initial date on output atcf records) = 2020090300
ATCF output frequency (in hours*100) = atcffreq =    300
  
 Values read in from trackerinfo namelist follow: 
 western boundary  = westbd  =    0.00
 eastern boundary  = eastbd  =    0.00
 northern boundary = northbd =    0.00
 southern boundary = southbd =    0.00
 tracker type = tracker
 mslp threshold = mslpthresh =  0.0015
 Flag for using backup mslp gradient check= use_backup_mslp_grad_check = y
 v850 threshold = v850thresh =  1.5000
 Flag for using backup 850 mb Vt check= use_backup_850_vt_check = y
 model grid type = global  
 Contour interval to be used =  100.00
 Flag for whether or not roci will be computed and written out for tracker-type case = T
 Flag for whether or not vitals will be written out = n
 Flag for whether or not a land mask will be used for tcgen candidate low filtering = n
 Flag for input data type (grib or netcdf) = netcdf
 Flag for which GRIB version (1 or 2) the input data will be in = 1
 Flag for input GRIB2 JPDTN (0 or 1) =  0
 Flag for input GRIB2 MSLP ID (1 or 192) =    0
 Flag for input GRIB1 MSLP ID (102 or 130) =    0
 Flag for input GRIB1 sfcwind level type (PDS Octet 10... should be 1 or 105) =    0
 Flag for input GRIB1 sfcwind level value  (PDS Octets 11 & 12... usually 0 or 10) =    0
  
  
 Values read in from netcdflist namelist: 
  
Total *possible* number of input NetCDF variables,
     including those that are included in the input
     NetCDF file and those that are not =  999
Input NetCDF filename = GKIM.tracker.2020090300.f %[FMIN]                                                                                                                                                   
 
List of NetCDF variables follows.  A value of X 
indicates the variable is not included in the 
input file and no attempt will be made to read in 
that variable: 
 
NetCDF variable name for 850 mb vort =       vorpot                        
NetCDF variable name for 700 mb vort =       vorpot                        
NetCDF variable name for 850 mb u-comp =     u                             
NetCDF variable name for 850 mb v-comp =     v                             
NetCDF variable name for 700 mb u-comp =     u                             
NetCDF variable name for 700 mb v-comp =     v                             
NetCDF variable name for 850 mb gp height =  hgt                           
NetCDF variable name for 700 mb gp height =  hgt                           
NetCDF variable name for MSLP =              psl                           
NetCDF variable name for near-sfc u-comp =   u                             
NetCDF variable name for near-sfc v-comp =   v                             
NetCDF variable name for 500 mb u-comp =     u                             
NetCDF variable name for 500 mb v-comp =     v                             
NetCDF variable name for 300-500 mb Mean T = X                             
NetCDF variable name for 500 mb gp height =  hgt                           
NetCDF variable name for 200 mb gp height =  hgt                           
NetCDF variable name for land-sea mask =     slmask                        
NetCDF variable name for 900 mb gp height =  hgt                           
NetCDF variable name for 800 mb gp height =  hgt                           
NetCDF variable name for 750 mb gp height =  hgt                           
NetCDF variable name for 650 mb gp height =  hgt                           
NetCDF variable name for 600 mb gp height =  hgt                           
NetCDF variable name for 550 mb gp height =  hgt                           
NetCDF variable name for 450 mb gp height =  hgt                           
NetCDF variable name for 400 mb gp height =  hgt                           
NetCDF variable name for 350 mb gp height =  hgt                           
NetCDF variable name for 300 mb gp height =  hgt                           
NetCDF variable name for time =              time                          
NetCDF variable name for longitudes =        lons                          
NetCDF variable name for latitudes =         lats                          
NetCDF time value (hours|days) =             hours                         
  
  
 Values read in from parmpreflist namelist: 
  
user_wants_to_track_zeta850=      y
user_wants_to_track_zeta700=      y
user_wants_to_track_wcirc850=     y
user_wants_to_track_wcirc700=     y
user_wants_to_track_gph850=       y
user_wants_to_track_gph700=       y
user_wants_to_track_mslp=         y
user_wants_to_track_wcircsfc=     y
user_wants_to_track_zetasfc=      y
user_wants_to_track_thick500850=  y
user_wants_to_track_thick200500=  y
user_wants_to_track_thick200850=  y
  
 Values read in from phaseinfo namelist: 
Storm phase flag = n  Phase scheme = n   
Storm phase, warm core depth (wcore_depth) =    1.00
  
 Values read in from structinfo namelist: 
Structure flag = n
IKE flag = n
  
 Values read in for grib file name from fnameinfo namelist: 
Model name description = gmodname = GKIM
Forecast run description = rundescr = tracker                                 
Optional ATCF / Storm name description = atcfdescr =                                         
  
 Value read in for verbose output for most output:
Value read in for verbose flag = verb =  1
  
 Value read in for verbose output for grib2 output:
Value read in for GRIB2 verbose flag = verb_g2 =  0
  
 Values read in from waitinfo namelist:
Flag for input file waiting = use_waitfor = n
min age (time in seconds since last mod) = wait_min_age =       10
min file size in bytes = wait_min_size =          100
max number of seconds to wait for each file = wait_max_wait =   3600
number of seconds to sleep between checks = wait_sleeptime =      5
command to run after every forecast time = "./deliver %[FHOUR] %[FMIN]"
   1     0
   2   360
   3   720
   4  1080
 in open_grib_files in multi else part....
   lugb=          201   lugi=         5201
 nlen1=           30  nlen2=           33
 gopen_g_file= 
 GKIM.tracker.2020090300.f00000                                                 
                                          
 gopen_i_file= 
 GKIM.tracker.2020090300.f00000.ix                                              
                                          
 tpm gopen_g_file= ...GKIM.tracker.2020090300.f00000...  gopen_i_file= ...GKIM.tracker.2020090300.f00000.ix...
 gopen_g_file= 
 GKIM.tracker.2020090300.f00000                                                 
                                          ....
 gopen_i_file= 
 GKIM.tracker.2020090300.f00000.ix                                              
                                          ....
 TEST open_grib_files, unit lugb=          201  is CLOSED
 TEST open_grib_files, unit lugi=         5201  is CLOSED
 TEST gname  open_grib_files, gfile= 
 GKIM.tracker.2020090300.f00000                                                 
                                           is CLOSED
 TEST iname  open_grib_files, ifile= 
 GKIM.tracker.2020090300.f00000.ix                                              
                                           is CLOSED
  
 !!! ERROR in sub open_grib_files opening grib file
 !!! or grib index file.  baopen return codes:
 !!! grib  file return code = igoret =            0
 !!! index file return code = iioret =          252
 !!! output file return code = iooret =            0

             !!! ERROR: from open_grib_files, rc=  113

 !!! Files after hour0 are missing, exiting normally
./gettrk.exe: symbol lookup error: ./gettrk.exe: undefined symbol: for_stop_core
 

The error messages show that the gettrk.exe tried to open grib files, were the ix files generated with netCDF setting?

There is a mismatch somewhere.

Thanks,

Linlin

Permalink

In reply to by linlin.pan

thank you for replying me. :)

I read that " Input GRIB or NetCDF model forecast and index files (for GRIB files)." at executables pages of user guide.

Isn't the grib file the only thing that needs index file?

I tried to run grbindex.exe with NetCDF file but I got error....

 

Yes, the ix file is only needed for grib file. However, the error messages showed that the ix files were opened.

Are you running the scripts on NOAA HPC (e.g., Jet, Hera) or NCAR machine (e.g., Cheyenne)? I have access to those machines and could check it.

Thanks,

Linlin 

Permalink

In reply to by linlin.pan

thank you for your suggestion. :)

unfortunately,  my machine is neither NOAA HPC nor NCAR machine...

isn't  the data (namelist or log message) enough to fix error?....

To fix error,  If you let me know what more information you need, I can attach it !

or.. do you have an example of namelist file for  NetCDF file?..

thank you :)

What's your email address? My email is Linlin.Pan@noaa.gov. I can email you a sample namelist/sample file.

Thanks,

Linlin