Hi everyone,
I want to restore model SST to observation by setting pt_interior_nml in a fully coupled run BHIST.
This is my pt_interior_nml:
&forcing_pt_interior_nml
pt_interior_data_inc = 24.
pt_interior_data_renorm = 1.
pt_interior_data_type = 'n-hour'
pt_interior_file_fmt = 'nc'
pt_interior_filename = '/public1/home/sc93477/xyw/CODE/forcing/pt_interior'
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 = 'bin'
pt_interior_restore_filename = 'unknown-pt_interior_restore'
pt_interior_restore_max_level = 0
pt_interior_restore_tau = 1
pt_interior_shr_stream_file = 'unknown-pt_interior_shr_stream'
pt_interior_shr_stream_year_align = 1
pt_interior_shr_stream_year_first = 1
pt_interior_shr_stream_year_last = 1
pt_interior_surface_restore = .true.
pt_interior_variable_restore = .false.
I pre-processed daily sst forcing data.The daily SST forcing is as follow:
filename : pt_interior.2010.001.12
path : ./pt_interior.2010.001.12
file global attributes:
dimensions:
time = 0 // unlimited
nlat = 384
nlon = 320
variables:
float TEMPERATURE ( nlat, nlon )
cell_methods : time: mean
grid_loc : 2110
coordinates : TLONG TLAT time
units : degC
long_name : Surface Potential Temperature
time : 734292
_FillValue : 9.96921e+36
valid_range : ( -3, 45 )
precision : 2
dataset : NOAA High-resolution Blended Analysis
var_desc : Sea Surface Temperature
level_desc : Surface
statistic : Mean
parent_stat : Individual Observations
actual_range : ( -1.8, 37.64 )
missing_value : 9.96921e+36
spatial_op : Bilinear remapping: 1st order: destarea: NCL: ./map_0.25x0.25_to_gx1v7_bilin_da_20220912.nc
integer nlat ( nlat )
long_name : Latitude index
integer nlon ( nlon )
long_name : Longitude index
I start running the case,but SST have no change in restoring run compared to the control run.
Is there any problem with my setting like this? Can you give me some advice?
Any help is appreciate!
Wang
I want to restore model SST to observation by setting pt_interior_nml in a fully coupled run BHIST.
This is my pt_interior_nml:
&forcing_pt_interior_nml
pt_interior_data_inc = 24.
pt_interior_data_renorm = 1.
pt_interior_data_type = 'n-hour'
pt_interior_file_fmt = 'nc'
pt_interior_filename = '/public1/home/sc93477/xyw/CODE/forcing/pt_interior'
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 = 'bin'
pt_interior_restore_filename = 'unknown-pt_interior_restore'
pt_interior_restore_max_level = 0
pt_interior_restore_tau = 1
pt_interior_shr_stream_file = 'unknown-pt_interior_shr_stream'
pt_interior_shr_stream_year_align = 1
pt_interior_shr_stream_year_first = 1
pt_interior_shr_stream_year_last = 1
pt_interior_surface_restore = .true.
pt_interior_variable_restore = .false.
I pre-processed daily sst forcing data.The daily SST forcing is as follow:
filename : pt_interior.2010.001.12
path : ./pt_interior.2010.001.12
file global attributes:
dimensions:
time = 0 // unlimited
nlat = 384
nlon = 320
variables:
float TEMPERATURE ( nlat, nlon )
cell_methods : time: mean
grid_loc : 2110
coordinates : TLONG TLAT time
units : degC
long_name : Surface Potential Temperature
time : 734292
_FillValue : 9.96921e+36
valid_range : ( -3, 45 )
precision : 2
dataset : NOAA High-resolution Blended Analysis
var_desc : Sea Surface Temperature
level_desc : Surface
statistic : Mean
parent_stat : Individual Observations
actual_range : ( -1.8, 37.64 )
missing_value : 9.96921e+36
spatial_op : Bilinear remapping: 1st order: destarea: NCL: ./map_0.25x0.25_to_gx1v7_bilin_da_20220912.nc
integer nlat ( nlat )
long_name : Latitude index
integer nlon ( nlon )
long_name : Longitude index
I start running the case,but SST have no change in restoring run compared to the control run.
Is there any problem with my setting like this? Can you give me some advice?
Any help is appreciate!
Wang
