MPAS | ERA5| Update SST and sea-ice fraction

Generating SST and sea-ice update files

Because MPAS-Atmosphere — at least, when run as a stand-alone model — does not contain prognostic equations for the SST and sea-ice fraction, these fields would remain constant if not updated from an external source; this is, of course, not realistic, and it will generally impact the quality of longer model simulations. Consequently, for MPAS-Atmosphere simulations longer than roughly a week, it is typically necessary to periodically update the sea-surface temperature (SST) field in the model. For real-time simulations, this is generally not an option, but it is feasible for retrospective simulations, where we have observed SST analyses available to us.

Prepare intermediate files containing SST and sea-ice analyses

For ERA5

Intermediate files

Reading intermediate files

1
2
3
4
5
6
#!/bin/bash
source /home/wpsze/WRF/WRFv440/wrf_env.sh
ln -sf /home/wpsze/WRF/WRFv440/wrf_install/WPS-4.4/bin/rd_intermediate.exe
ln -sf /home/wpsze/WRF/WRFv440/wrf_install/WPS-4.4/util/plotfmt.ncl
./rd_intermediate.exe $1
ncl plotfmt.ncl $1

and, on termianl,

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
$ ./read_int.sh ERA5:2022-06-30_00
......
================================================
FIELD = SEAICE
UNITS = fraction DESCRIPTION = Sea-Ice Fraction
DATE = 2022-06-30_00:00:00 FCST = 0.000000
SOURCE = ECMWF
LEVEL = 200100.000000
I,J DIMS = 1440, 721
IPROJ = 0  PROJECTION = LAT LON
  REF_X, REF_Y = 1.000000, 1.000000
  REF_LAT, REF_LON = 90.000008, 0.000000
  DLAT, DLON = -0.250000, 0.250000
  EARTH_RADIUS = 6367.470215
DATA(1,1)=0.927246

================================================
FIELD = SST
UNITS = K DESCRIPTION = Sea-Surface Temperature
DATE = 2022-06-30_00:00:00 FCST = 0.000000
SOURCE = ECMWF
LEVEL = 200100.000000
I,J DIMS = 1440, 721
IPROJ = 0  PROJECTION = LAT LON
  REF_X, REF_Y = 1.000000, 1.000000
  REF_LAT, REF_LON = 90.000008, 0.000000
  DLAT, DLON = -0.250000, 0.250000
  EARTH_RADIUS = 6367.470215
DATA(1,1)=271.459961
================================================
......

Plotting Intermediate Files: plotfmt.ncl

Plots the fields in the ungribbed intermediate files

$ ncl plotfmt.ncl 'filename="FNL:2007-09-15_00"'

See above shell script.

skintemp (skt) vs sst

Consider skt - sst form ERA5 grib file, they are different but why?

  • ERA5-0p25-SL-2022063000.grib
  • ERA5 grib: skt - sst
1
2
3
4
5
6
ncvis ERA5-0p25-SL-2022063000.nc
cdo selname,sst ERA5-0p25-SL-2022063000.nc sst.nc
cdo selname,skt ERA5-0p25-SL-2022063000.nc skt.nc
cdo chname,skt,sst skt.nc skt-sst.nc
ncdiff skt-sst.nc sst.nc diff.nc
ncvis/ncview diff.nc
ERA5 grib | skt - sst |
ERA5 grib | skt - sst |

