Ocean Modeling Discussion

turuncu

Thanks Hernan. Actually, i am using pyroms mostly to generate input files and i think that the python scripts i used are little bit old and did not add those specific attributes to the variables. Just for others, following shell scrip fix the netcdf files and also the problem. #!/bin/bash lst ...

turuncu

Hi, I am trying to run the latest version of the model (revision 923) but i am getting following error NETCDF_GET_TIME_1D - error while inquiring ID for variable: temp_time in input file: input/clim/AEGE16_clim_20040701_20040801_MEDINGV_fixed.nc call from: ROMS/Utility/get_cycle.F Found Error: 02 ...

turuncu

I was setting OBC as follows LBC(isFsur) == Cha Clo Clo Clo ! free-surface LBC(isUbar) == Fla Clo Clo Clo ! 2D U-momentum LBC(isVbar) == Fla Clo Clo Clo ! 2D V-momentum LBC(isUvel) == Rad Clo Clo Clo ! 3D U-momentum LBC(isVvel) == Rad Clo Clo Clo ! 3D V-momentum LBC(isMtke) == Cha Clo Clo Clo ...

turuncu

Thanks Kate. Actually, my domain is Mediterranean Sea and I am trying to define OBC for Atlantic Ocean. There is a buffer zone that extends to Atlantic Ocean from Strait of Gibraltar. So, only western boundary will be open. I am not sure what do you mean external values but for OBC, i will use ORAS4 ...

turuncu

Hi All, I am trying to setup a case with open boundary only in western side of the domain. The case without open boundary the simulation works without any problem but when i open western side of the domain the model produce very high ssh around 20 m at couple of time step and crash. First, i think ...

turuncu

By the way, if you want to add new atmospheric model component to RegESM modeling system. You just need to write a ESMF cap (in NUOPC terminology). Then rest of the modeling system will adjust to use new atmosphere model. To have an idea bout the cap, just look at a mod_esmf_atm.F90 in the RegESM ...

turuncu

Hi, ROMS is coupled with RegCM but the approach is little bit different when it is compared with COAWST modeling system. In this case, we are using driver approach and the the name of the modeling system is RegESM. You can find more information in the following link, https://github.com/uturuncoglu ...

turuncu

Hi, I applied ROMS to Caspian Sea. It is closed basin like lake. If you are interested, you could find more information in the following paper, Turuncoglu, U. U., Giuliani, G., Elguindi, N., and Giorgi, F.: Modelling the Caspian Sea and its catchment area using a coupled regional atmosphere-ocean ...

turuncu

Hi, I did a test with ROMS for Mediterranean domain as a part of PRACE preparatory phase project to find the optimal tiling combination. So, the roms_benchmark.pdf figure might help to define the tiling parameters. The 1/12° horizontal resolution ocean domain is used in this case. The grey and red ...

turuncu

Hi Kate, I am actively use pyroms. It is great tool. Thanks for it. I am using lpsolve to fix the bathymetry to prevent pressure gradient errors. As you suggest, it will be great to use ESMF (ESMPy) for interpolation instead of SCRIP. If you want we could collaborate about it, i have some experience ...

turuncu

okay. i implemented it. in this case i used ANA_SOURCE to define number of source point, their locations, direction, vertical shape (just 1 in surface and 0 in rest), qbar (as zero). Then, i calculate EmPmR and surface area in drag.F. I also modify qbar and qsrc in drag.F based on the values of ...

turuncu

I implemented code and tested with and without "factor of surface salinity" but code crashes after couple of time step. It seems that it is not stable and error grows. Qsrc option is related with rivers. Right? If i use Qsrc options then i have to define source points along the coastline and it won ...

turuncu

Hi Kate, I stated to implement the idea but i have one more minor question. Do i need to multiply surface flux with surface salinity (like t(i,j,N(ng),nrhs,salt) in SCORRECTION or SRELAXATION)? If i don't want to change the salinity in the modified region what is the best way? I also put my fresh ...

turuncu

Okay. Thanks. Sometimes, it is hard to find the correct way and place for these kind of custom setup. So, can i modify FORCES(ng)%stflx(salt) in set_vbc.F. The salt and temp correction or relaxation are applied in here. By the way, where is the correct place to calculate total stflx? diag.F is still ...

turuncu

Thanks Kate. I confirm that the volume is not changing in the original code unless applying rivers as a source points. So, your suggestion is the second one. Right? By the way, what do you think about zeta = zeta-dV. I think that it is not correct way to apply the correction. It seems that zeta ...

