WRF-CROCO run question

Hi dear croco users,

When running the WRF-CROCO coupled model, I encountered an issue. Part of the content of out_run.txt is as follows.

I have set the advection scheme to UV_HADV_UP3 , with the ocean model’s time step (dt ) set to 60 seconds. Based on the CFL condition, this should be sufficient to handle relatively high velocities. However, after running for 48 time steps, the simulation blows up. I configured the output to write a restart file (rst ) at each time step, and during this period, the maximum u and v values recorded in the restart files were 12.48 and 8.61, respectively, which are not very large. I am unsure why this issue occurs. What modifications can I make to resolve the problem?

Time splitting: ndtfast = 60 nfast = 82

Maximum grid stiffness ratios: rx0 = 0.325646355278143 rx1 =
5.16849074464534

GET_INITIAL – Processing data for time = 0.1850E+05 record = 1

GET_INITIAL - unable to find variable: hbl
in input NetCDF file: ./croco_ini.nc ==> Initialized to zero state.
>> CAUTION in case of #define EXACT_RESTART <<
If it is the case
- OK if it is a ‘cold start’ i.e coming from a 3rd-party initial file
- otherwise if it is a ‘hot start’ i.e from a restart file produced by this code:
=> problem: run is not restartable
=> check your initial file
CPL-CROCO: sent CROCO_SST 1
CPL-CROCO: sent CROCO_EOCE 2
CPL-CROCO: sent CROCO_NOCE 3
CPL-CROCO: received CROCO_EVPR 4
CPL-CROCO: received CROCO_SRFL 5
CPL-CROCO: received CROCO_STFL 6
CPL-CROCO: received CROCO_TAUM 7
CPL-CROCO: received CROCO_ETAU 8
CPL-CROCO: received CROCO_NTAU 9
CPL-CROCO: received CROCO_PSFC 10
CPL-CROCO: ENDDEF OK
CPL-CROCO: coupling_masks.nc not found
CPL-ROMS: OASIS_TIME 0
GET_BRY – Read all boundary data for time = 0.1850E+05 0
GET_BRY – Read all boundary data for time = 0.1850E+05 0
DEF_HIS/AVG - Created new netCDF file ‘./croco_his.nc’. mynode = 0
WRT_GRID – wrote grid data into file ‘./croco_his.nc’. mynode = 0
WRT_HIS – wrote history fields into time record = 1 / 1 mynode = 0

MAIN: started time-stepping