Solutions

  • SST input resolution and LANDMASK | Feb 2024
    • SST in ERA5 only has valid values over the ocean. When metgrid conducts horizontal interpolation, it may result in some discontinuity.
    • TSK is skin temperature predicted by land model (over the ocean without sst update, it remains as skintemp at the initial time. if sst_update =1, tsk is updated by SST.)
    • I am suspicious that SST is somehow not continuous over the coast, which affects tsk later.
    • To overcome this problem, please use skintemp as SST. You cam simply delete the line for SST in Vtable.ECMWF , then rerun ungrib and metgrid.
    • I would like to confirm that in ERA5, SKINTEMP and SST over the ocean are same. This will justify that we can use SKINTEMP to replace SST.
  • Metgrid SST interpolation artifacts | Nov 2023
    • One quick fix of the unrealistic SST along coastal areas is to replace SST by SKINTEMP. Please delete the line below from Vtable.ECMWF, then rerun ungrib.exe.
    • 34 | 1 | 0 | | SST | K | Sea-Surface Temperature |
    • In this case SST will not be extracted from ERA5 and SKINTEMP will be used as SST. We believe this is reasonable because SST and SKINTEMP is almost the same over ocean.
    • In the case you attempt to use external SST instead of ERA5 SST, you need to provide landmask specifically for the external SST data, then modify METGRID.TBL so that metgrid.exe will know what landmask is used for horizontal interpolation.
      • A detailed explanation can be found at https://www2.mmm.ucar.edu/wrf/users/tutorial/presentation_pdfs/201907/duda_wps_advanced.pdf, and the slides 47-58 give an example how to process external soil data. The same approach is also applied to external SST.
      • the default fill-values overland is 0K (as in -273.15ºC). If you're running in high res, the interpolation of SST leads to some wonky numbers around the landmask. Probably you can set fill_missing = 285. or some other value to remedy the problem.
    • does WRF needs TSK from global data | Nov 2020
    • Strong non-smooth temperature pattern over snow/ice/water land categories | Sep 2021
      • I see you're using sst_update. Are you using a separate SST input data than the ERA5 data? If so, is it possible that the landsea mask for the two data types is inconsistent? Otherwise, there is a particular process necessary for processing landsea for ERA5 data.
      • One suggestion is to use SKINTEMP only (instead of SST)
  • Issue with ERA5 SST fields | Feb 2023
    • This is an old issue. If you replace the SST by the skin temperature you will have unrealistic warm waters near the coast in summer.
  • WRF regional climate simulation | Sep 2021
    • This will potentially lead to problem in WRF simulation, especially for long term climate simulation with sst_update.
    • To solve this issue, please use skintemp to replace SST following the steps below:
        1. in Vtable.ERA-I, delete the line below:
        • 34 | 1 | 0 | | SST | K | Sea-Surface Temperature |
        1. rerun ungrib.exe and metgrid.exe
        1. run real.exe and wrf.exe
  • Merge Vtable.GFS and Vtable.SST | Jan 2021
    • It is perfectly fine that the SST field is derived from SKINTEMP; however, if you are not using high-resolution time-varying SST data input, there really is no need to use the sst_update option. If you are only running for 4 days, it shouldn't be necessary.
  • Ungrib error: not reading SST | Dec 2022
    • This means the code is not familiar with a level code of 160, and if you want to use that one, you'll need to modify the ungrib/src/rd_grib2.F file to allow it to work. You can search for that part of the code by looking for "Rd_grib2 does not know" in that file. After you make any modifications, you'll need to recompile ungrib.exe.

Summary

The ERA5.grib can be IC, LBC, FDDA and SST/Sea-Ice updated as well.

And, ERA5 corresponding interemediate file can be used.

Vtable

  • Vtable.GFS
    • Skin temperature (can use for SST also) 
      1
      2
      3
      4
      5
      6
      7
      8
      9
      GRIB1| Level| From |  To  | metgrid  | metgrid | metgrid                                 |GRIB2|GRIB2|GRIB2|GRIB2|
      Param| Type |Level1|Level2| Name     | Units   | Description                             |Discp|Catgy|Param|Level|
      -----+------+------+------+----------+---------+-----------------------------------------+-----------------------+
        91 |   1  |   0  |      | SEAICE   | proprtn | Ice flag                                | 10  |  2  |  0  |   1 |
        81 |   1  |   0  |      | LANDSEA  | proprtn | Land/Sea flag (1=land, 0 or 2=sea)      |  2  |  0  |  0  |   1 |
        81 |   1  |   0  |      | LANDN    | proprtn |                                         |  2  |  0  | 218 |   1 |
         7 |   1  |   0  |      | SOILHGT  | m       | Terrain field of source analysis        |  0  |  3  |  5  |   1 |
        11 |   1  |   0  |      | SKINTEMP | K       | Skin temperature                        |  0  |  0  |  0  |   1 |
      ......

  • Vtable.ERA-interim.pl 

    1
    2
    3
    4
    5
    6
    7
    GRIB | Level| Level| Level| metgrid  |  metgrid | metgrid                                  |
    Code | Code |   1  |   2  | Name     |  Units   | Description                              |
    -----+------+------+------+----------+----------+------------------------------------------+
     235 |  1   |   0  |      | SKINTEMP | K        | Sea-Surface Temperature                  |
      31 |  1   |   0  |      | SEAICE   | fraction | Sea-Ice Fraction                         |
      34 |  1   |   0  |      | SST      | K        | Sea-Surface Temperature                  |
    ......

  • Vtable.SST 

    1
    2
    3
    4
    5
    6
    7
    GRIB | Level| Level| Level| metgrid  |  metgrid | metgrid                                  |
    Code | Code |   1  |   2  | Name     |  Units   | Description                              |
    -----+------+------+------+----------+----------+------------------------------------------+
     172 |  1   |   0  |      | LANDSEA  | 0/1 Flag | Land/Sea flag                            |
     235 |  1   |   0  |      | SKINTEMP | K        | Sea-Surface Temperature                  |
      31 |  1   |   0  |      | SEAICE   | fraction | Sea-Ice Fraction                         |
    -----+------+------+------+----------+----------+------------------------------------------+

