Ocean Modeling Discussion

mengqingjun

Hi ROMS's users! I want to compute the disspiation of turbulent kinetic energy using ROMS outputs.I used Mellor-Yamada turbulence close schemes. I have ROMS output tke(oceantime,s_w,eta_rho,xi_rho) and gls(oceantime,s_w,eta_rho,xi_rho) . Now I have a question that if q^2=sqrt(2*tke) and l = gls/tke ...

mengqingjun

ROMS zeta is "absolute SSH" in the sense that it includes time mean and time varying SSH with respect to a datum (i.e. zeta=0) that is a geopotential surface, i.e. if zeta=constant there is no associated pressure gradient force. Altimetrists do not say "absolute" SSH so let me clarify: In altimetry ...

mengqingjun

Respected ROMS Users, In some earlier posts I found that ROMS Zeta is absolute SSH which can be considered equivalent to satellite derived MDT + SLA , where MDT is mean dynamic topography and SLA is Sea level anomaly. In AVISO site I found MADT, which is map of absolute dynamic topography is again ...

mengqingjun

Ok,thank you wilkin for your replies! I am modeling the sea over shallow sea shelf.I didn't define BULK_FLUXS,I have very shallow water and therefore very thin surface cells.So I need to define SOLAR_SOURCE option as your advives, no matter whether I use MY25_MIXING nor other vertival mixing schemes ...

mengqingjun

If you #define SOLAR_SOURCE then the specified net shortwave radiation is propagated vertically into the water column and absorbed internally according to a double exponential profile set by the Jerlov water type option WTYPE in ocean.in. This step is handled by routine lmd_swfrac even if the KPP ...

mengqingjun

