Ocean Modeling Discussion

kobl1201

Hi when i work out ROMS model, I try to use 2 cpu, but the cpu continues to work with 1. When i set paralle node =1, the model work out well. which option can I check?? --------------------------------- It is my run.sh file #!/bin/bash #PBS -V #PBS -q batch #PBS -l nodes=2:ppn=1 #PBS -N ROMS ulimit ...

kobl1201

Thanks Wilkin

Is it possible to perform the ROMS model by making only the grid data into lambert?

Or is there an option within the ROMS model that needs to be turned on?

kobl1201

Hi

We want to combine the atmospheric model based on the lambert grid with the ROMS model.

Can a ROMS model perform with lambert grid?

If there is a way, can you tell me how to do it?

thank you~

kobl1201

Hello.

I have a question about the mixing scheme.

What is different characteristics between harmonic mixing scheme(VIS_UV) and smagorinsky shceme??

and can i use the smagorinsky scheme in ROMS model??

Thank you!

kobl1201

Hello Kate. I tried to change the shorter time for DT, but it is not useful as well. And i also changed the bathmetry, but not useful. I printed the values as follows. TIME-STEP YYYY-MM-DD hh:mm:ss.ss u2v2 POTEN_ENRG my_avgpe NET_VOLUME C => (i,j,k) Cu Cv Cw Max Speed 1242 2019-02-07 10:21:00.00 0 ...

kobl1201

Hello i have some problem, when i worked the ROMS model. I applicated a boundary file every 1 day, and atmospheric forcing every 1 hour. Starting time is 2019-02-07 00:00:00 ROMS model worked well until 2019-02-07 10:32:00, and then the model status had blow-up. When i checked the log file, avgpe ...

kobl1201

Hello i recently had an experince which the model stopped after first step. I could find NaN from the log file. The model always stops in a certain area. Therefore, i checked my input files such as initial, boundary, atmospheric forcing and so on. I couldn't find erro file or missing file. I can't ...

kobl1201

Thank you Kate~

so, the time scale is demermined by myself???

kobl1201

Hi~ when i simulated model, which resolution is about 3km, i set visc2 ~ 1m^2/s. i heard the rule of thumb approximation like this. visc2 ~ (2*delta_x / pi)**2 / timescale According to this formaula, when i set the 3km resolution and time scale is 100s, i should set the VISC2 ~ (2*3000/3.14)**2/100s ...

kobl1201

Hi Kate I change from UV_VIS2, TS_DIF2 to UV_VIS4, TS_DIF4 at the cpp option, respectively, and the value was changed from the VISC2=1m^2/s to VISC4=1m^4/s. The model finished and then i could get the results. But i wonder that the VISC2=1m^2/s corresponds to VISC4=1m^4/s or can i need to calculated ...

kobl1201

Hello I have some questions. Until now, i have run the ROMS model using the harmonic mixing coefficients, which are 1 m2/s for tracer and momentum. But i wanna change biharmonic mixing coefficient. Therefore, i just change the cpp opotion. EX) UV_VIS2 -> UV_VIS4 TS_DIF2 -> TS_DIF4 Can i change only ...

kobl1201

When i run the mdoel, i have some problem from output file.

However, i can't understand the problem.

Which one can i change to slove the problem.

Here i uploaded *.in

kobl1201

Thank you Kate~

Do you know the paper of Kaitlin Alexander or other people about the this problem??

kobl1201

Hi everyone~ I have some question about the result I worked two experiments which are the same condition except tile number. For example, at the Exp 1, i used NtileI == 8, Ntilej == 10 and at the Exp 2 i set Ntile == 10, Ntilej == 16. And then when i compared two results, I found the different ...

kobl1201

Dear All., My model is blowing and I don't know why. I tried to change DT(10sec -> 5 sec) but I saw same things. I put boundary condition every day, and atmospheric forcing applied every 6 hour. I check Input data, but it is fine. I think if the input data is wrong, the model would blow up first ...

kobl1201

Hi~ i want to run hot start. So firstly, I worked model for 5 days. start time is 55927. And then I needed to run 5 days more. So i want to do restart. However, i saw the error like this ; SET_NGFLD - current model time exceeds ending value for variable: zeta_west TDAYS = 55932.0000 Data Tmin ...

kobl1201

Thank you Kate.

I understand the nudginf time scale.

kobl1201

Thanks Kate.

At that time, If I have that the external data is interval of 1 day, H should be updated every 1 day or every 30 days??

Have a good time

kobl1201

Hello~

I applied relaxation method like this equation

T = T + dt/k(H-T)

Here, T is calculated by model, dt is time-step , k is nudging time scale, and H is external data.

If i set k = 30 day * 86400 sec, dt = 30sec, H is applied each time-step, or it is applied every 30 days.

Please Help me

kobl1201

Thanks Kate.

Here, inflow and outflow mean current in the boundary ???

