Ocean Modeling Discussion

Hi everyone,

I really appreciate all the assistance and suggestions! Although I did not explicitly define both #define UV_LDRAG and #define UV_QDRAG, I decided to test my build using the BASIN example case this time, which already has #define UV_QDRAG set by default.

For some reason, the error did ...

Hi everyone,

Thanks so much for the assistance. I really appreciate this.
To ensure that I didn’t accidentally modify something that caused the rdrag2 (or rdrag) issue, I redownloaded the latest ROMS version (4.2) today and tested the upwelling example case with only two modifications:

1. I ...

I'm currently running a model in ROMS, building it based on the example upwelling cases. However, I encountered an error message:
NLM: CHECKVARS_NF90 - unable to find model variable: NO3
in file: /home/hyc006/ideal_model/clm_ellipse/ini_ellipse_0213.nc
Found Error: 2 Line: 317 Source: ROMS/Utility ...

Thank you so much for the clarification!

I have already run a git pull to ensure I am using the latest version. In fact, I recently installed the newest version on my new server, so I don’t believe the issue stems from an outdated version.

I’ve attached my set_vbc.F file for reference. set_vbc.F ...

Thanks for your response! I’ll try to fix the issue and see if any further errors occur.

I’ve attached my set_vbc.f90 file for reference. set_vbc.f90
I noticed that if I use ANA_DRAG instead of UV_QDRAG, the build_roms process completes without errors.
However, I would still prefer to use UV ...

Thank you so much for your response—I really appreciate it!

However, the error persists even after modifying the Ana_drag file. To address the issue, I updated my model (SHALLOW) to use the same drag coefficient definitions as the upwelling example case:
#if defined UPWELLING || SHALLOW
# if ...

Thanks for the help on my previous issue! That solution worked, but now I’m encountering a new error during compilation related to rdrag2:
set_vbc.f90:274:49:

274 | bustr(i,j)=0.5_r8*(rdrag2(i-1,j)+rdrag2(i,j))* &
| 1
Error: Function 'rdrag2' at (1) has no IMPLICIT type

set_vbc.f90:289:49 ...

Hi everyone,

I recently tried to build ROMS using build_roms.sh, but I'm encountering an error that I suspect is related to defining #define UV_QDRAG in my .h file. This option worked fine in an older version of ROMS, but I'm using a newer version now and am less familiar with its structure.

Here ...

I believe I've identified why BULKFLUX prevents me from defining my own sea surface heat flux. However, I need BULKFLUX to set ANA_PAIR. I'm attempting to create a linear pressure gradient from north to south. To prevent changes in Pair due to the periodic condition, I disable boundary changes and ...

Thanks for your reply! I still got a question for the net heat flux. I've make a forcing file for the shflux for the net heat flux just like you said.

The net heat flux is variable shflux . That is what you need to set.

However, I find out that in the log file the net heat flux forcing file is ...

I'm trying to run an idealized ROMS model and want to eliminate surface heat and salt fluxes. I've successfully turned off the shortwave and longwave radiation by setting `ANA_SRFLUX` and creating a longwave radiation forcing file.

Now, I'm trying to do the same for "sensible heat flux", "latent ...

Hello ROMS Community,

I am encountering an issue with maintaining ocean stratification over time in my idealized model. Despite using climatology data and applying nudging, the ocean becomes unstratified after approximately 30 days (My total model time is about 120 days). I am confident that the ...

Thank you for your reply. I will give this a try.

I haven't used nudging in climatology before. Do you have any recommended resources on this topic? I would greatly appreciate your help!

Hsin-Yi

Hello ROMS Community,

I’m currently working on an idealized ocean model, and I’m facing an issue with maintaining the ocean stratification. Over time, the stratification disappears. So I'm trying to clamping the temperature and salinity inflow boundary conditions on the west and north edges.
The ...

I'm currently attempting to implement a model with a sponge layer, but I'm encountering confusion regarding the specification of the sponge layer's width. Below is the relevant excerpt from my code in the ana_sponge.h file:

