Ocean Modeling Discussion

c.drinkorn

@John: Thanks! I am trying to make use of this now but can't seem to get around what av2 is supposed to be. Am I lacking a particular set of functions here? @Kate: We set up an Arctic/N.Atlantic sed-ocean-ice model for several sensitivity tests. E:g. our first paper: https://doi.org/10.1016/j.palaeo ...

c.drinkorn

I think I found it! Is it gridpak?

c.drinkorn

Thanks Kate! I am not much of a python expert (on my list!) but I worked my way through the installation tutorials for a local version of the grid tool. However, when I try to load A20 into the tool I get an "unloadable" error. Maybe I should also look into the "old old" grid generation tools ...

c.drinkorn

I didn't make the grid myself, it's the one from the former metno Arctic20km model. The original grid file has no psi coordinates. How can I extrapolate them from the rho coordinates? Also: I have experienced difficulties making a grid of my own with the North Pole in it (that's why I re-use ...

c.drinkorn

Hi all! I am currently trying to use conservative remapping of ROMS data instead of my previous bilinear method. For conservation the cell areas are needed and therefore the cell corners. I thought those were lon_psi and lat_psi, respectively. But the remapping routine complains about missing corner ...

c.drinkorn

HI Osean, our application has a vertical grid like the one you are aiming at. We achieved this with the vertical transformation formulation 2 and the vertical stretching formulation 3 with theta_b = 1.5 and theta_s = 1.0. Our vertical resolution ranges from 0.3m in places at the bottom and surface ...

c.drinkorn

Hi Neel! It seems like the string comparison never finds a match. So I guess there is something wrong with the reading of the constituents from the file. In line 486 of TPXO2ROMS_v5pt2.m you can see that the harmonicNames are read from the variable "con" in the elevation file. Can you check if the ...

c.drinkorn

Hi Neel, I think the error comes from the line 617 and following. I forgot to mention that you also have to change the path to your grid file here or remove mine if your grid file is in the Matlab working folder: elseif strcmp(TPXOfile_grid,'grid_tpxo9_atlas_v4.nc'); Then it should work. :) Cheers ...

c.drinkorn

Hi Neel and John the attached file is my Matlab function to read TPXO atlas data. I adapted it to the current TPXO9 version on the basis of John's TPXO2ROMS_v4 script. It works fine for me. It can read in several constituents after each other on the basis of individual TPXO files or just one, such ...

c.drinkorn

Hi Neel, I am very sorry I couldn't help you solve the problems. I suspect TPXO2ROMS_v5x is not compatible with TPXO9_atlas data but maybe I am wrong. TPXO9_atlas data has a very different structure than TPXO8 and even if John Luick writes in the comments of the function that there were a lot of ...

c.drinkorn

Hi Nee, yes, I know now I had to add the variable to the function readTPXOdata. Inside of the function file the code also calls subfunctions which need necessary variables passed on to them. readTPXOdata doesn't know ROMSnames in your file. You have to add it as a function variable: function TPXO ...

c.drinkorn

Hi Neel, I think I know what the issue is. The dimension ncread complains about is nH which is the number of constituents. You are only using one of them, Q1 as I can see. In your file the variables such as h have only two dimensions. I think in TPXO files with more than one constituent they have ...

c.drinkorn

Hi Neel, yes, it's basically the same problem as above, just now for the u-coordinate latitudes and longitudes and will be the same for the v-coordinate lat/lon, too. So I suggest you add the meshgrid lines there, too: Xi=ncread(TPXOfile.grid,'lon_u'); Yi=ncread(TPXOfile.grid,'lat_u'); [X,Y ...

c.drinkorn

Hi Neel, I really want to help you so I took another look at your function file. I know that John Luick is the absolute expert here but maybe this could help: Before the critical line (line 421 is the line where the depth is read from the TPXO grid file) the variables X and Y are determined on the ...

c.drinkorn

Hi Neel, the script is the Matlab program you should run from the command line or a bash script. There you can change the paths to your ROMS grid file, the output file name and location and the desired constituents as well as information about the time settings for your tides. This script then hands ...

