Custom Query (986 matches)

When the NEMURO model (nemuro_inp.h) is activated, ROMS reports the wrong parameter (AlphaZL, VO2S0, KD2N) values to standard output.

Many thanks to Paul Mattern for ​bringing this to our attention.

I also added a local CPP flag DEBUG_VOLCONS to obc_volcons.F to track the parallel global reductions:

# ifdef DEBUG_VOLCONS
          WRITE (150,10) MyRank, iic(ng), iif(ng), bc_area, buffer(1),  &
                         bc_flux, buffer(2), buffer(2)/buffer(1)
  10      FORMAT (i3,1x,i3.3,1x,i3.3,5(1x,1pe23.15))
          CALL my_flush (150)
# endif

This was done to track the round-off of global parallel reductions (big numbers) with different partitions. The default is to have this CPP off.

Added the Red Tide Ecosystem Model for Alexandrium fundyense. It can be activated with CPP option RED_TIDE. This model is primarily designed for the Gulf of Maine.

As code it, the user needs to provide the bottom cyst concentration (cyst_bottom; cell/m2; 2D array) and the observed dissolved inorganic nutrient (DIN; micromoles; 3D array). The cyst_bottom and DIN fields can be provided in a separated forcing NetCDF file. See varinfo.dat for the metadata for these input variables.

Many thanks to Dennis McGillicudy and Ruoying He for providing us their code. This model was rewritten from scratch.

Also, many thanks to Sasha Shchepetkin for providing us his updated semi-Lagrangian code sinking and ascending used for vertical swimming.

References:

Stock, C.A., D.J. McGillicudy, A.R. Solow, and D.M. Anderson, 2005: Evaluating hypotheses for the initiation and development of Alexandrium fundyense blooms in the western Gulf of Maine using a couple physical-biological model, Deep-Sea Research II, 52, 2715-2744.

He, R., D.J. McGillicuddy, B.A. Keafer, and D.M. Anderson, 2008: Historic 2005 toxic bloom of Alexandrium fundyence in the western Gulf of Maine: 2, Coupled biophysical modeling, J. Geophys. Res., 113, C07040, doi:10.1029/2007JC004602

In read_phypar.F, lines 2736 through 2750 check the point source methodologies

!
!  Check if both point sources methodologies are activated.  Only one
!  method is allowed for a particular grid.  Otherwise, the point
!  source will be applies twice.
!
       IF (LuvSrc(ng).and.LwSrc(ng)) THEN
         IF (Master) THEN
           WRITE (out,260) 'LuvSrc', LuvSrc(ng),                        &
     &           'Because LwSrc  is also ON; only one method is legal.'
           WRITE (out,260) 'LwSrc', LwSrc(ng),                          &
     &           'Because LuvSrc is also ON; only one method is legal.'
           exit_flag=4
           RETURN
         END IF
       END IF

This needs to be enclosed in a loop over ng:

      DO ng=1,Ngrids
       IF (LuvSrc(ng).and.LwSrc(ng)) THEN
      ...
       END IF
      END DO

Otherwise, ng has the value it acquired on exit from the previous such loop, which leads to an out-of-bound access in LuvSrc(ng) and LwSrc(ng).