I am trying and running RegCM 4.1 with the SST data from ERA-Interim at 1.5° resolution.

I have downloaded data directly from the ERA Interim website. The downloaded file has the following attributes:

ncdump -h sstERAIN.2004_2005.nc

netcdf sstERAIN.2004_2005 {

dimensions:

longitude = 240 ;

latitude = 121 ;

time = UNLIMITED ; // (2924 currently)

variables:

float longitude(longitude) ;

longitude:units = "degrees_east" ;

longitude:long_name = "longitude" ;

float latitude(latitude) ;

latitude:units = "degrees_north" ;

latitude:long_name = "latitude" ;

int time(time) ;

time:units = "hours since 1900-01-01 00:00:0.0" ;

time:long_name = "time" ;

short sstk(time, latitude, longitude) ;

sstk:scale_factor = 0.000600706197210905 ;

sstk:add_offset = 287.908701051615 ;

sstk:_FillValue = -32767s ;

sstk:missing_value = -32767s ;

sstk:units = "K" ;

sstk:long_name = "Sea surface temperature" ;

// global attributes:

:Conventions = "CF-1.0" ;

:history = "2012-03-26 08:46:58 GMT by mars2netcdf-0.92" ;

}

I have modified the subroutine mod_sst_ersst.f90:

- with regards to the name of the input file which covers the period of interest (2004-2005)

- to change varname from "sst" to "sstk"

But I am not sure whether the data file has the correct format as I get an error when executing sst:

The error I get looks like:

Cannot get sstk from file

NetCDF: Not a valid ID

ERROR READ SST

Does anyone have any idea why I get this error?

Many thanks by advance for your response.

Best regards,

Beatrice MOREL

Energy, Electronics and Process Lab
Unit EA 4079

15, avenue René Cassin • BP 7151
97 715 Saint-Denis cedex 9
Tél.: +262(0)262 93 86 71 • Fax: +262(0)262 93 86 73
mail: secretariat.le2p@univ-reunion.fr
internet: http://sciences.univ-reunion.fr/laboratoires/le2p