#elif defined SHALLOW
!
! Shallow Ocean With An Island, Island Wake be ...

Greetings ROMS Discussion Community,

I'm currently working on implementing latitude-varying air pressure in my ROMS model to create a pressure gradient for driving geostrophic currents. To achieve this, I've constructed the air pressure in the ana_pair.h file as follows:

#elif defined SHALLOW ...

Hello ROMS Discussion Community,

I'm encountering a persistent issue with the sea surface height (zeta) output from my ROMS model, particularly regarding the boundary conditions. My aim is to fix the sea surface height to zero at the north and south boundaries using the Chapman boundary condition ...

This works! Thanks a lot for your help.

Hello everyone,

I'm new to the ROMS community and currently attempting to implement analytical air pressure changes in space, specifically across latitudes. To achieve this, I attempted to incorporate the latitude variable into my pair functions. Here's what I've tried so far:
CALL ana_pair_tile ...

Running with both BODYFORCE and a vertical turbulent mixing closure scheme might be inconsistent. Using MY25_MIXING (or GLS_MIXING or LMD_MIXING) assumes the stress is applied as the surface boundary condition to vertical momentum mixing.


Thanks so much for your reply! For your comments on ...

Hello everyone,

I'm currently working on configuring ROMS to include the calculation of the surface boundary layer (Hsbl) and bottom boundary layer height (Hbbl) in the output files rom_his.nc and roms_avg.nc. However, despite enabling these options, I'm unable to see the results in the output ...

Thank you for providing clarification. I have an additional question that has arisen. I'm curious about where I could define my bustr and bvstr variables. It appears that it's not as straightforward as editing the ana_smflux.h file, unlike sustr and svstr.

Dear ROMS community,

I hope this message finds you well. I'm currently in the process of incorporating the force due to sea surface tilt as a body force in my ROMS model. I've initiated this process by defining 'BODYFORCE' in my .h file. Additionally, I've made adjustments to the parameter LEVSFRC ...

I truly appreciate the time and effort you took to assist me with my query. Your solution worked perfectly for me, and I was able to manage the force now. Thank you!

Hello everyone,

I'm a new user on this platform and currently in the process of building my own model based on the upwelling testing case. However, I've encountered an issue that I'm seeking assistance with.

Specifically, I want to disable the westward wind force in my model's forcing file. I ...

Hello ROMS community,

I'm encountering an issue with the output file generation in ROMS and would appreciate some assistance. After running the model, I've noticed that no output files are being generated in the designated folder. I've checked the roms.in file and ensured that the output file ...

Hello, I'm a newcomer to ROMS and currently learning the basics of setting up ROMS. I apologize for posing what might be considered a fundamental question. As I work on creating my roms_model.in file based on the roms_upwelling.in test cases, I've encountered a challenge with the lateral boundary ...

Got it! Thank you for your guidance!

Thanks a lot for your prompt reply! I truly appreciate your assistance.

Your suggestion to use c_grid for creating a new grid was helpful. Now, I have a follow-up question. If I make changes to specific parameters, such as f, mask_v, or lon_rho these kind of variables, within the grid structure ...

Hi ROMS community,

I'm a new user in ROMS and I'm currently working on setting up an idealized grid for the ROMS model using MATLAB and have run into a bit of a roadblock. I'm trying to save the grid structure to a NetCDF file, as required by ROMS.

On the ROMS wiki [https://www.myroms.org/wiki ...

I think I got what you mean. I've been trying to create a new gid. Is there a Matlab function for me to create it?

Additionally, I've seen that "write_contact.m" is also for a nested grid. So I was confused is there a function for helping me to write the grid I've made to save as the nc file I need ...

Hello, I'm encountering some issues with creating an idealized grid for my ROMS model using MATLAB scripts from [https://www.myroms.org/wiki/Grid_Processing_Scripts](https://www.myroms.org/wiki/Grid_Processing_Scripts). When attempting the `c_contact("grid_ideal.nc",1,1)` command, I encountered the ...