rivers example

[samad r]

hi alan
I am studying rivers in model
in your example (rivers_all.txt) what is the below numbers (red)?
 32153.0  -49.5000   -0.0500  Amazon A                Brazil          RIVDISP   RIVDISP   none
                                           24901.8   30546.8   35465.6   40194.7   43554.4   43605.2
                                           41303.5   36238.8   28495.0   21230.2   19398.9   20902.1
is this the case the below?
24901.8 m3/s   (Mean Jan Discharge)
30546.8  m3/s (Mean Feb Discharge)
35465.6 m3/s (Mean Mar Discharge)
and so on

Many Thanks
samad

[alan.wa...@hycom.org]

Yes, they are the monthly mean transports.

There is a README with this information:

COAPS 98> pwd
/home/awallcraft/hycom/ALL/bin
COAPS 99> ll *rivers.d RE*riv*
-rw-r-----. 1 awallcraft awallcraft 361295 May 17  2016 hycom_rivers.d
lrwxrwxrwx. 1 awallcraft awallcraft     35 Sep 28  2011 README.ALL.bin.rivers -> ../README.ALL/README.ALL.bin.rivers

The hycom_rivers program just removed the monthly temperature field from the entries in hycom_rivers.d that lie in the specified lon,lat box.

Alan.

[samad r]

Many Thanks alan

I am interesting to simulate the caspian sea. in the bin/README.ALL.bin.rivers and bin/hycom_rivers.d there is no information about rivers to caspian sea basin. and also data in RIVDISP is old (1984). what is your suggestion about newer river data?

and second can you send blkdat.input (configure file) that is used in this article
"Kara, A.B., Wallcraft Alan.J and Metzer E.J.,2010. Cunduz Murat, Impacts of freshwater on
the seasonal variations of surface salinity and circulation in the Caspian Sea. Continental Shelf
Research., No. 30"
 
and third:
the north caspian sea is covered by ice in winter. for use CICE what is the step to do?

Thanks alot
Best Regards
samad

[alan.wa...@hycom.org]

There is an old web link to NRL's Caspian Sea work: https://www7320.nrlssc.navy.mil/caspian/

I don't have the blkdat used in the paper, but I have an earlier one that is probably similar.  It is in the attached tar bundle with the rivers dataset we used.

As in this example, you can use any "rivers.d" file that is in the format of hycom_rivers.d.  It is also possible to use gridded runoff fields as a starting point instead.

It is tricky to setup CICE for use with HYCOM, you have to get an old version of ESMF working for example.  HYCOM comes with a built in thermodynamic sea ice model (i.e. the sea ice does not move) with is good enough for some purposes.  So I would start by turning it on (iceflg=1).

Finally, there are earlier posts about the Caspian Sea, particularly how to handle the fact that it is below sea level.  Simply search for Caspian on the forum to get these posts.

Alan.

[samad r]

Many Thanks alan

To study the Caspian Sea by ice Do I have enough ncep-cfsv2 data (TaqaQrQp.nc, strblk.nc, dlwsfc.nc, slwsfc.nc, ...)? or I need other data? Also, for interpolation, do I need to use the cice options in forcing scripts?

and for implement river, Should I apply an open boundary (ports.input and ...)?

Your advice is  helpful
Best Regards
samad

[alan.wa...@hycom.org]

Yes, CFSR forcing will work.

You don't need the cice options (i.e. don't include a subdirectory cice in the primary scratch directory).

Use HYCOM's "rivers as rain option", i.e. use force/rivers and set priver to 1 in blkdat.input.

Alan.

[samad r]

Thanks
Good Luck
samad

[samad r]

[alan.wa...@hycom.org]

I have no idea, beyond the usual advice about how to pick target densities.

If we were configuring a Caspian Sea model today, we would almost certainly use fixed sigma-Z layers (levels) down to the shelf break.  Based on nsigma=14, this is currently set at 54.76 m:

dp0k( 1) =   1.00 m    thkns =   1.00 m    depth =    1.00 m
dp0k( 2) =   1.19 m    thkns =   1.19 m    depth =    2.19 m
dp0k( 3) =   1.42 m    thkns =   1.42 m    depth =    3.61 m
dp0k( 4) =   1.68 m    thkns =   1.68 m    depth =    5.29 m
dp0k( 5) =   2.00 m    thkns =   2.00 m    depth =    7.29 m
dp0k( 6) =   2.38 m    thkns =   2.38 m    depth =    9.68 m
dp0k( 7) =   2.84 m    thkns =   2.84 m    depth =   12.52 m
dp0k( 6) =   2.38 m    thkns =   2.38 m    depth =    9.68 m
dp0k( 7) =   2.84 m    thkns =   2.84 m    depth =   12.52 m
dp0k( 8) =   3.38 m    thkns =   3.38 m    depth =   15.89 m
dp0k( 9) =   4.02 m    thkns =   4.02 m    depth =   19.91 m
dp0k(10) =   4.78 m    thkns =   4.78 m    depth =   24.69 m
dp0k(11) =   5.69 m    thkns =   5.69 m    depth =   30.37 m
dp0k(12) =   6.76 m    thkns =   6.76 m    depth =   37.14 m
dp0k(13) =   8.05 m    thkns =   8.05 m    depth =   45.18 m
dp0k(14) =   9.58 m    thkns =   9.58 m    depth =   54.76 m
dp0k(15) =  10.00 m    thkns =  10.00 m    depth =   64.76 m
dp0k(16) =  10.00 m    thkns =  10.00 m    depth =   74.76 m
dp0k(17) =  10.00 m    thkns =  10.00 m    depth =   84.76 m
dp0k(18) =  10.00 m    thkns =  10.00 m    depth =   94.76 m
dp0k(19) =  10.00 m    thkns =  10.00 m    depth =  104.76 m
dp0k(20) =  10.00 m    thkns =  10.00 m    depth =  114.76 m
dp0k(21) =  10.00 m    thkns =  10.00 m    depth =  124.76 m
dp0k(22) =  10.00 m    thkns =  10.00 m    depth =  134.76 m
dp0k(23) =  10.00 m    thkns =  10.00 m    depth =  144.76 m
dp0k(24) =  10.00 m    thkns =  10.00 m    depth =  154.76 m
dp0k(25) =  10.00 m    thkns =  10.00 m    depth =  164.76 m

You should check to see that ~55m is a good value for the shelf break (this is below the Caspian Sea surface, not below global mean sea level).

The Gulf of Mexico Reanalysis is an example of this aproach, although its target densities are sigma2.  See ftp://ftp.hycom.org/datasets/GOMu0.04/expt_50.1/conf/blkdat.input_501
This has nsigma=16 and the sill depth is about 83 m (note that isotop is set to 83.0).  also note that the terrain following dp00s's are all 1 m and the top deep layer, dp00(1), is 1 m.

Alan.

On Thursday, May 3, 2018 at 7:08:42 AM UTC-4, samad r wrote:

Thanks alan

This is work well
if i change the number of layer 25--> 35. What are the values for Sigma in blkdat.input?

   2.50   'sigma ' = layer  1 isopycnal target density (sigma units)
   3.25   'sigma ' = layer  2 isopycnal target density (sigma units)
   4.00   'sigma ' = layer  3 isopycnal target density (sigma units)
   4.75   'sigma ' = layer  4 isopycnal target density (sigma units)
   5.50   'sigma ' = layer  5 isopycnal target density (sigma units)
   6.00   'sigma ' = layer  6 isopycnal target density (sigma units)
   6.50   'sigma ' = layer  7 isopycnal target density (sigma units)
   7.00   'sigma ' = layer  8 isopycnal target density (sigma units)
   7.50   'sigma ' = layer  9 isopycnal target density (sigma units)
   7.90   'sigma ' = layer 10 isopycnal target density (sigma units)
   8.20   'sigma ' = layer 11 isopycnal target density (sigma units)
   8.50   'sigma ' = layer 12 isopycnal target density (sigma units)
   8.75   'sigma ' = layer 13 isopycnal target density (sigma units)
   9.00   'sigma ' = layer 14 isopycnal target density (sigma units)
   9.25   'sigma ' = layer 15 isopycnal target density (sigma units)
   9.50   'sigma ' = layer 16 isopycnal target density (sigma units)
   9.70   'sigma ' = layer 17 isopycnal target density (sigma units)
   9.92   'sigma ' = layer 18 isopycnal target density (sigma units)
  10.01   'sigma ' = layer 19 isopycnal target density (sigma units)
  10.10   'sigma ' = layer 20 isopycnal target density (sigma units)
  10.15   'sigma ' = layer 21 isopycnal target density (sigma units)
  10.18   'sigma ' = layer 22 isopycnal target density (sigma units)
  10.20   'sigma ' = layer 23 isopycnal target density (sigma units)
  10.22   'sigma ' = layer 24 isopycnal target density (sigma units)
  10.25   'sigma ' = layer 25 isopycnal target density (sigma units)

Thanks in advance
samad