Search found 60 matches

by simion1232006
Thu Feb 04, 2016 8:55 am
Forum: ROMS Installation
Topic: build fail
Replies: 1
Views: 3239

build fail

Hi there, I am trying to build a ROMS_SWAN coupling model and got some error messages.Did any one happen to see this before? I used to compile ROMS-swan on other machine and the model worked all right. The machine I am using is a Cray system. I guess it may be duo to the Cray compiler. Any suggestio...
by simion1232006
Thu Aug 20, 2015 2:47 am
Forum: ROMS Bugs
Topic: contact.m issue with fine grid on border of coarse grid
Replies: 2
Views: 4213

Re: contact.m issue with fine grid on border of coarse grid

I assume this is not a bug but rather a limited functionality. As long as your nested grid would not cross the border or within few grid of the coarse model border, there should not be any issue. When doing interpolation from donor grid to receiver grid, I think you'd better allow couple of grids of...
by simion1232006
Tue Jun 30, 2015 8:14 am
Forum: ROMS Adjoint
Topic: why the ajoint model need climatology input
Replies: 0
Views: 2574

why the ajoint model need climatology input

Hi everyOne, I am testing the ROMS strong constrain data assimilation model. I got a interesting error that asking me to prepare climatology inputs . This is only required for the updated ROMS version. For the version older than ROMS3.4 there is no such error. I am curious why it requires climatolog...
by simion1232006
Fri Apr 24, 2015 2:14 am
Forum: ROMS Bugs
Topic: A bug in nesting.F (2 way nesting with WET_DRY)
Replies: 2
Views: 3120

Re: A bug in nesting.F (2 way nesting with WET_DRY)

Good catch. I saw similar problems and did similar modifications. Still better let ROMS group to update an official version based on what we report.
by simion1232006
Tue Apr 21, 2015 7:31 am
Forum: ROMS Bugs
Topic: FINE2COARSE generate zero termperature value
Replies: 0
Views: 2372

FINE2COARSE generate zero termperature value

For two way nesting model with wet-Dry on, if the fine grid treat a area as ‘dry’, it will pass a zero value to the coarse grid, even this area is considered as 'wet' in coarse grid. This somehow did not affect too much the coarse model simulation, but will leave some unpleasant dots on temperature ...
by simion1232006
Mon Mar 09, 2015 6:55 am
Forum: ROMS Bugs
Topic: Possible bug in set_tides.F
Replies: 3
Views: 3482

Re: Possible bug in set_tides.F

The code should be correct. If my understanding is right, the phase angle should not have an a factor of period in it. Just make sure you have the correct phase angle input, and I think you will be fine.
by simion1232006
Thu Feb 12, 2015 2:05 am
Forum: ROMS Bugs
Topic: Bug in wet_dry- the lake_jersey case has the same issue
Replies: 8
Views: 5573

Re: Bug in wet_dry- the lake_jersey case has the same issue

Thanks a lot for your reply,Arango. I apologize for directly jumping into contact.m(Although there are couple of none essential typos mentioned by rock, the code works fine). The spherical grid also works great with the contact.m. I did not know the 659 ticket, but I do constantly update my ROMS at ...
by simion1232006
Wed Feb 11, 2015 9:47 am
Forum: ROMS Bugs
Topic: Bug in wet_dry- the lake_jersey case has the same issue
Replies: 8
Views: 5573

Re: Bug in wet_dry- the lake_jersey case has the same issue

Hi Warner, Thanks for your reply. Mine is svn #756, and I could not get the #659 you mentioned.I am quite curious of the ROMS version you use. I know you have developed your own branch of ROMS (forgot the name, sorry). Is that all right if I can use svn to download your version? I just looked into t...
by simion1232006
Wed Feb 11, 2015 1:19 am
Forum: ROMS Bugs
Topic: Bug in wet_dry- the lake_jersey case has the same issue
Replies: 8
Views: 5573

Re: There might be a bug in contact.m

