Search found 45 matches

by zhaoqian
Tue Jun 21, 2011 8:51 am
Forum: ROMS Discussion
Topic: nudge one kind of passive traces
Replies: 7
Views: 4018

Re: nudge one kind of passive traces

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:l...
by zhaoqian
Sun Jun 19, 2011 8:09 am
Forum: ROMS Discussion
Topic: nudge one kind of passive traces
Replies: 7
Views: 4018

Re: nudge one kind of passive traces

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,...
by zhaoqian
Fri Jun 17, 2011 1:04 pm
Forum: ROMS Discussion
Topic: nudge one kind of passive traces
Replies: 7
Views: 4018

Re: nudge one kind of passive traces

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, series' 'd...
by zhaoqian
Tue May 24, 2011 9:00 am
Forum: ROMS Discussion
Topic: nudge one kind of passive traces
Replies: 7
Views: 4018

nudge one kind of passive traces

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!
by zhaoqian
Sat Mar 05, 2011 3:28 am
Forum: ROMS Discussion
Topic: Can roms do the momentum and tracer modeling sepatately?
Replies: 3
Views: 2544

Re: Can roms do the momentum and tracer modeling sepatately?

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 trunk ...
by zhaoqian
Fri Mar 04, 2011 9:10 am
Forum: ROMS Discussion
Topic: Can roms do the momentum and tracer modeling sepatately?
Replies: 3
Views: 2544

Can roms do the momentum and tracer modeling sepatately?

Hi ,
I want to simulate the hydrodynamical fields first, and then use the output date to drive the tracer module.
Thank you!
by zhaoqian
Fri Sep 10, 2010 2:38 am
Forum: ROMS Sediment
Topic: question about transport equation of sediment
Replies: 4
Views: 4361

Re: question about transport equation of sediment

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 ?.
by zhaoqian
Fri Sep 10, 2010 1:34 am
Forum: ROMS Sediment
Topic: question about transport equation of sediment
Replies: 4
Views: 4361

Re: question about transport equation of sediment

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 ?
by zhaoqian
Fri Sep 10, 2010 1:00 am
Forum: ROMS Sediment
Topic: question about transport equation of sediment
Replies: 4
Views: 4361

question about transport equation of sediment

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 !
by zhaoqian
Tue Aug 31, 2010 3:21 am
Forum: ROMS Sediment
Topic: question about sediment.F
Replies: 2
Views: 4998

question about sediment.F

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. !---------------------------...
by zhaoqian
Tue Jun 15, 2010 1:41 pm
Forum: ROMS Problems
Topic: urgent problem on time ref and tide_start !!!!!!!
Replies: 11
Views: 7810

Re: urgent problem on time ref and tide_start !!!!!!!

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.
by zhaoqian
Fri Jun 04, 2010 9:10 am
Forum: ROMS Problems
Topic: How to create bry file for dye and sediment
Replies: 2
Views: 2049

Re: How to create bry file for dye and sediment

I have found them in varinfo.dat which are defined as dye_east ...,mud_east_ ...,sand_east_ ... .
by zhaoqian
Fri Jun 04, 2010 8:39 am
Forum: ROMS Problems
Topic: How to create bry file for dye and sediment
Replies: 2
Views: 2049

How to create bry file for dye and sediment

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!
by zhaoqian
Sun Mar 21, 2010 6:24 am
Forum: ROMS Problems
Topic: a question about NUDGING_SST || defined NUDGING_T
Replies: 3
Views: 3869

Re: a question about NUDGING_SST || defined NUDGING_T

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 nu...
by zhaoqian
Wed Mar 17, 2010 2:59 pm
Forum: ROMS Problems
Topic: a question about NUDGING_SST || defined NUDGING_T
Replies: 3
Views: 3869

