Ocean Modeling Discussion

hi
During the process of making obs files, how to write u_super_obs.nc and v_super_obs.nc into one file, easily stack？
When i set OBSNAME = u_super_obs.nc| v_super_obs.nc, ROMS only read in u, not both.
-Zongchen

hey!
Last time I came up this situation is "NRREC=-1". After making ini.nc, it's ok!

hi everyone,
I run IS4DVAR models. Adjust_wstress and adjust_boundary is actived, and compiling is ok. When outer =001, Inner =000, NLM proceeds well. While TLM runs, it says ''forrtl: severe(174): SIGSEGV, segmentation fault occurred'' after *_itl.nc, *_frc.nc and *_bry.nc have been read. Anyone ...

And also,
i prepare 24 gauge data along the coast, and restore it in form of grid.Then 'd_ssh_obs.m' is changed to process the grid data. The result is, only 7 diffetent gauge data writes to the obs.nc.Why?

hi,
i come with the same question. Gauge data is assimilated to the model and i check every point to make sure the spatial locations (obs_Xgrid, obs_Ygrid, obs_Zgrid) are in in grid fractional coordinates (0:L, 0:M, 1:N) and the time coordinate (survey_time, obs_time) is in days with the same ...

hi,
i do IS4D_var experiment by assimulating only gauge sea level to ROMS, the error shows that. Before, it prompts
[*]<<<<<<<Descent Algorithm>>>>>>>
[*]CG_DELTA not positive
why?

hi,
I want to restart form rst.nc. Error:
LBA_GETATT- error while reading global attribute: NLM_LBC
in restart file: ./output/ocean_rst.nc
...
Then,i check restart file and history file.

:NLM_LBC = "\n",
"EDGE: WEST SOUTH EAST NORTH \n",
"zeta: Cha Cha Clo Cha \n",
"ubar: Fla Fla Clo Fla \n ...

Thank U， Kate！
I see it in globaldefs.h. All needed cpp options can be found here. As a beginner, everyone can face 'inquire' problems. Share to all the experience.
Error prompt: INQUIRE - unable to find requested variable: ***
$cd trunk/
$grep -n *** -r ./ !find the variable *** in the folder and ...

hi,
I want to know which cpp option active the variable 'swrad'. I try to it, but failed. Anyone can tell me? Many thanks!
**********************************************************************
Activated C-preprocessing Options:

ECS East of China Sea
ANA_BTFLUX Analytical kinematic bottom ...

hi arango，
In which my 2D_ROMS run is a cluster, and i ask for a 100G space. I think one week simulation is able to satified.
U say 'I won't guarantee the shallow-water (only 2D) equations are running with 4D-Var nowadays.' . Do u mean if i want to do the way i want, i have to modify the code and ...

hi Kate
DEBUGGING and DIS4DVAR is actived in build.bash, and recompiling is ok. H-file follows:
********************************************************************************
*/
/*
**-----------------------------------------------------------------------------
** Nonlinear basic state settings ...

hi，
No prompt error！
All log information followes:
******************************************************************************
Model Input Parameters: ROMS/TOMS version 3.6
Tuesday - September 17, 2013 - 7:54:04 PM ...

1128 'IF'
1139 'END IF' should before 1138 '# endif'

hi kate,
After deleting '#define ADJUST_BOUNDARY' AND '#define ADJUST_STFLUX' in WC13/I4DVAR/wc13.h, this bug is prompted when compiling. By the way, thinking of computation cost and memory space limits, i want to do I4DVAR test by using ROMS_2D first. So every cpp option about solve3D is deleted ...

160 allocate (Lfout(Ngrids))
161 Lbout(1:Ngrids)=1

Do u mean if Netcdf-3 is used, only Matlab version ‘2008b’ can handle it? My Matlab version is '2009a' and Netcdf-3 is used. Does it mean i have to uninstall my version and install version '2008b'?
--zongchen

hi,
i come up with the same problem. Do U settle it and how?
--zongchen

i think i almost know sth. Actually, normalization factors have been computed successfully. As to wc13_adj, it does have some problems like arango said. It needs more computation technics.

hi arango
Thank u!
Normalization Factors for prior error correlations are computed for the unbalanced increments.I totally do the test according to what 4DVar_turorial_9 says, but failed.U say "We use different strategy".That is to say wc13_nrm_*.nc in Data file is not computed as the turorial_9's ...

hi
when i try to compute normalization factors by using standard deviations, it comes up,
AD_DEF_HIS - creating adjoint file :wc13_adj.nc
forrtl: severe(174):seigsegv, sementation fault occurred
Then i page up the log file, before that, it has had prompted,
unsupported option in adjoint_based ...

relates to ARPACK/LAPACK libraries which may need to be compiled and built (depending on your existing system set-up).

i sbatch workbooks in clusters, and i don't know whether ARPACK/PARPACK has been built. After searching the lib, it's sure there is nothing! So i try to compile and build the lib ...

hey arango,
Thank U for your guidance!
This is not spam information for fast reading.
i have really tasted it! What i want to do is assimulating gauge stations' data(now just sea level) to ROMS for improving forecasting results of storm surge.i have read so many papers about the field. At the same ...

hey Alan Berry,
Thank U for your reply!
Finding ARPACK\PARPACK libraries in the trunk/Lib, so i referened them as "pwd".
ifdef USE_ARPACK
ifdef USE_MPI
# PARPACK_LIBDIR ?= /opt/intelsoft/PARPACK
PARPACK_LIBDIR ?= /export/home/wangzc/swgfs/trunk/Lib/ARPACK/PARPACK
LIBS += -L$(PARPACK_LIBDIR ...

hi everyone：
My roms is $Id 645.
I try to test 4d_var with roms,so i have customized "build" script and provided exact values for ROOT directory, MY_ROOT_DIR, ..., as file(Readme) in the test/WC13/I4DVAR/ says
Then, i compiled the program, it prompts:

/export/home/wangzc/swgfs/projects/WC13 ...

Hi, everyone:
i download roms_tools from AGRIF. Make_grid is ok, make_tide is ok(an alone nc), as well as make_NCEP. while, when I submit a task, it say Unable to get value for attribute:cycle_length in variable sms_time. I check NCEP2_Y2012M1.nc, it is cycle_length=ncdouble(0). My question is what ...

Hi, Kate. Thank you very much!
I took the your suggestion and it did work.I think this is because ROMS will not read LBC values(ubar,vbar,zeta) from *_frc.nc. Either *_bry.nc or ana_fsobc(or ana_m2obc) is needed to ROMS, right? This is also why Lm(or Mn) in the *.in equals Lm(or Mn)-2 in the *_frc ...

hello:
I am trying to simulate tide by 2D-ROMS.
I only retain sustr and sms_time in the frc.nc. At the same time, tide is added to the frc.nc. ANA_INITIAL and ANA_BMFLUX is defined.Compiling is ok, but something is wrong when running :( .
Can anyone help me?

Output info:
INITIAL: Configuring ...

Hello everyone:
I am a new ROMS user and I wonder if 2D ROMS can simulate cyclone storm surge alone. If the answer is yes, where I can find some examples?

Thank everyone's help!
Zongchen