Change SST in met_em error

[kaignuel]

Hi guys. I am trying to do some sensitivity analysis by changing the sea surface temperature.
I downloaded FNL data from NCEP and ungribbed using Vtable.GFS. I successfully ran metgrib.exe, and have obtained a set of met_em files. The control run has been successful.
After that, I tried to amend the Sea Surface temperature by amending the met_em files directly, like make the values of SKINTEMP +1 on sea area considering LANDMASK. (There is no SST variable in met_em files)
However, after I amended the met_em files and ran real.exe, the following error occurs:
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 3110
grid%tsk unreasonable
-----------------------
Are there better ways to change the initial and boundary condition, in order to avoid the error?
Thank you

 

[meteoadriatic]

Hi,
Did you verify that met_em files were properly modified as you intended, like, compared the original and modified version using ncview or similar tool?

 

[Ming Chen]

Please clarify how you change SKINTEMP in met_em.

 

[kaignuel]

I change SKINTEMP by simple python script.
dsst = 1
initial = xr.open_dataset(met_em_path)
skintemp = initial['SKINTEMP']
landmask = initial['LANDMASK']
sea_mask = landmask == 0
new_skintemp = xr.where(sea_mask,skintemp + dsst ,skintemp)
initial['SKINTEMP'] = new_skintemp

And I did verify and confirm that the files are properly verified.

 

[Ming Chen]

I am not 100% sure whether the code works as expected. But the error message from real.exe implies that somehow there is inconsistency between landuse/landmask and skintemp.