Problems with stations.in

Report or discuss software problems and other woes

Moderators: arango, robertson

Post Reply
Message
Author
staalstrom
Posts: 31
Joined: Mon Feb 04, 2008 3:43 pm
Location: NIVA, OSLO, NORWAY

Problems with stations.in

#1 Unread post by staalstrom »

Hello

I want to write data from 2 positions to station output files.
I have defined stations in my cpp options.
The right staions.in file is spesified in ocean.in
But ROMS thinks ocean_sta.nc should already exist!
How do I make one?
Here's the output:

STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME

0 0 00:00:00 0.000000E+00 4.878751E+02 4.878751E+02 3.733559E+10
DEF_HIS - creating history file: ocean_his_s_0001.nc
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
DEF_STATION - inquiring stations file: ocean_sta.nc

NETCDF_OPEN - unable to open existing NetCDF file:
ocean_sta.nc
call from: def_station.F

best regards
Andre Staalstrom (ans@niva.no)

User avatar
kate
Posts: 4091
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: Problems with stations.in

#2 Unread post by kate »

Did you set LDEFOUT == T in the input file? Stations, floats, averages, and history files all share this flag. If it is F, ROMS will try to read existing files.

staalstrom
Posts: 31
Joined: Mon Feb 04, 2008 3:43 pm
Location: NIVA, OSLO, NORWAY

Re: Problems with stations.in

#3 Unread post by staalstrom »

Hi
Yes, according to ocean.out LDEFOUT=F.
But in ocean.in LDEFOUT is set to True.
There is no problems creating a new history file.
I see from some old runs without STATIONS that LDEFOUT=F even if it is set to T in ocean.in.

Is it set to F somewhere in the code?

Best regards
Andre Staalstrom (ans@niva.no)

staalstrom
Posts: 31
Joined: Mon Feb 04, 2008 3:43 pm
Location: NIVA, OSLO, NORWAY

Re: Problems with stations.in

#4 Unread post by staalstrom »

Hi again
I managed to make it work in a horrible way.
I went into ROMS/Utility/def_stations.F and changed

DEFINE : IF (ldef) THEN
CALL netcdf_create (ng, iNLM, TRIM(ncname), ncSTAid(ng))
to
DEFINE : IF (.not.ldef) THEN
CALL netcdf_create (ng, iNLM, TRIM(ncname), ncSTAid(ng))
and now the ocean_sta.nc file is created (and also ocean_his)

But still LDEFOUT is set to False even if it is set to True in ocean.in. Why?

Here are the CPPoptions I used:
#define UV_ADV /* turn ON or OFF advection terms */
#define UV_COR /* turn ON or OFF Coriolis term */
#define UV_SADVECTION
#define UV_VIS2
#define UV_QDRAG
#define TS_MPDATA
#define TS_DIF2
#define SALINITY /*
#define NONLIN_EOS
#define DJ_GRADPS
#define SOLVE3D
#define MASKING
#define SPLINES
#define ANA_BSFLUX /* analytical bottom salinity flux */
#define ANA_BTFLUX
#define ANA_SSFLUX /* analytical surface salinity flux */
#define ANA_STFLUX
#define MIX_GEO_UV
#define MIX_GEO_TS
#define MY25_MIXING
# define N2S2_HORAVG
# define KANTHA_CLAYSON /* Kantha and Clayson stability
#define SPONGE
#define RADIATION_2D
#define SOUTH_FSCHAPMAN
#define SOUTH_M2FLATHER
#define SOUTH_TRADIATION
#define SOUTH_TNUDGING
#define SOUTH_M3RADIATION
#define SOUTH_M3NUDGING
#define EASTERN_WALL
#define WESTERN_WALL
#define NORTHERN_WALL
#define TCLIMATOLOGY /* Processing of tracer climatology */
#define TCLM_NUDGING /* Nudging of tracer climatology */
#define STATIONS

brgds
Andre

User avatar
kate
Posts: 4091
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: Problems with stations.in

#5 Unread post by kate »

I wonder if there's a new bug? My floats stopped making floats files even with FLOATS defined. I'll investigate once I'm back in the office.

Post Reply