Thank you Professor Arango for your suggestions ! Firstly,about instability. physical instabilities (barotropic, baroclinic, mixed, etc) and they are usually associated with eddies .There are no obvious developing eddies in my domain,the truth is that there exist strange oscillations(the surface ...

mengqingjun

I have made a test forced only by the wind stress(3 hourly ecmwf) and output the results every 3 hour. ps. wind_stress_08th.gif is a gif animation of the wind stess in Aug; Total_Velocity_Surface_Aug.gif is a gif animation of the surface Total Velocity in Aug. The current looks like so unstable and ...

mengqingjun

Thank you LiangPeng for your reply! I did not test the oscillation because I think the oscillation is a little bit strange...
The inertial oscillation...um...The inertial oscillations fill all my domain? I have no idea...

mengqingjun

Hi redhat,I think if you are running a climatology case, you can use a cycle_length attribute, if you are running a real-time case, you can start your model from 2008-01-01_01:30:00. Whatever,try and see,good luck!

mengqingjun

Thank you Kate for reply ! The wind can excite inertial oscillations, yes, but I don't think the strange oscillations(the velocity is going round and round...) in my application are resonable. The velocity showed dramatic changes....After three years, I think the barotropic adjustments caused by the ...

mengqingjun

Recently I am similating a real-time ocean circulation case of East China Sea. I make initial condiations from HYCOM model on 2012.01.01 and the daily boundary condiations are also made from HYCOM model from 2012.01.01 to 2012.12.31. The forcing files are made from ECMWF datasets every 3 hours(using ...

mengqingjun

Hi d.kobashi~ :D I tried SOLAR_SOURCE and QCORRECTION to adjust heat flux (satellite-based SST). High SST/sub-surface temperature still persists. I am wondering how you gave the values of dQdSST to ROMS ? I meet a similar problem like yours. Without BULK_FLUXS,ROMS need only sstr,svstr,shflux,swflux ...

mengqingjun

Thank you kate for your replies ! I turn off the BULK_FLUXS, and as arango said, ROMS needs shflux, swflux, swrad (with some CPP options), sustr, and svstr as the atmospheric forcing. I am wondering we have already had shflux, then what swrad is used for ?

mengqingjun

I have a question that if I don't active BULK_FLUXS in .h file, as arango wrote in https://www.myroms.org/forum/viewtopic.php?f=30&t=3003&hilit=boundary+conditions that ROMS just needs shflux, swflux, swrad (with some CPP options), sustr, and svstr. .I want to apply the surface forcing without BULK ...

mengqingjun

When you say monthly, do you mean you have each individual month or do you have a generic January that is for all Januaries? If it's the first, you don't want to set cycle_length at all. If it's the second, I don't see that 365.25 is going to be so very terrible. It's not like you are picking up ...

mengqingjun

When you say monthly, do you mean you have each individual month or do you have a generic January that is for all Januaries? If it's the first, you don't want to set cycle_length at all. If it's the second, I don't see that 365.25 is going to be so very terrible. It's not like you are picking up ...

mengqingjun

Thank you Kate for reply ! in fact your nudging on outflow is stronger than my nudging on inflow, which could be trouble I see, I will try your's set. I'm not sure it's being used for anything since you don't have TCLM_NUDGING or similar turned on If I turn on TCLM_NUDGING or similar, ROMS will use ...

mengqingjun

Thank you rduran for reply! A monthly BC is not a good idea ,do you mean that a daily BC will do better for a ten year model ? Mybe I will try a daily BC... how come you have temperatures as low as 0.4 I think the temperature could reach as low as 0.4, may even be negtive one...I think it is ...

mengqingjun

Hi everyboday, I am running a real application without tide forcing on East China Sea. The model resultion is 1/12 degree. I run my model for 10 years in two cases: 1.The model is forced only by the boundary file *_bry.nc and it is interpolated from the HYCOM 1/12 degree results. The variables in ...

mengqingjun

Thank you Kate very much. It really helps

!

mengqingjun

If you use a Flather condition to specify ssh only with a propagating wave, you will only get half of the wave you expect. Using reduced physics might help this a bit, but you should do some experiments to be sure. I am running a real case which only has the tidal forcing with M2,S2,K1,O1;I provide ...

mengqingjun

Thank you John for your quick replies !
Except for SRELAXATION TCLM_NUDGING ,are there any other else CPP definations that will compile the file ana_nudgcoef.h?
Coulkd I think that even some files for examole ana_nudgcoef.h was compiled ,they may not be used when run the application ?
Thanks!

mengqingjun

Oh,you mean that the Tnudg value in ocean.in is first read by the file read_phypar.F and the Tnudg value is used the inversed value ! I see, thank you Kate! :D ps. Except for TCLIMATOLOGY TCLM_NUDGING that will compile ana_nudgcoef.h,what other else CPP definitations that will compile ana_nudgcoef.h ...

mengqingjun

ps. Are there some CPP definitations that could control wheather the file ana_nudgcoef.h will be compiled or not ?

mengqingjun

Sorry, I make a mistake ,I used LBC(): LBC(isFsur) == Clo Clo Clo Clo ! free-surface LBC(isUbar) == Clo Clo Clo Clo ! 2D U-momentum LBC(isVbar) == Clo Clo Clo Clo ! 2D V-momentum LBC(isUvel) == Clo Clo Clo Clo ! 3D U-momentum LBC(isVvel) == Clo Clo Clo Clo ! 3D V-momentum LBC(isMtke) == Clo Clo Clo ...

mengqingjun

ps. I mean that how the ana_nudgcoef.h was used ?

mengqingjun

It seems that the longer Nudging time scale (days) the weaker the nudging. But in ana_nudgcoef.h: IF (DOMAIN(ng)%SouthWest_Test(tile)) THEN Tobc_out(itrc,ng,iwest)=Tnudg(itrc,ng) Tobc_in (itrc,ng,iwest)=obcfac(ng)*Tnudg(itrc,ng) END IF As recommended in ocean.in, OBCFAC > 1 ; at the same time, Tnudg ...

mengqingjun

I am still confused about the nudging coefficient. In my application I set LBC as follows: ad_LBC(isFsur) == Clo Clo Clo Clo ! free-surface ad_LBC(isUbar) == Clo Clo Clo Clo ! 2D U-momentum ad_LBC(isVbar) == Clo Clo Clo Clo ! 2D U-momentum ad_LBC(isUvel) == Clo Clo Clo Clo ! 3D U-momentum ad_LBC ...

mengqingjun

Thank you Kate and I will have a try !

mengqingjun

Thank you for your replies Kate !Because I am a newer to ROMS, so I want get familiar with it through runing as mang test cases as possible. There are now global ocean models which we can use for RadNud boundary conditions and any of them should lead to better behavior Could I think that the *_clm ...

mengqingjun

I am testing the DAMEE_4 application in the test directory.There is no *_bry.nc inputing, and the Lateral Boundary Condition Type are all " CLO ". It is obvious that ROMS get the boundary condiations from the damee4_Lclm_b.nc. If so, I check the damee_4.h, I find that #define TCLIMATOLOGY and ...

mengqingjun

I am wandering how you set up your grid file?
Using Seagrid,or other else?
Could you kindly tell me ?
Thank you!

mengqingjun

I am preparing initial input files for ROMS using the matlab scripts d_initial.m . There is "OAname = fullfile(my_root, 'Data/OA', 'oa4_lev94_feb.nc');",and I don't know how to prepare my "oa_lev94_feb.nc".It seems like the file is produced in the old fortran code. I don't know how to produce my ...

mengqingjun

Thank you wilkin for your reply!
I am waiting for your update of the Latte input profiles, if you have arleady made it, just let me know!
Thank you!

mengqingjun

I have the same problem. I need hellp! Thanks! I had a problem when I am compiling my ROMS downloaded via the svn。I used Cygwin.When I input the commend:"./build.bash",it showed: "./buid.bash: line 397 cd: /home/mengqingjun/ocean/repository/trunk: No such file or directory meke:***没有规则可以创建目标“clean”.停 ...

Ocean Modeling Discussion

Search found 35 matches

About disspiation of turbulent kinetic energy

Re: 'zeta' in ROMS

Re: About ROMS Zeta

Re: Surface heat flux. Double the count?

Re: Surface heat flux. Double the count?

Re: Strange Oscillation

Re: Strange Oscillation

Re: Strange Oscillation

Re: forcing the ROMS with ERA interim

Re: Strange Oscillation

Strange Oscillation

Re: Surface heat flux. Double the count?

Re: about atmospheric forcing

about atmospheric forcing

cycle_length problem

Re: cycle_length for non-360days/yr monthly forcing files?

Re: Confusion about abnormal flow and great temperature grad

Re: Confusion about abnormal flow and great temperature grad

Confusion about abnormal flow and great temperature gradient

Re: do i have to add the UV_TIDES for a real case?

Re: do i have to add the UV_TIDES for a real case?

Re: Confusion about Nudging Coefficient

Re: Confusion about DAMEE4 test

Re: Confusion about Nudging Coefficient

Re: Confusion about Nudging Coefficient

Re: Confusion about Nudging Coefficient

Re: Confusion about DAMEE4 test

Re: Confusion about Nudging Coefficient

Re: Confusion about DAMEE4 test

Re: Confusion about DAMEE4 test

Confusion about DAMEE4 test

Re: Need help regarding free surface in South China Sea

d_initial.m problem

Re: Error running Latte example

Re: build.bash problem