turuncu

Dear All, I would like to modify surface averaged E-P-R in the model to keep the net freshwater contribution as zero in the domain. The idealised domain have two basin connected by a narrow strait (similar to dogbone example). In this case, the calculated E-P-R in a part (let's say part A, left of ...

turuncu

Hi Again, I am trying to generate 'dogbone_ngc_composite.nc' using Matlab toolbox using following Matlab commands input = {'roms/Data/ROMS/Grid/dogbone_grd_left.nc'; ... 'roms/Data/ROMS/Grid/dogbone_grd_right.nc'}; output = 'dogbone_ngc_composite.nc' [S,G] = contact(input, output); but it gives ...

turuncu

Hi,

Thanks for your help. The file is okay now but dogbone_ngc_composite.nc has also same issue. If you don't mind could you regenerate that too. By the way, i could not see any difference in the definition of CPP flags in ticket 673

https://www.myroms.org/projects/src/ticket/673

Regards,

--ufuk

turuncu

Hi,

I am having same problem. Is this file corrupted somehow? Any idea to run the doggone test case? It could be related with NetCDF version but i am not sure.

Thanks,

--ufuk

turuncu

You could also use Python pyroms module which can be found in following link,

https://github.com/kshedstrom/pyroms

it also has option for linear programming approach to smooth bathymetry and as i know, it also supports Vstretching = 4. I hope it helps.

turuncu

Hi, You can use NCL to extrapolate data. Look at the following page, https://www.ncl.ucar.edu/Applications/grid_fill.shtml You could also use ESMF regridding capability in NCL. In this case, you could perform two step interpolation. Look at the attached example. It interpolates low resolution ...

turuncu

Thanks Kate and Arango. Yes, KPP helps but in this case the surface currents are not well represented by model. It was better with k-omega. I'll try to do more test. By the way, i also tried to split advection scheme but in that case it also activates biharmonic horizontal diffusion and i could not ...

turuncu

Hi, I had a problem with ROMS. The vertical mixing seems too much and i could not generate cold intermadiate layer (CIL) in my domain (Black Sea). I heard that sigma layer models has a problem in vertical mixing when there is a sharp gradient in the bathymetry. So, i ma just looking for a solution ...

turuncu

Hi, In my case, RDRG2 is little bit high (0.02d0). Now, i am trying to run the model with original value. I hope it works because i stuck into this problem. The region is the Sicily Strait in the Mediterranean Sea and it is not so close to the boundary. The z layer model (MITgcm) does not show this ...

turuncu

Hi Kate, Actually this region is in the strait of sicily and it affects the stability of the model in decadal run. i think that it is not an observed process and just something goes wrong. The model circulation seems very wierd after 10-15 years of simulation. it accumulates and produces the problem ...

turuncu

Hi, I have a weird results in my simulation. When i check the results of the model, the lowest layer salinity has decreasing trend in a specified region (see attached file, plot_dig_lev_1.pdf). The rest of the domain seems okay and also the main circulation patter seems correct. So, i don't have any ...

turuncu

No. It is not COWAST. I wrote my own driver

RegESM - https://github.com/uturuncoglu/RegESM

to couple ROMS with ICTP's RegCM regional climate model. I just want to test the latest version of the ROMS with the driver, then i found this error.

--ufuk

turuncu

Hi Again, To fix hanging problem i also fix "ROMS/Utility/ntimestep.F". I just initialise my_Nsteps as 1. ELSE my_Nsteps(ig)=1 my_Nsteps(ig)=MAX(my_Nsteps(ig), WindowSteps(ig)+1) step_counter(ng)=step_counter(ng)+WindowSteps(ig)+1 END IF # ifdef ADJOINT The previous version had a problem because my ...

turuncu

Hi, I think, i found a bug in main3d. When i run the model (revision 737) in coupling mode (coupled with atmospheric model) the model is not run correctly. For example, if i set coupling interval as 1-hour the model runs 30 step (in my case dt = 120 s) for first 1-hour but after that it runs until ...

turuncu

Thanks Kate. Yes, for standalone case i am using bulk formula but i am also planing to use the model in fully coupled mode without bulk. So, it is better to find a way to suppress the effect outside of the bulk module.

turuncu

Hi, You did not give any detail about the average monthly discharge from your river but you might consider to represent river using more than one source point. By this way you could distribute the discharge to multiple source point and you could decrease the velocity at that tiny grid that has the ...

turuncu

Hi All, I am using the output of the regional climate model to drive the ROMS ocean model and i have a question related with the wind speed. I compared the modelled wind speed and direction with the available observational datasets such as Blended Sea Winds (http://www.ncdc.noaa.gov/oa/rsad/air-sea ...

turuncu

Hi, I confirmed that sward must be multiplied by -1 to convert it to instantaneous field and also multiplied by -1 to correct its direction. Otherwise, the model gets cooling and this is evident in all the seasons (except winter). Look at the attached plots. The first one is the original code and i ...

turuncu

I think that it is better to add support for time varying WTYPE. The WTYPE could change by the time especially in the costal regions that includes the river (high/low discharge rate) or lagoons.

--ufuk

turuncu

Yes, i know him. I'll look at the output to find the source of the problem. In any case, i could write to him about the problem. Thanks for your help.

turuncu

No, i did not activated WET_DRY is this case but i'll plan to activate it in the production run. This could be also effect of the storm surge but i am not sure. i could change the code to write the zeta value in this cell. it might be a negative value greater than the cell depth. According to my ...

turuncu

Hi, i tried to decrease the time step up to 180 s but i have still problem with the model and it is also running much more slower. It gives error related with ice frazil in the 27th day of the simulation. 13106 27 07:18:00 5.470609E-03 2.530466E+03 2.530471E+03 7.354573E+13 (37,088,32) 1.099318E-03 ...

turuncu

Actually, i did not try to reduce the time step because i did not change the horizontal resolution. Currently, i am using following parameters to define the vertical grid, Vtransform == 2 ! transformation equation Vstretching == 4 ! stretching function ! Vertical S-coordinates parameters (see below ...

turuncu

Hi, I am running the ROMS-ice model in Caspian Sea. The problem is that the model generally underestimates the ice thickness and coverage when i compare with the ARCLAKE dataset. Based on my previous experience, the ice formation is closely linked with the depth of the sea and in my case, i set the ...

turuncu

Thanks Kate. I start to work in this problem. I hope i'll fixed it soon. By the way, i am using ice branch which was retrived in 20-03-2012 (i think that it is based on ROMS 3.5) and i just wonder that did you do any major code change after that. I am getting ice frazzle problem somtimes when i use ...

turuncu

Hi, I have same issue and i just wonder that is this bug fixed or still open? I am using ROMS 3.5 with ice model and it seems that it has still problem. When i try to restart the model with perfect restart option, the wet-dry mask is not correct at the beginning. i think that the problem is in the ...

turuncu

Hi, I have a closed domain - Caspian Sea - which the minimum depth is set as 5 meter. In this case, i activated the wetting-drying scheme to have dynamic masking and i set the critical depth as 0.2 meter like most of the applications. I know that the cell assumed as dry when the following criteria ...

turuncu

Hi,

You can find installation steps for pyroms on a Mac OS system. In the following link,

http://houstonthereisaproblem.blogspot.com

I hope it helps,

--ufuk

turuncu

Hi, I am using nudging (only T and S) in my application using following CPP flags. #define TCLM_NUDGING #define TCLIMATOLOGY but in my case i want to apply nudging to only a specific sub-region (or i want to exclude a sub-region from nudging). Does it make sense in the ocean modelling point of view ...

turuncu

Hi Kate, Thanks for your help. Actually, i coupled ROMS with RegCM regional climate (RCM) model. In this case, the atmospheric component sends net heat flux to ROMS. To conserve the heat fluxes between the model components, i want to use exactly same heat fluxes in the ocean model because i plan to ...

turuncu

Hi, I just wonder that is there anybody who used the ICE model without BULK_FLUX option? In the source file "ice_mk.h", stflx(:,:,items) variable is recalculated base on the ice coverage and the stflx argument is defined as intent(out). Maybe, i am wrong but it seems that the net heat flux is not ...

turuncu

Hi, I have a question about calculation of shflux (or stflx) to provide surface net heat flux into the ROMS ocean model (no bulk option activated). Does it include the solar radiation component? or it is just a non-solar heat flux. Then, what is the actual formula to calculate shflux? Is it ...

turuncu

Thanks, Kate. Now, it makes sense. Sometimes, it might be very hard to keep the volume in the range (especially for the lakes). I just try different combination of configurations. Some of them increase volume and some others decreases. Anyway, thanks again.

Regards,

--ufuk

turuncu

Hi, I just want to ask you about the usage of RUNOFF cpp flag for the model configuration. In the varinfo.dat file you can find the definition of it like "Fresh water runoff from land" but what happened if i set both river runoff using UV_PSOURCE and TS_PSOURCE and RUNOFF. Are they same? What is the ...

turuncu

Hi, I am running ROMS-ICE and after 9 year simulation, i got the "trouble in ice_frazil" error than model crashed. The main thing is that i am also using wetting/drying with ice. I also checked the restart file (one step before blow-up) and some single grid points behave so differently than the grid ...

turuncu

Hi, I think, i found the problem. It is related with MPI itself. After hanging the model, i attched the Totalview into one of the processes to see the problem. The model always hanging in the mp_distributeXXX calls. So, it is obvious that the problem triggering by some limitation of MPI. Actually, i ...

turuncu

Thanks Kate. is it in the git repository or somewhere else? anyway, i will try the newer version and let you know.

regards,

--ufuk

turuncu

Hi, I am running ice code but after certain time (3 months of simulation) the model stop running and freezing. It is waiting without any error messages. I put some print statement to the code to find the place that causes the problem. In main3d the code hanging just after calling ice_frazil ...

turuncu

Hi Kate, I add the print statement inside the if statement. The code complaining about "index" like following, NetCDF: Index exceeds dimension bound WRT_RST - error while writing variable: wetdry_mask_rho into restart NetCDF file for time record: 4 I think that i found the problem in the ROMS ...

turuncu

Hi, I have same problem like, WRT_RST - error while writing variable: wetdry_mask_rho into restart NetCDF file for time record: 4 In my configuration following CPP options are defined. .... #define PERFECT_RESTART #define WET_DRY .... and also i defined to write restart file each day with following ...

turuncu

Hi Again, Just for your information, i add the major rivers (Volga, Kura and Ural) into the model. With rough estimation, i calculated that the river discharge compensates the volume change in annual scale but i have still problem in the total net volume change. The results is still same with my ...

turuncu

Hi, I am running ROMS in a closed basin (Caspian Sea) using perpetual forcing to reach steady state solution but i have some strange results in net volume and zeta variables. The net volume (extracted from ROMS stout file) has decreasing trend in time (see the attached plot) and same pattern exist ...

turuncu

Hi, I setup ice model for Caspian Sea but i have a problem about that. I started the model from ice free condition using ana_ice.h. The content of the file is, DO j=JstrR,JendR DO i=Istr,IendR ui(i,j,1) = 0._r8 uie(i,j,1) = 0._r8 ui(i,j,2) = ui(i,j,1) uie(i,j,2) = uie(i,j,1) ENDDO ENDDO DO j=Jstr ...

turuncu

Hi,

Thanks Kate. I tried your fix but it gives error in the compile time and i changed it like,

Code: Select all

# if defined ICE_MODEL && !defined ANA_ICE

and it seems working.

Thanks again,

--ufuk

turuncu

Hi, I try to run ice code with only ANA_ICE cpp option (no ANA_INITIAL cpp flag is used) but model try to find the ice related variables in the initial condition netcdf file and gives following error, NLM: CHECKVARS - unable to find model variable: uice I check the code and found that this error ...

turuncu

Hi, Actually i am using the combination of python and NCL scripts and i have already install pyroms into Mac OS 10.7.2 (It is really great tool, thanks again Kate). You can find the detailed installation procedure in the following link, http://houstonthereisaproblem.blogspot.com/2011/12/installation ...

turuncu

Hi,

I wrote a NCL script to create SODA grid file from SODA data file. It basically creates the the grid file which includes DEPTH_bnds, MASK_T and MASK_UV variables. It works in my case. You can get and try it

Regards,

--ufuk

turuncu

Hi, I try to create initial condition file using pyroms but get_nc_BGrid_SODA call gives me following error, In [22]: src_grd = pyroms_toolbox.BGrid_SODA.get_nc_BGrid_SODA(src_grid_file, name=src_grid_name, area='global ...

turuncu

Hi Kate, Thanks for your comments. I did not set ALBEDO in my application. Do you think that it is the problem? Anyway, i plan to do more test with the latest version of the ice code. Which one is the recent in git or roms repository? Because the version that i used does not have ice formation ...

turuncu

Hi, I try to produce ice in Caspian Sea but still i could not. In this case, i am using a regional climate model (RegCM4) output to drive ROMS. The one dimensional lake model which is coupled into RegCM model produces ice with correct spatial distribution but ROMS not. The lake surface temperature ...

turuncu

Hi, I just wonder about the grid that you used for the global simulation. I try to create same kind of grid (periodic in East-West direction) to simulate the small belt around equator but i have some difficulties to setup EW periodic grid. If it is possible could you share your grid with me. Best ...

turuncu

Hi, I just wonder that is there any ready to go code to force ROMS model to split history file as monthly. The NDEFHIS controls the time step to split file but it can't split file exactly end of each month. I think that such kind of logic can be added into ./ROMS/Nonlinear/output.F file somewhere ...

turuncu

Hi, First of all, thanks for your help. From Kate's post, i=0 is a copy of i=Lm and i=Lm+1 is a copy of i=Lm are you sure about that because i expect like following, i=0 is a copy of i=Lm and i=1 is a copy of i=Lm+1 0, 1 ......... Lm, Lm+1 0 , 1 .............. Lm, Lm+1 What do you think? Anyway ...

turuncu

Hi, I want to create a periodic grid in East-West boundary. The grid box is 0,20N --------- 360,20N | | | | 0,-20N --------- 360,-20N So, it is like a belt over the earth equator. So, can i use ROMS to simulate such domain? My question is that do i need to add grid point in Eastern or Western side ...

turuncu

Hi Kate, Thanks for your help. The heat fluxes and the atmospheric variables such as wind, pressure, humidity etc. are provided by regional atmosphere model. The atmospheric model has 20 km grid in horizontal. I also compare the results with observational dataset and the fluxes seems correct. The ...

turuncu

Hi, I ran the ROMS with ice model for Caspian Sea but the model did not produce any ice in the domian. Both aice and hice variables are empty. The model starts with zero ice conditions (ana_ice.h file is used) and i am using the standard ice.in file. Based on the observation the northern part of the ...

turuncu

Hi, I am still investigating the problem with different combinations of the CPP options (vertical mixing parametrization etc.). I also tried to run the model with SST and SSS correction but the results is same. The attached plot shows the vertical profile of the salinity data in three different ...

turuncu

Hi, I did a climatological run but i think that there is a problem in the vertical profile of salinity field. The attached figure belongs to MAM average salinity field in three different cross sections. As it can be seen from the figure, the saline water always in the surface. This is same in the ...

turuncu

Hi Kate, I think, i found the code section. The code is crashing in the following section. #ifdef DISTRIBUTE CALL mp_exchange3d (ng, tile, model, 4, & & LBi, UBi, LBj, UBj, 1, 2, & & NghostPoints, EWperiodic, NSperiodic, & & ui, uie, vi, vie) CALL mp_exchange3d (ng, tile, model, 4, & & LBi, UBi, LBj ...

turuncu

Hi Kate, Thanks for your help. I am running in parallel mode without debugging option. I don't know why but when i put some print statement into the ana_ice.h just after the line that starts with "Must ..." the code is runs without problem. Otherwise, it crashes. It is so strange !!! So, if i found ...

turuncu

Hi,

Actually i add some print statement into code to find the line that code crash. The place is the analytical ice subroutine. I want to start the model without ice and i am using following file.

--ufuk

turuncu

Hi, I prepare the forcing file using NCPE/NCAR data and the specific humidity unit is kg/kg in it. So, based on the bulk_flux.F file i think that i have to modify the varinfo.dat file as, 'Qair' ! Input 'surface air specific humidity' ! relative or specific 'g/kg' ! [percentage or g/kg ] 'Qair ...

turuncu

Hi, I have some problem when i try to run the ice code. The ice model included version of ROMS stops at the beginning of the simulation without giving any warning or error messages. The last output as follows, .... .... ndtfast, nfast = 20 29 nfast/ndtfast = 1.45000 Centers of gravity and integrals ...

turuncu

Yes. The ERA Interim is 3 hourly. So the conversion factor will be different in that case. I am using 4.62963e-5 for 6 hourly ERA40 data.

turuncu

Hi, I think that you should do following, W/m^2/s --> 1/(60*60*dt) W/m^2 In this case dt (in hour) should be the time interval of the input data. For example for 6 hourly data the conversion factor will be, 1.0/(60.0*60.0*6.0) = 1.0/21600.0 = 4.62963e-5 So, you can also check the results after ...

turuncu

Thanks for your suggestions and help !!! I am also planing to use the version that has ice model but first of all i need to be sure that model simulates the climatological behavior of the lake. The main thing is that the northern part of the lake is covered by ice in the winter season but the ...

turuncu

Yes. I have already done that test. Actually this results belong to it. The lake is much more warmer than this when i used the net longwave and net shortwave (without LONGWAVE_OUT) but i don't know why it is colder in this case.

--ufuk

turuncu

Hi, Thanks for the clarification. Now, it is much more clear for me. I also did some tests using NCEP/NCAR Reanalysis after fixing the calculation of the net fluxes but the model is still cold in almost each season after 10 year simulation. In this case, the same forcing is applied perpetually. You ...

turuncu

Hi, Thanks for your reply. I don't use the ROMS_Tools either. I wrote my own tool using NCL and i check the ROMS_Tools to be sure about calculation of the fluxes. If you suggest any other tools that can get inspiration, please let me know. The region that i actually interested in Caspian Sea and the ...

turuncu

Hi, It could be a basic question but i am little bit confused to calculate the net shortwave and longwave fluxes in forcing input. If i use the NCEP/NCAR as a forcing input, the surface upward and downward fluxes can be used to estimate the net fluxes. For the shortwave flux this must be, swrad ...

turuncu

Hi, I did it before using following steps, 1 - if you don't have actual numeric values of the data. You should get them from the plot. You can use the PlotDigitizer to do that and can be downloaded from following link, http://plotdigitizer.sourceforge.net/ 2 - After getting the data (it is values ...

turuncu

Sorry about that. I don't know why you didn't see the plots but i replaced them with their gif versions. So, i can see the plots now.

--ufuk

turuncu

Hi, I try to set-up ROMS model for Caspian Sea. The result is little bit strange. The unrealistic cold regions over deepest regions. It seems that there is a problem in horizontal and vertical diffusion but i am not sure. I plot the surface and vertical temperature distribution and they can be seen ...

turuncu

Hi Kate,

Now, i am able to access the ice code and also ice model related section. I totally understand you, it is not easy to manage all of these things. I will try to setup and install it tomorrow. Thanks,

Regards,

--ufuk

turuncu

Okay. The cycle_length must be defined as 366 for leap years. This is my fault. Sorry. Not it is working without problem. Anyway, thanks for your kindly help. Regards,

turuncu

Hi,

I wrote to him but i could not get any response yet. Thanks,

--ufuk

turuncu

Hi, I think that i figure out the problem. As i understand from the ROMS/Utility/get_cycle.F file. The problem is mainly related with the "cycle_length" NetCDF attribute. The attribute must be set also in the time variables (or only in time variables, wind_time etc.) but in my case, i set it as a ...

turuncu

I am not sure. The configuration is exactly same that i wrote in my previous post. Just an idea but maybe i can solve the problem using following procedure, 1 - run the model for first year (1996). 2 - modify the time axis of the same forcing files for next year simulation (1997). 3 - restart model ...

turuncu

Hi Kate, I follow your suggestions and correct the value of the time axis. Now the values are, 0, 0.25, 0.5, 0.75, 1, 1.25, ..., 364.75, 365, 365.25, 365.5, 365.75 after this modification i also set the DSTART as, DSTART = 0.0d0 ! days TIDE_START = 0.0d0 ! days TIME_REF = 19600101.0d0 ! yyyymmdd.dd ...

turuncu

It is 6 hourly forcing file and the times are, tair_time = 46020, 46020.25, 46020.5, 46020.75, .... , 46381.25, 46381.5, 46381.75, 46382, 46382.25, 46382.5, 46382.75, 46383, 46383.25, 46383.5, 46383.75, 46384, 46384.25, 46384.5, 46384.75, 46385, 46385.25, 46385.5, 46385.75 ; and in the configuration ...

turuncu

Hi,

I am using latest tag (ROMS version 3.4). How can i contact with him?

Thanks,

--ufuk

turuncu

In addition to my previous post, I also check the surface temperature and it is not increase after first year simulation finish. It seems that model does not use the actual data, it acts like only using first time step data for the rest of the simulation. There is a movie in the following link but ...

turuncu

Hi, How can i add the ice related variables (such as ice thickness, varinfo.dat has variable name "hice") into the history file. Is there any special preprocessor flag to enable sea ice model? I check the code but there is only one flag (ICE_SHELF) which seems not related with the output. Thanks ...

turuncu

Hi, It seems that the model works without problem in this solution but when i try to plot the total turbulent kinetic energy (sum over the entire lake volume) vs time (see attached plot), it seems little bit strange after end of the first year simulation. What do you think? I also check the stdout ...

turuncu

Hi, I have a problem with the setting cycle_length netcdf attribute for the forcing files. I plan to run the model with perpetual annual forcing until it reach the steady state. My aim is that using the results as an initial condition of the real simulation. I think that the cycle_length netcdf ...

turuncu

Hi All, I am also working on coupling ROMS and WRF using ESMF library. I have already created a prototype version of coupled system using WRF 3.1.1, ROMS rev 382 and ESMF 3.1.0rp2. The report about this study can be found in the attached progress report. It is not cover all issues but it can give a ...

turuncu

Hi, I try to interpolate heat fluxes from atmospheric model to ROMS and i want to use conservative type interpolation for heat and rain fluxes. To do that i need to calculate area of the each grid cell in ROMS model. I was using following code to calculate rho centered grid area, area(i,j) = (1.0d0 ...

turuncu

Hi All, When i try to couple ROMS and WRF using ESMF, i modify the ROMS/Utility/get_date.F file to get correct date value. Because ROMS used each gridded component start and end date to calculate the coupler simulation interval, start and end date. In my application, the ocean.in definition was ...

turuncu

Hi, I am little bit confused about units of the fields that are imported into ROMS. In mct_roms_wrf.h file the heat fluxes are converted from Watts/m2 to Celsius m/s by multiplying -1.0_r8/(rho0*Cp) ratio but in bulk_flux.F file the flux variables are multiplying with Hscale again. Why? What is the ...

turuncu

Hi, I try to couple WRF with ROMS but after getting data from atmospheric component via ESMF library but ROMS reads the forcing file and replace the atmospheric model data with forcing one because i am getting data before calling ROMS routine. Is there any way to prevent reading forcing data? I ...

turuncu

Hi, I had problem with setting TIME_REF parameter. When i set it into 19990701.0d0 i expect that the caldate method in ROMS/Utility/get_date.F must return model start date as 1999-07-01 but it returns 1999-06-30. I add some print statement into this method and i figure out that the INT must be ...

turuncu

I solve the problem. It is mainly related bug in Master/esmf_roms.F file. To solve it modify following lines in ROMS_GetImportData routine,

old one,
LBj=BOUNDS(ng)%LBi(MyRank)
UBj=BOUNDS(ng)%UBi(MyRank)

new one,
LBj=BOUNDS(ng)%LBj(MyRank)
UBj=BOUNDS(ng)%UBj(MyRank)

The i must be changed to j.

turuncu

Hi, I try to run ROMS in coupling mode. I implement the ESMF WRF+ROMS coupling code and i try to do some test with that using 2 CPU. Each model component (WRF and ROMS) uses one CPU (or 1 PET) after running coupled code and few time steps, i got following error, MP_EXCHANGE2D - communication buffer ...

turuncu

Hi, In coupling.dat file, to convert Kelvin to Celsius the offset value is defined as 273.16d0 for Tair variable 'Tair' ! at 2 m, usually 'Tair' 'surface air temperature' 'Celsius' ! [Kelvin] 'r2dvar' 273.16d0 ! Kelvin to Celsius 1.0d0 but i think that it must be -273.16d0. Am i correct? Because the ...

turuncu

sorry, the revison is 382 and i download it yesterday. I don't know why but may be it is my fault. So anyway, thank you look into this.

turuncu

I found a bug in ROMS rev 832 in ROMS/Utility/inp_par.F file which triggers the error. To solve it go line 736 and replace 'waves' with 'atmos' in if else statement. There is two waves definition out there.

The code works after fixing this file. Anyway, thanks

turuncu

Hi, I try to couple WRF and ROMS and i set the parameters but i got following error, MOD_COUPLER - Number assigned processors: 001 not equal to spawned MPI nodes: 004 MOD_COUPLER - Number assigned processors: 001 not equal to spawned MPI nodes: 004 MOD_COUPLER - Number assigned processors: 001 not ...

turuncu

Hi, I try to compile ROMS with ESMF library but it gives following error (after correcting some minor bug related with ESMF and AIX IBM XLF compiler), "esmf_roms.f90", line 1009.11: 1513-119 (S) Identifier roms_import2d appeared in a USE statement but is not publicly accessible from this module ...

turuncu

Hi, I want to calculate ubar and vbar using u and v. Is it calculated using following formula, ubar = ubar + u*dz vbar = vbar + v*dz but in this case the units of the 2d vertical integrated momentum component will be m2/s (m/s * m) because we multiply with the thickness of the layer. In climatology ...

turuncu

Okay, now it makes sense because, i interpolated wind components to u and v points. I will change my code and re try. I hop it will work in this time. Anyway, i think that this point is not clear in documentation but i will check it again. In ROMS forcing page sustr and svstr are in u and v points ...

turuncu

I check the grid dimensions and they are same. May be i am little bit confused, so please tolerate me but i think that the error is expected when i use the lat-lon rather than x-y because lat-lon is in Lambert Conformal projection. See the attached plots that represent the grid points location ...

turuncu

Hi, I try to run ROMS for Mediterranean region. I have already setup grid and input files but when i run the model, it gives, REGRID - input gridded data does not contain model grid: Gridded: LonMin = -12.83264334 LonMax = 44.14562002 LatMin = 27.69616640 LatMax = 49.65302851 Model: LonMin = -12 ...

turuncu

Hi, i have almost same problem. Model gives error as, NETCDF_CHECK_VAR - inconsistent value of variable: hc 1.000000E+01 1.000000E+00 in file: Data/roms_ini_SODA_Y1997M7.nc i check the all files and value is same in all input netcdf files (10). I found that it causes from ROMS/Modules/mod_netcdf.F ...

turuncu

I mean that, i want to use specific CCSM3 SRES scenario data such as 20C (b30.30e) and i need to use POP output that belongs to that CCSM output. Is it possible? I also check the SODA page (i think that http://www.atmos.umd.edu/~ocean/) and i found the base article about that but there is no ...

turuncu

Hi,

I want to use POP output as IC/BC for ROMS but i couldn't find any information about that. Is there any way to use it? Any suggestion will be helpful.

Thanks,

--ufuk

turuncu

Hi, I think the problem is related with differences between versions of the ESMF. The coupling code is created with older version of ESMF. So, I apply following patch into Master/esmf_roms.F 552,553c552 < ! CALL ESMF_ArrayGet (romsGrdArray(igrid,ng), & < CALL ESMF_ArrayGet (romsGrdArray(igrid,ng),0 ...

turuncu

Hi, I try to install ROMS using ESMF library but i get following error messages, cd Build; /SFS2/progs/intel/impi/3.0/bin64/mpiifort -c -ip -O3 -xW -I/SFS2/progs/atmos/esmf-3.1.0r/mod/modO/Linux.intel.64.intelmp i.default -I/SFS2/progs/atmos/esmf-3.1.0r/src/include esmf_roms.f90 fortcom: Error: esmf ...

turuncu

yes, i could import netCDF4 module without error. The version of it is 0.7.4.

best,

--ufuk

turuncu

Hi, I try to install pyroms but i get following error messages when importing pyroms, /Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/pyroms/Dataset.py:65: UserWarning: netCDF4 not found -- using pupynere. warnings.warn('netCDF4 not found -- using pupynere.') Traceback ...

turuncu

Hi, What about lakes? I try to generate grid for lake and it gives following results (see attached image). I also faced some problems defining X,Y and resol parameter. I know the size of the lake (122 km, 81 km) and i set the resol parameter as 1000 m. So in this case i expect that grid must be ...

turuncu

Hi, my intel compiler version is also 9.1.038. > icc -V Intel(R) C Compiler for Intel(R) EM64T-based applications, Version 9.1 Build 20060323 Package ID: l_cc_p_9.1.038 Copyright (C) 1985-2006 Intel Corporation. All rights reserved. i try to upgrade it. if it works for me i will let you know. ufuk

turuncu

intel mpi comes with intel cluster toolkit and it use generic environment for all types of interconnect using dapl drivers. It is also use mpd for connection between nodes. You can find detailed information in,

http://www.intel.com/cd/software/produc ... 308295.htm

turuncu

Hi, I try to install ROMS but it gives following error messages, cd Build; /SFS/progs/intel/ict/3.0/mpi/3.0/bin64/mpiifort -c -ip -O3 -xW white_noise.f90 fortcom: Severe: **Internal compiler error: segmentation violation signal raised** Please report this error along with the circumstances in which ...

Search found 128 matches