Ocean Modeling Discussion

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

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

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

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

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

Hi

In my remarks about inaccuracies there is an inaccuracy

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

Sorry

Best regards

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

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

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

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

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

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

Sorry, the big respect and for Paul

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

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

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

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

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

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

Hi

Thank you, Kate

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

Regards

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

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

Thank you, Kate.

We should repeat your way

Regards,
Boris

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

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

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

Kate,

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

Thanks,
Boris

Kate,

Ok, I need update ROMS.

Many thanks for the help.

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

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