a question about NUDGING_SST || defined NUDGING_T

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) t(i,j,k,n...
by zhaoqian
Mon Mar 15, 2010 11:28 am
Forum: ROMS Problems
Topic: problem when coupling roms and swan
Replies: 1
Views: 2774

problem when coupling roms and swan

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/mpir...
by zhaoqian
Thu Mar 11, 2010 7:56 am
Forum: ROMS Usage
Topic: About the TEST_HEAD case
Replies: 7
Views: 5494

Re: About the TEST_HEAD case

Hi liangliang:
Is the variable Tcline? In ocean_test_head.in, it is 0.1 while in *.nc file it is very large.
by zhaoqian
Fri Mar 05, 2010 6:32 am
Forum: ROMS Problems
Topic: a question about NUDGING_SSH and other NUDGING opptions
Replies: 2
Views: 2246

Re: a question about NUDGING_SSH and other NUDGING opptions

Dear arrango:
Thank you for your work. I have updated the cppdefs.h and globaldefs.h files.
by zhaoqian
Thu Mar 04, 2010 1:01 pm
Forum: ROMS Problems
Topic: a question about NUDGING_SSH and other NUDGING opptions
Replies: 2
Views: 2246

a question about NUDGING_SSH and other NUDGING opptions

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.
by zhaoqian
Wed Feb 24, 2010 10:45 am
Forum: ROMS Sediment
Topic: Can roms model the Sediment-bound Tracer Transport ?
Replies: 0
Views: 2138

Can roms model the Sediment-bound Tracer Transport ?

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...
by zhaoqian
Sun Jan 24, 2010 7:51 am
Forum: ROMS Discussion
Topic: trace the float trajectory back
Replies: 8
Views: 7887

trace the float trajectory back

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!
by zhaoqian
Wed Jan 20, 2010 4:32 pm
Forum: ROMS Discussion
Topic: about the coordinate conversion
Replies: 6
Views: 3503

Re: about the coordinate conversion

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!
by zhaoqian
Wed Jan 20, 2010 1:06 am
Forum: ROMS Discussion
Topic: how set a vertical discharge?
Replies: 8
Views: 9984

Re: how set a vertical discharge?

[img]d:\result.jpg[/img]
by zhaoqian
Mon Jan 18, 2010 9:24 am
Forum: ROMS Usage
Topic: How to model pollutant diffusion in a bay?
Replies: 3
Views: 3252

Re: How to model pollutant diffusion in a bay?

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? Thanks ...
by zhaoqian
Mon Jan 18, 2010 8:06 am
Forum: ROMS Discussion
Topic: how set a vertical discharge?
Replies: 8
Views: 9984

Re: how set a vertical discharge?

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.
by zhaoqian
Fri Jan 15, 2010 2:23 am
Forum: ROMS Discussion
Topic: about the coordinate conversion
Replies: 6
Views: 3503

Re: about the coordinate conversion

! 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=FOR...
by zhaoqian
Fri Jan 15, 2010 12:48 am
Forum: ROMS Discussion
Topic: about the coordinate conversion
Replies: 6
Views: 3503

Re: about the coordinate conversion

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 ROM...
by zhaoqian
Thu Jan 14, 2010 3:17 am
Forum: ROMS Discussion
Topic: about the coordinate conversion
Replies: 6
Views: 3503

about the coordinate conversion

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 ar...
by zhaoqian
Tue Jan 12, 2010 1:42 am
Forum: ROMS Usage
Topic: How to model pollutant diffusion in a bay?
Replies: 3
Views: 3252

Re: How to model pollutant diffusion in a bay?

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 while...
by zhaoqian
Mon Jan 11, 2010 4:15 pm
Forum: ROMS Usage
Topic: How to model pollutant diffusion in a bay?
Replies: 3
Views: 3252

How to model pollutant diffusion in a bay?

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 fol...
by zhaoqian
Sat Jan 09, 2010 1:36 am
Forum: ROMS Usage
Topic: about sea surface condition of T&S
Replies: 7
Views: 4347

Re: about sea surface condition of T&S

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 qual...
by zhaoqian
Fri Jan 08, 2010 5:18 pm
Forum: ROMS Usage
Topic: about sea surface condition of T&S
Replies: 7
Views: 4347

about sea surface condition of T&S

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 code...
by zhaoqian
Tue Dec 29, 2009 7:12 am
Forum: ROMS Information
Topic: Extracting Bathymetry and Coastline Data for SeaGrid
Replies: 16
Views: 30326

Re: Extracting Bathymetry and Coastline Data for SeaGrid

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 theOrthogonalit...
by zhaoqian
Tue Dec 22, 2009 3:13 pm
Forum: ROMS Usage
Topic: run mpi version on SGI altix 450 computer
Replies: 1
Views: 1958

Re: run mpi version on SGI altix 450 computer

I have solved this problem by defining INLINE_2DIO in the *.h file.
by zhaoqian
Sun Dec 20, 2009 4:13 am
Forum: ROMS Usage
Topic: run mpi version on SGI altix 450 computer
Replies: 1
Views: 1958

run mpi version on SGI altix 450 computer

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 way...
by zhaoqian
Mon Nov 23, 2009 12:41 pm
Forum: ROMS Usage
Topic: how many input files should I supply when simulating tide
Replies: 5
Views: 4074

Re: how many input files should I supply when simulating tide

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. Woul...
by zhaoqian
Sun Nov 22, 2009 2:54 am
Forum: ROMS Problems
Topic: Where does TIMELESS_DATA be difined?
Replies: 1
Views: 1485

Re: Where does TIMELESS_DATA be difined?

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
by zhaoqian
Sat Nov 21, 2009 7:32 am
Forum: ROMS Problems
Topic: Where does TIMELESS_DATA be difined?
Replies: 1
Views: 1485

Where does TIMELESS_DATA be difined?

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.
by zhaoqian
Sat Nov 21, 2009 7:20 am
Forum: ROMS Problems
Topic: about forcings
Replies: 3
Views: 2455

Re: about forcings

You can follow the instructions of the following page step by stey.
https://www.myroms.org/wiki/index.php/Tidal_Forcing
by zhaoqian
Sat Nov 21, 2009 6:59 am
Forum: ROMS Usage
Topic: Why the DZmin and DZmax become NaN?
Replies: 4
Views: 3137

Re: Why the DZmin and DZmax become NaN?

Thank you very much! I have solve the problem following your advice.
by zhaoqian
Sat Nov 21, 2009 6:57 am
Forum: ROMS Usage
Topic: how many input files should I supply when simulating tide
Replies: 5
Views: 4074

Re: how many input files should I supply when simulating tide

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 t...
by zhaoqian
Thu Nov 19, 2009 7:35 am
Forum: ROMS Usage
Topic: how many input files should I supply when simulating tide
Replies: 5
Views: 4074

how many input files should I supply when simulating tide

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 follow...
by zhaoqian
Thu Nov 19, 2009 7:27 am
Forum: ROMS Usage
Topic: Why the DZmin and DZmax become NaN?
Replies: 4
Views: 3137

Re: Why the DZmin and DZmax become NaN?

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 se...
by zhaoqian
Fri Nov 13, 2009 7:22 am
Forum: ROMS Usage
Topic: Why the DZmin and DZmax become NaN?
Replies: 4
Views: 3137

Why the DZmin and DZmax become NaN?

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 = 1.7...
by zhaoqian
Wed Oct 28, 2009 2:52 pm
Forum: ROMS Information
Topic: Problem when dealing with the boundary file of tidal forcing
Replies: 0
Views: 2175

Problem when dealing with the boundary file of tidal forcing

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. -------------------------------------------------...