Ocean Modeling Discussion

mashinde

What is the difference between eddy resolving and eddy permitting models ?
How to identify them ?
Is there any criterion ?

mashinde

HI,

It is better to read data FAQ from ECMWF.

http://old.ecmwf.int/products/data/arch ... a_faq.html

mashinde

Hi,

I am also interested in heat budget. Which quantities do I need to use from dia* file and which formula i have to use ?

Thanks

mashinde

Dear All,

what are the Advantages and DisAdvantages of GLS scheme over KPP ?
which scheme is preferable ?

Thanks

mashinde

Hi,

can we store EKE in roms output netcdf file ? if yes, how ?

Thanks.

mashinde

Hi,

How can I download ROMS3.4 ? can someone give me svn command ?

thanks

mashinde

Hi all, I am confused with following keys in cppdefs.h and it's physical role in ROMS. Also the possible good combination of them. solar_source, diurnal_srflx and atm_press. And what will be right choice of it for 6 hourly forcings. currently, i am using for 6 hrly forcings. #define BULK_FLUXES #if ...

mashinde

Dear all,

does anyone experiences slow performance of hpc cluster with netcdf with pnetcdf option ?
is it good option to build netcdf with pnetcdf option ? will it affect simulation ?

any hint/suggestion please ?

thanks

mashinde

Hi kate, (1) All forcings are 6 hourly. The short wave radiation is Net surface short wave radiation, I think it included cloud effect. Yes, there is albedo field in ERA. let me try without diurnal_srflux. (2) which short wave radiation field is required : Surface Net Short Wave Radiation OR ...

mashinde

Dear Roms Users, I am forcing the MODEL with 6hourly era intrim data for test run for 1998. I use following variables.... ROMS_uwind --> ERA_Intrim_u10 (10 metre U wind component) ROMS_vwind --> ERA_Intrim_v10 (10 metre V wind component) ROMS_tair --> ERA_Intrim_t2m (2 metre temperature) ROMS_pair ...

mashinde

Hi kate, I am following the USwest example for creating curvilinear grid with grid refinement. I want to learn that skill and apply in the Mediterranean sea, especially at the Gibraltar strait. But I am getting errors. I already copied files from USwest directory to Include. Now I am getting ...

mashinde

Hi kate,

I am sorry , i can't make it with curvilinear grid for USwest example.
can you please give steps how to proceed in ascending order ?

thanks

regards

mashinde

hi kate, It's works fine now. I successfully created bathymetry for gom example. I use editmask for handling land/sea mask, which i always use. Now creating spherical grid using Gridpak is done for GoM test area. Now I am trying to create curvilinear grid for USwest example. This is quite difficult ...

mashinde

Hi kate, it works fine... I successfully created GoM nc file with ETOPO5 which is a regular grid, but i don't understand some points, they are as follows (1) when i run command /bathtub ETOPO5 minimum depth (negative is land elevation) 10 maximum depth 8000 apply mask? (T/F) T im 4320 2161 0 ...

mashinde

Yes I downlaoded the etopo5.nc file. I set env variable as export etopo5=/data/rmed16_etp1_input/gridpak/etopo5.nc then I run ./bathtub etopo5 but it gives me error ... minimum depth (negative is land elevation) 10 maximum depth 4000 apply mask? (T/F) T EXTRACT - cannot find variable: topo_lon ...

mashinde

ok,

meanwhile I keep USwest aside and trying the example of Gulf of Mexico. for this example there grid_gom.nc file is present .
now how can I run bathtub ?

thanks

regards

mashinde

Hi kate, Ok.... now it is giving error while interpolating etopo5 into grid file. /bathtub minimum depth (negative is land elevation) 10 maximum depth 8000 apply mask? (T/F) T EXTRACT - cannot find variable: topo_lon EXTRACT - cannot find variable: topo_lat im 21601 10801 0.0000000000000000 EXTRACT ...

mashinde

Hi kate, Thanks for the new version. I successfully created all executables now. Now I am playing with giver USwest example. But it is not doing well. i am getting following errors for USwest example.. grid < grid.in rectangularity error in mapped contour at iteration 1 is 3.7708E-02 rectangularity ...

mashinde

