Detailed Description

Routines used to set up and use a set of (one for now) dynamically passive tracers in the ISOMIP configuration.

For now, just one passive tracer is injected in the sponge layer.

Data Types
type	isomip_tracer_cs
	ISOMIP tracer package control structure. More...

Functions/Subroutines
logical function, public	register_isomip_tracer (HI, GV, param_file, CS, tr_Reg, restart_CS)
	This subroutine is used to register tracer fields. More...

subroutine, public	initialize_isomip_tracer (restart, day, G, GV, h, diag, OBC, CS, ALE_sponge_CSp)
	Initializes the NTR tracer fields in tr(:,:,:,:) More...

subroutine, public	isomip_tracer_column_physics (h_old, h_new, ea, eb, fluxes, dt, G, GV, US, CS, evap_CFL_limit, minimum_forcing_depth)
	This subroutine applies diapycnal diffusion, including the surface boundary conditions and any other column tracer physics or chemistry to the tracers from this file. More...

subroutine, public	isomip_tracer_surface_state (sfc_state, h, G, CS)
	This subroutine extracts the surface fields from this tracer package that are to be shared with the atmosphere in coupled configurations. This particular tracer package does not report anything back to the coupler. More...

subroutine, public	isomip_tracer_end (CS)
	Deallocate any memory used by the ISOMIP tracer package. More...

Variables
integer, parameter	ntr = 1
	ntr is the number of tracers in this module.

Function/Subroutine Documentation

◆ initialize_isomip_tracer()

subroutine, public isomip_tracer::initialize_isomip_tracer	(	logical, intent(in)	restart,
		type(time_type), intent(in), target	day,
		type(ocean_grid_type), intent(in)	G,
		type(verticalgrid_type), intent(in)	GV,
		real, dimension( g %isd: g %ied, g %jsd: g %jed, g %ke), intent(in)	h,
		type(diag_ctrl), intent(in), target	diag,
		type(ocean_obc_type), pointer	OBC,
		type(isomip_tracer_cs), pointer	CS,
		type(ale_sponge_cs), pointer	ALE_sponge_CSp
	)

Initializes the NTR tracer fields in tr(:,:,:,:)

Parameters

[in]	g	Grid structure.
[in]	gv	The ocean's vertical grid structure.
[in]	restart	.true. if the fields have already been read from a restart file.
[in]	day	Time of the start of the run.
[in]	h	Layer thickness [H ~> m or kg m-2].
[in]	diag	A structure that is used to regulate diagnostic output.
	obc	This open boundary condition type specifies whether, where, and what open boundary conditions are used. This is not being used for now.
	cs	The control structure returned by a previous call to ISOMIP_register_tracer.
	ale_sponge_csp	A pointer to the control structure for the sponges, if they are in use. Otherwise this may be unassociated.