c.drinkorn

Hi everyone,

yes I think sonaljit.py would be a more suitable Username.

I think, if you chose the Option AVERAGES_DETIDED and set some of the detided variables True for output, ROMS will make a harmonics file ocean_har.nc which contains the constituents, am I right?

Merry Christmas!

c.drinkorn

Hi Alfredo and John the initial excessive removal of sediments in my model is physical: it is due to a deep sea water fall inducing strong erosion at its bottom (they correlate with high vertical local velocities). I could stabilize this erosion a bit by shortening consolidation times and increasing ...

c.drinkorn

So, I could trace back the reason for this error and find a rudimentary workaround for it: The problem in sed_bed_cohesive is that the catch for the case when the available sediment is not enough for the amount needed to add is not taking care of the variable Ksed. If all the sediments are not ...

c.drinkorn

Yes, of course.

But now that newlayer_thick is larger than zero the layers are all exactly that value.

c.drinkorn

That was a really good idea! The result is stunnning: Of all 30 layers à 30cm thickness, the lower 29 have a thickness of zero at the first time step already! The top layer maintains a thickness of several meters (EDIT: I just realize it now after I have read it another time - several meters!? OK, ...

c.drinkorn

UPDATE: Ok, so I ran into the same error with layer thicknesses of 30cm, i.e. 30 layers à 30cm which makes 9m of sediments. I doubt that it was possible to erode this within one time step... It seems as it is rather a numerical problem instead that may have been overlooked because it is ignored ...

c.drinkorn

HI, yes, the problem only occurs in Debug mode. In normal Mode the model runs into stability problems after some time. Unfortunately, I can not look into what happens with the bed layers because the error occurs already at the first time step so that both rst and dia are empty. In the first his file ...

c.drinkorn

Hi all, I am experiencing a strange error regarding the number of bed layers in a stratified cohesive sediment application: Subscript #3 of the array BED has value 31 which is greater than the upper bound of 30 The error is traced back to: coawstG 0000000001B1CFCF sed_bed_cohesive_ 581 sed_bed ...

c.drinkorn

Hi Sediment Modelers, I happen to be very confused about the units of suspended material in the water column, and deposited or eroded sediments. I want to capture how much material I lose via the boundaries. In order to do so I integrate deposition (positive change in bed thickness) and erosion ...

c.drinkorn

Hi Kate, thanks for the suggestion on adding a bad salinity alarm to diag! I did so and could therefore make the model stop running already when salt goes negative. Looking at the fields in rst then made it obvious that something suddenly goes wrong with the atmosphere-ocean saltflux. After looking ...

c.drinkorn

Hi, I am trying to solve an issue with salinity in my COAWST branch. It only appears with the ICE MODEL switched on. I get extreme (negative) salinity values. This also causes ice_frazil to prompt a trouble message because gamma gets negative: trouble in ice_frazil 316 197 31 0.927691919473818 -2326 ...

c.drinkorn

Hi all, I think I found a typo in sediment_inp.h where the sediment tracer diffusivity is read from sediment.in: After CASE ('MUD_TNU4') Npts=load_r(Nval, Rval, NCS, Ngrids, Rmud) DO ng=1,Ngrids DO itrc=1,NCS i=idsed(itrc) nl_tnu4(i,ng)=Rmud(itrc,ng) END DO END DO stores the value of MUD_TNU4 into ...

c.drinkorn

Hi folks,

I wonder where that spindown time formula from the documentation of the horizontal mixing comes from?

Thank you!

c.drinkorn

Okay, looking at the f90 code made it much easier to see what is happening. It also makes it very obvious that the code is not written to work with only ICE_MODEL (and ICE_BULK_FLUXES) switched on, e.g. if it comes to the calculation of the altered heat fluxes due to ice in the bulk_flux routine ...

c.drinkorn

Okay, since I want the same nudging time scales everywhere I just hard coded the number into ana_nudgecoef. Maybe for a more complex nudging time scales field I'd rather provide a file with the coefficients field. I can see that only switching on ICE_MODEL already takes the ice concentration into ...

c.drinkorn

Hi Kate,

thank you! Things are much clearer to me know. I am trying to nudge to a sea ice concentration climatology. I hope I got it right and ROMS will use the M2 nudging time scales for sic?

c.drinkorn

Thank you Kate and Hernan for your information on the sea ice model status! I have to clarify (since this was an old post of mine, it probably got lost up there...) that this was about the COAWST branch. But it would be good to hear respective information on this one, too. Is COAWST supposed to have ...

c.drinkorn

Hi all, it seems like I finally solved this one: In the makefile the link to the SeaIce folder in Nonlinear was (unintentionally?) missing. When I added this one the compiling error vanished. There was a warning that all old commands for the respective targets are being overridden by the new ...

c.drinkorn

I'm sorry, I meant nl_ocean not mod_ocean.

c.drinkorn

Hi John, thank you very much for the clarification! Yes, grep is my best friend. :) So I was somewhat familiar with the routines you pointed at. What's still unclear to me is from where t3dbc_im.F receives its Tobc_out which I suppose is already rescaled at that time. I guess, the thread goes back ...

