Ocean Modeling Discussion

Hi John, I think I added two river inflows, each with a discharge of 𝑄=5 m^3/s, primarily directed in the u-direction. However, I have noticed two different velocity directions in the model, and I am uncertain whether these indicate inflows as expected. I have attached the following files for your re ...

Dear John, Thank you for your complete reply. I am currently working on a model where I’ve added two freshwater rivers, each with a discharge of 5 m³/s. The model is set to run for 40 days with a time step of 40 seconds. Some relevant conditions include: The lake's salinity is set to 14 PSU. The su ...

Thank you John,
I attached my roms files. It is a simple test case. When I defend outflow, Salinity decreases and vice versa.

Thank you for your reply.
I did same settings for river definition. Salinity increased again. Is there any solution for increasing salinity?

I defined rivers both in lateral boundary and mask areas (mask = 0), but I still have the same problem with LuvSrc. I tried LwSrc (mask = 1), and the salinity increased (I attached a salinity vertical profile close to the river mouth). I changed my settings and defined one inflow and one outflow ...

Thank you for your reply. There is no place with a salinity of less than 14 PSU. The river salinity is 0 PSU, and the river temperature is 4 degrees C. Based on velocity fields, river directions are correct. I changed my f=0, but still, i have the same problem with salinity. I attached velocity ...

Dear all, I am working with an axisymmetric idealized lake. I added two freshwater river inflows with a constant temperature of 4 degrees centigrade (I activated LUVSrc and temperature and salinity tracers). The lake surface temperature is set to freezing point temperature, and the lake's salinity ...

Dear John, I want to run an idealized ice-covered lake with a homogenous salinity of 0 PSU (I did not add an ice model). The surface temperature was set to freezing point temperature, with an initial linear stratification within the water column to a temperature of maximum density at the bottom. My ...

Dear all, I want to run an idealized ice-covered lake with a homogenous salinity of 0 psu (I did not add an ice model). The surface temperature was set to freezing point temperature, with an initial linear stratification within the water column to a temperature of maximum density at the bottom ...

Thank you for your thoughtful response.
I added solar radiation effects to the model, but I couldn't keep the surface temperature at 0 degrees throughout the entire model run. Is there any way to achieve this?

Hi all, Thank you, John and Kate. I run a model with two river inflows. I added solar radiation effects to this model (#define SOLAR_SOURCE). I cannot run it because the model needs a forcing file (It works without river inflow). I defined #define ANA_SMFLUX #define ANA_STFLUX #define ANA_SSFLUX ...

Dear kate,
Thanks for getting back to me.
Are these phyton files on your GitHub?
I want to add some fresh water to my idealized lake. Adding several rivers is the best option. Can I add Qs? I cannot add a file like precipitation because it is an ice-covered lake.

Dear John and Kate, Thank you for your reply. I added two rivers (using create_roms_rivers.m in COAWST m files) to my idealized ice-covered lake with a homogenous salinity of 14 psu. The surface temperature was set to freezing point temperature, with an initial linear stratification within the water ...

Dear John,
I want to define a river in my idealized axisymmetric lake rather than as a point. I want to have distributed river input across the whole lake surface. Is it possible, or must I define it as source/sink points?

Best regards,
Fatima

Dear all, I want to add The decay law for the flux of solar radiation in my model. As you know, it reduces with depth exponentially. I = I0 * e^(-k*d) where I is the intensity of solar radiation at a given depth in the water, I0 is the intensity of radiation at the surface of the ice, k is the ...

Dear John,
Thank you for your answer. Here is my complete output for both fedora Linux and server.
Best regards,
Fatima

Dear Arango, I switched to the HLRN server and ran my model with export USE_MPI=on export USE_MPIF90=on export FORT=gfortran I still have this problem. I can't understand why after the first timestep my KEchar = NaN and PEchar =NaN? Could you please help me to fix this problem? NL ROMS/TOMS: started ...

Dear Arango,
Thank you for your answer. I am using gfortran and fedora (35) Linux. I did not change the definition. How can I solve the problem?
I have access to intel Fortran and a server. Do you think I should switch to the server and intel Fortran, or can I solve this problem?
Best regards,
Fatima

Dear all, I had a runtime error when I ran my model (after the first timestep, this error happened, but my model was not blown up.). Could you please help me to fix it? ran1.f90:150:51: runtime error: signed integer overflow: -1832521919 + -2138245891 cannot be represented in type 'integer(kind=4 ...

Dear all, When I run my model after first timestep, I had this error: ran1.f90:150:51: runtime error: signed integer overflow: -1832521919 + -2138245891 cannot be represented in type 'integer(kind=4)' But my model was run. After 20 timesteps were blow-up with this error Note: The following floating ...

Hi I try to run my model in fedora 30 and i used gfortran without mpi. I have this error Program received signal SIGSEGV: Segmentation fault - invalid memory reference. Backtrace for this error: #0 0x7f92562bfd51 in ??? #1 0x7f92562bef15 in ??? #2 0x7f92555fbf3f in ??? #3 0x7f92565b17f7 in ??? #4 ...

Dear
I have same problem, do you find a correct way to solve it?