CPL-CROCO: getvar first 0 0
CPL-CROCO: getvar first 0 0
CPL-CROCO: getvar first 0 0
CPL-CROCO: getvar first 0 0
CPL-CROCO: getvar first 0 0
CPL-CROCO: getvar first 0 0
CPL-CROCO: getvar first 0 0
STEP time[DAYS] KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME trd
0 18497.00000 0.000000000E+00 3.4364351E+01 3.4364351E+01 4.3742971E+15 0
1 18497.00069 1.567972486E-03 3.4369000E+01 3.4370568E+01 4.3745743E+15 0
DEF_RST - Created new netCDF file ‘./croco_rst.nc’. mynode = 0
WRT_GRID – wrote grid data into file ‘./croco_rst.nc’. mynode = 0
WRT_RST – wrote restart fields into time record = 1 / 1 mynode = 0
2 18497.00139 2.797860684E-03 3.4372950E+01 3.4375747E+01 4.3748045E+15 0
WRT_RST – wrote restart fields into time record = 2 / 2 mynode = 0
3 18497.00208 4.012026828E-03 3.4376847E+01 3.4380860E+01 4.3750315E+15 0
WRT_RST – wrote restart fields into time record = 3 / 3 mynode = 0
4 18497.00278 5.216414034E-03 3.4380711E+01 3.4385927E+01 4.3752565E+15 0
WRT_RST – wrote restart fields into time record = 4 / 4 mynode = 0
5 18497.00347 6.413504446E-03 3.4384541E+01 3.4390954E+01 4.3754801E+15 0
WRT_RST – wrote restart fields into time record = 5 / 5 mynode = 0
6 18497.00417 7.594097056E-03 3.4388342E+01 3.4395936E+01 4.3757021E+15 0
WRT_RST – wrote restart fields into time record = 6 / 6 mynode = 0
7 18497.00486 8.754860090E-03 3.4392138E+01 3.4400893E+01 4.3759227E+15 0
WRT_RST – wrote restart fields into time record = 7 / 7 mynode = 0
8 18497.00556 9.890312547E-03 3.4395931E+01 3.4405821E+01 4.3761419E+15 0
WRT_RST – wrote restart fields into time record = 8 / 8 mynode = 0
9 18497.00625 1.100146192E-02 3.4399698E+01 3.4410699E+01 4.3763597E+15 0
WRT_RST – wrote restart fields into time record = 9 / 9 mynode = 0
10 18497.00694 1.209296687E-02 3.4403419E+01 3.4415512E+01 4.3765762E+15 0
WRT_RST – wrote restart fields into time record = 10 / 10 mynode = 0
11 18497.00764 1.316330266E-02 3.4407083E+01 3.4420247E+01 4.3767913E+15 0
WRT_RST – wrote restart fields into time record = 11 / 11 mynode = 0
12 18497.00833 1.421262104E-02 3.4410671E+01 3.4424884E+01 4.3770049E+15 0
WRT_RST – wrote restart fields into time record = 12 / 12 mynode = 0
13 18497.00903 1.523799233E-02 3.4414179E+01 3.4429417E+01 4.3772170E+15 0
WRT_RST – wrote restart fields into time record = 13 / 13 mynode = 0
14 18497.00972 1.624013740E-02 3.4417594E+01 3.4433834E+01 4.3774276E+15 0
WRT_RST – wrote restart fields into time record = 14 / 14 mynode = 0
15 18497.01042 1.721824290E-02 3.4420907E+01 3.4438125E+01 4.3776365E+15 0
WRT_RST – wrote restart fields into time record = 15 / 15 mynode = 0
16 18497.01111 1.817342794E-02 3.4424130E+01 3.4442304E+01 4.3778438E+15 0
WRT_RST – wrote restart fields into time record = 16 / 16 mynode = 0
17 18497.01181 1.910034599E-02 3.4427287E+01 3.4446387E+01 4.3780493E+15 0
WRT_RST – wrote restart fields into time record = 17 / 17 mynode = 0
18 18497.01250 1.999532336E-02 3.4430389E+01 3.4450384E+01 4.3782533E+15 0
WRT_RST – wrote restart fields into time record = 18 / 18 mynode = 0
19 18497.01319 2.085455990E-02 3.4433469E+01 3.4454323E+01 4.3784558E+15 0
WRT_RST – wrote restart fields into time record = 19 / 19 mynode = 0
20 18497.01389 2.168467578E-02 3.4436550E+01 3.4458235E+01 4.3786568E+15 0
WRT_RST – wrote restart fields into time record = 20 / 20 mynode = 0
21 18497.01458 2.248564090E-02 3.4439636E+01 3.4462122E+01 4.3788563E+15 0
WRT_RST – wrote restart fields into time record = 21 / 21 mynode = 0
22 18497.01528 2.325863640E-02 3.4442711E+01 3.4465970E+01 4.3790542E+15 0
WRT_RST – wrote restart fields into time record = 22 / 22 mynode = 0
23 18497.01597 2.400446223E-02 3.4445774E+01 3.4469779E+01 4.3792505E+15 0
WRT_RST – wrote restart fields into time record = 23 / 23 mynode = 0
24 18497.01667 2.473288943E-02 3.4448835E+01 3.4473568E+01 4.3794451E+15 0
WRT_RST – wrote restart fields into time record = 24 / 24 mynode = 0
25 18497.01736 2.544410518E-02 3.4451908E+01 3.4477352E+01 4.3796381E+15 0
WRT_RST – wrote restart fields into time record = 25 / 25 mynode = 0
26 18497.01806 2.613507092E-02 3.4454989E+01 3.4481124E+01 4.3798293E+15 0
WRT_RST – wrote restart fields into time record = 26 / 26 mynode = 0
27 18497.01875 2.680649825E-02 3.4458051E+01 3.4484857E+01 4.3800189E+15 0
WRT_RST – wrote restart fields into time record = 27 / 27 mynode = 0
28 18497.01944 2.746680767E-02 3.4461059E+01 3.4488525E+01 4.3802067E+15 0
WRT_RST – wrote restart fields into time record = 28 / 28 mynode = 0
29 18497.02014 2.811183641E-02 3.4464008E+01 3.4492120E+01 4.3803928E+15 0
WRT_RST – wrote restart fields into time record = 29 / 29 mynode = 0
CPL-CROCO : 0 CROCO_SST sent at time 1740
CPL-CROCO : 0 CROCO_EOCE sent at time 1740
CPL-CROCO : 0 CROCO_NOCE sent at time 1740
CPL-CROCO: getvar first 1800 0
CPL-CROCO: getvar first 1800 0
CPL-CROCO: getvar first 1800 0
CPL-CROCO: getvar first 1800 0
CPL-CROCO: getvar first 1800 0
CPL-CROCO: getvar first 1800 0
CPL-CROCO: getvar first 1800 0
30 18497.02083 2.874440124E-02 3.4466896E+01 3.4495640E+01 4.3805771E+15 0
WRT_RST – wrote restart fields into time record = 30 / 30 mynode = 0
31 18497.02153 2.845697719E-02 3.4464265E+01 3.4492722E+01 4.3804822E+15 0
WRT_RST – wrote restart fields into time record = 31 / 31 mynode = 0
32 18497.02222 2.829027925E-02 3.4462447E+01 3.4490738E+01 4.3804327E+15 0
WRT_RST – wrote restart fields into time record = 32 / 32 mynode = 0
33 18497.02292 2.816638823E-02 3.4460657E+01 3.4488824E+01 4.3803844E+15 0
WRT_RST – wrote restart fields into time record = 33 / 33 mynode = 0
34 18497.02361 2.807117455E-02 3.4458888E+01 3.4486959E+01 4.3803361E+15 0
WRT_RST – wrote restart fields into time record = 34 / 34 mynode = 0
35 18497.02431 2.799010617E-02 3.4457152E+01 3.4485142E+01 4.3802873E+15 0
WRT_RST – wrote restart fields into time record = 35 / 35 mynode = 0
36 18497.02500 2.792813886E-02 3.4455430E+01 3.4483358E+01 4.3802378E+15 0
WRT_RST – wrote restart fields into time record = 36 / 36 mynode = 0
37 18497.02569 2.787962422E-02 3.4453698E+01 3.4481577E+01 4.3801876E+15 0
WRT_RST – wrote restart fields into time record = 37 / 37 mynode = 0
38 18497.02639 2.784333263E-02 3.4451962E+01 3.4479805E+01 4.3801367E+15 0
WRT_RST – wrote restart fields into time record = 38 / 38 mynode = 0
39 18497.02708 2.780277799E-02 3.4450270E+01 3.4478072E+01 4.3800853E+15 0
WRT_RST – wrote restart fields into time record = 39 / 39 mynode = 0
40 18497.02778 2.776701255E-02 3.4448633E+01 3.4476400E+01 4.3800332E+15 0
WRT_RST – wrote restart fields into time record = 40 / 40 mynode = 0
41 18497.02847 2.774634995E-02 3.4447054E+01 3.4474800E+01 4.3799806E+15 0
WRT_RST – wrote restart fields into time record = 41 / 41 mynode = 0
42 18497.02917 2.773250816E-02 3.4445558E+01 3.4473291E+01 4.3799273E+15 0
WRT_RST – wrote restart fields into time record = 42 / 42 mynode = 0
43 18497.02986 2.769461173E-02 3.4444193E+01 3.4471888E+01 4.3798735E+15 0
WRT_RST – wrote restart fields into time record = 43 / 43 mynode = 0
44 18497.03056 2.763508760E-02 3.4442969E+01 3.4470604E+01 4.3798192E+15 0
WRT_RST – wrote restart fields into time record = 44 / 44 mynode = 0
45 18497.03125 2.755509867E-02 3.4441883E+01 3.4469438E+01 4.3797644E+15 0
WRT_RST – wrote restart fields into time record = 45 / 45 mynode = 0
46 18497.03194 2.745308848E-02 3.4440920E+01 3.4468373E+01 4.3797095E+15 0
WRT_RST – wrote restart fields into time record = 46 / 46 mynode = 0
47 18497.03264 2.733311599E-02 3.4440058E+01 3.4467391E+01 4.3796549E+15 0
WRT_RST – wrote restart fields into time record = 47 / 47 mynode = 0
48 18497.03333 2.720839109E-02 3.4439262E+01 3.4466470E+01 4.3796005E+15 0
WRT_RST – wrote restart fields into time record = 48 / 48 mynode = 0

=======================================
= =
= STEP2D: ABNORMAL JOB END =
= BLOW UP =
= =

VMAX (M/S) = NaN
IMAX JMAX = 323 253
NODE = 19
IMAX JMAX = 80 39
IINT IEXT = 49 36

Hi, your maximum u and v values recorded in the restart files (12.48 and 8.61) are in fact quite large… Try deceasing the time step. Otherwise can you check if zeta is doing funny things? in coupled mode it can be due to the coupling of the surface pressure (what version of WRF are you using? do you pass the pressure as a coupled field?)

Dear sjullien

Thank you very much for your response! After reducing dt, the model can now run properly. I am using WRF-CROCO version 4.2, and I have enabled pressure transfer, i.e., WRF_d01_EXT_d01_PSFC → CROCO_PSFC.

However, I still have a question. I would like to use the CROCO restart file as the input file after spin-up. In the coupled workflow, it seems that I can either start a brand-new simulation or restart both WRF and CROCO together. What I want, though, is to restart CROCO while having WRF start as a new simulation. How can I achieve this? Should I rename the restart file as an ini file?

I look forward to your reply, and I wish you a pleasant day!

yes renaming your croco restart as ini file will work!