Search found 18 matches

by stlaur
Tue Mar 10, 2020 1:12 pm
Forum: ROMS Source
Topic: Question about HSIMT in step3d_t.F (trunk version)
Replies: 1
Views: 341

Question about HSIMT in step3d_t.F (trunk version)

Good morning, I apologize in advance if I'm wrong about what follows. I could be misinterpreting the code. Early on in step3d_t.F (lines 334-350), oHz is computed as: IF (Lmpdata) THEN DO k=1,N(ng) DO j=Jstrm2,Jendp2 DO i=Istrm2,Iendp2 oHz(i,j,k)=1.0_r8/Hz(i,j,k) END DO END DO END DO ELSE DO k=1,N(n...
by stlaur
Sun Mar 01, 2020 4:48 pm
Forum: ROMS Source
Topic: Typo in roms.in following overhaul of advection schemes
Replies: 4
Views: 525

Re: Typo in roms.in following overhaul of advection schemes

Just to complete the earlier post---here is the relevant part of ROMS/Utility/read_phypar.F (as of rev.1008): CASE ('Hadvection') IF (itracer.lt.(NAT+NPT)) THEN itracer=itracer+1 ELSE itracer=1 ! next nested grid END IF itrc=itracer Npts=load_tadv(Nval, Cval, line, nline, itrc, igrid, & & itracer, 1...
by stlaur
Fri Feb 28, 2020 9:38 pm
Forum: ROMS Source
Topic: Typo in roms.in following overhaul of advection schemes
Replies: 4
Views: 525

Typo in roms.in following overhaul of advection schemes

Good afternoon, I think there is a typo in the GLOSSARY section of roms.in. It says: ! or in nested applications ! ! Hadvection == U3 \ ! temperature, Grid 1 ! U3 \ ! temperature, Grid 2 ! U3 \ ! temperature, Grid 3 ! HSIMT \ ! salinity, Grid 1 ! HSIMT \ ! salinity, Grid 2 ! HSIMT ! salinity, Grid 3...
by stlaur
Wed Sep 11, 2019 7:11 pm
Forum: Job Opportunities
Topic: Postdoctoral scholar position in Estuarine Plastics Modeling
Replies: 0
Views: 593

Postdoctoral scholar position in Estuarine Plastics Modeling

A Post-Doctoral Scholar position is available at The Pennsylvania State University to study the transport and fate of plastics in the Chesapeake Bay using numerical models and to collaborate with experimentalists. The Post-Doctoral Scholar will work with an interdisciplinary team at two institutions...
by stlaur
Mon Jul 22, 2019 2:15 pm
Forum: ROMS Bugs
Topic: Typo in ad_set_vbc.F line 1165 rev 973
Replies: 1
Views: 285

Typo in ad_set_vbc.F line 1165 rev 973

Hi,
There is a typo in ad_set_vbc.F on line 1165:

Code: Select all

     &            ad_stfx(i,j,itemp)
should be

Code: Select all

     &            ad_stflx(i,j,itemp)
(needs a "l" between "f" and "x").
Pierre
by stlaur
Wed Aug 29, 2018 9:19 pm
Forum: ROMS Bugs
Topic: Typo in mod_scalars.F of rev. >=919
Replies: 0
Views: 812

Typo in mod_scalars.F of rev. >=919

Good afternoon, I think there's a tiny typo in the recent versions of mod_scalars.F (revision >=919): #ifdef T_PASSIVE IF (.not.allocated(inert)) THEN allocate ( inert(NPT) ) END IF #endif Since "inert" is a pointer, I think it should be "IF(.not.associated(inert)" rather than "IF(.not.allocated(ine...
by stlaur
Sat Dec 16, 2017 5:16 pm
Forum: ROMS Bugs
Topic: WTYPE=9 causes division by zero in lmd_swfrac.F
Replies: 2
Views: 1405

Re: WTYPE=9 causes division by zero in lmd_swfrac.F

Maybe you can figure out what the parameters should be - probably a very small but finite value not zero WTYPE=9 (dark coastal water) is absent from the original papers (e.g. Paulson & Simpson JPO 1977) but the following publication (attached) from: https://publications.europa.eu/en/publication-det...
by stlaur
Fri Dec 15, 2017 6:27 pm
Forum: ROMS Bugs
Topic: WTYPE=9 causes division by zero in lmd_swfrac.F
Replies: 2
Views: 1405

WTYPE=9 causes division by zero in lmd_swfrac.F

I would like to report a division by zero that occurs when WTYPE=9 and SOLAR_SOURCE are selected. WTYPE=9 is defined in mod_scalars.F as: real(r8), dimension(9) :: lmd_mu1 = & & (/ 0.35_r8, 0.6_r8, 1.0_r8, 1.5_r8, 1.4_r8, & & 0.42_r8, 0.37_r8, 0.33_r8, 0.0_r8 /) Note that lmd_mu1 is 0.0_r8 for WTYPE...
by stlaur
Thu Oct 13, 2016 1:06 pm
Forum: ROMS Discussion
Topic: Spurious supercooling in ROMS-CICE simulation
Replies: 33
Views: 7567

Re: Spurious supercooling in ROMS-CICE simulation

Just a suggestion---sensitivity to the choice of the advection scheme may possibly indicate that your grid stretching is more than it should be. Could be worth trying: (1) decreasing the number of vertical levels, (2) use conservative stretching parameters, (3) increase the horizontal resolution. Ch...
by stlaur
Thu Oct 23, 2014 3:11 pm
Forum: ROMS Tools and Techniques
Topic: Offline implementation of ROMS advection-diffusion schemes
Replies: 2
Views: 1734

Re: Offline implementation of ROMS advection-diffusion schem

Yes, that would work for my purpose. I see the option in my own version of the code (v.~3.3) and I will give it a try. Thanks for the reply.
Pierre
by stlaur
Tue Oct 21, 2014 7:13 pm
Forum: ROMS Tools and Techniques
Topic: Offline implementation of ROMS advection-diffusion schemes
Replies: 2
Views: 1734

Offline implementation of ROMS advection-diffusion schemes

Hi, I have high-resolution outputs from ROMS and I would like to advect+diffuse an arbitrary 3-D scalar field using an `offline' implementation of ROMS advection-diffusion schemes. In other words, I'm looking for a Fortran/C/whatever program that takes as input the netcdf files, interpolates them in...
by stlaur
Mon Aug 30, 2010 7:56 pm
Forum: ROMS Problems
Topic: Question on subroutine netcdf_get_fvar from mod_netcdf.F
Replies: 11
Views: 3205

Re: Question on subroutine netcdf_get_fvar from mod_netcdf.F

I updated my code and it works great!
Thank you very much for providing this fix.

Best regards,

Pierre
by stlaur
Sun Aug 29, 2010 6:23 pm
Forum: ROMS Problems
Topic: Question on subroutine netcdf_get_fvar from mod_netcdf.F
Replies: 11
Views: 3205

Re: Question on subroutine netcdf_get_fvar from mod_netcdf.F

I would be very grateful if anybody that does not encounter this issue could describe what is happening in his roms code during the reading of boundary input data such as zeta_north(:,:). It can be simply a matter of inserting a couple of write(*,*) in your code to see what roms does. Just (please) ...
by stlaur
Wed Aug 25, 2010 8:42 pm
Forum: ROMS Problems
Topic: Question on subroutine netcdf_get_fvar from mod_netcdf.F
Replies: 11
Views: 3205

Re: Question on subroutine netcdf_get_fvar from mod_netcdf.F

I get the runtime error with: gfortran (gcc version 4.3.2, from Linux Debian Lenny). I also get the error on a different machine that uses: f95: Sun Fortran 95 8.3 SunOS_i386 Patch 127002-04 2008/04/16 and it says ****** FORTRAN RUN-TIME SYSTEM ****** An array expression does not conform : dim 1 is ...
by stlaur
Wed Aug 25, 2010 6:49 pm
Forum: ROMS Problems
Topic: Question on subroutine netcdf_get_fvar from mod_netcdf.F
Replies: 11
Views: 3205

Re: Question on subroutine netcdf_get_fvar from mod_netcdf.F

Thanks for this reply! Unfortunately, I checked with ncdump -h, and by opening the file in Matlab, and zeta_north correctly has two dimensions (zeta_time and xi_rho), with one record in time. Numerous users of ROMS use boundary files and this is the first time that this is reported. I agree that thi...
by stlaur
Wed Aug 25, 2010 2:15 pm
Forum: ROMS Problems
Topic: Question on subroutine netcdf_get_fvar from mod_netcdf.F
Replies: 11
Views: 3205

Re: Question on subroutine netcdf_get_fvar from mod_netcdf.F

Perhaps a minimal testcase would make things clearer: program testcase implicit none integer, parameter :: UBi = 403, LBi = -2, UBj = 1 ! Dummy values. integer :: ng, model, ncid, Trec, Iend, Istr character(400) :: ncfile, Vname real, dimension((UBi-LBi+1)*UBj) :: A ! Dummy values for testcase. ng =...
by stlaur
Sat Aug 21, 2010 3:05 pm
Forum: ROMS Problems
Topic: Question on subroutine netcdf_get_fvar from mod_netcdf.F
Replies: 11
Views: 3205

Re: Question on subroutine netcdf_get_fvar from mod_netcdf.F

Thanks for your help! ROMS is trying to get zeta_north(:,:) from my ocean_bry.nc file. My zeta_north(:,:) array has only one record in time, and the other dimension is space (size Lm+2). What happens then is that net_cdf_get_fvar is called over line 400 of get_ngfld: CALL netcdf_get_fvar (ng, mode, ...
by stlaur
Fri Aug 20, 2010 3:53 pm
Forum: ROMS Problems
Topic: Question on subroutine netcdf_get_fvar from mod_netcdf.F
Replies: 11
Views: 3205

Question on subroutine netcdf_get_fvar from mod_netcdf.F

Hello, I have a small question regarding subroutine netcdf_get_fvar. This routine is called within get_ngfld.F around line 400 (SVN rev.491M, very recent) to obtain a one-dimensional array called 'A(:)', yet the 'start' and 'total' arguments each have two elements (as if A(:) was two-dimensional). I...