For sst = skintemp

  • Vtable.ERA-interim.pl is modified as Vtable.ERA5
  • replace the SST by the skin temperature !!
  • Remove the row : 34 | 1 | 0 | | SST | K | Sea-Surface Temperature |

For using ERA5 as SST source

  • use ERA5's skintemp as SST source instead of sst itself
  • modify Vtable.SST with replace 11 by 235

  • try to keep three varibles (not only SST and SEAICE)

  • Vtable.ERA5.SST 

    1
    2
    3
    4
    5
    6
    GRIB | Level| Level| Level| metgrid  |  metgrid | metgrid                                  |
    Code | Code |   1  |   2  | Name     |  Units   | Description                              |
    -----+------+------+------+----------+----------+------------------------------------------+
     235 |  1   |   0  |      | SKINTEMP | K        | Sea-Surface Temperature                  |
      31 |  1   |   0  |      | SEAICE   | fraction | Sea-Ice Fraction                         |
     235 |  1   |   0  |      | SST      | K        | Sea-Surface Temperature                  |

ln -sf ${TEMPLATE_DIR}/Variable_Tables/Vtable.ERA5.SST Vtable

For MPAS

code

mpas_stream_manager.F

  • MPAS_stream_mgr_get_property
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
interface MPAS_stream_mgr_get_property
    module procedure MPAS_stream_mgr_get_property_int
    module procedure MPAS_stream_mgr_get_property_char
    module procedure MPAS_stream_mgr_get_property_logical
end interface

!-----------------------------------------------------------------------
!  routine MPAS_stream_mgr_get_property_int
!
!> \brief Gets a property of a stream in an MPAS stream manager.
!> \author Michael Duda, Doug Jacobsen
!> \date   13 June 2014
!> \details
!>  Retrieves the value of a stream property within an MPAS stream manager.
!>  NOTE: This routine does not support streamID regular expressions
!
!-----------------------------------------------------------------------
subroutine MPAS_stream_mgr_get_property_int(manager, streamID, propertyName, propertyValue, direction, ierr)

MPAS init.nc SST issue

  • sst = skintemp that sst values is from skintemp.
  • src/core_init_atmosphere/mpas_init_atm_cases.F
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
if (config_met_interp) then

!ldf (2011-11-19): added initialization of the sea-surface temperature, seaice fraction, and
!seaice flag:
  sst = 0.0
  xice = 0.0
  seaice = 0.0
!ldf end.

! Set SST based on SKINTEMP field if it wasn't found in input data
if (minval(sst) == 0.0 .and. maxval(sst) == 0.0) then
    call mpas_log_write('Setting SST from SKINTEMP')
    !where (landmask == 0) sst = skintemp
    sst = skintemp
end if

But, ERA5's skintemp and sst are different as above.

Left) init.nc, Right) sst.nc init.nc - sst.nc, and same as above figure of ERA5 grib file It is init.nc: skintemp - sst that has shown sst is same as skintemp in init_atmosphere step in MPAS.

