Hi Jin,Your pt_interior_restore_filename file should have two fields (both 2D):PT_RESTORE_MAX_LEVEL* This should be an integer corresponding to the deepest level you want to mix at (for each point). For a gx3v7 grid, if you want to mix everywhere, you can set it to 60 at every point. PT_RESTORE_RTAU* This is the inverse time scale that you want to restore over (in units of day^-1). For example, if you want to restore on a timescale of 5 days at a given point, you would set this value to 0.2 at the specific (i,j) coordinate. Good luck!~Mike
Hi, Mike, I am very sory to bother you! But I have a same problem with Jin and don't know how to fix it. I really need your help. My namelist is as follows:
&forcing_pt_interior_nml
pt_interior_data_inc = 1.
pt_interior_data_renorm = 1.
pt_interior_data_type = 'monthly'
pt_interior_file_fmt = 'nc'
pt_interior_filename = '/ncsfs/home/trifle/prescribed_ice_sst/ncl_scripts/model_input_test/sst/model_input/TEMPERATURE.gx1.2018.hh_60.nc'
pt_interior_formulation = 'restoring'
pt_interior_interp_freq = 'every-timestep'
pt_interior_interp_inc = 72.
pt_interior_interp_type = 'linear'
pt_interior_restore_file_fmt = 'nc'
pt_interior_restore_filename = '/ncsfs/home/trifle/prescribed_ice_sst/ncl_scripts/model_input_test/sst/model_input/PT.gx1.hh_60.nc'
pt_interior_restore_max_level = 1
pt_interior_restore_tau = 0.1
pt_interior_variable_restore = .false.
The TEMPERATURE.gx1.2018.hh_60.nc includs 12 fields(60x384x320): TEMPERATURE01 TEMPERATURE02 ...TEMPERATURE12. The PT.gx1.hh_60.nc includs 2 fields(384x320): PT_RESTORE_MAX_LEVEL and PT_RESTORE_RTAU.
But at begining of the run, I got the error:
POP aborting...
Error in getting varid for netCDF field
------------------------------------------------------------------------
application called MPI_Abort(MPI_COMM_WORLD, 0) - process 1
------------------------------------------------------------------------
Is there any wrong with my nc files? The following is some details of these nc files:
filename: TEMPERATURE.gx1.2018.hh_60
path: TEMPERATURE.gx1.2018.hh_60.nc
file global attributes:
CDI : Climate Data Interface version 1.6.8 (http://mpimet.mpg.de/cdi)
Conventions : CF-1.4
history : Sun May 10 13:38:31 2020: cdo -remapbil,grid.txt ./org/TEMPERATURE_1x1_plus_amoice_hh_2018_60_12.nc ./TEMPERATURE.gx1.2018.hh_60.nc
CDO : Climate Data Operators version 1.6.8 (http://mpimet.mpg.de/cdo)
dimensions:
x = 320
y = 384
z_t = 60
time = 12 // unlimited
variables:
double TLON ( y, x )
standard_name : longitude
long_name : T grid center longitude
units : degrees_east
_CoordinateAxisType : Lon
double TLAT ( y, x )
standard_name : latitude
long_name : T grid center latitude
units : degrees_north
_CoordinateAxisType : Lat
float z_t ( z_t )
long_name : depth from surface to midpoint of layer
units : centimeters
positive : down
axis : Z
double time ( time )
standard_name : time
long_name : model time
units : days since 0-01-01 00:00:00
calendar : 365_day
float TEMPERATURE01 ( z_t, y, x )
long_name : Potential Temperature
units : degC
coordinates : TLON TLAT
_FillValue : 9.96921e+36
missing_value : 9.96921e+36
cell_methods : time: mean
grid_loc : 3111
float TEMPERATURE02 ( z_t, y, x )
long_name : Potential Temperature
units : degC
coordinates : TLON TLAT
_FillValue : 9.96921e+36
missing_value : 9.96921e+36
grid_loc : 3111
cell_methods : time: mean
.............
float TEMPERATURE12 ( z_t, y, x )
long_name : Potential Temperature
units : degC
coordinates : TLON TLAT
_FillValue : 9.96921e+36
missing_value : 9.96921e+36
grid_loc : 3111
cell_methods : time: mean
integer datesec ( time )
long_name : current seconds of current date
integer date ( time )
long_name : current date (YYYYMMDD)
filename: PT.gx1.hh_60
path: PT.gx1.hh_60.nc
file global attributes:
dimensions:
y = 384
x = 320
variables:
integer PT_RESTORE_MAX_LEVEL ( y, x )
lon2d : <ARRAY of 122880 elements>
lat2d : <ARRAY of 122880 elements>
missing_value : -2147483647
missing_value_original : 9.96921e+36
_FillValue : -2147483647
_FillValue_original : 9.96921e+36
grid_loc : 3111
cell_methods : time: mean
coordinates : TLON TLAT
units : degC
long_name : Potential Temperature
z_t : 500
float PT_RESTORE_RTAU ( y, x )
lon2d : <ARRAY of 122880 elements>
lat2d : <ARRAY of 122880 elements>
z_t : 500
long_name : Potential Temperature
units : degC
coordinates : TLON TLAT
_FillValue : 9.96921e+36
missing_value : 9.96921e+36
cell_methods : time: mean
grid_loc : 3111
double TLAT ( y, x )
standard_name : latitude
long_name : T grid center latitude
units : degrees_north
_CoordinateAxisType : Lat
double TLON ( y, x )
standard_name : longitude
long_name : T grid center longitude
units : degrees_east
_CoordinateAxisType : Lon
During the process of generating nc files, I mainly pay attention to the variable name and resolution. Is tere anything else I haven't noticed? Could you give me some advice? Thank you very much!
Best wishes!
Zhuzhu