Ocean Modeling Discussion

cloudy4568

Hi, there

The manual of Coawst version 3.7 has no presentations of sea ice.
Does any one know the progress of Coawst on sea ice?
Is Coawst version 3.7 the latest version?
Is Coawst still going to couple CICE?

Thank you! Any answers will be appreciated.

cloudy4568

Thank you for replying, Kate!

Yes, I referred most settings to your regional Beaufort domain. I have tried the RadNud option, the blowing up of high speed happens at the north-east corner (north and east are both open boundaries) and middle in depth, the speed of surrounding points are also high ...

cloudy4568

Hello everyone:

I am simulating the Beaufort Sea with three open boundaries and ICE_MODEL activated, but I found the water velocities terrible.
There is an eddy in the middle and a stronger current along the east boundary, could it be the problem of boundary conditions? I have tried other boundary ...

cloudy4568

Thank you, John! I have got the grid I expected with your method. Just make a rectangle grid in the Lambert Projection and then convert it back to lon/lat. Next I will try this curvilinear grid in ROMS. Thank you again for your help!

cloudy4568

Thank you, Kate! I think maybe I should put aside the issue of angles and make the grid first.

cloudy4568

Hello, everyone!

I am trying to make a curvilinear grid for ROMS, this grid is at equal distance at polar regions so it is curvilinear in longitude and latitude, just like the attached picture. Does anyone have the experience of making such a curvilinear grid for ROMS? Like how to calculate the ...

cloudy4568

Hi, the two problems were solved for my case.

