Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

sst_update and ocean coupling

Dear All,
my atmnospheric and oceanic domains are the same as like the computational grids. Hence I do not need to merge SST data from other sources, so that CPLMASK = 1 everywhere.
I want to try the sst_skin module, but to do that I must activate the sst_update.
Real was run with sst_update=1 .
So before testing sst_skin I want to check the sst_update does not interfere with the SST coming from the coupler.
I did 2 tests:
1) with sst_update=0
2) with sst_update=1
the other part of the namelist is:

Code:
&time_control
start_year                   =  2021
start_month                  =  01
start_day                    =  31
start_hour                   =  00
run_days                     =  1
run_hours                    =  00
run_minutes                  =  00
interval_seconds             =  10800
input_from_file              = .false.
history_outname              = './SEAS_COUPLED_ATM_<date>.nc'
history_interval             = 180
frames_per_outfile           = 8
restart                      =  .false.
rst_outname                  = './restart/atm'
rst_inname                   = './restart/atm'
restart_interval             = -1
nocolons                     = .true.
io_form_history              = 11
input_inname                 = '/work/csp/fm05020/wrf/met_prep/ecmwf_fcd2/2021new_monalb/20210131/wrfinput_d01'
io_form_input                = 2
bdy_inname                   = '/work/csp/fm05020/wrf/met_prep/ecmwf_fcd2/2021new_monalb/20210131/wrfbdy_d01'
io_form_boundary             = 2
auxinput4_inname             = '/work/csp/fm05020/wrf/met_prep/ecmwf_fcd2/2021new_monalb/20210131/wrflowinp_d01'
auxinput4_interval           = 180
io_form_auxinput4            = 2
iofields_filename            = 'non_li_voglio.txt'
/
&domains
time_step                           = 20
time_step_fract_num                 = 0
time_step_fract_den                 = 1
max_dom                             = 1
e_we                                = 1658
e_sn                                = 752
e_vert                              = 40
p_top_requested                     = 5000
num_metgrid_levels                  = 26
num_metgrid_soil_levels             = 4
dx                                  = 4632.395
dy                                  = 4632.395
grid_id                             = 1
parent_id                           = 0
i_parent_start                      = 1
j_parent_start                      = 1
parent_grid_ratio                   = 1
parent_time_step_ratio              = 1
feedback                            = 1
smooth_option                       = 0
num_ext_model_couple_dom            = 1
nproc_x                             = 16
nproc_y                             = 64
/
&physics
mp_physics                          = 8
cu_physics                          = 6
cudt                                = 15
ra_lw_physics                       = 4
ra_sw_physics                       = 4
radt                                = 15
bl_pbl_physics                      = 2
sf_sfclay_physics                   = 2
bldt                                = 15
sf_surface_physics                  = 2
sf_lake_physics                     = 0
lakedepth_default                   = 15
use_lakedepth                       = 0
sst_update                          = 0
prec_acc_dt                         = 180
icloud                              = 1
num_land_cat                        = 21
sf_urban_physics                    = 0
usemonalb                           = .true.
rdmaxalb                            = .true.
rdlai2d                             = .true.
/
&dynamics
hybrid_opt                          = 2
w_damping                           = 0
diff_opt                            = 1
km_opt                              = 4
diff_6th_opt                        = 0
diff_6th_factor                     = 0.12
base_temp                           = 290.
damp_opt                            = 3
zdamp                               = 5000.
dampcoef                            = 0.2
khdif                               = 0
kvdif                               = 0
non_hydrostatic                     = .true.
moist_adv_opt                       = 1
scalar_adv_opt                      = 1
gwd_opt                             = 1
/
&fdda
/
&bdy_control
spec_bdy_width                      = 5
specified                           = .true.
/
&grib2
/

I would have expected the results being identical because the CPLMASK=1 everywhere should basically get rid of the SST_INPUT.
Unfortunately the results are different, no matter of which surface variable I select.
What am I missing?
Thank you for your help
 
Hi,
Apologies for the delay in response. I don't have experience with running coupled simulations, but I can try my best to help (and/or ask a colleague who may be able to assist, if needed). First, is this still an issue for you? If so, can you let me know what version of WRF you are using? Thanks!
 
Top