I retested the lake_jersey case. With wet_dry on, the test case have the similar problem. I guess now I can have a rest (too tired). It confirms that the wet_dry option is still buggy in the newest ROMS. See below: lake_jersey_wetdry.jpg Anyone can try this out by switch on the wet_dry in lake_jerse...
by simion1232006
Tue Feb 10, 2015 1:26 pm
Forum: ROMS Bugs
Topic: Bug in wet_dry- the lake_jersey case has the same issue
Replies: 8
Views: 5573

Re: There might be a bug in contact.m

Thank you. That is helpful. I did a couple of tests, and the run without wet_dry looks fine, the wired results only appear with wet_dry on. I will keep updating my findings. Below is another nesting domain - twoway nesting, wet_dry. Notably there is no dry land in the domain. two_way_wetdry.jpg I as...
by simion1232006
Tue Feb 10, 2015 3:06 am
Forum: ROMS Bugs
Topic: Bug in wet_dry- the lake_jersey case has the same issue
Replies: 8
Views: 5573

Bug in wet_dry- the lake_jersey case has the same issue

Dear ROMS users, I had a problem with the information exchange between donor and receiver grid months ago, while after the update at the end of 2014, the problem disappear. But now it comes back. I hope someone can give me a hint. Simulation with one_way nesting. Donor_grid.jpg receiver_grid.jpg The...
by simion1232006
Thu Dec 18, 2014 1:56 am
Forum: ROMS Bugs
Topic: nesting_bug
Replies: 3
Views: 3596

Re: nesting_bug

Oh, I did not notice I need to regenerate the refined grid. That could be the problems. I actually modified the refined the batyhymetry to improve the depth, and matched them at the coarse fine boundary. I will need to recreate this grid and see how it goes. Many thanks.
by simion1232006
Wed Dec 17, 2014 7:59 am
Forum: ROMS Bugs
Topic: nesting_bug
Replies: 3
Views: 3596

Re: nesting_bug

This is my previous model output. It looks all fine. Note that I have also updated the matlab repository, created an new contact input file.
by simion1232006
Wed Dec 17, 2014 7:55 am
Forum: ROMS Bugs
Topic: nesting_bug
Replies: 3
Views: 3596

nesting_bug

I have no problem to use the previous ROMS version(nesting, no wet_dry), while after I updated the code yesterday,I could not run the model any more. It blows up after ~50 time step. see the attachment (SST of restart file): the SST changes a lot on the coarse fine boundary. I guess there is a bug i...
by simion1232006
Tue Dec 09, 2014 8:42 am
Forum: ROMS Bugs
Topic: WET_DRY and NESTING (refined) possible bug
Replies: 6
Views: 5053

Re: WET_DRY and NESTING (refined) possible bug

HI Prof arango. This is very good news to me. Thanks a lot for you hard work on ROMS.
by simion1232006
Tue Dec 09, 2014 2:28 am
Forum: ROMS Bugs
Topic: WET_DRY and NESTING (refined) possible bug
Replies: 6
Views: 5053

Re: WET_DRY and NESTING (refined) possible bug

I am facing the similar issue. I attempted to modified the nesting.F subroutine, while it is quite hard. I noticed that Lweight is directly taken from the contact file (created by matlab),and is constantly used to interpolate value from donor to receiver or from receiver to donor. Given that contact...
by simion1232006
Tue Nov 25, 2014 8:35 am
Forum: ROMS Discussion
Topic: Surface heat retention, year-round fake "upwelling"
Replies: 21
Views: 20560

Re: Surface heat retention, year-round fake "upwelling"

I noticed that the ERA-interim data set has a little cooling problem near the coast. It appears that the land side gets less heat than the ocean, while the land-ocean boundary is too coarse for some coastal circulation model. Thus some areas of ocean are gridded as land area receiving less heat.
by simion1232006
Tue Nov 25, 2014 8:24 am
Forum: ROMS Discussion
Topic: Surface heat retention, year-round fake "upwelling"
Replies: 21
Views: 20560

Re: Surface heat retention, year-round fake "upwelling"

That is just the description, the actual name are 'eewss' and 'nsss'.
by simion1232006
Sat Sep 20, 2014 6:34 am
Forum: ROMS Discussion
Topic: the ssh is abnormal
Replies: 5
Views: 2887

