Ocean Modeling Discussion

Thanks, John.
Yes. I checked and plotted his and rst file, looking normal.
When I turned on the wet_and dry flag, the model just passes that blow-up points.
May I know the reasons behind that?
Thanks

Hi all, I met a problem for model blowing-up in the form of outputting NaN in the log file with ROMS3.7. My model is used in temperate shallow-water estuarine systems in East US, the grid resolution is around 400m*400m, mean-depth is less than 9 meters. The model starts in Jan 1st-XXXX, and It blows ...

Thanks John,
Problem solved by changing in
mat2roms_mw.m
line 93

projection = s.projection;
m_proj(s.projection);
switch lower(projection)

to
projection = s.projection;
m_proj(s.projection.name);
switch lower(projection.name)

Hi, I met a proble when hoping to reproduce a COAWST inlet case gird, by the script: COAWST-master\Tools\mfiles\mtools: create_roms_xygrid.m, but met the error: >> create_roms_xygrid ## Defining Global Attributes... ## Defining Dimensions... ## Defining Variables and Attributes... ## Filling ...

May I ask if ROMS could process this: Assuming I am going to model a lake and a river is in the upstream connecting the lake. But the river discharge was controlled by a water sluice, some times discharge is zero. So the water level in the late is influenced by that discharge. Could ROMS process ...

Hi I am testing the lake Jersey nesting case. However, results did show same as in the. It is a constant temperature all the time.It might due to surface wind setting. May I know how to fix that by the analytical setting? Here is the log.file and also the my results and roms website results. Thanks ...

Hi All, I set to run the sediment transport example in COAWST/Projects/Sed_floc_toy/ It contains 20 layers in the sea-bed as in the *.in file. But when I finished the run then looked the results e.g, The output variable: bed_porosity seemed unchanged, it is the same as the input file which is 0.5 ...

Hi All,
May I ask the possible and simple way for obtaining the ROMS mixed layer depth under GLS mixing scheme.
The ROMS only has surface boundary depth for application under LMD scheme, but that scheme mixes almost whole water column. Any suggestion?

Thanks

Here are small ROMS bugs for DIAGNOSTICS_BIO terms in ROMS. In def_diag.F # ifdef BIO_FENNEL DO ivar=1,NDbio2d ifield=iDbio2(ivar) IF (TRIM(var_name(i)).eq.TRIM(Vname(1,ifield))) THEN got_var(ifield)=.TRUE. DIA(ng)%Vid(ifield)=var_id(i) END IF END DO DO ivar=1,NDbio3d ifield=iDbio3(ivar) IF (TRIM ...

Hello All, I am going through the ROMS code. What I found is in the ROMS user manual, it said it will use Crank_niskson for calculating the vertical diffusion term time-stepping. And it will consistent with ROMS has a pre3d_t.F routine, then it finally gets to step3d_uv.F and step3d_t.F. But when I ...

Hello All, I am going through this ROMS Lake_Jersey example. I have a very basic question, is this lake a real lake, since I see a lot of strong surface gradient in salinity and temperature? Could lake have so high-density gradient? Or this is just for a test case. Since I want to use ROMS for some ...

Hi All, I just want to ask a question related to ROMS boundary code for processing zeta. As we know the zeta=zeta_tide+zeta_notide at the boundary and it is how rooms calculated in set_tide.F IF (LBC(ieast,isFsur,ng)%acquire.or. & & LBC(ieast,isUbar,ng)%acquire.or. & & LBC(ieast,isVbar,ng)%acquire ...

Hi Kate,
I am using TS_MPDATA advection scheme and try both K-w and GLS turbulent scheme, parameter just as Warner(2005).

Hello, I am modeling a vertical profile of a salinity in a estuary, depth around 25-30 meters. My case is that wind blow down-estuary, so the the pycnocline might be deeper, so the surface mixed layer might be about 12 meters according to measurement. But my model choosing a fix station did not ...

Hi all, I met some problem and confusion for the bustr/bvstr units. In the ROMS code: Modules/mod_forces.F: ! Bottom momentum stresses. ! ! ! ! bustr Kinematic bottom momentum flux (bottom stress) in ! ! the XI-direction (m2/s2) at horizontal U-points. ! ! bvstr Kinematic bottom momentum flux ...

Hi all, I am try to run dissolved oxygen simulation in real bay, with ROMS+fennel.h BGC part. I am try to add the NO3 and NH4 deposition from ATMOS in the model to put more nutrients to draw down the bottom dissolved oxygen. General thought, more nutrient , low DO. But strangely, what I see is high ...

Hi all, I am using ROMS with TS_MPDATA scheme to calculate the vertical salinity profile in a real estuary. The platform I choose is Mac Pro(single working station, 8 processes), OS. When I choose openMP without TS_MPDATA(not define that in CPP), it works. However, when I choose TS_MPDATA (define TS ...