Definition at line 150 of file ISOMIP_tracer.F90.

 
   type(ocean_grid_type),                 intent(in) :: g !< Grid structure.
   type(verticalgrid_type),               intent(in) :: gv !< The ocean's vertical grid structure.
   logical,                               intent(in) :: restart !< .true. if the fields have already
                                                        !! been read from a restart file.
   type(time_type), target,               intent(in) :: day !< Time of the start of the run.
   real, dimension(SZI_(G),SZJ_(G),SZK_(G)), intent(in) :: h !< Layer thickness [H ~> m or kg m-2].
   type(diag_ctrl),               target, intent(in) :: diag !< A structure that is used to regulate
                                                        !! diagnostic output.
   type(ocean_obc_type),                  pointer    :: obc !< This open boundary condition type specifies
                                                        !! whether, where, and what open boundary conditions
                                                        !! are used. This is not being used for now.
   type(isomip_tracer_cs),                pointer    :: cs !< The control structure returned by a previous call
                                                        !! to ISOMIP_register_tracer.
   type(ale_sponge_cs),                   pointer    :: ale_sponge_csp !< A pointer to the control structure for
                                                        !! the sponges, if they are in use.  Otherwise this
                                                        !! may be unassociated.
 
   real, allocatable :: temp(:,:,:)
   real, pointer, dimension(:,:,:) :: &
     obc_tr1_u => null(), & ! These arrays should be allocated and set to
     obc_tr1_v => null()    ! specify the values of tracer 1 that should come
                            ! in through u- and v- points through the open
                            ! boundary conditions, in the same units as tr.
   character(len=16) :: name     ! A variable's name in a NetCDF file.
   character(len=72) :: longname ! The long name of that variable.
   character(len=48) :: units    ! The dimensions of the variable.
   character(len=48) :: flux_units ! The units for tracer fluxes, usually
                             ! kg(tracer) kg(water)-1 m3 s-1 or kg(tracer) s-1.
   real, pointer :: tr_ptr(:,:,:) => null()
   real :: h_neglect         ! A thickness that is so small it is usually lost
                             ! in roundoff and can be neglected [H ~> m or kg m-2].
   real :: e(szk_(g)+1), e_top, e_bot, d_tr
   integer :: i, j, k, is, ie, js, je, isd, ied, jsd, jed, nz, m
   integer :: isdb, iedb, jsdb, jedb
 
   if (.not.associated(cs)) return
   is = g%isc ; ie = g%iec ; js = g%jsc ; je = g%jec ; nz = g%ke
   isd = g%isd ; ied = g%ied ; jsd = g%jsd ; jed = g%jed
   isdb = g%IsdB ; iedb = g%IedB ; jsdb = g%JsdB ; jedb = g%JedB
   h_neglect = gv%H_subroundoff
 
   cs%Time => day
   cs%diag => diag
 
   if (.not.restart) then
     if (len_trim(cs%tracer_IC_file) >= 1) then
       !  Read the tracer concentrations from a netcdf file.
       if (.not.file_exists(cs%tracer_IC_file, g%Domain)) &
         call mom_error(fatal, "ISOMIP_initialize_tracer: Unable to open "// &
                         cs%tracer_IC_file)
       do m=1,ntr
         call query_vardesc(cs%tr_desc(m), name, caller="initialize_ISOMIP_tracer")
         call mom_read_data(cs%tracer_IC_file, trim(name), cs%tr(:,:,:,m), g%Domain)
       enddo
     else
       do m=1,ntr
         do k=1,nz ; do j=js,je ; do i=is,ie
           cs%tr(i,j,k,m) = 0.0
         enddo ; enddo ; enddo
       enddo
     endif
   endif ! restart
 
 ! the following does not work in layer mode yet
 !!  if ( CS%use_sponge ) then
   !   If sponges are used, this example damps tracers in sponges in the
   ! northern half of the domain to 1 and tracers in the southern half
   ! to 0.  For any tracers that are not damped in the sponge, the call
   ! to set_up_sponge_field can simply be omitted.
 !    if (.not.associated(ALE_sponge_CSp)) &
 !      call MOM_error(FATAL, "ISOMIP_initialize_tracer: "// &
 !        "The pointer to ALEsponge_CSp must be associated if SPONGE is defined.")
 
 !    allocate(temp(G%isd:G%ied,G%jsd:G%jed,nz))
 
 !    do j=js,je ; do i=is,ie
 !      if (G%geoLonT(i,j) >= 790.0 .AND. G%geoLonT(i,j) <= 800.0) then
 !        temp(i,j,:) = 1.0
 !      else
 !        temp(i,j,:) = 0.0
 !      endif
 !    enddo ; enddo
 
       !   do m=1,NTR
 !    do m=1,1
       ! This is needed to force the compiler not to do a copy in the sponge
       ! calls.  Curses on the designers and implementers of Fortran90.
 !      tr_ptr => CS%tr(:,:,:,m)
 !      call set_up_ALE_sponge_field(temp, G, tr_ptr, ALE_sponge_CSp)
 !    enddo
 !    deallocate(temp)
 !  endif
 

◆ isomip_tracer_column_physics()

subroutine, public isomip_tracer::isomip_tracer_column_physics	(	real, dimension(szi_(g),szj_(g),szk_(g)), intent(in)	h_old,
		real, dimension(szi_(g),szj_(g),szk_(g)), intent(in)	h_new,
		real, dimension(szi_(g),szj_(g),szk_(g)), intent(in)	ea,
		real, dimension(szi_(g),szj_(g),szk_(g)), intent(in)	eb,
		type(forcing), intent(in)	fluxes,
		real, intent(in)	dt,
		type(ocean_grid_type), intent(in)	G,
		type(verticalgrid_type), intent(in)	GV,
		type(unit_scale_type), intent(in)	US,
		type(isomip_tracer_cs), pointer	CS,
		real, intent(in), optional	evap_CFL_limit,
		real, intent(in), optional	minimum_forcing_depth
	)

This subroutine applies diapycnal diffusion, including the surface boundary conditions and any other column tracer physics or chemistry to the tracers from this file.

Parameters

