Search found 14 matches

by winterbellgirl
Mon Feb 18, 2013 7:48 pm
Forum: ROMS Discussion
Topic: 3-D stress
Replies: 0
Views: 780

3-D stress

Hi, I wanted to look at the vertical structure of stresses, and so I computed rho*AKv*du/dz (or dv/dz). The vertical gradient is computed using the finite difference method, not splines. I compared the values at the surface and bottom boundaries to sustr/bustr, which are output into history files, a...
by winterbellgirl
Fri Oct 12, 2012 7:38 pm
Forum: ROMS Discussion
Topic: pressure gradient in the diagnostic output
Replies: 4
Views: 2687

Re: pressure gradient in the diagnostic output

How about just calculate the barotropic PGF from the standard zeta output? I think it is better than changing the ROMS codes. you can get baroclinic PGF by subtracting barotropic PGF from total PGF. Thank you for your reply. Actually I did the opposite recently. I computed the baroclinic pressure g...
by winterbellgirl
Fri Jun 22, 2012 2:33 am
Forum: ROMS Discussion
Topic: pressure gradient in the diagnostic output
Replies: 4
Views: 2687

pressure gradient in the diagnostic output

Hi, Is there any way to separate the barotropic and baroclinic pressure gradient in the 2D diagnostic output? There is only one term (ubar_prsgrd/vbar_prsgrd) in the diagnostic output. I learned that this term contains both the barotropic and baroclinic parts. I found where the baroclinic pressure g...
by winterbellgirl
Thu May 10, 2012 3:04 am
Forum: ROMS Discussion
Topic: diagnostics
Replies: 7
Views: 2686

Re: diagnostics

Hi Dr. Wilkins, yes your suggestion works! This time I computed 'ubar_sustr+ubar_bustr+ubar_prsgrd+ubar_cor+ubar_hadv+ubar_hvisc-ubar_accel' and the result is zero. Thank you!
by winterbellgirl
Wed May 09, 2012 7:43 pm
Forum: ROMS Discussion
Topic: diagnostics
Replies: 7
Views: 2686

diagnostics

Hi, I'm trying to do some momentum balance analysis using ROMS. I output the 2D momentum terms into the diagnostic output files, but I'm not quite sure how the momentum balance is completed. I chose a point in the model, and added up all the terms (ubar_sustr+ubar_bustr+ubar_prsgrd+ubar_cor+ubar_had...
by winterbellgirl
Thu Jan 12, 2012 10:11 pm
Forum: ROMS Discussion
Topic: script for computing the momentum terms
Replies: 0
Views: 1085

script for computing the momentum terms

Hi, I'm wondering if anyone has a script for computing the 2-D momentum equation terms for ROMS. I know this can be realized by turning on the diagnostic output option, but my model run really takes time and I don't want to repeat it to just produce the diagnostic files. So a well-written code for m...
by winterbellgirl
Thu Oct 13, 2011 6:21 pm
Forum: ROMS Discussion
Topic: outflow at the open boundary condition
Replies: 6
Views: 2000

Re: outflow at the open boundary condition

Thanks arango, and then I can try with the volcon.
by winterbellgirl
Tue Oct 11, 2011 6:57 pm
Forum: ROMS Discussion
Topic: outflow at the open boundary condition
Replies: 6
Views: 2000

Re: outflow at the open boundary condition

Hi rduran,

I did see some discussions on VOLCON before. But i remember that somebody mentioned using VOLCON would increase the computational time by several times. This could be an alternative if there are no other solutions. Thank you for your hint!
by winterbellgirl
Tue Oct 11, 2011 6:53 pm
Forum: ROMS Discussion
Topic: outflow at the open boundary condition
Replies: 6
Views: 2000

Re: outflow at the open boundary condition

Hi Simon, Thank you for your suggestion for using the radiation boundary for M2. Actually when I increased the integration time, the surface elevation began to stabilize in the 2nd year, although the values are still high compared to observations (but the variability matches the obs well).So I guess...
by winterbellgirl
Sat Oct 08, 2011 8:57 pm
Forum: ROMS Discussion
Topic: outflow at the open boundary condition
Replies: 6
Views: 2000

outflow at the open boundary condition

Hi all, I'm running a model for the Gulf of Mexico. Now I'm trying to set up the open boundary conditions. The model has two open boundaries, the southern boundary opened at Yucatan and the eastern boundary opened at Florida strait. Since now I'm mainly focusing on the circulation of northwestern Gu...
by winterbellgirl
Fri Jul 08, 2011 1:54 am
Forum: ROMS Problems
Topic: fortran runtime error with inp_par.f90
Replies: 4
Views: 2807

Re: fortran runtime error with inp_par.f90

Hi Kate, Arango and Shchepet, Thanks for your replies! Previously the model was working fine, and after I changed the values of some parameters in the ocean.in file, such as NTIMES, NHIS, and DSTART, that error occurred. I'm using gfortran for the compilation and openmpi to do the parallel computing...
by winterbellgirl
Tue Jul 05, 2011 8:09 pm
Forum: ROMS Problems
Topic: fortran runtime error with inp_par.f90
Replies: 4
Views: 2807

fortran runtime error with inp_par.f90

Hi, I tried to submit a job and found the following error in the .err file. At line 2033 of file inp_par.f90 Fortran runtime error: Bad real number in item 1 of list input I looked at line 2033 of the inp_par.f90 file and the statements are as follows: 2031 string=Vstring(is:ie) 2032 LenS=LEN_TRIM(s...
by winterbellgirl
Thu Mar 24, 2011 7:01 pm
Forum: ROMS Discussion
Topic: How can make sure proper river point position in ROMS
Replies: 15
Views: 7307

Re: How can make sure proper river point position in ROMS

Thanks, mathiu, and I'm going to doublecheck it.
by winterbellgirl
Sat Mar 19, 2011 6:07 pm
Forum: ROMS Discussion
Topic: How can make sure proper river point position in ROMS
Replies: 15
Views: 7307

Re: How can make sure proper river point position in ROMS

The river inflow in ROMS is quite complicated. In the netcdf file of rivers you have the following arrays: double river_Xposition(river) ; river_Xposition:long_name = "river ETA-position at RHO-points" ; river_Xposition:units = "nondimensional" ; river_Xposition:field = "river_Xposition, scalar" ; ...