Hi kate sorry about that, now i fixed it.But no i am getting following error .... make -f Makefile.Gnu gfortran -c -g -O2 -I/usr/local/include grid.F gfortran -c -g -O2 -I/usr/local/include checkdefs.F gfortran -c -g -O2 -I/usr/local/include def_grid.F gfortran -c -g -O2 -I/usr/local/include get ...

mashinde

Hi kate, I changed it to -lnetcdf , but still the same error. my flags in the Makefile_Gnu are... CFT = g77 CPPFLAGS = CPP = /lib/cpp -P $(CPPFLAGS) # GLDR = ncargf77 GLDR = g77 LDR = g77 RANLIB = ranlib LDFLAGS = # FFLAGS = -g $(CPPFLAGS) FFLAGS = -g -O2 $(CPPFLAGS) # Add -I /usr/local/include for ...

mashinde

Hi, I am using Fortran based gridpak to create curvilinear grid with refinement at the straits. But when I am trying to make to create grid exe file, then it is giving me following error, ofcourse the error is related with netcdf configuration, but i did not understand how to fix it. [mshinde ...

mashinde

(1) Yes, I have the SSS monthly climatology file from WOA. (2) I still did not understand clearly the physical role of Tnudge and Relaxation and timescale and how to do it in ROMS code, only enabling SRelaxation in cppdef.h is not sufficient, I think. Suppose I want to do SRelax SSS for 10 days for ...

mashinde

Thanks for reply. I have another doubts (1) So to imply SRelaxation and SCorrection, ROMS needs salinity forcings files irrespective of bulk forcings. (2) If I am using 6 hourly atmospheric forcing then what will be Tnudg in days ? how to determine this. (3) What will happen if I don't apply ...

mashinde

Hi All, I am using SRELAXATION ON in cppdefs.h for salinity relaxation and defining in ocean.in file as .. ! Nudging/relaxation time scales, inverse scales will be computed ! internally, [1:Ngrids]. TNUDG == 2*30.0d0 ! days. (1) Will it relax model salinity towards climatology in entire domain or at ...

mashinde

Hi, I am forcing ROMS with era interim 6hrly data for 1998. for one year simulation it gives cold ssts and temp is not varying. I don't know why this is happening. my cppdefs.h is as below... /* Basic physics options */ #define UV_ADV #define UV_COR #define UV_VIS4 #define MIX_S_UV #define TS_DIF4 ...

mashinde

Hi Kate,

Thanks for the comment. I am also getting cold SSTs throughput the year. No change in temperature. Do you think, it is not maintaining the relation between short wave radiation and ocean temperature. I did not have such problems in case of climatological runs.

mashinde

Hi, I am forcing the model with era40 data. I am also saving net swrad and lwrad to his file. But it swrad gives zero values. here I am attaching ncdump of swrad from his file and cppdefs.h used. The lward is givng values in the his file. I checked input file of swrad and it is not with zero values ...

mashinde

Hi,

I am using era interim 6hrly data with bulk forcing, SRelaxation. Now I want to OFF only rainfall and keep SRelaxation ON. How can I do it ?
my bulk options are ...

#define BULK_FLUXES
#ifdef BULK_FLUXES
#define EMINUSP
#define LONGWAVE_OUT
#endif

any comment/hint ?

thanks.

mashinde

Hi,

Is it necessary to apply climatology file in realistic applications ? what will happen if clim file is not used ?

Thanks

mashinde

Hi, After doing some search along with a paper by Large et. al.(1994), I found that the mixed layer depth computation by KPP is based on the critical value of turbulent processes parametrized by a bulk Richardson number. So it is not based on temp or density criterion ???? For more and short ...

mashinde

Hi Kate, Thanks for the kind reply. At present I am using GLS with epsilon option because with KPP model blows up after some days. I am applying ROMS in the Mediterranean sea. Therefore, I want to check inflow i.e. incoming Atlantic freshwater flux and outgoing i.e. outgoing Mediterranean freshwater ...

mashinde

Hi,

I have following questions

(1) can ROMS store mixed layer depth in a output file? Is depth of surface boundary layer (Hsbl) is same as mixed layer depth?
(2) How to compute freshwater transport(inflow, outflow and net) and it's time series ?

Thanking you

mashinde

Dear Prof. Arango, Thank you for kind suggestion. I am unable to get rainfall and radiation fields at every 6hrs. Other fields except rainfall and radiation are available at every 6hr for analysis time [00,06,12,18] and at every step for forecasting time [0,3,6,9,12]. So I extracted these fields ...

mashinde

Dear Julian,

Yes, I tried this data. But I am not getting total rainfall and radiation fluxes data at every 6 hours.
Do you have any idea/suggestion how to get it ?

Thanks

mashinde

Dear Friends,

I am getting problem in finding 6 hourly data for radiation fluxes and rainfall for ERA Interim. Is there anyone who already use 6hourly data from ERA Interim, then kindly suggest/give idea how to obtain these data ?

Thanking you

mashinde

Hi Shou,

I replace "Surface thermal radiation" with "Surface thermal radiation downwards" from ERA40 data set.

regards

mashinde

Dear julian,

I fix the problem. The problem was in input radiation flux. Model runs fine and no more negative temperatures.

Thanks

mashinde

Dear ROMS users,

I am running model with ERA40 daily data(00hrs)for 1995-2000. I am getting negative temperature. I don't know why it is giving such -ve temperature.
I attached here one figure from 249 day and also log file.

any suggestions/comments ?

Thanking you

mashinde

Hi Kate, Thanks for the kind reply. I am still not understood the time attributes in forcing file. let's take example of monthly forcing. I am forcing the model with ERA-INTERIM monthly forcing for 1990-2010. So I download the data from 1990-2010 which has 252 months. I had made forcing file for ...

mashinde

Hi kate, Thanks for the kind help. now model is running. But I have following doubts about "cycle_length" (1) I am using era40 monthly data from 1990 to 2010(total 252 months). Following are my ncdump -h of Tair. ncdump -h eraint_tair_mn_1990_2010.nc netcdf eraint_tair_mn_1990_2010 { dimensions: xi ...

mashinde

Hi Kate, Thanks for your kind suggestions. The problem solved and problem is running, but it blowing-up at 56 day. I am forcing the model with ERA-Interim monthly data from 1990-2010. if restart file save the records for latest 2 records why it is giving the error as .... Blowing-up: Saving latest ...

mashinde

Hi, I am trying to run ROMD with ERA-Interim forcing for 1990-2010. I am getting floating invalid error.... -------------------------------------------------------------------------- forrtl: error (65): floating invalid Image PC Routine Line Source oceanG 000000000078495E nf_fread2d_mod_mp 234 nf ...

mashinde

Dear friend, I am trying ROMS with era-interim monthly forcings from 1990-2010. I have following errors. (1) In log file it always shows the NaN values in forcing file(log file attached here). (2) Model starts at STEP=864 and writing the NaN values and model blows up at 1st step. for details please ...

mashinde

Dear friends, I am running the script 'd_core2_frc_md.m' for medi sea. But I am getting following error. d_core2_frc_md ??? Undefined variable "netcdf" or class "netcdf.getConstant". Error in ==> d_core2_frc_md at 120 mode = netcdf.getConstant('CLOBBER'); % overwrite my matlab version is R2008a. any ...

mashinde

Hi

I am using ERA40 from 00hrs, 06hrs, 12hrs and 18hrs. so in that case what should be the multiplication factor for ssr ?

any suggestions ?

thanks

mashinde

Dear friends, I have following questions.... (1) I have climatological runs for 10 years. Now I want to force the model from ERA40 dataset. For example, the starting date is 01-01-1995 and model run time is for next 5 years. So, How can use my previous climatological run ? e.g., I want to start ...

mashinde

Dear Friends,

I have error while restarting the model. When I restart the model it looks for the avg_rmed16_etp1_0000.nc. This files is not created before restarting the model. I am using ROMS3.4. A log file attached here for details.

any comments/suggestions.

mashinde

Dear Mathieu,

Thanks for the reply. I use GRID_LaplacianSelectSmooth_rx0 of LP bathy to smooth the bathymetry. How can I compute Horizontal Pressure Gradient term ?

-Mahesh

mashinde

Hi,

Model blows up at day 344, when I use bulk forcing. I use ROMS 3.4
I attached here log file.

Thanks

mashinde

Hi kate,

Yes,I put Vtransform=2 and Vstretching=2 in the script. The version of script is 586. I also attached here the script for details.
I am using ROMS 3.4. My matlab version is R2008a.

mashinde

Dear friends,

I have a problem in initial forcing file. I created initial file using arango's matlab script, i.e. d_initial.m . The file is created, but it gives error while running the model. The log file is attached here.

any suggestion/solution how I can overcome this error.

thanking you

mashinde

Hi Andres Sepulveda, I am working on ROMS application in the Mediterranean sea. I need to modify manually the bathymetry after smoothing to reach the required sill depth at Gibraltar straits. I had seen this post and tried the script. But I don't know how to use it properly. Is there any user manual ...

mashinde

Dear Friends , I am building a grid with seagrid. I have following errors while reading bathymetry via read_srtm30plus.m script, [xbathy,ybathy,zbathy]=read_srtm30plus([-20 37],[30 49],60); ??? Error using ==> rtifc /tmp/tp8377d624_0366_4942_b7c9_623a47471fcf: Not a TIFF file, bad magic number 16188 ...

mashinde

Dear Friends,

I want to save vertical transect of ROMS fields in output nc file. How can I do this ?

Thanks

regards
Mahesh

mashinde

Dear Friends, I am running model with bulk forcings. I have following doubts. what is the relation between SRELAXATION and EMINUSP ? Is it like this .... case(1) when I have to switch ON SRELAXATION in cppdef.h then I have to switch OFF EMINUSP. case(2) when I have to switch OFF SRELAXATION in ...

mashinde

Dear Friends, After successful completion of climatological run(without bulk forcing), now I am trying to run the model with Bulk forcing. Initially, I am using monthly climatological Bulk forcing from COADS. The model blows up at timestep=19889 and at 47th day. here I am attaching cppdefs.h and ...

mashinde

Dear kkobchul and John,

I am also interested in ROMS coupling with AGCM. Kindly share your suggestions/comments.
My office email id is : samurai@tropmet.res.in

Thanking you

-Mahesh

mashinde

Good questions all, let's see what I can answer: Thanks Kate for the valuable suggestions. (1) I tried ROMS with climatological forcing in the Mediterranean sea. And it works well after initial settings. The large scale surface circulation features seems to be comparable when validated with the ...

mashinde

Dear Friends, (1) I tried ROMS with climatological forcing in the Mediterranean sea. And it works well after initial settings. The large scale surface circulation features seems to be comparable when validated with the MEDATLAS data with some biases. It further needs some fine tunning as the current ...

mashinde

Dear kate,

yes....i did it for upwelling test.
it is creating multiple file onwards where it stops..

in this case i just put his file in ininame and restart the model and works

thanks !

mashinde

Dear ROMS frds,

I have a very serious problem about how to restart the model ....i need urgent and helpful solution....

for details..please see the link....

viewtopic.php?f=14&t=2111

Thanks

mashinde

Dear Kate,

Yes I replace ininame by rst file while making restart.
I am also trying with upwelling test.

thanks

mashinde

Dear Kate, Yes, I got the same error. The only thing I discovered is that I have to store every records in rst file which is not a good idea as it will make rst file too large. currently, in my application rst file stored recent two records(by default) and I need to know how can I use this latest ...

mashinde

Dear Kate, Can you restart say the upwelling case after a day, saving at its default frequency? yes i did it for upwelling test.....in this case i stored all the records in rst file by changing LcycleRST == F, so model restart ok for upweling case. but in my application.....ROMS already completed ...

mashinde

dear kate, well i did some experiments with NHIS/NAVG/NDIA as follows........ ---------------------------------------------------------------------------- ! Input/Output parameters. NRREC == -1 !for new job LcycleRST == T NRST == 1 NSTA == 1 NFLT == 1 NINFO == 1 ! Output history, average, diagnostic ...

mashinde

Hi Kate, Thanks for the reply. I am still checking it..... is there any connection with mixing schemes...i am using KPP.? perfect_restart option is undefined. I am still unable to restart the model to create multiple avg/his/dia files . It is very important to know how can I restart the model in my ...

mashinde

Dear friends, I have error during restarting the model. -------------------------------------------------------------------------- 0 1 00:02:00 6.266689E-08 1.830794E+04 1.830794E+04 1.331604E+16 DEF_AVG - inquiring average file: /share/scratch/odc/mahesh/roms_files/output/opti_test/avg_rmed16_grd_e ...

mashinde

Hi frds, Now I want to force the model with realistic forcing from say. ECMWF 2009 monthly mean data. so my questions are : (1) I want to force the model with realistic data, for example, say Jan 2009(monthly mean here only one month) ECMWF data. so what should be the time in the nc file. I do not ...

mashinde

Thanks Tony,

do you any idea where I can get these scripts?

Mahesh

mashinde

Hi frds , I want to create river forcing for my application i.e. for Mediterranean sea. Does anybody have a idea about the river runoff data and how to create the river forcing files. OR is it possible to use a river runoff data from hydrological model. Even I also seen runoff data in ECMWF dataset ...

mashinde

Hi Kate,

I think I solved the problem........model run with bulk forcing...

Thanks for the advices.

mashinde

hi kate, Thanks ! these are the monthly climatology from COADS data. when I print time from the nc file it gives me like..... nc=netcdf('rmed16_grd_etopo5_e_uwnd_n.nc'); >> time=nc{'wind_time'}(:); >> time time = 15 45 75 105 135 165 195 225 255 285 315 345 Total 12 time elements so these are the ...

mashinde

Hi kate, thanks for the reply. I am doing climatological run. my u-wind nc file looks like. dimensions: xi_rho = 940 ; eta_rho = 424 ; xi_psi = 939 ; eta_psi = 423 ; xi_u = 939 ; eta_u = 424 ; xi_v = 940 ; eta_v = 423 ; wind_time = UNLIMITED ; // (12 currently) variables: double wind_time(wind_time ...

mashinde

Dear frds, I have following error.... This error comes when I use bulk forcing......... forrtl: error (65): floating invalid Image PC Routine Line Source oceanG 000000000069E78C time_string_ 490 get_date.f90 oceanG 00000000006947EF get_2dfld_ 456 get_2dfld.f90 oceanG 00000000005B200C get_data_ 73 ...

mashinde

Hi Kate, Thanks for the suggestion. But I am still confuse about the ntimes while restarting the model. Are you also writing to a restart file? Yes I typically need to change three things from the initial run to the restart runs: * Change nrrec from 0 to -1 to read newest record. YES * Change ...

mashinde

Hi frd, I have following questions regarding the restart the model..... suppose I submit a job for 10 years(long term) on HPC.I am saving his and avg files for each month. In between , due to some technical problem in HPC, model stops at certain timestep and after certain number of avg file. suppose ...

mashinde

Dear frds, I have error while reading the input forcing file.......the error msg is as follows: (I copy/past error msg only...) T Dout(M3vvis) Write out 3D momentum vertical viscosity. Output/Input Files: Output Restart File: /share/scratch/udic/mahesh/roms_files/output/rmed16_grd_e/rst_rmed16_grd ...

mashinde

Hi Dear frds, I have very general question : (1)about tiling: How tiling settings will affect the model and blow up ?(mostly the blow up...) What should be the ideal number of tiling for parallel computing ? Is there any relation between horizontal/vertical resolution and number of tiling ...

mashinde

Hi, I am still not get satisfactory time step to run the model for long simulation. I need the time step of the order of 100 to 200 seconds.I am using ETOPO5. I am smoothing bathymetry according to smth_bath.m from the ROMS matlab package. I have following questions about the r factor. (1)what ...

mashinde

Hi frds, I have doubts of time steping. I created 1/16 deg resolution with 72 levels over Mediterranean sea. when I use time step say dt=180 and ndfast=30, the it blows up after few time steps. If I put very small time step say dt=10 and ndfast=5 then model works fine. it is not good idea to keep ...

mashinde

Hi kate, here I am pasting the forcing files from input file. NFFILES == 8 ! number of forcing files FRCNAME == /share/scratch/mahesh/roms_files/input/COADS/rmed16_lwrd.nc \ /share/scratch/mahesh/roms_files/input/COADS/rmed16_pair.nc \ /share/scratch/mahesh/roms_files/input/COADS/rmed16_rhum.nc ...

mashinde

Hi Kate,

YES. This is correctly spelled and file path is correct. I don't know why it is reading only 1st forcing file.

any suggestion.

Thanks

Mahesh

mashinde

Dear Frds, I have an error while reading the second forcing file, here i am pasting error msg. Input Forcing File 01: /share/scratch/mahesh/roms_files/input/COADS/rmed16_lwrd.nc READ_PHYPAR - could not find input file: /share/scratch/mahesh/roms_files/input/COADS/rmed16_pair.nc Elapsed CPU time ...

mashinde

Dear frds,

I am running the model with dt=10s...it runs without blowing...but when i check the avg.nc file in matlab. then it can't store the fields.
so i am confused why it is not storing any field in avg file.

thanks a lot !

regards
mahesh

mashinde

hi frds, i am getting this error recently....when i submit the job to cluster /home/mshinde/ROMS/src/Projects/rmed16/oceanG: symbol lookup error: /home/mshinde/ROMS/src/Projects/rmed16/oceanG: undefined symbol: message_catalog is this error linked with cluster or it's model error ?? thanks mahesh

mashinde

Dear kate,

I am also running the model with dt=10s...it runs without blowing...but when i check the avg.nc file in matlab. then it can't store the fields.
so i am confused why it is not storing any field. Here i am attaching input file.

thanks a lot !

regards
mahesh

mashinde

Dear kate, here I am sending a log file. It blows up when i keep time step(dt) greater than 10s. If this is the case then it will take very long time to complete the spin up of say 20 years. You also find the time stepping and time stamps in this file. I want to store daily avgs for the months ...

mashinde

Hello Kate,

I really do not understand.....how can i set this time parameters in ocean.in file. It will be better if you extrapolate with examples.

Thanks

Mahesh

mashinde

Dear frds, I have following doubts.... (1)How to store hourly data ? (2)I want to run ROMS for 10 years to spin up..how to do that..where i set the time parameters in ocean.in file..any example ? (3)How can I store monthly data for 10 years of spin up ...daily avgs for each month for 10 years ? (4 ...

mashinde

Hi frds, I am running ROMS on newly installed SUN cluster with 32 processors.....i am getting following error. I am not geting the message from the error..is it coming from the model or is it coming from the netcdf/hdf lib installed in the cluster. I know there is something wrong with I/O. Below ...

mashinde

Hi ROMS users, I am RNT toolbox. I have following error ... >> MAIN_createClima Current Settings: setng = nameit: 'medi16' outputdir: './' Return to continue ... (CTRL-C to stop) Creating clim .. - Time index is > 12, please make sure to put in the times vclm_time tclm_time sclm_time ssh_time uclm ...

mashinde

Hi, I have some doubts of following input/output files Input NetCDF file names, [1:Ngrids]. for e.g., !Input files ITLNAME == myapp_itl.nc IRPNAME == myapp_irp.nc IADNAME == myapp_iad.nc FWDNAME == myapp_fwd.nc ADSNAME == myapp_ads.nc so what is this file and how to create it..is it mandatory to use ...

mashinde

Hi...

problem solved......

Thanks !

mashinde

Hi, I am using mpi to build exe file...but it gives following error..... /usr/bin/gfortran -frepack-arrays -g -fbounds-check /home/mahesh/ROMS/src/Projects/rmed16/Build/esmf_roms.o /home/mahesh/ROMS/src/Projects/rmed16/Build/master.o /home/mahesh/ROMS/src/Projects/rmed16/Build/ocean_control.o /home ...

mashinde

Hi frds, I update my svn to 471M, Now i am compiling building bash again but it gives me following error..... .. ... .. /usr/bin/cpp -P -traditional -DMPI -DLINUX -DX86_64 -DGFORTRAN -D'ROOT_DIR="/home/mahesh/ROMS/src"' -DRMED16 -D'HEADER="rmed16.h"' -D'ROMS_HEADER="/home/mahesh/ROMS/src/Projects ...

mashinde

Hi frds, I am RNT toolbox. I have following error ... >> MAIN_createClima Current Settings: setng = nameit: 'medi16' outputdir: './' Return to continue ... (CTRL-C to stop) Creating clim .. - Time index is > 12, please make sure to put in the times vclm_time tclm_time sclm_time ssh_time uclm_time ...

mashinde

Dear frds, I have following two doubts... (A)I am using d_climatology.m and c_climatology.m from roms matlab tools to create clim file. I have created a grid with 1/16 spatial resolution with 72 vertical levels. The geographical area of grid is -20W to 37E and 29N to 47N. I am using R2009b version ...

mashinde

Dear frds, while running the d_climatology.m in matlab r2009b...it gives me following error d_climatology myException = MException Properties: identifier: 'MATLAB:netcdf:endDef:formatSizeConstraintsViolated' message: 'Enddef operation failed, library error message "NetCDF: One or more variable sizes ...

mashinde

Dear frds, I am using all_oa.m and oa_cat.m it appends bottom annual levels to monthly oa fields. please look the attachments. when I run all_oa.m it gives me following error...... ??? Error using ==> NetCDF file: m16_oa4_lev_jan.nc nc_inq at 49 mexnc: ncinquire - error while inquiring file: m16_oa4 ...

mashinde

Dear RubenDiez-Lazaro, I have questions about OA package. I downloaded OA pkg from ROMS website. I have global levitus data in hdat format. Now i want to process OA package and create nc files to my application grid. Do you know how to process this OA pkg, to create monthly and annual clim files ...

mashinde

Hello Kate, Thanks for reply. (1) Actually I am using Levitus98 monthly global files in hdat formats as given in roms packages on roms website. It also contains global coads clim data in hdat format. Now I want to map these data to my application grid. So I can use output nc files that mapped ...

mashinde

dear frds, when I compiled Makefile_GNU from oa pakcag, it creats 'oa' exe file in ${HOME}/bin directory....and when i use this 'oa' exe to extract the global hdat data to my subdomain then it gives me following error..... here is the part of shell script 'xhdat_job_medi ...

mashinde

Dear frds,

the error is gone......but i did not find the exe file 'extr_hdat' which is in the shell script 'xhdat_job'

if anybody knows where i have to look/get this exe file. ?

thank you in advance

mashinde

Dear frds, when i compile make file as make -f Makefile_GNU...........it gives me following error make -f Makefile.Gnu gfortran -DNO_EXIT -g -O2 -o /home/sadmin/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 dynht.o exitus.o filter.o ...

mashinde

Dear All, I use sun ultra 24 workstation...... I have following error after running 'make -f MakeFile.Sun' f77 -c -fast -O3 -u -U -I/usr/local/include oamain.F MAIN oamain: Error on line 41 of <stdin>: Statement order error: declaration after DATA Warning on line 316 of <stdin>: nf_put_varm_int ...

mashinde

Dear kate,

I am going to use ecmwf data for forcing....

regards
mahesh

mashinde

Dear frds.

I want to use OA pkg for creating OA.nc files......
i have following questions..

(1) How to run/execute OA pkg to create nc files?
(2) How 'extr_hdat' exe file created?

Thanks !

mashinde

Dear frds,

does anybody knows how the two boundary files viz. boundary1.cdl and boundary2.cdl is created.

Thanks

Mahesh

mashinde

Dear frds,

how to create boundary file ?

Thanks

regards
Mahesh

mashinde

Dear Frds,

I want to create forcing files for the roms. I read the 'readme' given in the forcing package on wiki roms page. But i do not understand how to make forcings using scripts/programs there.

Thanks

Mahesh

mashinde

Dear frds,

where i get preprocessing and post processing tools in matlab for ROMS

mashinde

Dear frds,

I am running now roms for the spin up for 10 years.
when i checked output of monthly avg.nc files they seems to be identical.

is this correct.?

thanks

mashinde

Dear Kate,

Yes. I also need tidal forcing. where i can get roms tools for matlab for generation of these files.

Thanks

mashinde

Dear frds,

i have general question. I created a grid. Now i want to create clim.nc,forcing.nc,blk.nc,ini.nc files....
can anybody explain how to make it.

thanks in advance

mashinde

Thank you Kate !

It works...

Mahesh

mashinde

i am running build.bash and i have following error cd /home/mahesh/ROMS/src/Projects/medi8/Build; /usr/bin/mpif90 -c -frepack-arrays -O3 -ffast-math get_data.f90 get_data.f90:114.41: & FORCES(ng) % lrflxG(LBi,LBj,1)) 1 Error: 'lrflxg' at (1) is not a member of the 't_forces' structure make ...

mashinde

Dear frds, I have one doubt regarding the make_clim. I have the problem with make_clim. I created 1/16 deg grid over medi. sea having 70 vertical levels. The errors are as below >> make_clim Making the clim: /home/sadmin/Desktop/ROMS/Roms_tools/Run/ROMS_FILES/roms_clm.nc Title: Mediterrenean Test ...

mashinde

I want to the plot bathymetry nc file which i created with seagrid. I am using roms_wilkin's matlab plotting tools for the same. I created a simple script as below... clear; % grd='roms_grd.nc'; file=grd; cmap=zebra(2,64); clev=[100 250 500 1000:1000:4000]; var = nc_varget(grd,'h'); % han=roms_plot ...

mashinde

please see the attached doc file..........

mahesh

mashinde

Dear perezro,

here is my coastline extractor file in mat fomrat

mahesh

mashinde

Dear perezro,

ok.....i have other question...
i want to make a grid of 1/16 deg x 1/16 deg having 70 levels ......so how to make it.....

-mahesh

mashinde

when i run the test application Adria02, then it gives following error

READ_PHYPAR - could not find input file: adria02_grid2.nc

where i get adria02_grid2.nc file....ç

mahesh

mashinde

ok here is my coastline data file.....

mashinde

Dear perezro,

i did as you reply me...here i am sending you figure...please check the arrows...

please see the figure.....

mashinde

hello dear perezro,

i did as your suggestions but it won't work fine. where i have to add the points...

mahesh

mashinde

please see the attached doc file with figure....... mahesh fillseg.doc ------------------------------------------------------------------------ x=axis; >> lon_range=ax(1:2), lon_range = -20 50 >> lat_range=ax(3:4); >> lon_range=ax(1:2); >> lat_range=ax(3:4); >> lon_range=ax(1:2); >> lat_range=ax(3:4 ...

mashinde

i am new to ROMS... (1)I want to make a grid of 1/16 deg x 1/16 deg over Mediterranean sea. how to make it. I tried with seagrid in matlab, but it take so long time and some grids are missing at the coastline. (2) I want to change the parameters theta_s,theta_b and Tcline in benchmark applications ...

mashinde

when i run make file i got the following errors..... ---------------- ./makefile ./makefile: line 31: NEED_VERSION: command not found ./makefile: command substitution: line 34: syntax error: unexpected end of file ./makefile: line 39: sources: command not found ./makefile: line 40: libraries ...

mashinde

when i keep netcdf=ON then following erros coming /usr/bin/gfortran -frepack-arrays -g -fbounds-check Build/esmf_roms.o Build/master.o Build/ocean_control.o Build/ocean_coupler.o Build/propagator.o Build/roms_export.o Build/roms_import.o -o oceanG Build/libNLM.a Build/libANA.a Build/libUTIL.a Build ...

mashinde

this error still persist when i run make -j 4 cd Build; /usr/local//mpich/mpich2-1.1/bin/mpif90 -c -frepack-arrays -g -fbounds-check master.f90 /usr/local//mpich/mpich2-1.1/bin/mpif90 -frepack-arrays -g -fbounds-check Build/esmf_roms.o Build/master.o Build/ocean_control.o Build/ocean_coupler.o Build ...

mashinde

when i run make file ....following error occurs.....even i solve these errors but still this error persists even after i write make clean and make ______________________________________________________________________________________ ROMS/Utility/set_scoord.F > Build/set_scoord.f90 ROMS/Bin/cpp ...

mashinde

i compiled netcdf4 in parallel with hdf5 library.....still i have the mpirun error as follows $mpirun -np 4 ./oceanM ROMS/External/ocean_upwelling.in rank 3 in job 73 ubuntu_50363 caused collective abort of all ranks exit status of rank 3: killed by signal 11 rank 2 in job 73 ubuntu_50363 caused ...

mashinde

i installed new HDF5 and Netcdf4 I also compile it in parallel I/O mode. In serial it work fine. Also some test example with mpich2 runs okay. now i am getting the following error. $mpirun -l -n 4 ./oceanM ROMS/External/ocean_upwelling.in rank 3 in job 18 ubuntu_50363 caused collective abort of all ...

mashinde

i run oceanM file with mpirun -np 4 ./oceanM ocean_upwelling.in, which gives me following error ******************************************************************++++++*******************************++++++++ $ mpirun -np 4 ./oceanM ROMS/External/ocean_upwelling.in [ubuntu:15999] *** Process received ...

mashinde

when i running make second time it gives me following errors......... cd Build; /usr/bin/gfortran -c -frepack-arrays -g -fbounds-check master.f90 /usr/bin/gfortran -frepack-arrays -g -fbounds-check Build/esmf_roms.o Build/master.o Build/ocean_control.o Build/ocean_coupler.o Build/propagator.o Build ...

Search found 135 matches