c.drinkorn

HI all, to get back at this and be fair to the ROMS code I want to share that I investigated the *nudg and *obc_out / _in values provided in the output and they are very much consistent with the expected internal rescaling from days to seconds and the obc factor, respectively. :oops: However, I ...

c.drinkorn

DEar ROMS community, I am pondering over this for quite some time now. Maybe the answer is fairly easy but I couldn't quite figure it on my own: How does ROMS calculate the RadNud boundary nudging time scales? Does it check the time unit of the boundary file in order to determine whether the ...

c.drinkorn

Thank you, Kate and John! Yes, at some point I also want to calculate the transport across a transect and I already wrote a script for that. I will compare the procedure carried out in Oceanspy with mine! For now I am using the volume transport to depict a mean vertically integrated flow vector ...

c.drinkorn

Hi John, these are really helpful explanations! Thank you! Just one more confusion: When I calculate the vertically integrated volume transport on the basis of the barotropic flow components, I can make use of the already interpolated and shifted eastward and northward components. However, then ...

c.drinkorn

Hi John, thank you for the very comprehensive explanation! To be honest, many of those average flow variables confused me. I was looking for a variable to create a vertically integrated flow field, preferably in m3/s or Sv, respectively. I post-processed Huon and Hvom using NCL but I don't trust the ...

c.drinkorn

Hi Kate, thank you for your helpful advice! I want Huon and Hvom at rho-point locations (that is in the 3D sense) and then rotate them towards East or North, respectively. I am still a bit confused as to why on_u is used for the calculation of Huon and not on_p. Wouldn't the distance of the psi ...

c.drinkorn

Hi ROMS cummunity! I am currently trying to make ROMS output the volume transport components as eastward and northward, too, as is being done for the 2D and 3D velocity fields. I was wondering, since the volume transports are Average out variables, why the rotation of u and v is only being conducted ...

c.drinkorn