Is it right??

kobl1201

Hi

OBCFAC is the radio of inflow and outflow nudging time scale.

Here, I wonder the meaning of inflow and outflow for nudgiinh time scale.

Inflow means data of climatology and outflow means the data calculated in model???

IF i set obcfac =5 to obcfac = 10, what does it mean??

kobl1201

Hello ~

i removed things such as !free-surface and then modified M2NUDG == 0.0d0 ! days, M3NUDG == 0.0d0 ! days. However, the results still showed same things.

Do you have other ideas??

Please help me

kobl1201

Hi Kate

I set up the model like this.

ocean_ub_3k_new.in: (117.57 KiB) Downloaded 454 times

I ,also applied relaxation method and sponge layer.

However, this methods didn't slove this problems

kobl1201

Sorry Kate.

I used HYCOM data for boundary condition and initial data.

The boundary condition applied every 1 day.

Usually, currents in the north region flow easternward. However, my result showed converse.

kobl1201

Hello~ i have some problems. When i run the ROMS model, the currets flow along the boundary. As time passed, the results varied as follows : 사본 -그림1.gif It seems the current can't go out the northern and eastern boundaris. I hope that outflow is much stronger than now. The model region is east/japan ...

kobl1201

Hi kate

when i turn off LuvSponge, the model run.

but if i turn on the LuvSponge, the model stop above.

the problem is also same reason, which you told me???

thank you~

kobl1201

Hi Kate~ i compile Ana_sponge. The model didn't give message. just stop the model like this. ** Omit *** ndtfast, nfast = 40 56 nfast/ndtfast = 1.40000 Centers of gravity and integrals (values must be 1, 1, approx 1/2, 1, 1): 1.000000000000 1.039265571007 0.519632785503 1.000000000000 1.000000000000 ...

kobl1201

Hello~

i want to apply sponge layer. Therefore, i change ana_sponge.h and then applied hear file.

and then i think i will turn on switch about sponge layer.

Which one should i turn on??

When i turned on the switch, which is LuvSponge and LtracerSponge, but the model didn't run.

Please help me~

kobl1201

Hi Kate~ i changed ININAME to restart file. however, i have some problems. the model didn't run. the model show like this. ======================================================================================= SET_NGFLD - current model time exceeds ending value for variable: zeta_west TDAYS = 55928 ...

kobl1201

when i set DSTART = 56109, the ocean_time in restart file is 4832179200. But i set boundary file from 55927 to 56292. if i devide 4832179200/86400, it is 55928. you mean that i just change ININAME to the name of the restart file, even though the ocean_time in restart file is 4832179200. Is it right??

kobl1201

Hello~ i want to know the method of restart i just know that LcycleRST set T. 1) How can I do if i do restart. Please give some information in detail. and also 2) when i run model. the start time is 56658 and then i run model during 1 day. but, the ocean_time in restart file is 484890400. I think it ...

kobl1201

Hi When I run the ROMS model, the model, blew up after running about 1 hours. The resolution is 1/72 and 41 vertical layers. Now i have simulated East sea area. However, when i set VISC2 = 5000, the model ran until the end. I thought the VISC2(=5000) is too big. It is not reasonable number. How can ...

kobl1201

hi i have some questions i saw ROMS code about sponge layer #if defined ADRIA02 ! ! Adriatic Sea southern sponge areas. ! fac=4.0_r8 # if defined UV_VIS2 DO i=IstrT,IendT DO j=JstrT,MIN(6,JendT) cff=visc2(ng)+REAL(6-j,r8)*(fac*visc2(ng)-visc2(ng))/6.0_r8 MIXING(ng) % visc2_r(i,j)=cff MIXING(ng ...

kobl1201

Hello I have some questions When i run ROMS model, the model stops without error. Here is my log file. --------------------------------------------------------------------------------- Model Input Parameters: ROMS/TOMS version 3.7 Thursday - July 28, 2016 - 9:26:36 PM ...

kobl1201

Hello guys i need your help when i built ROMS model, i saw the some error. but i could't understand why the error came out. ' ' >> ./build bash /usr/bin/cpp -P -traditional -DMPI -DLINUX -DX86_64 -DPGI -D'ROOT_DIR="/home/Class/ocean/kobl/ROMS"' -Dupwelling -D'HEADER="upwelling.h"' -D'ROMS_HEADER ...

kobl1201

Hello~~

i want to your help.

i want to make model grid file. However i can't run make_grid.m used in matlab.

when i ran make_grid.m, i saw 'undefined function romstool_param.m .

how can i run make_grid.m.

Please tell me some example.

kobl1201

Hello

what is effect for U and v 3-D velocity if i set Luvsponge and Ltracersponge are T or F ??

kobl1201

i have a question one more.

cff is nudging coefficient or not??

what is the nudging coeficient ??

kobl1201

