- Log in to post comments
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