Search found 128 matches

by turuncu
Mon Oct 01, 2018 7:58 pm
Forum: ROMS Discussion
Topic: missing temp_time
Replies: 2
Views: 2153

Re: missing temp_time

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=`ncdum...
by turuncu
Mon Oct 01, 2018 12:00 pm
Forum: ROMS Discussion
Topic: missing temp_time
Replies: 2
Views: 2153

missing temp_time

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 Li...
by turuncu
Tue Jul 17, 2018 7:33 pm
Forum: ROMS Discussion
Topic: problem in setting open boundary ...
Replies: 5
Views: 3199

Re: problem in setting open boundary ...

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 ! mixi...
by turuncu
Tue Jul 17, 2018 7:12 pm
Forum: ROMS Discussion
Topic: problem in setting open boundary ...
Replies: 5
Views: 3199

Re: problem in setting open boundary ...

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...
by turuncu
Tue Jul 17, 2018 9:52 am
Forum: ROMS Discussion
Topic: problem in setting open boundary ...
Replies: 5
Views: 3199

problem in setting open boundary ...

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 t...
by turuncu
Mon Apr 23, 2018 8:24 pm
Forum: ROMS Discussion
Topic: Coupling any Atmospheric model to ROMS
Replies: 3
Views: 2756

Re: Coupling any Atmospheric model to ROMS

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 mo...
by turuncu
Mon Apr 23, 2018 8:21 pm
Forum: ROMS Discussion
Topic: Coupling any Atmospheric model to ROMS
Replies: 3
Views: 2756

Re: Coupling any Atmospheric model to ROMS

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/R...
by turuncu
Sun May 07, 2017 12:41 pm
Forum: ROMS Discussion
Topic: ROMS application in Lake
Replies: 4
Views: 3213

Re: ROMS application in Lake

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 mo...
by turuncu
Sat Feb 20, 2016 8:27 am
Forum: ROMS Discussion
Topic: Optimal tilling once again
Replies: 4
Views: 4388

Re: Optimal tilling once again

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 a...
by turuncu
Thu Nov 05, 2015 6:47 am
Forum: ROMS Tools and Techniques
Topic: PYROMS installation
Replies: 45
Views: 70289

Re: PYROMS installation

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...
by turuncu
Fri Oct 30, 2015 4:52 pm
Forum: ROMS Discussion
Topic: artificially generated current along the strait ...
Replies: 11
Views: 4867

Re: artificially generated current along the strait ...

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 EmPm...
by turuncu
Wed Oct 28, 2015 8:24 pm
Forum: ROMS Discussion
Topic: artificially generated current along the strait ...
Replies: 11
Views: 4867

Re: artificially generated current along the strait ...

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 a...
by turuncu
Wed Oct 28, 2015 9:03 am
Forum: ROMS Discussion
Topic: artificially generated current along the strait ...
Replies: 11
Views: 4867

Re: artificially generated current along the strait ...

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 wa...
by turuncu
Mon Oct 26, 2015 6:12 pm
Forum: ROMS Discussion
Topic: artificially generated current along the strait ...
Replies: 11
Views: 4867

Re: artificially generated current along the strait ...

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...
by turuncu
Mon Oct 26, 2015 5:49 pm
Forum: ROMS Discussion
Topic: artificially generated current along the strait ...
Replies: 11
Views: 4867

Re: artificially generated current along the strait ...

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 = -d...
by turuncu
Mon Oct 26, 2015 4:33 pm
Forum: ROMS Discussion
Topic: artificially generated current along the strait ...
Replies: 11
Views: 4867

artificially generated current along the strait ...

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 t...
by turuncu
Mon Jul 27, 2015 8:56 am
Forum: ROMS Problems
Topic: Problem with netcdf file for Dogbone example
Replies: 6
Views: 5566

Re: Problem with netcdf file for Dogbone example

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 erro...
by turuncu
Thu Jul 16, 2015 6:59 am
Forum: ROMS Problems
Topic: Problem with netcdf file for Dogbone example
Replies: 6
Views: 5566

Re: Problem with netcdf file for Dogbone example

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
by turuncu
Wed Jul 15, 2015 11:55 am
Forum: ROMS Problems
Topic: Problem with netcdf file for Dogbone example
Replies: 6
Views: 5566

Re: Problem with netcdf file for Dogbone example

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
by turuncu
Sat May 16, 2015 10:24 am
Forum: ROMS Tools and Techniques
Topic: Bathymetry smoothing schemes
Replies: 4
Views: 6269

Re: Bathymetry smoothing schemes

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.
by turuncu
Wed Mar 04, 2015 6:28 am
Forum: ROMS Discussion
Topic: extrapolation
Replies: 5
Views: 4179

Re: extrapolation

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 atmosp...
by turuncu
Sat Nov 15, 2014 12:33 pm
Forum: ROMS Discussion
Topic: way of reducing vertical mixing ...
Replies: 3
Views: 2837

Re: way of reducing vertical mixing ...

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 ...
by turuncu
Thu Nov 13, 2014 8:12 am
Forum: ROMS Discussion
Topic: way of reducing vertical mixing ...
Replies: 3
Views: 2837

way of reducing vertical mixing ...

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 t...
by turuncu
Mon Oct 27, 2014 2:37 pm
Forum: ROMS Discussion
Topic: decreasing trend in bottom layer salinity ...
Replies: 4
Views: 3316

Re: decreasing trend in bottom layer salinity ...

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 ...
by turuncu
Mon Oct 20, 2014 7:45 pm
Forum: ROMS Discussion
Topic: decreasing trend in bottom layer salinity ...
Replies: 4
Views: 3316

Re: decreasing trend in bottom layer salinity ...

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...
by turuncu
Sat Oct 18, 2014 2:38 pm
Forum: ROMS Discussion
Topic: decreasing trend in bottom layer salinity ...
Replies: 4
Views: 3316

decreasing trend in bottom layer salinity ...

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...
by turuncu
Mon Sep 22, 2014 5:49 pm
Forum: ROMS Bugs
Topic: bug in main3d
Replies: 4
Views: 3512

Re: bug in main3d

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
by turuncu
Mon Sep 22, 2014 12:18 pm
Forum: ROMS Bugs
Topic: bug in main3d
Replies: 4
Views: 3512

Re: bug in main3d

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 b...
by turuncu
Mon Sep 22, 2014 10:58 am
Forum: ROMS Bugs
Topic: bug in main3d
Replies: 4
Views: 3512

bug in main3d

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 e...
by turuncu
Tue Sep 02, 2014 6:18 pm
Forum: ROMS Discussion
Topic: suppress the effect of the strong wind speed
Replies: 2
Views: 2256

Re: suppress the effect of the strong wind speed

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.
by turuncu
Tue Sep 02, 2014 6:08 pm
Forum: ROMS Discussion
Topic: river point in roms
Replies: 7
Views: 4536

Re: river point in roms

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 r...
by turuncu
Tue Sep 02, 2014 6:00 pm
Forum: ROMS Discussion
Topic: suppress the effect of the strong wind speed
Replies: 2
Views: 2256

suppress the effect of the strong wind speed

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/...
by turuncu
Mon Mar 03, 2014 7:39 am
Forum: ROMS Discussion
Topic: ECMWF ERA Interim Products related doubts
Replies: 8
Views: 8900

Re: ECMWF ERA Interim Products related doubts

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 ...
by turuncu
Sun Oct 06, 2013 2:12 pm
Forum: ROMS Discussion
Topic: spatially varying Jwtype in a forcing file
Replies: 5
Views: 4372

Re: spatially varying Jwtype in a forcing file

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
by turuncu
Tue Jul 02, 2013 5:25 pm
Forum: ROMS Ice
Topic: tuning ice model ...
Replies: 9
Views: 12321

Re: tuning ice model ...

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.
by turuncu
Sun Jun 30, 2013 9:04 am
Forum: ROMS Ice
Topic: tuning ice model ...
Replies: 9
Views: 12321

Re: tuning ice model ...

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 pre...
by turuncu
Sat Jun 29, 2013 12:30 pm
Forum: ROMS Ice
Topic: tuning ice model ...
Replies: 9
Views: 12321

Re: tuning ice model ...

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 ...
by turuncu
Fri Jun 28, 2013 7:34 am
Forum: ROMS Ice
Topic: tuning ice model ...
Replies: 9
Views: 12321

Re: tuning ice model ...

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 ...
by turuncu
Thu Jun 27, 2013 7:52 pm
Forum: ROMS Ice
Topic: tuning ice model ...
Replies: 9
Views: 12321

tuning ice model ...

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 ...
by turuncu
Fri May 24, 2013 5:53 pm
Forum: ROMS Discussion
Topic: ROMS wet/dry mask
Replies: 6
Views: 4109

Re: ROMS wet/dry mask

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 ...
by turuncu
Fri May 24, 2013 8:28 am
Forum: ROMS Discussion
Topic: ROMS wet/dry mask
Replies: 6
Views: 4109

Re: ROMS wet/dry mask

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 N...
by turuncu
Fri May 17, 2013 12:39 pm
Forum: ROMS Discussion
Topic: strange behaviour of wetting-drying ...
Replies: 2
Views: 3077

strange behaviour of wetting-drying ...

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 v...
by turuncu
Sat Nov 17, 2012 8:46 am
Forum: ROMS Tools and Techniques
Topic: Pyroms or Octant and Enthought
Replies: 2
Views: 3438

Re: Pyroms or Octant and Enthought

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
by turuncu
Mon Sep 03, 2012 10:26 am
Forum: ROMS Usage
Topic: applying nudging to the results in a sub region?
Replies: 1
Views: 2243

applying nudging to the results in a sub region?

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?...
by turuncu
Mon Jun 25, 2012 6:21 am
Forum: ROMS Ice
Topic: question about ice_mk.h and net heat flux
Replies: 3
Views: 5411

Re: question about ice_mk.h and net heat flux

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 ...
by turuncu
Sun Jun 24, 2012 8:51 pm
Forum: ROMS Ice
Topic: question about ice_mk.h and net heat flux
Replies: 3
Views: 5411

question about ice_mk.h and net heat flux

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...
by turuncu
Mon Jun 11, 2012 9:50 am
Forum: ROMS Usage
Topic: question about calculation of shflux ...
Replies: 2
Views: 3635

question about calculation of shflux ...

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 necessar...
by turuncu
Mon May 21, 2012 7:09 pm
Forum: ROMS Usage
Topic: question about RUNOFF cpp flag ...
Replies: 5
Views: 4548

Re: question about RUNOFF cpp flag ...

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
by turuncu
Mon May 21, 2012 7:44 am
Forum: ROMS Usage
Topic: question about RUNOFF cpp flag ...
Replies: 5
Views: 4548

question about RUNOFF cpp flag ...

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? W...
by turuncu
Wed May 09, 2012 9:58 am
Forum: ROMS Ice
Topic: question about ice frazil ...
Replies: 1
Views: 3913

question about ice frazil ...

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 tha...
by turuncu
Sat May 05, 2012 7:53 am
Forum: ROMS Ice
Topic: ice code hanging / freezing ...
Replies: 7
Views: 8946

Re: ice code hanging / freezing ...

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...
by turuncu
Wed May 02, 2012 7:26 am
Forum: ROMS Ice
Topic: ice code hanging / freezing ...
Replies: 7
Views: 8946

Re: ice code hanging / freezing ...

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

regards,

--ufuk
by turuncu
Tue May 01, 2012 11:22 am
Forum: ROMS Ice
Topic: ice code hanging / freezing ...
Replies: 7
Views: 8946

ice code hanging / freezing ...

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 subrouti...
by turuncu
Fri Apr 27, 2012 2:22 pm
Forum: ROMS Results
Topic: problem with perfect restart
Replies: 13
Views: 14801

Re: problem with perfect restart

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 R...
by turuncu
Sat Apr 21, 2012 5:37 pm
Forum: ROMS Results
Topic: problem with perfect restart
Replies: 13
Views: 14801

Re: problem with perfect restart

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 ...
by turuncu
Wed Apr 11, 2012 10:29 am
Forum: ROMS Discussion
Topic: conservation of volume ?
Replies: 1
Views: 1934

Re: conservation of volume ?

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 pre...
by turuncu
Tue Apr 10, 2012 10:18 am
Forum: ROMS Discussion
Topic: conservation of volume ?
Replies: 1
Views: 1934

conservation of volume ?

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 i...
by turuncu
Mon Mar 26, 2012 2:16 pm
Forum: ROMS Ice
Topic: some strange results from ice code ...
Replies: 1
Views: 3842

some strange results from ice code ...

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,Je...
by turuncu
Tue Mar 20, 2012 3:41 pm
Forum: ROMS Ice
Topic: using ANA_ICE and ANA_INITIAL
Replies: 3
Views: 5314

Re: using ANA_ICE and ANA_INITIAL

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
by turuncu
Mon Mar 19, 2012 1:54 pm
Forum: ROMS Ice
Topic: using ANA_ICE and ANA_INITIAL
Replies: 3
Views: 5314

using ANA_ICE and ANA_INITIAL

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 com...
by turuncu
Fri Mar 02, 2012 8:45 am
Forum: ROMS Tools and Techniques
Topic: ROMS Pre- and Post-Processing Packages on MAC?
Replies: 6
Views: 4295

Re: ROMS Pre- and Post-Processing Packages on MAC?

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...
by turuncu
Sat Feb 11, 2012 5:08 pm
Forum: ROMS Tools and Techniques
Topic: SODA grid problem in pyroms ...
Replies: 1
Views: 2810

Re: SODA grid problem in pyroms ...

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
by turuncu
Thu Feb 09, 2012 3:49 pm
Forum: ROMS Tools and Techniques
Topic: SODA grid problem in pyroms ...
Replies: 1
Views: 2810

SODA grid problem in pyroms ...

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')--------------------------------------------------------------------------- Key...
by turuncu
Thu Jan 12, 2012 8:05 am
Forum: ROMS Ice
Topic: running ice code ...
Replies: 22
Views: 24191

Re: running ice code ...

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 varia...
by turuncu
Wed Jan 11, 2012 8:38 am
Forum: ROMS Ice
Topic: running ice code ...
Replies: 22
Views: 24191

Re: running ice code ...

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 i...
by turuncu
Sun Nov 13, 2011 5:41 pm
Forum: ROMS Problems
Topic: ROMS kinetic energy and other problems
Replies: 3
Views: 5907

Re: ROMS kinetic energy and other problems

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 re...
by turuncu
Fri Nov 11, 2011 11:45 am
Forum: ROMS Usage
Topic: splitting the history file in monthly period ???
Replies: 1
Views: 2152

splitting the history file in monthly period ???

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 li...
by turuncu
Thu Nov 10, 2011 8:12 pm
Forum: ROMS Usage
Topic: creating periodic grid ...
Replies: 4
Views: 5115

Re: creating periodic grid ...

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, tha...
by turuncu
Fri Nov 04, 2011 3:23 pm
Forum: ROMS Usage
Topic: creating periodic grid ...
Replies: 4
Views: 5115

creating periodic grid ...

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 t...
by turuncu
Wed Aug 03, 2011 7:38 am
Forum: ROMS Ice
Topic: running ice code ...
Replies: 22
Views: 24191

Re: running ice code ...

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 on...
by turuncu
Tue Aug 02, 2011 11:05 am
Forum: ROMS Ice
Topic: running ice code ...
Replies: 22
Views: 24191

Re: running ice code ...

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...
by turuncu
Tue Aug 02, 2011 10:01 am
Forum: ROMS Discussion
Topic: possible problem in salinity profile
Replies: 1
Views: 1819

Re: possible problem in salinity profile

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 loca...
by turuncu
Tue Jul 19, 2011 7:03 pm
Forum: ROMS Discussion
Topic: possible problem in salinity profile
Replies: 1
Views: 1819

possible problem in salinity profile

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 ot...
by turuncu
Fri Jul 15, 2011 12:10 pm
Forum: ROMS Ice
Topic: running ice code ...
Replies: 22
Views: 24191

Re: running ice code ...

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, mod...
by turuncu
Fri Jul 15, 2011 10:17 am
Forum: ROMS Ice
Topic: running ice code ...
Replies: 22
Views: 24191

Re: running ice code ...

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, i...
by turuncu
Thu Jul 14, 2011 1:20 pm
Forum: ROMS Ice
Topic: running ice code ...
Replies: 22
Views: 24191

Re: running ice code ...

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
by turuncu
Thu Jul 14, 2011 9:15 am
Forum: ROMS Information
Topic: units of specific humidity
Replies: 10
Views: 12895

units of specific humidity

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, scal...
by turuncu
Thu Jul 14, 2011 8:08 am
Forum: ROMS Ice
Topic: running ice code ...
Replies: 22
Views: 24191

running ice code ...

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 ...
by turuncu
Thu Jun 16, 2011 7:39 am
Forum: ROMS Tools and Techniques
Topic: Simple question about conversion factor
Replies: 6
Views: 6641

Re: Simple question about conversion factor

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.
by turuncu
Wed Jun 15, 2011 8:53 am
Forum: ROMS Tools and Techniques
Topic: Simple question about conversion factor
Replies: 6
Views: 6641

Re: Simple question about conversion factor

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 conve...
by turuncu
Wed Jun 08, 2011 1:32 pm
Forum: ROMS Discussion
Topic: Which heat fluxes to input when using bulk formulation?
Replies: 24
Views: 16307

Re: Which heat fluxes to input when using bulk formulation?

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 surfac...
by turuncu
Wed Jun 08, 2011 12:20 pm
Forum: ROMS Discussion
Topic: Which heat fluxes to input when using bulk formulation?
Replies: 24
Views: 16307

Re: Which heat fluxes to input when using bulk formulation?

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
by turuncu
Wed Jun 08, 2011 8:41 am
Forum: ROMS Discussion
Topic: Which heat fluxes to input when using bulk formulation?
Replies: 24
Views: 16307

Re: Which heat fluxes to input when using bulk formulation?

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 ...
by turuncu
Mon May 30, 2011 9:05 am
Forum: ROMS Discussion
Topic: Which heat fluxes to input when using bulk formulation?
Replies: 24
Views: 16307

Re: Which heat fluxes to input when using bulk formulation?

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...
by turuncu
Sun May 29, 2011 9:08 am
Forum: ROMS Discussion
Topic: Which heat fluxes to input when using bulk formulation?
Replies: 24
Views: 16307

Re: Which heat fluxes to input when using bulk formulation?

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 = ds...
by turuncu
Thu May 26, 2011 7:41 am
Forum: ROMS Discussion
Topic: initial and boundary conditions from measured data
Replies: 6
Views: 3429

Re: initial and boundary conditions from measured data

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 an...
by turuncu
Sat May 21, 2011 9:27 am
Forum: ROMS Usage
Topic: unexpected temperature distribution ...
Replies: 3
Views: 3953

Re: unexpected temperature distribution ...

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
by turuncu
Fri May 20, 2011 9:31 pm
Forum: ROMS Usage
Topic: unexpected temperature distribution ...
Replies: 3
Views: 3953

unexpected temperature distribution ...

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 ...
by turuncu
Wed May 11, 2011 7:28 pm
Forum: ROMS Usage
Topic: activating or enable sea-ice model ???
Replies: 6
Views: 4267

Re: activating or enable sea-ice model ???

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
by turuncu
Wed May 11, 2011 10:37 am
Forum: ROMS Usage
Topic: set cycle_length attribute when using multiple forcing file
Replies: 11
Views: 8814

Re: set cycle_length attribute when using multiple forcing f

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,
by turuncu
Wed May 11, 2011 9:39 am
Forum: ROMS Usage
Topic: activating or enable sea-ice model ???
Replies: 6
Views: 4267

Re: activating or enable sea-ice model ???

Hi,

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

--ufuk
by turuncu
Wed May 11, 2011 9:37 am
Forum: ROMS Usage
Topic: set cycle_length attribute when using multiple forcing file
Replies: 11
Views: 8814

Re: set cycle_length attribute when using multiple forcing f

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 ...
by turuncu
Tue May 10, 2011 6:36 pm
Forum: ROMS Usage
Topic: set cycle_length attribute when using multiple forcing file
Replies: 11
Views: 8814

Re: set cycle_length attribute when using multiple forcing f

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 ...
by turuncu
Tue May 10, 2011 3:39 pm
Forum: ROMS Usage
Topic: set cycle_length attribute when using multiple forcing file
Replies: 11
Views: 8814

Re: set cycle_length attribute when using multiple forcing f

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 ...
by turuncu
Mon May 09, 2011 6:17 pm
Forum: ROMS Usage
Topic: set cycle_length attribute when using multiple forcing file
Replies: 11
Views: 8814

Re: set cycle_length attribute when using multiple forcing f

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...
by turuncu
Mon May 09, 2011 6:12 pm
Forum: ROMS Usage
Topic: activating or enable sea-ice model ???
Replies: 6
Views: 4267

Re: activating or enable sea-ice model ???

Hi,

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

Thanks,

--ufuk
by turuncu
Mon May 09, 2011 4:16 pm
Forum: ROMS Usage
Topic: set cycle_length attribute when using multiple forcing file
Replies: 11
Views: 8814

Re: set cycle_length attribute when using multiple forcing f

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 t...
by turuncu
Mon May 09, 2011 12:20 pm
Forum: ROMS Usage
Topic: activating or enable sea-ice model ???
Replies: 6
Views: 4267

activating or enable sea-ice model ???

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. ...
by turuncu
Mon May 09, 2011 10:03 am
Forum: ROMS Usage
Topic: set cycle_length attribute when using multiple forcing file
Replies: 11
Views: 8814

Re: set cycle_length attribute when using multiple forcing f

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 o...
by turuncu
Wed May 04, 2011 10:27 am
Forum: ROMS Usage
Topic: set cycle_length attribute when using multiple forcing file
Replies: 11
Views: 8814

set cycle_length attribute when using multiple forcing file

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 attr...
by turuncu
Tue Mar 30, 2010 5:43 am
Forum: ROMS Discussion
Topic: WRF_ROMS coupling compiling error
Replies: 6
Views: 2936

Re: WRF_ROMS coupling compiling error

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 ...
by turuncu
Mon Nov 02, 2009 2:41 pm
Forum: ROMS Discussion
Topic: conservative interpolation
Replies: 1
Views: 1914

conservative interpolation

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/...
by turuncu
Thu Oct 08, 2009 1:03 am
Forum: ROMS Discussion
Topic: Model Time Clocks
Replies: 10
Views: 24399

Re: Model Time Clocks

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, DST...
by turuncu
Mon Sep 07, 2009 9:24 pm
Forum: ROMS Discussion
Topic: coupling.dat definitions
Replies: 0
Views: 1302

coupling.dat definitions

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 ...
by turuncu
Fri Sep 04, 2009 4:02 am
Forum: ROMS Discussion
Topic: FRC_FILE definition for coupling
Replies: 0
Views: 1308

FRC_FILE definition for coupling

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 thin...
by turuncu
Wed Sep 02, 2009 7:03 am
Forum: ROMS Problems
Topic: problem in setting TIME_REF parameter
Replies: 0
Views: 1372

problem in setting TIME_REF parameter

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 repla...
by turuncu
Tue Sep 01, 2009 1:47 am
Forum: ROMS Problems
Topic: MP_EXCHANGE2D - communication buffer too small
Replies: 6
Views: 6141

Re: MP_EXCHANGE2D - communication buffer too small

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.
by turuncu
Mon Aug 31, 2009 9:36 pm
Forum: ROMS Problems
Topic: MP_EXCHANGE2D - communication buffer too small
Replies: 6
Views: 6141

MP_EXCHANGE2D - communication buffer too small

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 ...
by turuncu
Mon Aug 31, 2009 3:27 am
Forum: ROMS Discussion
Topic: converting kelvin to celsius coupling.dat file
Replies: 0
Views: 1411

converting kelvin to celsius coupling.dat file

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...
by turuncu
Fri Aug 14, 2009 6:29 pm
Forum: ROMS Problems
Topic: coupling parameters
Replies: 3
Views: 2393

Re: coupling parameters

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.
by turuncu
Fri Aug 14, 2009 4:18 pm
Forum: ROMS Problems
Topic: coupling parameters
Replies: 3
Views: 2393

Re: coupling parameters

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 :)
by turuncu
Fri Aug 14, 2009 3:53 pm
Forum: ROMS Problems
Topic: coupling parameters
Replies: 3
Views: 2393

coupling parameters

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 e...
by turuncu
Tue Aug 04, 2009 4:27 am
Forum: ROMS Installation
Topic: Compiling ROMS with ESMF
Replies: 5
Views: 5933

Re: Compiling ROMS with ESMF

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 m...
by turuncu
Fri Jul 31, 2009 5:30 pm
Forum: ROMS Discussion
Topic: calculating ubar and vbar
Replies: 10
Views: 6267

calculating ubar and vbar

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 ...
by turuncu
Tue Jul 28, 2009 4:46 pm
Forum: ROMS Problems
Topic: REGRID - input gridded data does not contain model grid:
Replies: 5
Views: 3651

Re: REGRID - input gridded data does not contain model grid:

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 (...
by turuncu
Tue Jul 28, 2009 7:18 am
Forum: ROMS Problems
Topic: REGRID - input gridded data does not contain model grid:
Replies: 5
Views: 3651

Re: REGRID - input gridded data does not contain model grid:

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 (lowe...
by turuncu
Mon Jul 27, 2009 10:33 pm
Forum: ROMS Problems
Topic: REGRID - input gridded data does not contain model grid:
Replies: 5
Views: 3651

REGRID - input gridded data does not contain model grid:

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.916...
by turuncu
Sun Jul 26, 2009 1:28 am
Forum: ROMS Problems
Topic: odd problem with hc
Replies: 4
Views: 7108

Re: odd problem with hc

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 f...
by turuncu
Mon Jul 06, 2009 5:42 pm
Forum: ROMS Usage
Topic: using POP fields as input
Replies: 3
Views: 2595

Re: using POP fields as input

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 informa...
by turuncu
Wed Jul 01, 2009 2:54 am
Forum: ROMS Usage
Topic: using POP fields as input
Replies: 3
Views: 2595

using POP fields as input

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
by turuncu
Tue Aug 12, 2008 8:51 am
Forum: ROMS Installation
Topic: Compiling ROMS with ESMF
Replies: 5
Views: 5933

Re: Compiling ROMS with ESMF

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...
by turuncu
Thu Jul 31, 2008 2:00 pm
Forum: ROMS Installation
Topic: Compiling ROMS with ESMF
Replies: 5
Views: 5933

Compiling ROMS with ESMF

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...
by turuncu
Fri Jul 11, 2008 10:32 am
Forum: ROMS Tools and Techniques
Topic: pyroms -- a suite of tools for working with ROMS in Python
Replies: 5
Views: 7742

Re: pyroms -- a suite of tools for working with ROMS in Python

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

best,

--ufuk
by turuncu
Thu Jul 10, 2008 8:17 pm
Forum: ROMS Tools and Techniques
Topic: pyroms -- a suite of tools for working with ROMS in Python
Replies: 5
Views: 7742

Re: pyroms -- a suite of tools for working with ROMS in Python

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 ...
by turuncu
Wed Jun 25, 2008 2:05 pm
Forum: ROMS Tools and Techniques
Topic: EASYGRID: Grid generation for ROMS
Replies: 24
Views: 30344

Re: EASYGRID: Grid generation for ROMS

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 cove...
by turuncu
Mon Jun 11, 2007 6:26 am
Forum: ROMS Installation
Topic: Internal compiler error when installing ROMS
Replies: 8
Views: 6133

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
by turuncu
Sun May 27, 2007 7:55 am
Forum: ROMS Installation
Topic: Internal compiler error when installing ROMS
Replies: 8
Views: 6133

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
by turuncu
Sat May 26, 2007 6:55 pm
Forum: ROMS Installation
Topic: Internal compiler error when installing ROMS
Replies: 8
Views: 6133

Internal compiler error when installing ROMS

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 ...