Ocean Modeling Discussion

Thanks arango and klyons.I just use the RNT to set up my ncfiles.
Thanks for your suggestion,I will try
Thanks!

Thanks klyons. Do you have any land in your model domain? yes,I have land in my model domain,sometimes,I set up the grid.nc file once again,this question is not exit,but sometimes it still exit.In the programe rnt_2s.m ,I think the -9999999.0 just take place of the deficit value of the variable. one...

Thanks klyons. GET_2DFLD - surface net heat flux sensitivity to SST, t = 45 00:00:00 (Rec=0002, Index=2, File: nwpacific-forc.nc) (Tmin= 15.0000 Tmax= 345.0000) (Min = 1.00000000E+37 Max = -1.00000000E+37) This is really a question.I set up the ncfiles(MAIN_ncfiles.m)once again and get the dQdSST. I...

Thanks kate,christianmohn and susonic,you are all so kind. I would try to run the model on 1 CPU as a quick check to make sure that there is no issue with MPI parallelization. If memory is a constraint, you could do this test with a coarser horizontal grid. Yes,I just test it with 1 CPU,but it blow ...

Thanks kate.I redefine TS_FIXED,and test again,but unfortunatly it blow up again and I close all 3-D and traces（alos test include open that）,just run the 2-D(I don't know whether it is right?) Also,test with different dt,dt=2,1800,... I really have no ideal for that. What's the matter? Thanks Physic...

Thank you lfowler. I just check my initial nc files ,the temp,salt value are all normal,the same with zeta,u,v,ubar,vbar,and some u,v are 0,but others are not 0. Is it dt<sqrt(((dx^2+dy^2)/(9.8*h))?If it is,then dt<4132? Minimum X-grid spacing, DXmin = 4.15066970E+01 km Maximum X-grid spacing, DXmax...

Hi kate,thanks for your reply,you are so kind. I just try more shorter timestp,and I set dt=1800;900;100;20;4;0.5 ,but unfortunately,all the result is the same: -------------------------------------------------------- dt=20 half a year STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME ...

hello everyone,I just want to set up my own application ,my domin is northwest pacific(9.08°N-46.80°N，105.39°E-179°E)，but the model blowing up just 1 step. ---------------------------------------------------------- NORTHWEST PACIFIC OCEAN Operating system : Linux CPU/hardware : x86_64 Compiler syste...

ok,thanks!
Now it is ok,but i still don't know the really reason for this question.

You should define 'mode' when you want to open NetCDF file like follows; nc=netcdf('ncfile', 'nowrite') Thank you hyy!I know you mean.But the question is the same >nc=netcdf('NCEP_landmask.nc','nowrite') Index exceeds the matrix dimensions and even >nc=netcdf('NCEP_landmask.nc','w') Index exceeds t...

Hello everyone,recently I have encoutered a strange question about netcdf; in matlab; >nc=netcdf('ncfile');(for example nc=netcdf('NECP_landmask.nc')) Index exceeds the matrix dimensions I just want to open the netcdf file ,and even the file is not so big. I do no know anything abot the question . t...

I think RNT(http://eros.eas.gatech.edu/RNT/) is useful
And the ROMS tutorials(http://eros.eas.gatech.edu/ROMS-Tutorial/tutorials.html)

Could anyone can help me?This is my first application,and I have struggled it for some days.
Thanke you!

Hello, I just follow the ROMS Tutorials(http://eros.eas.gatech.edu/ROMS-Tutorial/tutorials.html)(MAIN_grid.m)and generate the grid file 'nwpacific-grid.nc' --------------------------------------------------------------- nwpacific-grid.nc gridindi.id = 'nwpacific' gridindo.name = 'northwestern pacifi...

Thanks ,now i konw how to resolve this trouble,in file:rnt/rnt_gridinfo.m ,I add
gridindo.cstfile=which('World_coast.mat')to that m file.
World_coast.mat is the world coast file I have set up;
Then :
>>rnt_plcm(.1/grd.pm/1000,grd)
the figure have the coast line
thanks !

Hello everyone ,when I follow the ROMS Tutorials (http://eros.eas.gatech.edu/ROMS-Tutorial/tutorials.html)(thanks for this tutorials)to se up my own grid file(the northwest pacific ocean:105°-179°E,9.6°-47.8°N),follow the MAIN_grid.m %========================================================= % Prepa...

[quote][1)>> [lon, lat, z] = read_srtm30plus([10 20], [40 46], 30.1);
Error in ==> read_srtm30plus at 60
pix = imread (url, 'tif');
/quote]
hello mistral,I use the first method to extract the bathymetry ,and I have the proble the same with you.
Have you resolved it?
help!

Thank you kate,I am sorry for my late reply.
I kown your mean and I hanvn't the arpack and parpack
I have download the arpack and parpack file,but I can not install them successful.And I need to try it again!

Maybe your Roms need to update.
svn update

hello,everyone.I want to test the application of A4DVAR_TOY . But first,the make is not accessed. slvester@OceanStar:~/roms> make ...... ...... .... /opt/mpi/mvapich/1.1/gcc.pgf90/bin/mpif90 -O3 -tp k8-64 Build/esmf_roms.o Build/master.o Build/ocean_control.o Build/ocean_coupler.o Build/propagator.o...

OK ,now it is ok,thank you robertson.

Yes,Thank you roberston,you are right Do you have a copy of default.cnt in the directory where you are executing the : Code: csec < csec_upwelling.in yes. ------------------------------------------- I set the csec_upwelling.in /dcfs2/users/slvester/roms/roms_plot/Data/varid.dat /dcfs2/users/slvester...

hello everyone.I donwload the roms plot package follow instruction: svn checkout https://www.myroms.org/svn/src/plot roms_plot after some change,I install it successful . And there are ccnt,cnt,et al in the PLT_BINDIR. But when I want to plot the upwelling output,follwing the Cross-section Plots : s...

Hello stone.I will do some flow simulation with ROMS IS4DVAR in the Northenwest Pacific domain(It is also very large) and assimilate Argo data et al. I am the beginner with ROMS. OI is more easy ,but if OI is still in the ROMS?I can not find it. there are some advice for data assimilation,I think it...

I noticed that you just joined ROMS a few days ago. The first rule of ocean modeling and data assimilation is to know your application well before you even think what data is needed to assimilate. From the scientific point of view, data assimalation cannot be treated as a black box because it requi...

ok,thank you!

arango wrote:Your posts are telling me that you are not reading any of the information that it is available on WikiROMS.

thanks for your suggestion!I will read it more seriously

Hello everyone. I finished the make operator yesterday and got the file oceanM .But I find that there are so much confusion code in the file,is it normal? so when I enter: ./oceanM < ROMS/External/ocean_upwelling.in there are error:PMGR_COLLECTIVE ERROR: unitialized MPI task: Missing required enviro...

you are right .It is ok now.Thank you very much.
the error is so small that it is hard to find.
Thanks!

It is /oceanM
I sets the BINDIR for
BINDIR ?=
I also try to set the BINDIR?= /opt/defs2/users/slvester/roms/.but the result is the same,and then I set the BINDIR ?= 'pwd',the result is also:
/usr/bin/ld: cannot open output file /oceanM: Permission denied
make: *** [/oceanM] Error 2
I don't know why

thanks for you ,you are very patient.I am sorry for what I said,but I really can not find the oceanM file,I should try it again. some results are there: /usr/bin/cpp -P -traditional -DMPI -DLINUX -DX86_64 -DPGI -D'ROOT_DIR="/dcfs2/users/slvester/roms"' -DUPWELLING -D'HEADER="upwellin...

Thanks! But there are two question.first,I can't find the file :oceanM,maybey there are not this file.Second,/usr/bin/ld,what's this mean?

I get the following error message at the last two lines of the "make" run :
/usr/bin/ld:cannot open output file/oceanM :Permission denied
make: ***[/oceanM] Error 2

How can I resolve it?
Thanks