core_atmosphere

  • src/core_atmosphere/physics/mpas_atmphys_manager.F
  • read sfc input 
    1
    2
    3
    type(mpas_pool_type),pointer:: sfc_input
    block => domain % blocklist
    call mpas_pool_get_subpool(block%structs,'sfc_input'   ,sfc_input   )

  • physics_update_sst 

    1
    2
    3
    4
    5
    6
    7
    !update surface boundary conditions with input sea-surface temperatures and fractional
    !sea-ice coverage:
    if(mpas_is_alarm_ringing(clock,sfcbdyAlarmID,ierr=ierr)) then
        call mpas_reset_clock_alarm(clock,sfcbdyAlarmID,ierr=ierr)
        if(config_sst_update) &
          call physics_update_sst(domain%dminfo,config_frac_seaice,mesh,sfc_input,diag_physics)
    endif

  • MPAS_stream_mgr_get_property read sst.nc files

1
2
3
4
5
6
7
8
9
10
!set alarm for updating the surface boundary conditions:
 if (config_sst_update) then
    call MPAS_stream_mgr_get_property(stream_manager, 'surface', MPAS_STREAM_PROPERTY_RECORD_INTV, stream_interval, &
                                      direction=MPAS_STREAM_INPUT, ierr=ierr)
    call mpas_set_timeInterval(alarmTimeStep,timeString=stream_interval,ierr=ierr)
    alarmStartTime = startTime
    call mpas_add_clock_alarm(clock,sfcbdyAlarmID,alarmStartTime,alarmTimeStep,ierr=ierr)
    if(ierr /= 0) &
       call physics_error_fatal('subroutine physics_init: error creating alarm sfcbdy')
 endif
  • src/core_atmosphere/Registry.xml 
    1
    <var_struct name="sfc_input" time_levs="1">

Update SST on MPAS time-step,

  • For uptated interval = 3 hour,
    • T=0, from SST(T=0) of init.nc
    • T=0 to T=3, are same as SST(T=0)
    • T>3 to T=6, are same as SST(T=3)
    • T>6 to T=9, are same as SST(T=6)
  • For uptated interval = 6 hour,
    • T=0, from SST(T=0) of init.nc
    • T=0 to T=6, are same as SST(T=0)
    • T>6 to T=12, are same as SST(T=6)
    • T>12 to T=18, are same as SST(T=12)

src/core_init_atmosphere/Registry.xml

1
2
3
4
5
6
7
8
9
10
11
12
13
<package name="sfc_update" description="Used by test cases that produce surface updates."/>

<nml_option name="config_input_sst"             type="logical"       default_value="false"
      units="-"
      description="Whether to re-compute SST and sea-ice fields from surface input data set; should be set to .true. when running case 8"
      possible_values="true or false"/>

<var name="sst" type="real" dimensions="nCells Time" units="K"
        description="sea-surface temperature"
        packages="met_stage_out;sfc_update"/>
<var name="xice" type="real" dimensions="nCells Time" units="unitless"
        description="fractional area coverage of sea-ice"
        packages="met_stage_out;sfc_update"/>

namelist.init_atmosphere

The key namelist options that must be set are shown below; other options can be ignored.

Note in particular that we have set the config_init_case variable to 8! This is the initialization case used to create surface update files, instead of real-data initial conditions files.

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
&nhyd_model
    config_init_case = 8
    config_start_time = '2014-09-10_00:00:00'
    config_stop_time = '2014-09-20_00:00:00'
/

&data_sources
    config_sfc_prefix = 'SST' # SST or ERA5: the prefix of the intermediate data files containing SST and sea-ice
    config_fg_interval = 86400
/

&preproc_stages
    config_static_interp = false
    config_native_gwd_static = false
    config_vertical_grid = false
    config_met_interp = false
    config_input_sst = true
    config_frac_seaice = true
/
&decomposition
    config_block_decomp_file_prefix = 'x1.10242.graph.info.part.'
/

streams.init_atmosphere

We have set both the filename_template and output_interval attributes of the "surface" stream. The output_interval should match the interval specified in the namelist for the config_fg_interval variable. The other streams ("input", "output", and "lbc") can remain unchanged — the input file should still be set to the name of the static file.

  • Noted: can set filename_template="sst.$Y-$M-$D_$h.$m.$s.nc"
1
2
3
4
5
6
<immutable_stream name="surface"
                  type="output"
                  filename_template="x1.10242.sfc_update.nc"
                  filename_interval="none"
                  packages="sfc_update"
                  output_interval="86400"/>

