Search found 30 matches

by Pysh
Mon Jun 09, 2014 5:21 pm
Forum: ROMS Usage
Topic: how about nudging SST to observation?
Replies: 8
Views: 5164

Re: how about nudging SST to observation?

And my question. If in the upper mixed layer temperature of the water can be considered equal SST, in the case of LMD parameterization of vertical mixing (and determination in a model of the depth of the upper mixed layer Hsbl) can I use nudging the temperature of the water to the measured SST in th...
by Pysh
Wed Apr 23, 2014 6:35 am
Forum: ROMS Tools and Techniques
Topic: Circumpolar grid
Replies: 5
Views: 1697

Re: Circumpolar grid

Hi May be this script? Sometimes it works normally :) (Gurvan Madec, Maurice Imbard; A global ocean mesh to overcome the North Pole singularity; Climate Dynamics (1996) 12:381-388) Cheers Boris function glob_grid global Jeq JM dd % ---- net parameters ------------ IM = 280; % number of longitudes (w...
by Pysh
Mon Dec 30, 2013 8:51 pm
Forum: ROMS Usage
Topic: climate
Replies: 4
Views: 2112

Re: climate

Dear redhat007 We don't use agrif, because of absence ice model. Russia is cold country :wink: , all seas near Russia (North Caspian also) covered with ice in winter. Some of seas - year round :( . For data preparing (from raw GFS data) we use wgrib2 and our fortran utility. You can contact with me ...
by Pysh
Thu Dec 26, 2013 8:31 am
Forum: ROMS Usage
Topic: climate
Replies: 4
Views: 2112

Re: climate

Hello In the Russian Hydrometcenter attempts modeling Caspian Sea. For forcing (with bulk-formula) we use GFS data (Global Forecast System, http://www.nco.ncep.noaa.gov/pmb/products/gfs/#GFS ), for river runoff (Volga, Kura, Ural) data from Russian sources (Государственный водный кадастр, 1980-1987,...
by Pysh
Fri Nov 08, 2013 4:59 pm
Forum: ROMS Discussion
Topic: ANALYTICAL flag
Replies: 0
Views: 656

ANALYTICAL flag

Hi If I don't use any ANA_ module I should include #define ANALYTICAL flag into header file, because of unconditional call of ana_nudgcoef routine from initial.f90 routine, that ( ana_nudgcoef routine) will be included into analytical.f90 only when ANALYTICAL flag is defined. Is it true or I can inc...
by Pysh
Fri Aug 02, 2013 11:07 am
Forum: ROMS Ice
Topic: WET_DRY and PERFECT_RESTART in Ice branch
Replies: 6
Views: 3956

Re: WET_DRY and PERFECT_RESTART in Ice branch

Hi

In my remarks about inaccuracies there is an inaccuracy :D

Neither in the file set_state.F, nor in the file get_state.F. In the file set_masks.F

Sorry

Best regards
by Pysh
Wed Jul 31, 2013 11:37 pm
Forum: ROMS Ice
Topic: WET_DRY and PERFECT_RESTART in Ice branch
Replies: 6
Views: 3956

Re: WET_DRY and PERFECT_RESTART in Ice branch

Hi After last code correction some inaccuracies are appeared, I believe in set_state.F string 78 # if defined PERFECT_RESTART && defined WET_DRY string 128 # if defined PERFECT_RESTART && defined WET_DRY string 182 # if defined PERFECT_RESTART && defined WET_DRY string 226 # if defined PERFECT_RESTA...
by Pysh
Mon Jul 15, 2013 10:07 pm
Forum: ROMS Ice
Topic: WET_DRY and PERFECT_RESTART in Ice branch
Replies: 6
Views: 3956

Re: WET_DRY and PERFECT_RESTART in Ice branch

And in file set_vbc.F !*********************************************************************** SUBROUTINE set_vbc (ng, tile) !*********************************************************************** ! USE mod_param USE mod_grid USE mod_forces USE mod_ocean USE mod_stepping ..............................
by Pysh
Mon Jul 15, 2013 8:13 pm
Forum: ROMS Ice
Topic: WET_DRY and PERFECT_RESTART in Ice branch
Replies: 6
Views: 3956

WET_DRY and PERFECT_RESTART in Ice branch

I run model with flags #define PERFECT_RESTART #define WET_DRY #define SOLVE3D ...................................... #define ICE_MODEL #ifdef SOLVE3D # ifdef ICE_MODEL # define ICE_THERMO # define ICE_MK # undef ICE_ALB_EC92 # define ICE_SMOOTH # define ICE_MOMENTUM # define ICE_MOM_BULK # define I...
by Pysh
Sat Nov 24, 2012 9:05 pm
Forum: ROMS Problems
Topic: Wet/dry stiffness
Replies: 4
Views: 1898

Re: Wet/dry stiffness

In the module stiffness.F values of rx0 and rx1 are computed through DO j=Jstr,Jend DO i=IstrU,Iend IF (umask(i,j).gt.0.0_r8) THEN my_rx0=MAX(my_rx0,ABS((z_w(i,j,0)-z_w(i-1,j,0))/ & & (z_w(i,j,0)+z_w(i-1,j,0)))) DO k=1,N(ng) my_rx1=MAX(my_rx1,ABS((z_w(i,j,k )-z_w(i-1,j,k )+ & & z_w(i,j,k-1)-z_w(i-1,...
by Pysh
Sun Nov 18, 2012 10:26 pm
Forum: ROMS Information
Topic: units of specific humidity
Replies: 10
Views: 5112

Re: units of specific humidity

I use the module ccsm_bulk (not bulk_flux). I'm not sure, but in this algorithm is not provided for by entering RH, only SH in kg/kg. Coefficient which transform the input variable is governed by the parameters of the Qair in file varinfo.dat By default, this is, as I understand, 0.01 (for RH in per...
by Pysh
Sat Nov 17, 2012 7:30 am
Forum: ROMS Information
Topic: units of specific humidity
Replies: 10
Views: 5112

Re: units of specific humidity

As I understand it is correct to take humidity SH but as a percentage, i.e., multiply the NECP/NCAR specific humidity by 100. In my experiments I do not multiply NECP/NCAR specific humidity by 100 (1, kg/kg) and multiply it by 1,000 (2, g/kg). Coincidence of computing the of latent heat flux in the ...
by Pysh
Tue Aug 28, 2012 7:25 am
Forum: ROMS Ice
Topic: Sea ice in shallow water (anchor ice)
Replies: 4
Views: 3071

Re: Sea ice in shallow water (anchor ice)

Sorry, the big respect and for Paul
by Pysh
Sat Aug 25, 2012 7:05 am
Forum: ROMS Ice
Topic: Sea ice in shallow water (anchor ice)
Replies: 4
Views: 3071

Re: Sea ice in shallow water (anchor ice)

Thank you, Kate I've calculated Caspian model in your code with ICE_SHOREFAST now only on a month. Results are normal. As I understand in your code thermal flux at the bottom starts to be limited with depth of 'clear' water (under ice) of 1.5 meter and is equal to zero at depth of 'clear' water less...
by Pysh
Sat Aug 25, 2012 6:57 am
Forum: ROMS Ice
Topic: Double time in avg file
Replies: 6
Views: 3629

Re: Double time in avg file

Thank you, Kate

Bug in time variable in avg-files is now absent.

Coding of a part of the nature will soon demand all natural resources :)

Good luck
by Pysh
Thu Aug 23, 2012 10:24 pm
Forum: ROMS Ice
Topic: Sea ice in shallow water (anchor ice)
Replies: 4
Views: 3071

Sea ice in shallow water (anchor ice)

Hi In my experiments (Caspian sea with ice model) there was a following situation: ice thickness was more then depth in the same points. (7 meters ice thickness in points with 2 meters depth). In north part of Caspian sea large area has depth ~ 2 meters and in winter this part becomes covered by an ...
by Pysh
Thu Aug 23, 2012 5:28 am
Forum: ROMS Ice
Topic: Double time in avg file
Replies: 6
Views: 3629

Re: Double time in avg file

Not with NAVG = 120, NHIS = 120, but with NAVG = 180, NHIS = 180. Before substituting of new vertion of set_avg.F netcdf casp_avg_20111125 { dimensions: xi_rho = 200 ; xi_u = 199 ; xi_v = 200 ; xi_psi = 199 ; eta_rho = 500 ; eta_u = 500 ; eta_v = 499 ; eta_psi = 499 ; s_rho = 20 ; s_w = 21 ; tracer ...
by Pysh
Wed Aug 22, 2012 5:29 am
Forum: ROMS Ice
Topic: Double time in avg file
Replies: 6
Views: 3629

Re: Double time in avg file

Thank you for reply, duran. I'v looked messages you recommended, but I believe the problem has remained. My values of NAVG, NHIS, NDIA are not optimal (as is not optimal to run the model step by step with one day period with restarts). It is caused by local circumstances and debugging version of a p...
by Pysh
Tue Aug 21, 2012 11:06 pm
Forum: ROMS Ice
Topic: Double time in avg file
Replies: 6
Views: 3629

Double time in avg file

Hi Strange result for ocean_time variable in avg file. My parameters for output files: 1440 ntimes Number of timesteps for 3-D equations. 60.000 dt Timestep size (s) for 3-D equations. 60 ndtfast Number of timesteps for 2-D equations between each 3D timestep. 1 ERstr Starting ensemble/perturbation r...
by Pysh
Tue Aug 14, 2012 12:28 pm
Forum: ROMS Ice
Topic: Question about set_avg.F
Replies: 2
Views: 2030

Re: Question about set_avg.F

Hi

Thank you, Kate

I have replaced set_avg.F with last your version. Program work normally. My corrections were not absolutely correct.

Regards
by Pysh
Thu Aug 09, 2012 7:26 pm
Forum: ROMS Ice
Topic: Question about set_avg.F
Replies: 2
Views: 2030

Question about set_avg.F

Hi

I try to work with a seaice code (omlab ver 1488). Tanslation process show error in set_avg.F. May be I haven't last vertion? In vertion 1488, as it seems to me, some corrections should be made (in attached files, left - original, right - corrected) I am not sure.

Thanks
by Pysh
Sat Jul 14, 2012 12:01 pm
Forum: ROMS Results
Topic: problem with perfect restart
Replies: 13
Views: 6504

Re: problem with perfect restart

Hi, all I know that it isn't new subject in the forum, but I can't understand is it real to use correctly PERFECT_RESTART with WET_DRY options now. In my case (changeset_624) restart file include fields WETDRY_MASK_RHO = 1 in all points for both records (I use LcycleRST == T), WETDRY_MASK_U = 1, WET...
by Pysh
Mon Jul 02, 2012 6:45 pm
Forum: ROMS Usage
Topic: question about RUNOFF cpp flag ...
Replies: 5
Views: 2346

Re: question about RUNOFF cpp flag ...

Thank you, Kate.

We should repeat your way :)

Regards,
Boris
by Pysh
Mon Jul 02, 2012 2:51 pm
Forum: ROMS Usage
Topic: question about RUNOFF cpp flag ...
Replies: 5
Views: 2346

Re: question about RUNOFF cpp flag ...

Hi, We have some problems with salinity distribution from the river runoff turuncu some times ago (21 May 2012) wrote I just want to ask you about the usage of RUNOFF cpp flag for the model configuration and Kate answered The purpose of RUNOFF is to distribute fresh water over more points - some far...
by Pysh
Sun Mar 18, 2012 4:20 pm
Forum: ROMS Discussion
Topic: Weights in VSHAPE and layers thickness?
Replies: 6
Views: 1873

Re: Weights in VSHAPE and layers thickness?

Hi, John

It's clear now. I'am working at small program to prepare river input. And thank you for your matlab tools. It is very usefull.

Thank you very much,
Boris
by Pysh
Sat Mar 17, 2012 5:56 pm
Forum: ROMS Discussion
Topic: Weights in VSHAPE and layers thickness?
Replies: 6
Views: 1873

Weights in VSHAPE and layers thickness?

Dear all, If I need distribute river input among the layers, are the weights of VSHAPE (for example 0.5,0.3,0.2) take into account layers thickness or these are the specific intensities of input, t.e real disstribution among the layers are (0.5 * h1, 0.3 * h2, 0.2 * h3) / (h1 + h2 + h3), where h1, h...
by Pysh
Fri Mar 09, 2012 8:20 pm
Forum: ROMS Problems
Topic: Run error: Subscript #2 of the array AKT_BAK has value 1
Replies: 5
Views: 2454

Re: Run error: Subscript #2 of the array AKT_BAK has value 1

Kate,

I'v updated the version of ROMS. The model run without problems.

Thanks,
Boris
by Pysh
Fri Mar 09, 2012 6:42 pm
Forum: ROMS Problems
Topic: Run error: Subscript #2 of the array AKT_BAK has value 1
Replies: 5
Views: 2454

Re: Run error: Subscript #2 of the array AKT_BAK has value 1

Kate,

Ok, I need update ROMS.

Many thanks for the help.
by Pysh
Fri Mar 09, 2012 6:26 pm
Forum: ROMS Problems
Topic: Run error: Subscript #2 of the array AKT_BAK has value 1
Replies: 5
Views: 2454

Re: Run error: Subscript #2 of the array AKT_BAK has value 1

Thank you, Kate

In casp.in and in ocea.in value NGRIDS = 1. How it can be equal zero I'll try to undestand. Now I don't see it in makefile
I have the version of ROMS 580 (2011-11-28) The ocean.in is also from this version.

Regards,
Boris
by Pysh
Fri Mar 09, 2012 5:09 pm
Forum: ROMS Problems
Topic: Run error: Subscript #2 of the array AKT_BAK has value 1
Replies: 5
Views: 2454

Run error: Subscript #2 of the array AKT_BAK has value 1

Hi! I tried to run ROMS application for Caspian sea. I am beginner. The model not runs with the message forrtl: severe (408): fort: (2): Subscript #2 of the array AKT_BAK has value 1 which is greater than the upper bound of 0 This message arises when in configuration file (casp.h) there are #undef M...