Ocean Modeling Discussion

bibi951

Hello,Kate I do tide simulation,and I think bottom friction really matters a lot.I use RDRG2 but it seems not good for my case,since my area is too large(the whole China Sea). I notice that you use spatial distributed bottom drag. What's your rationale for a grid size dependent (if that's what it is ...

bibi951

Hello,everyone. I tried another initial condition derived from hycom, and used nudging mentioned above (20 days of nudging time scale for outflow , 0.5 day for inflow).It partly works(I don't know which contributes more to this, intialization or nudging),however not good enough. ssh and barotropic ...

bibi951

Thanks for your suggustion.
My OB data is given daily to the interior domain.
I have not tried Clamped OBC yet,maybe that would help.
I am also thinking about nudging ,would 20 days of nudging time scale be OK for outflow , 0.5 day for inflow ?

bibi951

Yes,parameter f is correct and it varies along latitude.
Maybe I should try another initial condition since it may conflict someway with the open boundary information,consindering they come from different model and differ in time.
I will report the result after this is done.

bibi951

Thank you for teaching me how to diagnose where the problem may be,Kate. I really appreciate it. I should have shown the CPP flags before. This is what I use NANHAI NANHAI SHUZHISHIYAN ANA_BTFLUX Analytical kinematic bottom temperature flux. ANA_SMFLUX Analytical kinematic surface momentum flux. ANA ...

bibi951

Sorry, I didn't mention my initial condition is taken from SODA-1.4.2 and is climatology yearly mean.Both my soda and hycom data are on z-level surfaces. Maybe the initial condition is the problem? I will show pictures of the process from intial condition to the day 87;and the time varing East open ...

bibi951

Thank you,Kate. You are always so kind. Sorry,I didn't realize the login problem. Here is the new address http://photo.blog.sina.com.cn/photo/46f52cc946f32f1f7e4be http://photo.blog.sina.com.cn/photo/46f52cc946f32f1f81612#pic and I will send you an email for pictures. My grid is North-South-West ...

bibi951

Hi,all I do a simulation using Hycom's Global simulation product for my open boundary. I have checked the Open Boundary section field after interpolating them to Roms's s- coordinate.It seems fine although not exactly the same. But I have some doubt about the field consistent.I just use the same ...

bibi951

does the RNT_tools(http://www.o3d.org/RNT/) only run on MATLAB higher than 2007a? i encountered many problems when compile and run on MATLAB (R14)2004 I learn from this forum that RNT is highly OS and ARCH dependent.So someone use this TOOL please show me the detailed OS and ARCH requirement.I ...

bibi951

I change from Winxp to Ubuntu. But I failed again to recompile the mex-files.I followed the complie.txt said in the RNT website,changing integer*8 to integer since my 32bit. Here is the error message: rnt_compile Warning: You are using gcc version "4.3.3-5ubuntu4)". The earliest gcc version ...

bibi951

Hi,thanks for your advice.I use win32,and I always encounter this error during running. http://photo.xiaonei.com/getphoto.do?id=1793430397&owner=89107472&curpage=0&t=#1793430397 mainly says forrtl:severe (38):error during write,unit -1,file CONOUT$ rnt_oa2d.mex.dll(etc) unknown unknown what is that ...

bibi951

Hi,everyone. I recently learn to use RNT matlab toolbox released by E. Di Lorenzo (Thanks for his excellent job).I find it really wonderful. However,when I mex the Mex-files,a warning occurs,saying: rnt_compile rnt_fill3ab_mex.F(371) : Warning: In the call to INVERSE, actual argument #5 does not ...

bibi951

Thank you for your kind reply, wilkin.
I checked my netcdf,finally I found that my latitude orientation is wrong.
Now Roms reads the right data.
Thanks

bibi951

Hi,everyone I took data from Ncep (2.5*2.5) and oaflux(1*1),And I wrote cdl file as Kate adviced.However I encounter this error NLM: GET_STATE - Read state initial conditions, t = 0 00:00:00 (File: initial_lJ.nc, Rec=0001, Index=1) .......................... GET_2DFLD - surface u-wind component, t ...

bibi951

Hi,Kate Thank you for your help.It is really kind of you. I use it and I guess Roms does this work,and I don't have to worry about the Missing_Value problem? The second question is about setting the initial condition.I think ROMS will take it as in the xi_eta_s coordinates.(ROMS will not do this ...

bibi951

Hi,Mauricio
Do you check the _REDUCED CASE , and see why it blowup at the open boundry.
I am quite confused about how to debug it as I mentioned above.

Regards!

bibi951

Hi everyone I am preparing ROMS Atmosphere foring. And I quite confused about this: 1)Suppose I have interpolate data to the computing grid coordinate(xi,eta),and I have done the rotaion. But how about the scale factor during the coordinate transformation, I just see the rotation, lengh scale doesn ...

bibi951

I just think of what is the balance in the momentum equition in my simple case? Since my Open Boundry is very close to the equator, will there be a prolem considering f ? I don't think the horizontal pressure gradient is the problem since the depth gradient and my only one vertical layer? But why ...

bibi951

I had another try.There is no zeta problem mentioned above(seems reasonable-all around 0.6 when ubar goes to the sky).The blowup begines here(red),and spread up-right. The depth is :where 10 means land,and I have checked my mask to be consistent. 10 10 874.75 929.36 833.89 713.79 10 489.22 805.24 10 ...

bibi951

Hi,everyone. My case underestimate the tide amplitude,so I suppose that my ubar vbar are not consistent with my zeta.So I use reduced physics by defining FSOBC_REDUCED and ****_M2REDUCED.However it blowup after 8.5days simulation. I checked the ubar and the zeta.I find that there is only one point ...

bibi951

I am doing a tidal forcing only 2D case.I am quite confused about my test cases : case 1: minimum depth 10m,RDRG2 == 3.0d-03,tide elevation amplitude 30-100cm lower than public quite recognized results.And lower(above30cm) than the tpxo_data I used for forcing. case 2: minimum depth 10m,RDRG2 == 2 ...

bibi951

Hello,I have a simple question to ask. I don't quite understand this : In set_tides.F line 434 and below, ! To avoid having two boundary ! points for these arrays, the values of "zeta_west" and "zeta_east" ! are averaged at u-points. # ifdef ADD_FSOBC zeta_west(j)=zeta_west(j)+ & & 0.5_r8*(Etide ...

bibi951

Hi,everyone. I want to use only real tides on the open boundry.So I write this in my .h file /* Basic physics options */ #define UV_ADV #define UV_COR #define UV_VIS2 #define UV_QDRAG #define MIX_S_UV #define SOLVE3D #define NONLIN_EOS #undef SALINITY /* Basic numerics options */ /*#define UV ...

bibi951

Thanks for the above replies,they are very clear. I get this error message: Minimum barotropic Courant Number = 9.63854617E-03 Maximum barotropic Courant Number = 2.99427149E-01 Maximum Coriolis Courant Number = 1.12416582E-02 GET_NGFLD - tidal period (Min = 4.32000000E+04 Max = 9.29494800E+04) GET ...

bibi951

Hi,everyone I created my domain grid with seagrid.The first corner(*) is southwest(in real world),so the first edge is south(in real world).But what is the west boundary in Roms Model? Is it the first edge(south in real world)? Or ...? If one boundary is half open,(there is a gap ),I should use OBC ...

bibi951

Hi,Kate. I tried Gfortran,it also failed to do it. I have rewrite this section to a explicit do loop. As follows: line 25 : character*2000 fmt1,fmt2,ftmp aroud 176: c output format fmt1='(f9.3,f9.3)' if(APRI)then fmt2='(f8.3,f8.1)' else fmt2='(f8.3,f8.3)' endif aroud 227: if(APRI)then write(11,trim ...

bibi951

Thank you,Kate.
I really appreciate your guide.Thank you.

bibi951

Hi,everyone. I am sorry to bother anyone here.I meet the annoying error again which I met a couple of days ago. I can't find what's wrong, so I hope someone can help me out. This is when I use OTPS for tide simulation. g95 -fno-second-underscore -fleading-underscore -L/usr/local/lib -lnetcdf -I/usr ...

bibi951

Hello,everyone. Back again. I now use OTPS to make tide forcing.I use Cygwin for windows,and G95,Gfortran are my fortran compiler.Now it seem that G95 don't recognize this.Any advice? $ make extract_HC g95 -o extract_HC extract_HC.f subs.f -L/usr/local/lib -lnetcdf -I/usr/local/include In file ...

bibi951

Hi.Thanks for your replies.
I have modified the bathymetre,and it works.
Very nice.
Thank you all

bibi951

Hello,everyone.I am new to Roms. I just have build up my first application program.As I ween,there is something wrong. It is just an imatation of windbasin.I use my grid-The area range is 3-41N,100-137E.I didn't do much preprocess on bathymethy.And I use four wall boundry conditions,surface and ...

bibi951

[quote="kate"] Paul Budgell has this code and it is now in my branch of ROMS. [/quote]
Is this shared? If it is,Would you please give me a link to this?

bibi951

If you have a large enough domain to require the body force potential tide, then that's another story. Paul Budgell has coded it up and it is in my svn branch. You then have to provide the tidal potential everywhere as more input fields. Do you mean that tide can be generated by ROMS by introducing ...

bibi951

Hi,Kate.I also use matlab . But how can I transform my realistic Forcing data onto my computation grid(curvilinear coordinates ).Do I need to do this myself or ROMS does this job? I have just found a matlab ROMS toolbox.But it seems have few discription files.Can this toolbox give me the right input ...

bibi951

Hi,I use Field-preparation package released by Hernan G. Arango. I find some small mistakes in Hydro,OA,Forcing,Initialization package. I fix some but I am not sure whether I did the right modification. for examples: OA : def_oa.f 240: status=nf_put_att_text(ncoutid,varid,'units',5,'meter',varid) or ...

bibi951

Thanks for your guide. I had a try. When use G95 or Gfortran : g95 -DNO_EXIT -g -O2 -I /usr/local/include -L/usr/local/lib -lnetcdf -o /h oamain.o: In function `MAIN_': /home/Sun/roms/oa/oamain.F:268: undefined reference to `nf_put_var1_double__' ... /home/Sun/roms/oa/oamain.F:79: undefined ...

bibi951

I use Cygwin. When I build OA package.This error occur：detail in the attached file Any suggestion would be appreciated. g77 -c -DNO_EXIT -g -O2 def_oa.F g77 -DNO_EXIT -g -O2 -o /home/Sun/bin/oa oamain.o acor.o asm.o atg.o bcintrp.o bcucof.o caldate.o day_code.o def_oa.o dhlev.o diagn.o down_shift.o d ...

bibi951

Oh ,
I use
mpirun -np 4 ocean < ocean.in

not
mpirun -np 4 ocean ocean.in

A silly mistake

bibi951

Hi. I now use SGI to run this example cblast in parallel. Its compilation ends with no error.BUT when execute oceanM: MPI: On host a3700, Program /lun/home/lujiuyou/roms/Projects/cblast/Forward/oceanM, Rank 0, Process 32711 received signal SIGSEGV(11) MPI: --------stack traceback------- sh: line 1 ...

bibi951

Thank you for your replies.
I've just realized having made a simple mistake
qsub a job in a incorrect way on my SGI.Not the problem of this model.
Many Thanks and Best Regards

bibi951

lu@a330:/> ls -l /usr/local/netcdf3.6.2/include/netcdf.mod -rw-r--r-- 1 root root 263985 2007-08-03 11:21 /usr/local/netcdf3.6.2/include/netcdf.mod "cp /usr/local/netcdf3.6.2/include/netcdf.mod" Sorry, I did not complete it with its destinition(Build directory for scratch). My SGI OS is unix. But ...

bibi951

I now use SGI to run this model(Serial).But I am not the administrator,and I canot install Cygwin as I wish.So I use SSH (installed) and netcdf(installed). What can I only do is to modify my build_cblast.sh to be adapted to the installed system.So In build_cblast.sh setenv NETCDF_INCDIR /usr/local ...

bibi951

First,I would like to thank all your replies. I now use SGI to run this model example cblast in parallel.This error occur: cp: cannot stat `/usr/local/netcdf3.6.2/include/netcdf.mod': No such file or directory make: *** [/lun/home/lujiuyou/roms/Projects/cblast/Forward/Build/netcdf.mod] Error 1 But I ...

bibi951

$ ./oceanS < cblast_stations.in ocean_cblast.in Error: Resolution, Grid 01: 0000x0000x000, Parallel Threads: 1, Tiling: 000x000 ROMS/TOMS: Wrong choice of domain 001 partition or number of parallel threads. NtileI*NtileJ must be a positive multiple of the number of threads. Change number of threads ...

bibi951

HI,I have use : export NETCDF_INCDIR=/usr/local/include export NETCDF_LIBDIR=/usr/local/lib before cd ${MY_ROMS_SRC} in build_cblast.sh script. It also leads to this error mentioned above. But if I use NETCDF_INCDIR ?=/usr/local/include NETCDF_LIBDIR ?=/usr/local/lib in the CYGWIN-g95.mk It can ...

bibi951

I am a very beginner of ROMS.When I run the application example of Calast I get this error makefile:233: INCLUDING FILE Compilers/make_macros.mk WHICH CONTAINS APPLICATION-DEPENDENT MAKE DEFINITIONS cp -f /usr/g95soft/netcdf/include/netcdf.mod Build cp: cannot stat `/usr/g95soft/netcdf/include ...

bibi951

thank you very much for the tutorial! i'm new in both ROMS and Unix,so i have to follow step by step the istall instruction. but why this tutorial is different with the one on the wikiROMS at https://www.myroms.org/wiki/index.php/ROMS_Cygwin#Install_Cygwin ? because on that one there are Install ...

Search found 47 matches