Re: the ssh is abnormal

It is very hard to understand what you mean. Although SSH is more relevant to your netcdf driving file than to .in file, you might want to check your driving force.
by simion1232006
Tue Aug 19, 2014 8:19 am
Forum: ROMS Bugs
Topic: bug in obs_roms2roms.m
Replies: 0
Views: 1854

bug in obs_roms2roms.m

I found a bug in obs_roms2roms.m. see code below % Receiver grid. if (isvec && D.vector_rotation ), if ( D.spherical ), Xname = 'lon_rho'; % If requested, interpolate Yname = 'lat_rho'; % U- and V-points variables else % to RHO-points instead to Xname = 'x_rho'; % facilitate the curvilinear ...
by simion1232006
Tue Jul 08, 2014 9:53 am
Forum: ROMS Discussion
Topic: Missing PAR variable in the roms_metadata.m
Replies: 1
Views: 1711

Re: Missing PAR variable in the roms_metadata.m

As far as I know, roms_metadata doesn't include that variable. You need to make a little modification for your own purpose.
by simion1232006
Thu Jun 26, 2014 1:46 am
Forum: ROMS Discussion
Topic: Problem about yellow_sea case
Replies: 12
Views: 5574

Re: Problem about yellow_sea case

Best choice would be using the ROMS official input file generating scripts (Matlab). You can modified it for your own purpose. Those scripts guarantee an standard input and usually would not make any mistake. I personally don't recommend change anything in varinfo.dat, as it might cause confusion la...
by simion1232006
Fri Jun 20, 2014 4:15 am
Forum: ROMS Discussion
Topic: Problem about yellow_sea case
Replies: 12
Views: 5574

Re: Problem about yellow_sea case

#define ANA_SMFLUX and #undef BULK_FLUX. Otherwise it needs to read in the heat flux and wind forcing.
by simion1232006
Fri Jun 13, 2014 2:04 am
Forum: ROMS Discussion
Topic: nesting blow up
Replies: 6
Views: 4287

Re: nesting blow up

On behalf of ROMS users, I need to say my poster here might be misleading. Prof Wilkin sent me a message about the potential volume conservative risk of one way nesting. So it might not be the best choice. The nesting scheme is still under active development, so will wait for the newest progress.
by simion1232006
Thu Jun 12, 2014 8:23 am
Forum: ROMS Discussion
Topic: relative humidity
Replies: 0
Views: 1448

relative humidity