hello when i run ROMS model using nuding scheme, i want to change nudging coefficient. At that time, i have a question. how much is the range of the coefficient?? now, the roms model was set like this DO j=JstrV,Jend DO i=Istr,Iend cff=0.25_r8*(CLIMA(ng)%M3nudgcof(i,j-1,k)+ & & CLIMA(ng)%M3nudgcof(i ...

kobl1201

when i checked zeta_west, this time range is 55927 to 56292.

But this model don't run.

How can i do?? please help me.

kobl1201

i set my model ! Input variable information file name. This file needs to be processed ! first so all information arrays can be initialized properly. VARNAME = ../varinfo_cfsr.dat ! VARNAME = ROMS/External/varinfo.dat ! Number of nested grids. Ngrids = 1 ! Number of grid nesting layers. This ...

kobl1201

When i change initial time like t = 0, will i change the other kind of times?? For example, ocean time, temp time and salt time and so on.

kobl1201

the dimension of this variable, zeta west,is (zeta_time, erho).

And i run the model for 1 year(366 days)

This picture shows that zeta time.

and i use nudging scheme for 3D velocity. This time interval is 1 day.

kobl1201

hello

When i run roms model, i have some errors.

the log file show that like this.

'current model time exceeds ending value for variable : zeta west'

Even if, i check time of zeta west, i could't understand why this error represented.

kobl1201

hello

when i want to use nudging method about 3d-field, shoud i need u and v values in climate_file??

if i need u,v values in clim_file, i can't slove this problem.

how can i do??

kobl1201

Hi~

How can i can the depth for one points from ROMS output???

i knew ROMS model uses sigma coordinate.

But i want to know the depth (Z- coordinate) for one point or all grids.

how can i do for getting the depth?? Which value will i check for this work???

kobl1201

hi~

i don't understand T0(background potential temperature) and theta_S.

my model was set as T0 = 14d0, theta_S = 1.0d0

what is effect this parameter in roms model, if i set T0 = 10.0d0 and theta_s = 2.0d0???

kobl1201

i use open boundary in ROMS model.

i was OBCFAC. it mean that ratio of inflow and outflow nudging time scale.

if i set OBCFAC = 1.0d0, what different effect between OBCFAC = 366.0d0 and OBCFAC = 1.0d0 for running model over 1 year.

my output average is 1 day, and M3NUDG = 1.0d0.

kobl1201

i want to do data assimilation of the observational data.

How can i do??? Actually, i found Roms wiki. But i couldn't understand it.

i just changed M3NUDGE=366.0d0 -> M3NUDGE=1.0d0 or can i need other work??

my climatology data was combined only temp and salinty.

kobl1201

i will make new initial data. This data was made a method. Firstly, i was running the Roms model set output from HYCOM model. And then, I compared the results, which averaged daily, with observation. However, the results showed too bad. So i thought that the results influenced used obsercational ...

kobl1201

Thank you. i will tell you my purpose. I was running the model set HYCOM model for initail data. And then i checked the results. The results for currents was bad comapred with observation data. So i changed initial data. and the initial data apply each 1 day for 1 year. and i will be running model ...

kobl1201

thank your reply. I read your answer ,and suddenly i have other question. What does Nudging mean?? I have understood when the model is running, it mean chganging initial data. i mean if the model calculate 2012.01.01, intial data apply data of 2012.01.01 time. And if model calculate 2012.01.02, the ...

kobl1201

I want to be nuding for initail data. When i run the model, i need to the climatology file, boundary file and initial file. Among three data, if i wanto to nudge for 3-D fields, which files should i change?? I think i need to change initial data by making 1 year data (interval 1 day), and M3NUDG ...

kobl1201

Hello when i use ROMS Model, I want to use Nudging method. Actually, i want to output data for u,v velocity. Therefore, i made initial data over 1 year, and then want to nudge interval 1 day. How can i do?? I saw TNUG and M3NUDG. How can i write down this parameter for nudging interval 1 day. thank ...

kobl1201

Hello

i have a very simple quations.

what does mean ETA-term??

thank you

kobl1201

Thank you very much.

Have a good day~

kobl1201

Generally, u or v velocities are computed through governing equations. However, i want to convert the method. Therefore, i need what program are finally computed u or v velocities for getting NC files before. I think that the velocities are computed in the directorty, which is Nonlinear or Utility ...

kobl1201

i can't recognize that what programs are finally computed u or v velocities.

Could you tell me the programs' name ??

Thank you

kobl1201

I mean when Roms model creat output file, there has been computed temperature, salinity, u and v velocities and so on based on governing equation. However, i want to convert the governing equations which are computed u or v velocities. Because i want to use u and v velocities for new initial data ...

kobl1201

Hi

I have some questions

I want to convert the method of computing u or v velocities. However, i can't recognize what programs are computed for u or v velocity.

Could you tell me the programs?

Thank you for your help.

I am looking forward to your reply.

Have a good day~

Search found 60 matches