[in]	g	The ocean's grid structure
[in]	gv	The ocean's vertical grid structure
[in]	h_old	Layer thickness before entrainment [H ~> m or kg m-2].
[in]	h_new	Layer thickness after entrainment [H ~> m or kg m-2].
[in]	ea	an array to which the amount of fluid entrained
[in]	eb	an array to which the amount of fluid entrained
[in]	fluxes	A structure containing pointers to thermodynamic and tracer forcing fields. Unused fields have NULL ptrs.
[in]	dt	The amount of time covered by this call [T ~> s]
[in]	us	A dimensional unit scaling type
	cs	The control structure returned by a previous call to ISOMIP_register_tracer.
[in]	evap_cfl_limit	Limit on the fraction of the water that can be fluxed out of the top layer in a timestep [nondim]
[in]	minimum_forcing_depth	The smallest depth over which fluxes can be applied [H ~> m or kg m-2]

Definition at line 250 of file ISOMIP_tracer.F90.

   type(ocean_grid_type),   intent(in) :: g    !< The ocean's grid structure
   type(verticalgrid_type), intent(in) :: gv   !< The ocean's vertical grid structure
   real, dimension(SZI_(G),SZJ_(G),SZK_(G)), &
                            intent(in) :: h_old !< Layer thickness before entrainment [H ~> m or kg m-2].
   real, dimension(SZI_(G),SZJ_(G),SZK_(G)), &
                            intent(in) :: h_new !< Layer thickness after entrainment [H ~> m or kg m-2].
   real, dimension(SZI_(G),SZJ_(G),SZK_(G)), &
                            intent(in) :: ea   !< an array to which the amount of fluid entrained
                                               !! from the layer above during this call will be
                                               !! added [H ~> m or kg m-2].
   real, dimension(SZI_(G),SZJ_(G),SZK_(G)), &
                            intent(in) :: eb   !< an array to which the amount of fluid entrained
                                               !! from the layer below during this call will be
                                               !! added [H ~> m or kg m-2].
   type(forcing),           intent(in) :: fluxes !< A structure containing pointers to thermodynamic
                                               !! and tracer forcing fields.  Unused fields have NULL ptrs.
   real,                    intent(in) :: dt   !< The amount of time covered by this call [T ~> s]
   type(unit_scale_type),   intent(in) :: us   !< A dimensional unit scaling type
   type(isomip_tracer_cs),  pointer    :: cs !< The control structure returned by a previous
                                               !! call to ISOMIP_register_tracer.
   real,          optional, intent(in) :: evap_cfl_limit !< Limit on the fraction of the water that can
                                               !! be fluxed out of the top layer in a timestep [nondim]
   real,          optional, intent(in) :: minimum_forcing_depth !< The smallest depth over which
                                               !! fluxes can be applied [H ~> m or kg m-2]
 
 ! The arguments to this subroutine are redundant in that
 !     h_new(k) = h_old(k) + ea(k) - eb(k-1) + eb(k) - ea(k+1)
 
   ! Local variables
   real :: b1(szi_(g))          ! b1 and c1 are variables used by the
   real :: c1(szi_(g),szk_(g))  ! tridiagonal solver.
   real, dimension(SZI_(G),SZJ_(G),SZK_(G)) :: h_work ! Used so that h can be modified
   real :: melt(szi_(g),szj_(g)) ! melt water (positive for melting, negative for freezing) [R Z T-1 ~> kg m-2 s-1]
   real :: mmax                ! The global maximum melting rate [R Z T-1 ~> kg m-2 s-1]
   character(len=256) :: mesg  ! The text of an error message
   integer :: i, j, k, is, ie, js, je, nz, m
   is = g%isc ; ie = g%iec ; js = g%jsc ; je = g%jec ; nz = g%ke
 
   if (.not.associated(cs)) return
 
   melt(:,:) = fluxes%iceshelf_melt(:,:)
 
   ! max. melt
   mmax = maxval(melt(is:ie,js:je))
   call max_across_pes(mmax)
   ! write(mesg,*) 'max melt = ', mmax
   ! call MOM_mesg(mesg, 5)
   ! dye melt water (m=1), dye = 1 if melt=max(melt)
   do m=1,ntr
     do j=js,je ; do i=is,ie
       if (melt(i,j) > 0.0) then ! melting
         cs%tr(i,j,1:2,m) = melt(i,j)/mmax ! inject dye in the ML
       else ! freezing
         cs%tr(i,j,1:2,m) = 0.0
       endif
     enddo ; enddo
   enddo
 
   if (present(evap_cfl_limit) .and. present(minimum_forcing_depth)) then
     do m=1,ntr
       do k=1,nz ;do j=js,je ; do i=is,ie
         h_work(i,j,k) = h_old(i,j,k)
       enddo ; enddo ; enddo
       call applytracerboundaryfluxesinout(g, gv, cs%tr(:,:,:,m), dt, fluxes, h_work, &
                                           evap_cfl_limit, minimum_forcing_depth)
       call tracer_vertdiff(h_work, ea, eb, dt, cs%tr(:,:,:,m), g, gv)
     enddo
   else
     do m=1,ntr
       call tracer_vertdiff(h_old, ea, eb, dt, cs%tr(:,:,:,m), g, gv)
     enddo
   endif
 