I know this has been posted before. While I still have a little confusion about the humidity used in bulk flux script. IF (RH.lt.2.0_r8) THEN !RH fraction cff=cff*RH !Vapor pres (mb) Q(i)=0.62197_r8*(cff/(PairM-0.378_r8*cff)) !Spec hum (kg/kg) ELSE !RH input was actually specific humidity in g/kg Q(...
by simion1232006
Thu Jun 12, 2014 4:13 am
Forum: ROMS Discussion
Topic: nesting blow up
Replies: 6
Views: 4287

Re: nesting blow up

I have checked everything. Turning on the one_way nesting is so far the only solution.
by simion1232006
Wed Jun 11, 2014 3:18 pm
Forum: ROMS Discussion
Topic: nesting blow up
Replies: 6
Views: 4287

Re: nesting blow up

For the interests of anyone who working on ROMS nesting scheme, I feel the one_way nesting is more stable than the two_way nesting. I will dig into the codes to see their differences . while if someone know the reason, I really appreciate you can post them here.
by simion1232006
Tue Jun 10, 2014 2:22 pm
Forum: ROMS Discussion
Topic: nesting blow up
Replies: 6
Views: 4287

Re: nesting blow up

I just found the radiation_2D flag should be switched on. Another test with this flag is running fine for 3000 time steps already. It could be the problem.
by simion1232006
Tue Jun 10, 2014 2:09 pm
Forum: ROMS Discussion
Topic: nesting blow up
Replies: 6
Views: 4287

nesting blow up

I have experienced blowing up problem with ROMS nesting scheme for my application. Since the real case is complicated and hard to tell , I cleaned up all unnecessary issues by turning off the surface heat flux and wind, making the refined grid exactly the same as the master grid (the same masking an...
by simion1232006
Sun Mar 16, 2014 5:11 am
Forum: ROMS Adjoint
Topic: zero value in tangent linear initial file
Replies: 2
Views: 3546

Re: zero value in tangent linear initial file

I ran the WC13 case (the strong constrained case) and got the similar log information. In the next loop, the model updated the itl file, and the value are not zero any more. In the weak constrained case, everything is normal.

There should be nothing wrong with it. Never mind.
by simion1232006
Sat Mar 15, 2014 2:43 pm
Forum: ROMS Adjoint
Topic: zero value in tangent linear initial file
Replies: 2
Views: 3546

zero value in tangent linear initial file

I cut part of the log file, see below: The nonlinear model runs all right, while when the model switch to the tangent model, it creates a initial file with all zero values (NWS_200911_itl.nc). Does anyone have any experience on this? I am quite new in regional model data assimilation. Many thanks. ...
by simion1232006
Tue Feb 18, 2014 9:04 am
Forum: ROMS Discussion
Topic: online nesting test problem
Replies: 2
Views: 2373

Re: online nesting test problem

Problem solved. I should not define the uv_logdrag and MB_BBL at the same time.
by simion1232006
Tue Feb 18, 2014 1:41 am
Forum: ROMS Discussion
Topic: online nesting test problem
Replies: 2
Views: 2373

online nesting test problem

I am working on testing the newest nesting scheme for my application. It seems that the model start to ran, while stop at the initialization and configuration step. Does any of you have experiences on that? Also I found a problem to run the nesting model:The tiles specified for the coarse grid shoul...
by simion1232006
Tue May 14, 2013 2:48 am
Forum: User Applications
Topic: LTRANS: Larval TRANSport Lagrangian model
Replies: 40
Views: 64872

Re: LTRANS: Larval TRANSport Lagrangian model

Thanks, Kate. I will have a look. There is another issue coming with the code. It is not a bug but very annoying(see the figure). somehow the particle jump over something and the model stop to run. I mean is there a way to turn this thing off or just let the particle stop if it passed through certai...
by simion1232006
Mon May 13, 2013 5:10 am
Forum: User Applications
Topic: LTRANS: Larval TRANSport Lagrangian model
Replies: 40
Views: 64872

Re: LTRANS: Larval TRANSport Lagrangian model

Thanks for your reply. The good side is it actually solved my problem. The trajectories of these particles look pretty smooth and realistic. So without further debugging of the codes, I would just use it as a temporary solution.
by simion1232006
Mon May 13, 2013 3:35 am
Forum: User Applications
Topic: LTRANS: Larval TRANSport Lagrangian model
Replies: 40
Views: 64872

Re: LTRANS: Larval TRANSport Lagrangian model

I don't know exactly what happens here. Seemingly, adding commands to print out variables on the screen can totally change the model results. The original code, in my experience, gave me unrealistic vertical motion of particles. If somebody has the same issue, I hope my post would help.
by simion1232006
Mon May 13, 2013 3:29 am
Forum: User Applications
Topic: LTRANS: Larval TRANSport Lagrangian model
Replies: 40
Views: 64872

Re: LTRANS: Larval TRANSport Lagrangian model

I added the test command because I found something wrong in the outputs. The interesting thing is I solved the problem by adding another command (see attachments).
by simion1232006
Mon May 13, 2013 1:13 am
Forum: User Applications
Topic: LTRANS: Larval TRANSport Lagrangian model
Replies: 40
Views: 64872

Re: LTRANS: Larval TRANSport Lagrangian model

Dear all, I tried to use the LTRANs version 2 for my ROMS application. It works great if I release the particles on the surface. However, I found something weird of the model outputs if I released the particle in the subsurface layer. In the attached figure, the depth of the first particle moves imm...
by simion1232006
Tue Feb 07, 2012 6:18 am
Forum: ROMS Problems
Topic: Compilation Problem
Replies: 21
Views: 11920

Re: Compilation Problem

I had the similar problems. I spent a whole day but still could not figure out the problems, until I saw Hadfield's reply.
Coping the uv_rotate.f to Utility folder really works. Maybe the ROMS team could do this tiny change.
by simion1232006
Tue Oct 18, 2011 7:47 am
Forum: ROMS Usage
Topic: periodic boundary condition
Replies: 16
Views: 13890

Re: periodic boundary condition

Sincerely appreciate your favor. Yes, this is the right one. I will have a look. -------------------------------------------------------------------------------- I get a reply from Dr zhang. It works for me now.I can see the atmosphere pressure driving flow now even with periodic boundary condition....
by simion1232006
Tue Oct 18, 2011 4:07 am
Forum: ROMS Usage
Topic: periodic boundary condition
Replies: 16
Views: 13890

Re: periodic boundary condition

Thank you.Wilkin. Also did you post the right link? I did not see any things related to periodic boundary with along-shelf pressure gradient in this paper. Also it is in Hawaiian , not Mid Atlantic Bight. I really appreciate if you can send me a good link. Currently , I am working on the instability...
by simion1232006
Mon Oct 10, 2011 6:55 am
Forum: ROMS Discussion
Topic: outflow at the open boundary condition
Replies: 6
Views: 4991

Re: outflow at the open boundary condition

First thing is why you want to nudge your baroclinic velocity to zero. I think radiation for u , v should be enough. Also, i had the similar problems when I ran a model with open boundary with no boundary value provided(you put it as zero).The flow can hardly cross the boundary when i use flather fo...
by simion1232006
Sun Oct 09, 2011 6:23 am
Forum: ROMS Usage
Topic: periodic boundary condition
Replies: 16
Views: 13890

Re: periodic boundary condition

Thank you very much. Kate.
I will reconsider it.
by simion1232006
Fri Oct 07, 2011 8:14 am
Forum: ROMS Discussion
Topic: Can we have only passive tracers in rivers/point sources?
Replies: 40
Views: 37273

Re: Can we have only passive tracers in rivers/point sources

I used MIX_GEO_TS and have not found any problems yet. Maybe it is case by case.
by simion1232006
Fri Oct 07, 2011 7:33 am
Forum: ROMS Usage
Topic: periodic boundary condition
Replies: 16
Views: 13890

Re: periodic boundary condition

Oh ,this is tricky. I simply want a channel flow driven by both pressure gradient and wind. For simplicity , it could be infinite in long-shore direction. However , seems it is very hard to put into practice
by simion1232006
Thu Oct 06, 2011 8:03 am
Forum: ROMS Usage
Topic: periodic boundary condition
Replies: 16
Views: 13890

Re: periodic boundary condition

Periodic boundary conditions are a set of boundary conditions that are often used to simulate a large system by modelling a small part that is far from its edge. Basically on these boundaries the dynamic properties are exactly the same. It is suitable for simulating the system that is infinite in on...
by simion1232006
Thu Oct 06, 2011 2:56 am
Forum: ROMS Usage
Topic: periodic boundary condition
Replies: 16
Views: 13890

Re: periodic boundary condition

Thank you. Kate. I think your methods will work. (flather with specified boundary at north, radiation at south, solid both east and west). Also it is a problem to decide the ubar corresponding to certain pressure gradients. Actually my case is more about pressure gradient via sloping isopycnals,so I...
by simion1232006
Wed Oct 05, 2011 2:02 am
Forum: ROMS Usage
Topic: periodic boundary condition
Replies: 16
Views: 13890

Re: periodic boundary condition

Hi Prof Arango, Thank you very much. However , is there a way to apply a negative pressure gradient (opposite to wind )in idealized up-welling cases. I did some hind cast simulations fully driven by global model,but it is too complicated to analyse. I would like to simplify it by using a idealized m...
by simion1232006
Tue Oct 04, 2011 7:57 am
Forum: ROMS Usage
Topic: periodic boundary condition
Replies: 16
Views: 13890

periodic boundary condition

Hi everyone. I have some questions related to the periodic boundary conditions. Does anyone have the experience in specifying a periodic boundary condition in idealized upwelling case with a pressure gradient (barotropic or baroclinic). So basically only the velocity fields are periodic, and both th...
by simion1232006
Wed Jun 01, 2011 2:56 am
Forum: ROMS Problems
Topic: mpirun problem
Replies: 7
Views: 6394

Re: mpirun problem

Try mpiexec command see whether it works. As if you use shared memory machine and use the gforker manager, the mpirun command might not work.
by simion1232006
Tue May 31, 2011 2:02 am
Forum: ROMS Problems
Topic: mpirun problem
Replies: 7
Views: 6394

Re: mpirun problem

I used to have the similar problem. You might not install the openmpi in the right way. I use mpich2, so not sure whether it is the same.
by simion1232006
Fri May 27, 2011 2:29 am
Forum: ROMS Installation
Topic: Error installing mpi on cygwin
Replies: 8
Views: 6756

Re: Error installing mpi on cygwin

You need install gfortran gcc g++ first in ubuntu using synaptic package manager(fantastic feature in ubuntu). Then you probably need to install netcdf and mpich2 library. For this part you really need to figure it out by yourself as it is worth a exploration. Be careful about the flag as all the li...
by simion1232006
Thu May 26, 2011 1:11 am
Forum: ROMS Installation
Topic: Error installing mpi on cygwin
Replies: 8
Views: 6756

Re: Error installing mpi on cygwin

Download it from the official website ,then use virtualbox to install it.You do not need a real CD, a ISO file is fine.
See this: http://www.ubuntu.com/download. Also Linux system is free of charge. 8)
by simion1232006
Wed May 25, 2011 2:51 am
Forum: ROMS Installation
Topic: Error installing mpi on cygwin
Replies: 8
Views: 6756

Re: Error installing mpi on cygwin

I agree with that. As long as you do not bother the multi processors, cygwin is enough to use. Virtualbox things is actually running a real linux system within window . So if you have a workstation with more than 8 core. you definitely can run the real simulation. Another thing: I am ubuntu fan~~~ :...
by simion1232006
Wed May 25, 2011 1:38 am
Forum: ROMS Installation
Topic: Error installing mpi on cygwin
Replies: 8
Views: 6756

Re: Error installing mpi on cygwin

I tried to use cygwin long time ago and found it is really tricky and also not efficient to run. I have a suggestion you can install a software called virtulbox, and install a ubuntu OS in the virtualbox. Then you can run this linux system within the virtualbox under you host computer. All the other...
by simion1232006
Wed May 25, 2011 1:25 am
Forum: ROMS Discussion
Topic: bottom vertical velocity
Replies: 2
Views: 2956

Re: bottom vertical velocity

Thank you.Dr Warner. My bathymetry is not flat.That should be the reason. Can i ask another question: If I calculate the heaflux from ROMS history file, how could I make sure the hearflux is conservative. For each grid, i can compute the hearflux to the x-direvtion, y-direction,but for vertical dire...
by simion1232006
Tue May 24, 2011 2:41 am
Forum: ROMS Discussion
Topic: bottom vertical velocity
Replies: 2
Views: 2956

bottom vertical velocity

Hi everyone: Is there anyone notice in ROMS output the vertical velocity at the bottom is not zero. In my study, I tried to look at the tracer transports. I integrated the vertical heat flux through the depth and found it is not zero, so I checked the vertical velocity from the history file, and fou...
by simion1232006
Fri Oct 29, 2010 2:18 am
Forum: ROMS Adjoint
Topic: observation data rejected!
Replies: 9
Views: 10981

observation data rejected!

Number of State Observations Processed: ObsTime = 40124.0000 Variable IstrObs IendObs Count Rejected temp 1 33 33 32 Total 33 32 Obs Tally 33 32 Wrote NLM state at observation locations, datum = 0000001 - 0000033 Almost all of my data were rejected.Does anyone have the similar problems?
by simion1232006
Fri Apr 30, 2010 3:29 am
Forum: ROMS Usage
Topic: nest
Replies: 0
Views: 1701

nest

I heard the nesting feature has been added into ROMs. I plan to use the one way nest ,but I can not find introduction about this. I want to nest a 250m resolution small model near the coast into a 2~4 km resolution big model. Please help ! Thanks ----It seems that I can not do nesting in the officia...