Ocean Modeling Discussion

Rocky_Wang

Hi everyone, I am using the ROMS in the COAWST model to simulate the temperature of some parts of the South China Sea, but I find that the temperature of the model is dropping rapidly. Figure 1 is the sea surface temperature of 195 days in a year of the HYCOM model, and Figure 2 is the simulation ...

Rocky_Wang

Hello everyone, Recently, I wanted to use Pyroms to more easily download and handle HYCOM data, but I encountered problems during installation and could not install Pyroms properly. My installation steps are as follows: step1 create a python environment by conda: $conda create -n model2roms python=3 ...

Rocky_Wang

to the Compilers/Linux-ifort.mk, you will need to add ifdef USE_SWAN FFLAGS += -assume byterecl endif i would suggest that if you intend on doing a lot of research with swan+roms, that you switch to coawst. We have updated swan to more recent versions, include vortex force wave effects on currents ...

Rocky_Wang

Hi John Thanks for your quick reply, It's a good idea to test the Sandy case with initial file. I run the sandy case with three data files : Sandy_roms_grid.nc, Sandy_ini.nc and tide_forc_Sandy.nc, both of them were from COAWST/Projects/Sandy and just the ROMS model was turned on. At this condition ...

Rocky_Wang

Hello everyone, I usually use the COAWST model to couple the ROMS and SWAN and simulate sediment transport. In recent days, I try to install COAWST model in a new server. I download the recent version of the COAWST and successfully run the UPWELLING case, the coupled Inlet_Test case. However, when I ...

Rocky_Wang

Your stdout says " 17974 62 09:50:00 9.874450E-02 2.791750E+02 2.792738E+02 2.775544E+13 (095,245,20) 0.000000E+00 1.108514E-02 8.481048E+00 6.956400E+01 DIAG speed trouble 69.5640000538570 Blowing-up: Saving latest model state into RESTART file DIAG speed ijk 98.3770376595520 95 246 1 1" So you ...

Rocky_Wang

Thanks John. Sorry for the late reply, I attended a field survey in recent days. My model is calculated from 2007 to 2016, It always Blowing-up in 2007/3/04. So I have made some test to find the reason. For example, use the 2016 wind file instead of the 2007 wind file for swan Model. I find only in ...

Rocky_Wang

Hello everyone, I used to run the ROMS model with SWAN output files(wave-orbital velocity ub,wave period T, and wave-propagation direction,) to investigate sediment transported processes dominated by currents and waves. now I want to use COAWST model to couple them together. The C-preprocessing ...

Rocky_Wang

It runs successfully, but I think the result may be not what I want, Because I changed so much of my model. the model can't output u and v, just ubar and vbar. There is no need to run in barotropic way. I want to know if there are app opinions that can keep the temp and salt constant, even though I ...

Rocky_Wang

Thank you Kate, that is to say, when I turn off SOLVE3D, I can't use SEDIMENT and BULK_FLUXES anymore. If I want to add wind, I must change wind speed into wind stress in my forcing files?

Rocky_Wang

Hello everyone, I am going to make some tests that how some factors such as tide, wind, temp and salt influence the residual current. I know how to run the model without tidal and wind forcing. But I don't know how to run the model with uniform and unchanged temp and salt that without baroclinic ...

Rocky_Wang

Hello Kate, I use the BULK_FLUXES in my model. My CPPDEFS options are : # define BULK_FLUXES # define ANA_BTFLUX # define ANA_BSFLUX # define ATM_PRESS # define EMINUSP # define SOLAR_SOURCE and I download data from ERA-interim to make my forcing files. My variables are: lwrad_down % downward ...

Rocky_Wang

Thank you Kate! It works!

Rocky_Wang

Hello everyone, I usually use the Mellor–Yamada Level 2.5 scheme in my model. Bnd when I try to use the GLS turbulence closure today, It always blowing-up after hundreds of time-steps. And I can not find the reason. I think the grid that I use is normal(I use this grid file for a longtime). I try to ...

Rocky_Wang

Ok, I get the point.Thanks John

Rocky_Wang

Thanks John, That is to say, If I want ROMS one way coupled with swan (make the swan output hwave, tmbot wdir into a wave forcing file) , the CPP option define SSW_CALC_UB can be activated. And if I want to ROMS two way coupled with SWAN, There is no need to activate it. Right? Yes, I want to look ...

