Scheduled Downtime
On Tuesday 24 October 2023 @ 5pm MT the forums will be in read only mode in preparation for the downtime. On Wednesday 25 October 2023 @ 5am MT, this website will be down for maintenance and expected to return online later in the morning.
Normal Operations
The forums are back online with normal operations. If you notice any issues or errors related to the forums, please reach out to help@ucar.edu

how to write a restore nc file in forcing_pt_interior_nml

J

jinyishuai@126_com

New Member
Dear professors:
I am so sorry to ask for your help.i have asked many questions....But it is really important to me.
This is my forcing namelist in pop
forcing_pt_interior_nml
pt_interior_data_type = 'monthly-calendar'
pt_interior_data_inc = 24.
pt_interior_interp_freq = 'every-timestep'
pt_interior_interp_type = 'linear'
pt_interior_interp_inc = 72.
pt_interior_restore_tau = 365.
pt_interior_filename = '/home/jys/work/cesm/archive/PB30/ocn/climate_nc_dou/temp_1501.nc'
pt_interior_file_fmt = 'nc'
pt_interior_restore_max_level = 60
pt_interior_formulation = 'restoring'
pt_interior_data_renorm(1) = 1.
pt_interior_variable_restore = .true.
pt_interior_restore_filename = '/home/jys/work/cesm/archive/PB30/ocn/climate_nc_dou/pt_restore.nc'
pt_interior_restore_file_fmt = 'nc'
/
My pt_restore file contains two variables:rs_depth and rs_time,they are all nlat*nlon,nlat=116,nlon=100,but in ccsm.log*,
POP aborting...
Error in getting varid for netCDF field
Please tell me how to write this nc file.
Many thanks for your helping!
Jin
 
M

mlevy

Michael Levy
CSEG and Liaisons
Staff member
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  
 
Z

zhuouc@163_com

New Member
Hi Mike,Sorry to bother you, but I am confused when managing pt_interior_restore. Although the pt_interior_restore_filename file I specify contains two fields as you suggested, the case failed at the very begining. The error information in ocn.log file readsPOP Exiting...
POP_SolversChronGear: solver not converged
POP_SolverRun: error in ChronGear
POP_BarotropicDriver: error in solver
Step: error in barotropicLook forward for any advice. Thank you!Chenyu
 
H

happy

happy
New Member
mlevy said:
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
Click to expand...
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
 
You must log in or register to reply here.
Top