Hi Charles, I attached the grid file and the topography (remapped to the grid file's grid). I included your fix into getDefaultBathymetry but unfortunately this didn't solve the problem so I still had to add the lines for axlim(1) and axlim(2). I appreciate if you try and find a solution for the ...

c.drinkorn

Okay, I managed to load my grid file into gridBuilder by adding these two lines to the getDefaultBathymetry function: if axlim(1)<-180;axlim(1)=-180;end if axlim(2)>180;axlim(2)=180;end By remapping my new topography to the old grid file grid and writing a simple script to load lon_rho, lat_rho and ...

c.drinkorn

Hi Charles,

I sent you an email with the files. In the meanwhile, I will try myself to find a solution, too.

Cheers,
Cate

c.drinkorn

Probably due to case-insensitivity, there is a compiling problem with the variable dP in the pressure gradient routines for quartic pressure Jacobians. The variable conflicts with the parameter DP for double precision from the mod_kinds module. Renaming dP to say dPr solved the problem. I know that ...

c.drinkorn

Hi all, I am trying to "replace" the topography in an existing ROMS grid file. But when I load the grid file into GridBuilder I get the error: Error using griddedInterpolant Grid is improperly defined. Error in setROMSgrid (line 59) BathyInterpolant=griddedInterpolant(bathymetry.xbathy,bathymetry ...

c.drinkorn

HI all, I have been looking into this issue for quite some time now and wrote about pre-steps of it in other threads. At this point I am not getting any further and desperately hope for help by the community. Here is a short outline: I am running a high-resolution (delta x ~ 20km) Arctic application ...

c.drinkorn

I got another step further: Only undefining SED_MORPH makes the restart blow up disappear. But since I want a time-evolving bathymetry I will investigate this further and report. SED_MORPH is meant to work with restart, right?

c.drinkorn

Ok, I got one step further but still didn't succeed. I figured the problem with the additional restart time step was nonsense. The additional time step was simply the one written after blow up. :roll: However, the blowing up still remains and without the sediment option restart works fine. I now try ...

c.drinkorn

Hi Kate, all of my forcing and boundary files range from Jan-01 00:00:00 this to Jan-01 00:00:00 of the next year. So this should be fine. I realized that all sediment related variables do have values in the additional rst fields. Now that I switched Lcycle back to false I get 367 rst fields of ...

c.drinkorn

As for MIX_S_TS, it's a really bad idea over a sloping bottom and can lead to unphysical cross-isopycnal mixing. I hear you! But there was no way to solve the issue of punctual alternating high/low temperatures at the slope when mixing along z-planes which would even cause the model to blow up. The ...

c.drinkorn

Then the only idea I have is to enable TS_MPDATA which will "force" the fields to be in 3d consistency...

c.drinkorn

Hi, I found this old issue while I am having had the same. I realized it is connected to the usage of PERFECT RESTART. In wrt_rst.F, when choosing PERFECT_RESTART, t2dgrd is becoming an extra value, i.e. DimIDs(31). Later in def_var.F this causes an inconsistency while nvd3 has the value 3 and Vdim ...

c.drinkorn

Hi Xavier,

I switched from MIX_GEO_TS to MIX_S_TS in order to solve a similar problem at the bottom. Maybe it also helps you.

//Cate

c.drinkorn

Hi Azam, just an idea: Did you try switching between Laplacian and biharmonic horizontal diffusion? I.e., define either TS_DIF2 or TS_DIF4. I don't know whether the issue is too less horizontal or vertical mixing but I guess if your subsurface temperature AND salinity is too low you need more ...

c.drinkorn

Hi all, I have a simple question: Why does ROMS store the last time step of the restart file one hour after the model time domain instead of the very last time step so I can use it as a restart initial? My model time step is 600s and I am simulating 52560 time steps (one year). NRST is set to 144 in ...

c.drinkorn

Hi Van! I think, you need to specify the path to your MPI library. I marked the line with a long arrow: gfortran ) if [ "${which_MPI}" = "mpich2" ]; then export PATH=/usr/bin:$PATH elif [ "${which_MPI}" = "openmpi" ]; then ------> export PATH=/opt/openmpi:$PATH fi

c.drinkorn

Hi all, I recently switched to the COAWST because I want to run both ice and sediments in an Arctic domain but couldn't solve some inconsistencies in the Kate branch code... Howver, now I am unable to solve this compiling error which seems top be about some ice routines: /scratch/b/b380636 ...

c.drinkorn

Unfortunately, I am reproducing the error with one of your ocean.in files. I took my time to look into the referred code lines but there is nothing suspicious or inconsistent. I guess, I'll switch to the COAWST code, then, unless anyone has an idea what produces this problem.

c.drinkorn

If you are using my code, you really, really need to start from one of my ocean.in files. In particular, you need not only NFFILES for the number of forcing files, but also NBCFILES and NCLMFILES for the number of boundary and climatology files as well. Maybe even NSSFFILES for the rivers if the ...

c.drinkorn

I am trying to set up Kate's branch for an Arctic application including ice and sediments. When I first used the code as I received it from github, there was a compilation error about the dimension of Rval in inp_par.F Obviously, there is a major difference between the SVN version and Kate's branch ...

c.drinkorn

Maybe look at the rho, temp and sal values right before blow up.

c.drinkorn

Yessss, that was it! I didn't look in the Qout section.

Thank you!

c.drinkorn

In my restart files I get mudmass written out so in principle this works. Just not in the history file...

c.drinkorn

Ok, here is my sediment.in

c.drinkorn

It does define mud as sediment class but not mention bed mass as Hout output: Sediment Parameters, Grid: 01 ============================= Size Sd50 Csed Srho Wsed Erate poros Class (mm) (kg/m3) (kg/m3) (mm/s) (kg/m2/s) (nondim) 1 0.0000E+00 0.0000E+00 1.8000E+03 1.0000E-01 5.0000E-04 9.0000E-01 tau ...

c.drinkorn

MUd, so NCS equals 1. I thought I found the reason but still no writing of bed mass: In the Write-out section of sediment_inp.h it says ... DO itrc=1,NST i=idfrac(itrc) IF (Hout(i,ng)) WRITE (out,160) Hout(i,ng), & & 'Hout(idfrac)', & & 'Write out bed fraction, sediment ', itrc, & & TRIM(Vname(1,i ...

c.drinkorn

Mysterious. I think, I have everything set as it should: ... Hout(idmud) == T ! mud_01, ... suspended concentration Hout(iMfrac) == T ! mudfrac_01, ... bed layer fraction Hout(iMmass) == T ! mudmass_01, ... bed layer mass Hout(iMUbld) == F ! bedload_Umud_01, ... bed load at U-points Hout(iMVbld ...

c.drinkorn

Something else: Is there anything I need to switch on in order to write out bed mass? Setting the true flag in sediment.in didn't do the trick.

c.drinkorn

Oh and by the way, I also needed to undefine the vertical tracer advection scheme, of course. Just in case someone is reading this later...

c.drinkorn

MPDATA did really solve the problem which does of course make sense because it just restricts the local simulation results. Even though I am always a bit concerned with too much correction terms in the calculation, I think switching to MPDATA was a really good idea. Thank you for this great advice ...

c.drinkorn

Thanks for your promt answer! Speaking of sharp horizontal gradients: I do indeed discover erroneous negative rho values (in fact of the same order as the susp. load concentration) next to the river discharge, so maybe this isn't exactly a sediment issue. T and S look fine, though. The errors occur ...

c.drinkorn

Hi, I get very low negative suspended load concentration values (down to several hundreds of kg/m^3) near river runoff. I have no sediment input via rivers. But the initial concentration in the water column is zero so I wonder what exactly a negative concentration should mean. Is there a way to ...

c.drinkorn

Thank you for your explanation, Kate. The error actually vanished with the recent updates so this might have just been a bug... I still don't use parallel I/O because I keep receiving an error about the scratch arrays in the global periodic patch. However, I don't mind. So far using a serial I/O is ...

c.drinkorn

The user has an Arctic domain. I added a GLOBAL_PERIODIC flag to ROMS to handle just this case, interpolating across the dateline. I've attached the patch file. Hi Kate, I'm sorry I have to get back at this again but I just wanted to switch to Parallel IO and got an error while compiling about too ...

c.drinkorn

It seems to be possible by using a bottom boundary layer model...

I will try something like that.
I'd be happy about any further hints.

c.drinkorn

The sediment part in roms can be run with just roms. The swan coupling provides wave information that can be used by the ocean model to: - compute enhanced surface mixing due to tke wave breaking - computed wave driven flow using vortex force approach - compute enhanced bottom stresses due to ...

c.drinkorn

Hi all,

I wonder if and how it would be possible to start a sedimentation simulation with initial sediment layers coming with e.g. thickness and densities, by running just the sediment module in ROMS. Or will I have to run a coupled SWAN simulation in order to do that?

c.drinkorn

Hi Scarlett, I had sort of the same problem (my Pair time record had switched 12h and 18h steps - it was what the conversion from grib to netcdf resulted in while my grib file had both analytical and forecast). Since I didn't retrieve the data myself I wrote a workaround to fix the field myself ...

c.drinkorn

Hi all, this isn't a big deal but when running a simulation using time evolving bathymetry, I received an error about a variable "ibath" from get_state. I am convinced it should be "idbath" in that line. Fixing this solved the error for me and from the code's logic this also makes sense. :wink:

c.drinkorn

try using this in the project.h file # define ANA_SMFLUX # define ANA_BPFLUX # define ANA_BSFLUX # define ANA_BTFLUX # define ANA_SPFLUX # define ANA_SRFLUX # define ANA_SSFLUX # define ANA_STFLUX Thanks! Switching on ANA_SPFLUX was the solution. I have an open vertical boundary so I didn't even ...

c.drinkorn

I am trying to set up a very simple first test case for the implementation of the sediment tool into my ROMS model. I want a 1m layer of sand with the typical sand properties (high porosity, no cohesion etc.) evenly distributed throughout the domain with closed borders and and no suspension in the ...

c.drinkorn

I realized lately that global_periodic gives a seam at lon 180/-180 if the coordinates are not in the order of -180 to 180 instead of 0 to -0 (which would be the outcome of the d_ecmwf2roms.m). Therefore, I added two lines to the script: In the section where lat and lon are read, right after the ...

c.drinkorn

Yes, my inputs files (the files that I downloaded from ECMWF) are monthly, and it was easier for me to modify the script to read them and create them monthly. I think that the better solution is make an annual file and run the code for the entire year once, like you suggest. If you could give me ...

c.drinkorn

HI Scarlett, I worked with the ecmwf2roms script a lot recently. How do you loop through the monthly input files? Maybe you accidentally call the same file all the time (i.e. January)? Just a thought... ADD: I looked into the ecmwf2roms file you shared in the other post and can see that you run the ...

c.drinkorn

If it is your first application in ROMS, I suggest that you gain experience with a smaller regional grid, so you can transfer all that knowledge when building and setting up your global application. Dear Hernan, your are right, I am a bloody beginner. :oops: But thanks to the comprehensive ...

c.drinkorn

Dear Kate, million thanks for the patch! I needed to adjust my Matlab processing a bit so that my lon/lat variables are 1d. In the Matlab package one needs to change the metadata (in roms_metadata.m) to just one variable dimension entry and in the main program comment out the repmat procedure ...

c.drinkorn

Hey all, after your hints about double-entries in the .in file I finally looked really really close once again and, alas, I found an entire block of entries at the very end of my file (after the Glossary and coupling section). My guess how it ended up there: an accidental mousewheel click while ...

c.drinkorn

Oh yes, tabs! I encountered them already in the .in file. :D I use Vi on a UNIX. And .e.g tabs before the forcing file links are a bad idea! :twisted: However, so far I couldn't find them to be the reason of double-allocation. At the moment I am still busy setting up the forcing fields to the ...

c.drinkorn

Considering the recent repository update underlines even more that the IF-condition wrap is a workaround but shouldn't be necessary. In fact, I did not encounter the problem when using another ocean.in. I tried hard to figure which setting causes the obsolete allocation but was not successful. Which ...

c.drinkorn

Thanks for your reply, Kate. No, there is no double entry for NtileJ. What I did now is to wrap the allocation of IO_BOUNDS inside an IF-condition, too. I guess, this is a nasty workaround and there sure is an issue in my model setup somewhere but it solved the problem for now. Let's see how the ...

c.drinkorn

After setting up my model once again, I re-encounter this problem.

No suggestions on how to solve this?

c.drinkorn

Hello everyone, while trying to run ROMS in parallel mode with analytical initial conditions I receive the following error output from debugging mode (excerpt from the whole error file): forrtl: severe (151): allocatable array is already allocated Image PC Routine Line Source oceanG 000000000220D767 ...

c.drinkorn

rm -f -r core *.ipo /home/yaduvendra/projects/upwelling/Build /home/yaduvendra/make_macros.mk ./build.bash: line 409: syntax error: unexpected end of file plz help This is some time ago but I just had the same error and fixed it. Reason and solution were rather trivial: There was a non-closed if ...

c.drinkorn

kate wrote:Paul Budgell is doing a global ROMS. On the small end, I think it has been run to simulate a laboratory rotating tank to the same scale. Four meters across? Something like that.

I know, this is old but is the global set up still available somewhere (GitHub etc.)? Thanks!

c.drinkorn

Hi Shauryasingh!

Do you use NetCDf4? (Usually, it is Kate's part to ask that question...)
If not, you need to set the flags to the different libraries manually (e.g. in the compiler file).
If yes, there might be something wrong with the path to your netcdf library.

c.drinkorn

It works! :D After re-compiling and re-running and re-compiling and re-running using Netcdf4 or not using it... I finally figured what the problem was: The path to the grid nc (and apparently all the other input nc files) was wrong. It had to be the absolute path, no home directory shortcuts (or at ...

c.drinkorn

Ok, I found the section in mod_netcdf.F and it looks like this now: IF (FoundError(status, nf90_noerr, __LINE__, & & __FILE__ // & & ", netcdf_open")) THEN WRITE (stdout,10) TRIM(ncname), TRIM(SourceFile) PRINT *, trim(nf90_strerror(status)) exit_flag=3 ioerror=status END IF There was no "IF (Master ...

c.drinkorn

Thank you, Kate! But where will I have to include it?

c.drinkorn

Yes, it is there and I can open it with ncview. All variables are in it and it seems alright.

c.drinkorn

Hi everyone! I am sure this is trivial but I am lacking experience: When trying to run the DAMEE_4 test case I receive the following error in the output file INITIAL: Configuring and initializing forward nonlinear model ... ******* Found Error: 02 Line: 6565 Source: ROMS/Modules/mod_netcdf.F, netcdf ...

c.drinkorn

Yes, I know. But I did comment out USE_DEBUG=on...
However, I found that in the makefile it was not. When I did it I got an oceanM. Happy end.

Thank you for your help, Kate!

c.drinkorn

They are only warnings. Did it build an executable? Can you try running it? Also, you want to make sure netcdf-fortran was compiled with the same compiler you are using to compile ROMS. Always. I missed this one... Ok, well it built an oceanG. But I was expecting an oceanM. :| Regarding the netcdf ...

c.drinkorn

I just solved it myself. I saw that you have a much newer openmpi so I looked again and found openmpi-1.10.7_hpcx-intel14. I thought it would be a wrong one since it comes with the hpcx suffix but the compiling worked without error. So thank you again! :D HOWEVER, I only have an oceanG and no oceanM ...

c.drinkorn

Hi Kate! Your hint was very valuable, I managed to find the module which provides the openmpi mpif90 environment for intel! Thank you! :D (For users having the same issue: typing something like "module whatis" gives additonal descriptions of the modules and makes it easier to identify the right one ...

c.drinkorn

Hi everyone, after having tried to find the source of error myself during the whole day, I resign and ask for help here: I want to compile a test case for parallel processing and use openmpi and ifort. I already switched on USE_MPI, of course, and also pointed to mpiifort instead of mpif90 in the ...

c.drinkorn

Hi, if you do not want to use NetCDF4, it may be necessary to change the library path settings in the compiler file. In the respective file (according to your system and fortran compiler) in the 'Library locations' section you may paste the output of "path/to/your/nf-config -flibs" into the "ifdef ...

c.drinkorn

Thank you very much for your reply, Kate! And sorry for the late answer; I just found the time to take a closer look at your suggested domain very recently. I also found a GitHub (yours?) with ARC-NATL header files so I'll take a look into them in order to examine the resolution and extent. From ...

c.drinkorn

Hi everyone!

I am new to ROMS and this forum.
While I am digging into how to set up a high-resolution model for the above mentioned domain, I thought it was worthwhile to ask around who might have done a similar set up already. I would be very grateful for exchange or hints!

Cate

Search found 110 matches