◆ isomip_tracer_end()

subroutine, public isomip_tracer::isomip_tracer_end ( type(isomip_tracer_cs), pointer CS )

Deallocate any memory used by the ISOMIP tracer package.

Parameters

cs	The control structure returned by a previous call to ISOMIP_register_tracer.

Definition at line 360 of file ISOMIP_tracer.F90.

   type(isomip_tracer_cs), pointer :: cs !< The control structure returned by a previous
                                         !! call to ISOMIP_register_tracer.
   integer :: m
 
   if (associated(cs)) then
     if (associated(cs%tr)) deallocate(cs%tr)
     deallocate(cs)
   endif

◆ isomip_tracer_surface_state()

subroutine, public isomip_tracer::isomip_tracer_surface_state	(	type(surface), intent(inout)	sfc_state,
		real, dimension( g %isd: g %ied, g %jsd: g %jed, g %ke), intent(in)	h,
		type(ocean_grid_type), intent(in)	G,
		type(isomip_tracer_cs), pointer	CS
	)

This subroutine extracts the surface fields from this tracer package that are to be shared with the atmosphere in coupled configurations. This particular tracer package does not report anything back to the coupler.

Parameters

[in]	g	The ocean's grid structure.
[in,out]	sfc_state	A structure containing fields that describe the surface state of the ocean.
[in]	h	Layer thickness [H ~> m or kg m-2].
	cs	The control structure returned by a previous call to ISOMIP_register_tracer.

Definition at line 329 of file ISOMIP_tracer.F90.

   type(ocean_grid_type),  intent(in)    :: g  !< The ocean's grid structure.
   type(surface),          intent(inout) :: sfc_state !< A structure containing fields that
                                               !! describe the surface state of the ocean.
   real, dimension(SZI_(G),SZJ_(G),SZK_(G)), &
                           intent(in)    :: h  !< Layer thickness [H ~> m or kg m-2].
   type(isomip_tracer_cs), pointer       :: cs !< The control structure returned by a previous
                                               !! call to ISOMIP_register_tracer.
 
   ! This particular tracer package does not report anything back to the coupler.
   ! The code that is here is just a rough guide for packages that would.
 
   integer :: m, is, ie, js, je, isd, ied, jsd, jed
   is = g%isc ; ie = g%iec ; js = g%jsc ; je = g%jec
   isd = g%isd ; ied = g%ied ; jsd = g%jsd ; jed = g%jed
 
   if (.not.associated(cs)) return
 
   if (cs%coupled_tracers) then
     do m=1,ntr
       !   This call loads the surface values into the appropriate array in the
       ! coupler-type structure.
       call coupler_type_set_data(cs%tr(:,:,1,m), cs%ind_tr(m), ind_csurf, &
                    sfc_state%tr_fields, idim=(/isd, is, ie, ied/), &
                    jdim=(/jsd, js, je, jed/) )
     enddo
   endif
 

◆ register_isomip_tracer()

logical function, public isomip_tracer::register_isomip_tracer	(	type(hor_index_type), intent(in)	HI,
		type(verticalgrid_type), intent(in)	GV,
		type(param_file_type), intent(in)	param_file,
		type(isomip_tracer_cs), pointer	CS,
		type(tracer_registry_type), pointer	tr_Reg,
		type(mom_restart_cs), pointer	restart_CS
	)

This subroutine is used to register tracer fields.

Parameters

[in]	hi	A horizontal index type structure.
[in]	gv	The ocean's vertical grid structure.
[in]	param_file	A structure indicating the open file to parse for model parameter values.
	cs	A pointer that is set to point to the control structure for this module (in/out).
	tr_reg	A pointer to the tracer registry.
	restart_cs	A pointer to the restart control structure.