Rocky_Wang

Hello, I use the COAWST model to couple SWAN and ROMS, and they exchange variables throuth MCT every 3 hours. I use the SSW_bbl for the wave-current interactions.the CPP options are as following: #define SSW_BBL /* Sherwood et al. BBL closure */ #ifdef SSW_BBL # define SSW_CALC_ZNOT /* Computing ...

Rocky_Wang

If the dmix_offset dmix_slope and dmix_time in initial file are used to describe the initial critical shear stress profile; and the MUD_TAUCR_SLOPE MUD_TAUCR_OFF and MUD_TAUCR_TIME in sediment.in are used to describe critical shear stress equilibrium profile?

Rocky_Wang

Hello Everyone， I have some small questions when use COAWST model to design a mixed bed. First,I want to konw the what is the unit of the parameter MUD_TAUCR_TIME. Is it days or time steps? MUD_TAUCR_TIME == 28800.0d0 ! Tau_crit consolidation rate and the Swelling time scale is order 100 times s ...

Rocky_Wang

Thanks， I use the 'set_depth.m' to get all depths of the grid and extract the point data according to I J positions of the station.

Rocky_Wang

hello everyone;
I want to know if I can use the station.in to get the grid depths of every layers at the station. Thanks

Rocky_Wang

One way could be to decrease the value of Zob in ocean.in. Zos is the surface roughness length, and you wouldn't need to tweak that here. Another would be to decrease rdrg2, also in the ocean.in you show. And another way could be to prescribe a constant, lower, bottom drag coefficient to increase ...

Rocky_Wang

Thank you Kate! that is to say the latent heat flux is calculated by variables that I provided above, and there is also no need to creat a file that contain latent heat flux. Right?

Rocky_Wang

Dear roms users, I'm new user of ROMS model and I'm having some problems about the use of BULK_FLUXES. My CPPDEFS options are : # define BULK_FLUXES # define ANA_BTFLUX # define ANA_BSFLUX # define ATM_PRESS # define EMINUSP # define SOLAR_SOURCE and I download data from ERA-interim to make my ...

Rocky_Wang

Hello everyone, I have a confusion about TIDE_START, I am going to develop my case followling 'create_sandy_application.m'. In ocean_sandy.in file, some time parameters setting are : DSTART = 56228.5d0 56228.5d0 ! days TIDE_START = 52866.0d0 ! days TIME_REF = 18581117.0d0 ! yyyymmdd.dd Sandy case is ...

Rocky_Wang

Hello every, I try to couple swan and roms using COAWST model, to simulate how wave and tide dominate the sediment transport. I do this work folloeing the Inlet_Test. I made a 3 hour wind file to force the swan model, at the same time roms model only forced by tide. When I run them separately, both ...

Rocky_Wang

Hello, I run the roms model with only tide forcing. We have put an ADCP into the modeled region to record the tide current velocity. When I compare the computed current velocity with observed data,I find the simulated result is smaller. which method or parameter I should adjust to increse the ...

Rocky_Wang

if you are using roms+swan coupling, you really should use our coawst system. wee have updated the swan version, implemented grid refinement in swan, have vortex force, etc. we have just been updating the coupling in the coawst system, and staying very close to the roms distribution. send me an ...

Rocky_Wang

Hi everyone:
When I couple ROMS with SWAN ,they exchange data every six hours, but I found my wave bottom orbital velocity that ROMS imported are always zero, Is there something wrong with my swan.in? The attachment is my .h file,ocean.in file and my swan.in file?
Thanks very much!!!

Rocky_Wang

Thanks very much John, I sloved the problem. change nothing and just add what you said. I sent an application letter about COAWST to your E-mail Address, MY name is Rocky. Thanks again!

Rocky_Wang

Sorry, I forget to say that the ROMS I use is Rutgers version of swan coupling.

Rocky_Wang

Thanks John, I use the Linux-ifort.mk, Could you mind show me how to add the flags?

Rocky_Wang

Hi everyone, Recently, in order to couple the SWAN model, I have run the inlet-test case successfully, but when I use matlab to check the output files of SWAN, such as hsig.mat, fric.mat and so on , I found cannot open them. The attachment are my .in , log ,print file and hsig.mat. Can you tell me ...

