Ocean Modeling Discussion

zhaoqian

Hi Kate and all :
the varinfo of dye in clm_ts.cdl is
float dye_01(dye_time, s_rho, eta_rho, xi_rho) ;
dye_01:long_name = "COD" ;
dye_01:time = "dye_time" ;
dye_01:coordinates = "lon_rho lat_rho s_rho dye_time" ;
float dye_02(dye_time, s_rho, eta_rho, xi_rho) ;
dye_02:long_name = "N" ;
dye ...

zhaoqian

Hi Kate and all, I have known the reason now,but I can't solve it yet.
In get_data.F, the code getting the climatology data is
# if !defined ANA_TCLIMA && defined TCLIMATOLOGY
DO i=1,NAT
CALL get_3dfld (ng, iNLM, idTclm(i), ncCLMid(ng), 1, &
& CLMname(ng), update(1), &
& LBi, UBi, LBj, UBj, 1 ...

zhaoqian

Hi all:
In order to nudge the tracers（active and passive） , first I fill the climatology date into the file of clm_ts.nc which is created by clm_ts.cdl file, then I define the Climatotogy variable of dye in varinfo.dat as follows:
'dye_'
'dye climatology'
'kilogram meter-3'
'dye_, scalar ...

zhaoqian

Hi all:
If I want to nudge one kind of passive traces during simulating, then first I should make the climatology file of this trace. Should I define the structure of this climatology variable in varinfo.dat file ?
Thank you!

zhaoqian

Hi Kate:
So you want to run with momentum and active tracers, then later with passive/bio tracers?

Yes，I want to run with momentum and active tracers(T、S),store the data, and then run with passive tracers(N,cod) using these data.

I think my branch can do that while it would be harder with the ...

zhaoqian

Hi ,
I want to simulate the hydrodynamical fields first, and then use the output date to drive the tracer module.
Thank you!

zhaoqian

jcwarner wrote:i am not quite clear about what u mean 'lower boundary of sediment'

Sorry. I means lower boundary condition of sediment. I don't know where the code is, in set_tide.F ( which calls ana_btflux or the forcing given ) or somewhere in sediment.F ?.

zhaoqian

Thank you for solving my first question.
I am confused with the lower boundary of sediment transport. Where is it defined , in set_tide.F ( which calls ana_btflux or the forcing given ) or somewhere in sediment.F ?

zhaoqian

Hi all,
When dealing with transport equation of sediment, (1) Should I define SUSPLOAD and TS_MPDATA as well as SEDIMENT? (2)What's the function of subroutine sediment.F and step3d_t.F?
Thanks in advance !

zhaoqian

Hi all:
Where are the following two cpps(marked in red) in sediment.F defined?
I search them if cppdefs.h and globaldefs.h, but can't find them.
Thanks for help!

!-----------------------------------------------------------------------
! Compute sediment bed stratigraphy ...

zhaoqian

Hi viviocean and all:
What is the program of your doing tidal harmonic analysis, t-tide or others? Do they take the parameters discussing above into account. I think maybe different program will lead to different results.

zhaoqian

I have found them in varinfo.dat which are defined as dye_east ...,mud_east_ ...,sand_east_ ... .

zhaoqian

Hi all:
I want to create bry file for dye and sediment, but I can't find the boundary variables information of dye and sediment in varinfo.dat. So I dont't konw how to do it! Please give me some instructions!
Thanks a lot!

zhaoqian

First, Thanks mathieu a lot for your reply. I have read set_data.F and understood something.
According to your answer, when the NUDGING_SST and NUDGING_T flags are actived, we should set the assimilation.in file and include it when executing the program. Is that true?
If the sst and t(3D) used in ...

zhaoqian

Hi all:
In trunk/ROMS/Nonlinear,
# if defined NUDGING_SST || defined NUDGING_T
!
! Assimilate tracer observations via nudging.
!
IF (update_T(itrc,ng)) THEN
DO k=1,N(ng)
DO j=JstrR,JendR
DO i=IstrR,IendR
cff=MIN(1.0_r8,MAX(0.0_r8,EobsT(i,j,k,itrc)))
cff=dt(ng)*Tnudass(itrc,ng)*(1.0_r8-cff ...

zhaoqian

Hi all:
I want to run the couple ROMS and SWAN, and I run the INLET_TEST case at first, but I encounter some problem when executing oceanM. The following is the related informations with the errors marked in red, would you please help me resolve it?

qzhao@altix450:~/trunk> /opt/mpich2-ifort/bin ...

zhaoqian

Hi liangliang:
Is the variable Tcline? In ocean_test_head.in, it is 0.1 while in *.nc file it is very large.

zhaoqian

Dear arrango:
Thank you for your work. I have updated the cppdefs.h and globaldefs.h files.

zhaoqian

Hi all:
If the NUDGING_SSH flag is actived, where is the nudging of ssh performed in ROMS? I didn't find the code related.
I think ROMS only has the interface of nudging ssh, and we can't use this function directly. However, we can nudge M2,M3 and T directly. Am I right?
Thank for replying.

zhaoqian

Dear all:
Can roms model the Sediment-bound Tracer Transport, which is analogous to a hydrophobic contaminant, e.g., organic chemical, heavy metal or radionuclide,that adsorbs to fine-grained sediment particles? The contaminants mentioned above are partially soluble and exist in both particulate and ...

zhaoqian

Hi all:
Can Roms be used to trace the float trajectory back in the past few days?
If can, which program?
Thanks for your replying!

zhaoqian

zhaoqian wrote:
In set_data.F, the transfromation about winds are listed above. According to the code, the values of sustr are positive if the wind is from north. Is that true?

If the wind is from north, the svstr is negitive!

zhaoqian

[img]d:\result.jpg[/img]

zhaoqian

Hi all:
When I use Q_PSOURCE modeling transport and diffusion of pollutant discharging into the sea, and set Qbar and Qsrc positive, the concentration of the pollutant is negative, and when I set Qbar and Qsrc negative, the concentration becomes posivie. Why? Can anybody explain this for me ...

zhaoqian

I think you can active Q_PSOURCE and TS_PSOURCE and set Qbar、Qsrc and Tsrc in ana_psource.h file.

The ana_psource.h file can not only sets the imformation for T&S, but can set the imformation for other tracers.

zhaoqian

! If input point wind stress, rotate to curvilinear grid. Notice
! that rotation is done at RHO-points. It does not matter.
!
IF (.not.Linfo(1,idUsms,ng).or. &
& (Iinfo(5,idUsms,ng).ne.Lm(ng)+1).or. &
& (Iinfo(6,idUsms,ng).ne.Mm(ng)+2)) THEN
DO j=JstrR,JendR
DO i=IstrR,IendR
cff1=FORCES(ng ...

zhaoqian

If you are giving ROMS a global wind field and expecting ROMS to do the interpolation, ROMS will also do the rotation. Same for a single value of wind data that gets used over the whole domain.
Dear Kate and other users:
Would you please tell where to find the interpolation and rotation part in ...

zhaoqian

Hi all users:
We notice that in orthogonal curvilinear grids the velocity in nc file are in eta and xi directions. So if we want to plot current fields using the matlab command quiver, should we change the velocity fields into the cartesian coordinate first? And if the external forcing like winds ...

zhaoqian

Dear Kate and all:
I notice in your explaination and the cppdefs.h that there have three options (TS_PSOURCE, UV_PSOURCE, Q_PSOURCE) about sources.
I know that TS_PSOURCE is giving the value of T/S or the concentration of tracers, and notice that TS_PSOURCE and Q_PSOURCE are used in step3d_t.F ...

zhaoqian

Hi all:
I want to model pollutants diffusion in a bay with one point source discharging about 1000 tons pollutants in one year into the sea. The initial concentration of the tracer is 1mg/L and the boundary concentration is set to be a constant(for example 1mg/L). How should I set the model? The ...

zhaoqian

You are hoping to replace the topmost values of T & S directly into the model?
Dear Kate and all:
Thanks for replying, and I know how the upper boundary conditions of T &S set in ROMS .
Yes! I want to use the use the values of T & S directly instead of using flux because of its bad quality in ...

zhaoqian

Dear all roms user:
I found in set_vbc.F that the flux condition(plus sst or sss for the purpose of flux correction sometimes) are used as the sea surface condition.
I want to konw that can SST and SSS be used directly as the sea surface condition, which is the same as POM. If can, where is the ...

zhaoqian

Dear arango and other roms users:

I also meet the problem metioned above when calling boundary file using seagrid. The following is the matlab executing information.

>> seagrid

thePoints =

117.5000 41.0000 1.0000
117.5000 36.5000 1.0000
127.0000 36.5000 1.0000
127.0000 41.0000 1.0000 ...

zhaoqian

I have solved this problem by defining INLINE_2DIO in the *.h file.

zhaoqian

Hi all:
I run mpi roms in my SGI altix 450 computer. But the program stopped when finishing initialization and starting time-stepping. The following is the details of the screen output.It seems that the program stopped when creating output file. Would you please give me some instructions? By the ...

zhaoqian

Dear kate and all:
Thank you for your instruction! I have changed my cpp definitions and the model can run for nearly one day, and blow up.
The following is the cpp flags of my realistic application: tidal driven,southern open boudary, t and s fixed. I am not sure if they are right or suitable ...

zhaoqian

I have found it in globaldefs.h.
/*
** Check if processing timeless data.
*/

#if (!defined ANA_PSOURCE && \
(defined UV_PSOURCE || defined TS_PSOURCE || \
defined Q_PSOURCE)) || \
(defined SSH_TIDES || defined UV_TIDES)
# define TIMELESS_DATA
#endif

zhaoqian

In \Nonliear\initial.F, it write
! If applicable, read in input data.
!
# ifdef TIMELESS_DATA
!
CALL get_idata (ng)
# endif
CALL get_data (ng)

but I can't understand the meaning of TIMELESS_DATA and can't find where(in which file) it is difined.

Thanks a lot for your instruction.

zhaoqian

You can follow the instructions of the following page step by stey.
https://www.myroms.org/wiki/index.php/Tidal_Forcing

zhaoqian

Thank you very much! I have solve the problem following your advice.

zhaoqian

I set ANA_INITIAL flags and give the t & s of the computing sea area analytically. But the program stop when executing the get_data.f90 which is in initial.f90. I read the get_data.f90 and find that it reads the bry file.

My question is: if i have the tidal constituents and want to use them as the ...

zhaoqian

If simulating only tide (t/s fixed )in some sea area, how many input files should I supply? I only supply the ocean_grd.nc and ocean_frc.nc(contain tidal imformation). The initial condition and wind forcing is set analytically.
Am i right? Plsease give me some suggestion! Thanks for help.
The ...

zhaoqian

Thanks for your reply!
Do you mean that the depth should be nonzero in the water points and zero in the land points?
I find that in some water points,the depth is zero. This may be the reason why the error metioned occur. In order to avoiding this problem, i set the minimum depth to 5m when using ...

zhaoqian

When executing oceanM, some error occur and the program stop. The details are as follows:
--------------------------------------------------------------------------
Minimum X-grid spacing, DXmin = 1.39883582E+00 km
Maximum X-grid spacing, DXmax = 1.48958137E+00 km
Minimum Y-grid spacing, DYmin ...

zhaoqian

When I use the matlab code "otps2frc_v3.m" to generate the nc file of tidal forcing, some mistake occurs. I am so confused :( .Please help me!
The OS of my computer is Suselinux 10.1, the matlab R14(sp2) and the netcdf for it are installed ...

Search found 45 matches