Definition at line 69 of file ISOMIP_tracer.F90.

   type(hor_index_type),      intent(in) :: hi    !<A horizontal index type structure.
   type(verticalgrid_type),    intent(in) :: gv   !< The ocean's vertical grid structure.
   type(param_file_type),      intent(in) :: param_file !< A structure indicating the open file
                                                        !! to parse for model parameter values.
   type(isomip_tracer_cs),     pointer    :: cs !<A pointer that is set to point to the control
                                                        !! structure for this module (in/out).
   type(tracer_registry_type), pointer    :: tr_reg !<A pointer to the tracer registry.
   type(mom_restart_cs),       pointer    :: restart_cs !<A pointer to the restart control structure.
 
   character(len=80)  :: name, longname
 ! This include declares and sets the variable "version".
 #include "version_variable.h"
   character(len=40)  :: mdl = "ISOMIP_tracer" ! This module's name.
   character(len=200) :: inputdir
   character(len=48)  :: flux_units ! The units for tracer fluxes, usually
                             ! kg(tracer) kg(water)-1 m3 s-1 or kg(tracer) s-1.
   real, pointer :: tr_ptr(:,:,:) => null()
   logical :: register_isomip_tracer
   integer :: isd, ied, jsd, jed, nz, m
   isd = hi%isd ; ied = hi%ied ; jsd = hi%jsd ; jed = hi%jed ; nz = gv%ke
 
   if (associated(cs)) then
     call mom_error(warning, "ISOMIP_register_tracer called with an "// &
                             "associated control structure.")
     return
   endif
   allocate(cs)
 
   ! Read all relevant parameters and write them to the model log.
   call log_version(param_file, mdl, version, "")
   call get_param(param_file, mdl, "ISOMIP_TRACER_IC_FILE", cs%tracer_IC_file, &
                  "The name of a file from which to read the initial "//&
                  "conditions for the ISOMIP tracers, or blank to initialize "//&
                  "them internally.", default=" ")
   if (len_trim(cs%tracer_IC_file) >= 1) then
     call get_param(param_file, mdl, "INPUTDIR", inputdir, default=".")
     inputdir = slasher(inputdir)
     cs%tracer_IC_file = trim(inputdir)//trim(cs%tracer_IC_file)
     call log_param(param_file, mdl, "INPUTDIR/ISOMIP_TRACER_IC_FILE", &
                    cs%tracer_IC_file)
   endif
   call get_param(param_file, mdl, "SPONGE", cs%use_sponge, &
                  "If true, sponges may be applied anywhere in the domain. "//&
                  "The exact location and properties of those sponges are "//&
                  "specified from MOM_initialization.F90.", default=.false.)
 
   allocate(cs%tr(isd:ied,jsd:jed,nz,ntr)) ; cs%tr(:,:,:,:) = 0.0
 
   do m=1,ntr
     if (m < 10) then ; write(name,'("tr_D",I1.1)') m
     else ; write(name,'("tr_D",I2.2)') m ; endif
     write(longname,'("Concentration of ISOMIP Tracer ",I2.2)') m
     cs%tr_desc(m) = var_desc(name, units="kg kg-1", longname=longname, caller=mdl)
     if (gv%Boussinesq) then ; flux_units = "kg kg-1 m3 s-1"
     else ; flux_units = "kg s-1" ; endif
 
     ! This is needed to force the compiler not to do a copy in the registration
     ! calls.  Curses on the designers and implementers of Fortran90.
     tr_ptr => cs%tr(:,:,:,m)
     ! Register the tracer for horizontal advection, diffusion, and restarts.
     call register_tracer(tr_ptr, tr_reg, param_file, hi, gv, &
                          name=name, longname=longname, units="kg kg-1", &
                          registry_diags=.true., flux_units=flux_units, &
                          restart_cs=restart_cs)
 
     !   Set coupled_tracers to be true (hard-coded above) to provide the surface
     ! values to the coupler (if any).  This is meta-code and its arguments will
     ! currently (deliberately) give fatal errors if it is used.
     if (cs%coupled_tracers) &
       cs%ind_tr(m) = aof_set_coupler_flux(trim(name)//'_flux', &
           flux_type=' ', implementation=' ', caller="register_ISOMIP_tracer")
   enddo
 
   cs%tr_Reg => tr_reg
   register_isomip_tracer = .true.

Detailed Description

Data Types

Functions/Subroutines

Variables

Function/Subroutine Documentation

◆ initialize_isomip_tracer()

◆ isomip_tracer_column_physics()

◆ isomip_tracer_end()

◆ isomip_tracer_surface_state()

◆ register_isomip_tracer()