results

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 $ ncdump -h sst.2022-06-30_00.00.00.nc
netcdf sst.2022-06-30_00.00.00 {
dimensions:
	StrLen = 64 ;
	Time = UNLIMITED ; // (1 currently)
	nCells = 40962 ;
variables:
	char xtime(Time, StrLen) ;
		xtime:units = "YYYY-MM-DD_hh:mm:ss" ;
		xtime:long_name = "Model valid time" ;
	float sst(Time, nCells) ;
		sst:units = "K" ;
		sst:long_name = "sea-surface temperature" ;
	float xice(Time, nCells) ;
		xice:units = "unitless" ;
		xice:long_name = "fractional area coverage of sea-ice" ;

Plot

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
#!/usr/bin/env python

import Ngl
import numpy as np
from netCDF4 import Dataset
import math
import sys


if __name__ == "__main__":

    #
    # Get the name of the file containing the static information
    #
    if len(sys.argv) < 2 or len(sys.argv) > 3:
        print('')
        print('Usage: '+sys.argv[0]+' [mesh filename] <field filename>')
        print('')
        exit(0)

    #
    # The (lat,lon) the plot is to be centered over
    #
    cenLat   = 0.0
    cenLon   = 0.0

    #
    # Projection to use for plot
    #
    projection = 'CylindricalEquidistant'


    r2d = 180.0 / math.pi             # radians to degrees

    if len(sys.argv) == 2:
        f = Dataset(sys.argv[1])
        g = f
    elif len(sys.argv) == 3:
        f = Dataset(sys.argv[2])
        g = Dataset(sys.argv[1])

    rlist = Ngl.Resources()
#    rlist.wkWidth = 1200
#    rlist.wkHeight = 1200

    rlist.wkColorMap = 'gui_default'

    wks_type = 'png'
    wks = Ngl.open_wks(wks_type, 'delta_sst', rlist)

    lonCell = g.variables['lonCell'][:] * r2d
    latCell = g.variables['latCell'][:] * r2d

    res = Ngl.Resources()

    res.sfXArray             = lonCell
    res.sfYArray             = latCell

    res.cnFillMode           = 'AreaFill'

    res.cnFillOn             = True
    res.cnLinesOn            = False
    res.cnLineLabelsOn       = False

    res.cnInfoLabelOn        = True

    res.lbLabelAutoStride    = True
    res.lbBoxLinesOn         = False

    res.mpProjection      = projection
    res.mpDataBaseVersion = 'MediumRes'
    res.mpCenterLatF      = cenLat
    res.mpCenterLonF      = cenLon
    res.mpGridAndLimbOn   = True
    res.mpGridAndLimbDrawOrder = 'PreDraw'
    res.mpGridLineColor   = 'Background'
    res.mpOutlineOn       = True
    res.mpDataBaseVersion = 'Ncarg4_1'
    res.mpDataSetName     = 'Earth..3'
    res.mpOutlineBoundarySets = 'Geophysical'
    res.mpPerimOn         = True
    res.mpLimitMode = 'LatLon'
    res.mpMinLonF = -180.0
    res.mpMaxLonF = 180.0
    res.mpMinLatF = -90.0
    res.mpMaxLatF = 90.0

    res.cnLevelSelectionMode = 'ManualLevels'
    res.cnMinLevelValF = -2.0
    res.cnMaxLevelValF = 2.0
    res.cnLevelSpacingF = 0.10
    res.lbAutoManage = False
    res.lbOrientation = 'Horizontal'
    res.lbBoxEndCapStyle = 'TriangleBothEnds'
    res.lbLabelAngleF = 90.0
    res.lbLabelFontHeightF = 0.01

    res.mpFillOn              = True            # Turn on map fill.
    res.mpFillAreaSpecifiers  = ['Land']
    res.mpSpecifiedFillColors = [0]
    res.mpAreaMaskingOn       = True            # Indicate we want to 
    res.mpMaskAreaSpecifiers  = ['Water']
    res.cnFillDrawOrder       = 'Predraw'       # Draw contours first.

    fld = f.variables['sst'][10,:] - f.variables['sst'][0,:]
    map = Ngl.contour_map(wks, fld, res)

    Ngl.end()

  • or $ ncdiff sst.nc init.nc diff-sst.nc

src/core_atmosphere/Registry.xml

1
2
3
4
<nml_option name="config_sst_update" type="logical" default_value="false"
      units="-"
      description="logical for configuration of sea-surface temperature"
      possible_values=".true. for time-varying sea-surface temperatures; .false., otherwise"/>

streams.atmosphere

  • stream_list.atmosphere.surface
1
2
3
4
5
6
7
8
9
<stream name="surface"
        type="input"
        filename_template="x1.10242.sfc_update.nc"
        filename_interval="none"
        input_interval="86400">

    <file name="stream_list.atmosphere.surface"/>

</stream>

namelist.atmosphere

1
2
&physics
    config_sst_update = true

Plot diff

For WRF

The WRF model physics do not predict sea-surface temperature, vegetation fraction, albedo or sea ice. For long simulations, the model provides an alternative to read-in the time-varying data and update these fields (the sst_update option). In order to use this option, one must have access to time-varying SST and sea ice fields. Twelve monthly values of vegetation fraction and albedo are available from the geogrid program. Once these fields are processed via WPS, one may activate the sst_update option within the namelist.input file.

You will need to input your SST data during the WPS process. Take a look at the SST Example from our WRF Online Tutorial for information on how to do this. Because your SST data are in netcdf format, however, you will need to write them to intermediate format in order for metgrid to incorporate the data. To do so, you can follow the instructions in Chapter 3 of the WRF Users' Guide. You can also find an example script for writing a netcdf file to intermediate format here. Just note that this script was provided to us from a user who wrote it specifically for their application, so you will need to make modifications for your own use.

  • SST's are typically added to the model:
      1. Use the SST at the initial time as a constant field for all time periods (this is good for short runs, like real-time runs, where SST is not updated during the WRF model run)
      1. As an extra input at each model input time (this is good for long-months-model runs)
      1. SST:* intermediate files generated from the SST input data (if you do not have these files, first run ungrib.exe)
      • ln -sf ungrib/Variable_Tables/Vtable.SST Vtable

Interpolate the input data onto our model domain (metgrid.exe)

1
2
3
4
start_date = '2016-10-06_00:00:00',
end_date = '2016-10-08_00:00:00',
interval_seconds = 21600,
fg_name = 'FILE', 'SST',
  • ./metgrid.exe
  • generate met_em.d01.xxx.nc

Run model

  • Run real.exe
1
2
3
4
5
6
7
&time_control
auxinput4_inname = "wrflowinp_d<domain>",
auxinput4_interval = 360,
io_form_auxinput4 = 2
 
&physics
sst_update = 1,

Note: Do not change the syntax "wrflowinp_d<domain>", to "wrflowinp_d01". The syntax should be left exactly as above. Input interval is in minutes.

  • generate one more file as wrflowinp_d01 that contains SST, VEGFRA, ALBBCK and SEAICE (if available) for each input time.
  • Run wrf.exe
  • To check if SST is updated in the model, look at field TSK over water.

WRF&MPAS-A Support Forum

  • MPAS tutorial questions about ungrid and SST data | Jun 2023
      1. Starting from WPSV4.3, a configure option, --nowrf, has been added to allow for configuration of the WPS without reference to a compiled WRF model. if you only need the ungrib component of the WPS, you can compile WPS with this option.
      1. Most of the reanalysis products like GFS, ERA5 etc include SST. When you run ungrib to extract variables from the reanalysis data, SST will be included in the intermediate file.
  • SST update on the lakes | Jul 2023
    • SST_UPDATE only works over water points, i.e., those grid points where landmask =0. Those points with lakemask =1 should correspond to points with landmask=0, and thus sst_update should apply to such points.
  • MPAS SST update file | Feb 2020
  • Updating SST not working | Sep 2018
    • forrtl: severe (64): input conversion error, unit -5, file Internal List-Directed Read
    • The issue seems to be that parts of the timestamp string ('xtime') that are read from input files contain garbage, which confuses the parsing code in the timekeeping module.
MPAS
#NWP #WRF #MPAS #ERA5 #PGW #GFS #GRIB #NetCDF #CMIP6
MPAS | ERA5| Update SST and sea-ice fraction
https://waipangsze.github.io/2025/02/14/MPAS-ERA5-Update-SST-and-sea-ice-fraction/
Author
wpsze
Posted on
February 14, 2025
Updated on
March 18, 2025
Licensed under