COHERENS forum

Hi there,

This is a rather general request, I think we need a little more info on what you are actually trying to accomplish before we can help. Did you have a stroll around the website and look at the manual or the COHERENS for dummies document?

Cheers,

Katrijn

Hi there,
Your question is about nesting, so try to read the chapter in the manual and look up the test cases that use nesting (such as the optos, nos test cases).

Hope this helps you get started,

KAtrijn

Dear Sou, I don't use the modfiles to read the discharges in my file, I open the file directly. Here's my code, the data are stored in rivers.loc file. don't forget the USE statement USE inout_routines, ONLY: close_file, close_filepars, open_file, open_filepars CALL open_file(iunit,'rivers.loc','IN'...

Hi Sou, Normally windspeed is listed as a standard output so just adding iarr_uwindatc in the standard 'SUBROUTINE usrdef_tsr_params' should be sufficient: !---varids and ranks tsrvars%ivarid = (/0,iarr_uwindatc,iarr_vwindatc,iarr_windatc,iarr_uvel,iarr_vvel,iarr_temp/) tsrvars%nrank = (/2,2,2,2,3,3...

Hi Sou, there is not a real alternative in V2, each proces has to be processed individually, hope this helps

Hi Sou,

iopt_grid_htype is the grid of the model, which is different than the surface grid. If you want the two to be the same you should enter nhtype = 4 (and not iopt_grid_htyp = 4)

Good luck,

Katrijn

Hi Sou,

You can check the dimensions of your surface grid by 'print *, size(surdata,1), size(surdata,2).

By default (nhtype = 0 in Usrdef_Model) the surface grid has 1 single grid point, if you want it to be equal to the model grid then set nhtype = 4

Good luck

Hi Sue,

Its not completely clear to me what you mean, I sent you a private message.

Cheers,

Kat

Hi Sou,

What exactly do you mean with periodic? What kind of boundary condition did you have in mind? Because inducing a periodic condition is not that difficult if you include time and/or a wave function in the program setup?

Hi Sou,

I do notice a few things:
east boundary is commented (probably you testing the code)
At the north boundary you reset jj to 0, this should not be done (iobv(1:20)), if you reset you don't fill the whole array

Hope this helps,

KAt