1. Extreme salt values using Coawst when ICE_MODEL is defined.
The equation of salt flux in set_vbc.F is different:in Coawst:
# if defined ICE_MODEL || defined CICE_MODEL
# ifdef NO_SCORRECTION_ICE
fac_ice = 1.0_r8
# ifdef ICE_MODEL
IF(ai(i,j,liold ...

cloudy4568

Yes, I defined CORE_FORCING so
Q(i)=RH
RH=Hair(i,j)
and Hair(Qair) are from about 0.0001 to 0.01 (kg/kg) in my input data.

cloudy4568

Thank you, John!

Pair seems good, there is no 0 value and the minimum value is 99817 (mb) in data.

But maybe rhoAir does have problems because when CCSM_FLUX was tried (instead of BULK_FLUX), line 468 also has error:
rhoAir(i) = Pair(i,j)*100.0_r8/(blk_Rgas*TairK(i)* &
& (1.0_r8+0.61_r8*Hair(i,j ...

cloudy4568

Thank you, Kate!

The problem of line 222 of get_data.f90 seems to be solved. Then there are other errors:
Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.

Backtrace for this error:
#0 0x2b3cdab703af in ???
#0 0x2ae96a6e93af in ???
#1 0x1854904 in bulk ...

cloudy4568

Dear Kate and John, thank you so much for your replies!

When ICE_MODEL and CCSM_FLUXES are turned off, Coawst run ends at 120th step without any statements:
120 2018-10-02 01:00:00.00 1.017956E-02 1.425615E+02 1.425717E+02 2.320034E+12
(65,29,11) 2.889733E-03 8.124888E-04 1.192085E-01 8.955182E ...

cloudy4568

Thank you, Kate!

I have tried other advection and vertical diffusion schemes, checked input files(boundary, initial and forcing) which are all normal, and the error is still:
1 2018-10-02 00:00:30.00 5.450352E-03 1.430441E+02 1.430496E+02 2.326514E+12
(74,33,09) 4.527786E-03 2.774973E-04 5 ...

cloudy4568

Hi Roms users,

I am using Kate's version Roms with ice module, it could run well as Roms. But the same case is not running in coawst.

In Coawst, if define Roms only, with the same cpp files and input files, it has 'trouble in frazil' error in first steps:
trouble in ice_frazil 40 31 9 21 ...

cloudy4568

Hi all,

I encountered the same question, while compiling coawst with ICE_MODEL:

./Build/libNLM.a(initial.o): In function `initial_':
initial.f90:(.text+0x152b): undefined reference to `__ini_strengthcoef_mod_MOD_ini_strengthcoef'
initial.f90:(.text+0x1547): undefined reference to `__ini ...

cloudy4568

Thank you, John. You are right, I changed the following code in Compilers/Linux-gfortran.mk:
FFLAGS += -ftree-vectorize -ftree-loop-linear -funroll-loops -w -ffree-form -ffree-line-length-none -frecord-marker=4 -fconvert=big-endian
to
FFLAGS += -ftree-vectorize -funroll-loops -w -ffree-form ...

cloudy4568

Thank you, John. The full build is:

rm -f -r core *.ipo ./Build /mnt/home/cyue/make_macros.mk
cp -f /opt/ohpc/pub/libs/gnu8/openmpi3/netcdf-fortran/4.5.2/include/netcdf.mod ./Build
cp -f /opt/ohpc/pub/libs/gnu8/openmpi3/netcdf-fortran/4.5.2/include/typesizes.mod ./Build
cp -p /mnt/home/cyue/make ...

cloudy4568

Dear all,

I encountered the same error while compiling coawst(v3.6) project inlet_test and installing SCRIP_COAWST :

f951: sorry, unimplemented: Graphite loop optimizations cannot be used (isl is not available) (-fgraphite, -fgraphite-identity, -floop-nest-optimize, -floop-parallelize-all)
make ...

cloudy4568

Problem solved :D . I changed unit of specific humidity from percentage to kg/kg, and unit of air pressure from millibar to Pascals, then the result of water temperature is reasonable. Seems that different versions of ROMS wants different units for the same variable.
Thank you very much, Kate, I ...

cloudy4568

Thank you Kate, I mean from where we can know what unit of each input that Roms wants?

cloudy4568

Yes, the compiling of ice models was failed when both BULK_FLUXES and CCSM_FLUXES are turned off.

May I confirm the units of inputs? Should all the units of input the same with varinfo.dat? But I see two units under variable.

'Pair' ! Input
'surface air pressure'
'Pascal' ! [ROMS wants ...

cloudy4568

Thank you, Kate.

I have turned off ICE_ALB_EC92 and tried BULK_FLUXES. But the latent heat flux calculation of BULK_FLUXES is unstable and the model blew up. And I have checked the units of all inputs, they are reasonable.

I think I'd better try directly input surface total heat flux instead of ...

cloudy4568

Ok, thank you. In this version, both of BULK_FLUXES and CCSM_FLUXES need to be activated to activate CCSM_FLUXES.

# ifdef BULK_FLUXES
# ifdef CCSM_FLUXES
USE ccsm_flux_mod, ONLY : ccsm_flux
# else
USE bulk_flux_mod, ONLY : bulk_flux
# endif
# endif
-----------------
# ifdef BULK_FLUXES
# if ...

cloudy4568

Thank you, Hernan.
I turned off SPLINES_VDIFF and SPLINES_VVISC but the result is almost unchanged. I think there are something wrong with latent heat and sensible heat calculation. The CCSM_FLUXES calculates heat flux over open water with iteration which I can't understand

cloudy4568

we have been advancing the coupling to swan a lot more, it will be a more recent version of swan, we also added more wave current interactions, more sediment algorithms, vegetation/marsh, wave watch 3 coupling, wrf coupling, infragravity wave model, new wave-current asymmetry bedload, ... bunch of ...

cloudy4568

Hi all,

Thank you for your help in advance. There are some problems with my result of water temperature.

My roms result of water temperature is much higher than data from ECMWF or HYCOM from attached figure 1. And I think there is something wrong with heat flux so I compared each component in heat ...

cloudy4568

That's the reason! Thank you, Kate!

cloudy4568

Hi all,

I have a confused question. The boundary conditions for T and S are set (RadNud) in .in file:

! W S E N
! e o a o
! s u s r
! t t t t
! h h
!
! 1 2 3 4

LBC(isFsur) == Clo Clo Che Clo ! free-surface
LBC(isUbar) == Clo Clo Shc Clo ! 2D U-momentum
LBC(isVbar) == Clo Clo Shc Clo ! 2D V ...

cloudy4568

Problem solved. Tide_period in Variables is not the only one model is looking for, but also tide_period in Dimensions, like 'tide_period = 8'. I was using the other name('period=8'), so this error occurred. Thank you for your help.

cloudy4568

Hi, I have met the same problem with yours, but can't solve it by your method.

When .in file is this:

TIDENAME == /home/cyue/cyue/bohai/data/bohaitide5.nc
NFFILES == 9 ! number of forcing files
FRCNAME == /home/cyue/cyue/bohai/data/gom_Pair_era.nc \
/home/cyue/cyue/bohai/data/gom_Qair_era.nc ...

cloudy4568

Thank you, but I think the path is right. And I have modified the grid NC file already, last time the error is that it can't find 'xl' in the grid NC file, then I added it to the NC file.

cloudy4568

Thank you, Kate! I have replaced varinfo.dat with the one in your Apps repository. But another error happened:

INITIAL: Configuring and initializing forward nonlinear model ...
*******

GET_GRID - unable to open grid NetCDF file: wtype_grid

GET_GRID - unable to open grid NetCDF file: /home ...

cloudy4568

Thank you, Kate. Then there is another error,

Tile partition information for Grid 01: 81x59x11 tiling: 1x1

tile Istr Iend Jstr Jend Npts

Number of tracers: 2
0 1 81 1 59 52569

Tile minimum and maximum fractional coordinates for Grid 01:
(interior points only)

tile Xmin Xmax Ymin Ymax ...

cloudy4568

It is:

max_Ffiles=MAXVAL(nSSFfiles)
allocate ( SSF(max_Ffiles,Ngrids) )
allocate ( SSFids(max_Ffiles,Ngrids) )
allocate ( NSSFcount(max_Ffiles,Ngrids) )
SSFids(1:max_Ffiles,1:Ngrids)=-1
NSSFcount(1:max_Ffiles,1:Ngrids)=0
CASE ('SSFNAME')
label='SSF - Sources/Sinks forcing fields'
Npts ...

cloudy4568

Then there is another error:

forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image PC Routine Line Source
romsS 00000000009D6963 Unknown Unknown Unknown
libpthread-2.27.s 00007F83AFE81890 Unknown Unknown Unknown
romsS 00000000007DF325 inp_decode_mod_mp 1662 inp_decode.f90
romsS ...

cloudy4568

Thank you Kate, this problem has been solved.

cloudy4568

Hi all,

I am running my customized example for the first time. The source code I use is Kate's trunk (https://github.com/kshedstrom/roms.git), and the grid file, forcing files, boundary and initial files are prepared. The error is:

READ_PHYPAR - Variable index not yet loaded, idpthR
Found Error ...

cloudy4568

kate wrote: Thu Oct 24, 2019 7:14 pm Since my Perl is dated 2015 and yours is rather older, it might help if you could try a newer Perl.

Thank you, Kate, that's the problem.
I updated Perl by re-installing Linux on my desktop, then the compilation is successful now.

cloudy4568

This is a problem with your Perl. What version of Perl is this? My Perl works:

chinook03.rcs.alaska.edu 435% perl -v

This is perl 5, version 22, subversion 1 (v5.22.1) built for x86_64-linux-thread-multi

The version of perl on the cluster: 'This is perl, v5.10.1 (*) built for x86_64-linux ...

cloudy4568

Hi all,

This is an error while I compiling the 'upwelling' test on our school's cluster.
When everything was done and typed './build_roms.bash', the error occured

'rm -f -r core *.ipo /home/cyue/roms/upwelling/Build_roms /home/cyue/make_macros.mk
cp -f /home/cyue/nc/nc/include/netcdf.mod /home ...

Search found 39 matches