Rocky_Wang

Thanks very much, I have found the reason that I didn't turn on the #export USE_MY_LIBS=on # use my library paths below,and did't change the path in my compiler file.

Rocky_Wang

this is my log file, and permission.

Rocky_Wang

Thanks your answer, and I checked my path, you mean I did't compile the MCT rightly?

Rocky_Wang

Hi, I tryed to compile the inlet (swan-roms coupling) test. first I compile the Model Coupling Toolkit which generated lib and include files. second I add the lib and include paths in my build.bash export MCT_INCDIR=/home/chenghao/soft/include export MCT_LIBDIR=/home/chenghao/soft/lib third I run ...

Rocky_Wang

Thanks your detailed description, I got the key point. Thanks!!!

Rocky_Wang

Thank you wilkin,Here is the situation, I didn't prepare boundary conditions files, I have a grid.nc and a tidal-forceing.nc, just run the model with only tidal force. In this case, I can't change the *h file into this? /* %%%%%%%%%%%%%%%%%% tide %%%%%%%%%%%%%%%%%*/ # define SSH_TIDES # define UV ...

Rocky_Wang

I change "define ADD_M2OBC" into "define ANA_M2OBC" in my *h file and the model can run successfully

Rocky_Wang

Thanks all your kind answer and suggestion, I got a feedback about CHECK_MULTIFILE, could you mind tell me how to deal with the initialization time.Attachment are my *.h and *.in files, if my cpp options are wrong.

Rocky_Wang

I set '#undef SOLVE3D' in my *.h file, and set Time-Stepping parameters 'NDTFAST == 0' in my *.in file, but it doesn't work. Then I recompile with USE_DEBUG turned on and rerun.

Rocky_Wang

Thank you Kate, but I don't know which Lateral Boundary Condition Types that I should change 'Fla' to, after the generation of the grid, I just want to add the tidal forcing file to see the variation of water level, and draw the cotidal lines. I didn't want to think about the changes of temperature ...

Rocky_Wang

HELLO, I am going to run roms in a domain with one open boundary(in south) in tidal region for first time.I have made a ROMS netcdf forcing file, using OSU Tidal Data Prediction Software(OTPS). when I closed every bounary, the model run successfully, but when I open the south boundary, a problem had ...

Rocky_Wang

hello kate In order to avoid building C and fortran libraries respectively， I installed netcdf in version before 4.2 and run the upwelling case successfullly. With the purpose of getting the Cross-section plots，I installed a binary distribution of NCL6.1.2, which compiled against NETCDF-4.2. So I re ...

Rocky_Wang

[root@centos upwelling]# ./oceanS < ocean_upwelling.in Model Input Parameters: ROMS/TOMS version 3.7 Monday - March 28, 2016 - 6:15:09 AM ----------------------------------------------------------------------------- Wind-Driven Upwelling/Downwelling over a Periodic Channel Operating system : Linux ...

Rocky_Wang

I changed the path in ocean_upwelling.in , tried it agian, lt looks like a success but a failure,oh..... NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00001440) STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME C => (i,j,k) Cu Cv Cw Max Speed 0 0 00:00:00 0 ...

Rocky_Wang

Thank you very much Kate,I reinstalled my Netcdf, the executable oceanS file is created successfully. But when I was going to run ROMS, another error happened : [root@centos upwelling]# ./oceanS < ocean_upwelling.in Model Input Parameters: ROMS/TOMS version 3.7 Monday - March 28, 2016 - 5:28:07 AM ...

Rocky_Wang

I Customize the Build Script like that: #export USE_MPI=on # distributed-memory parallelism #export USE_MPIF90=on # compile with mpif90 script #export which_MPI=mpich # compile with MPICH library #export which_MPI=mpich2 # compile with MPICH2 library #export which_MPI=openmpi # compile with OpenMPI ...

Rocky_Wang

I am a first time user to install and run ROMS, when I attempt to compile the Upwelling case, I meet this and don't know what's going wrong.I set the build.bash just as the Tutorial(https://www.myroms.org/wiki/ROMS_UNSW2008),I have installed the intel compiler and use ifort. could you mind tell me ...

Search found 50 matches