Changes between Initial Version and Version 1 of Ticket #870


Ignore:
Timestamp:
10/16/20 00:06:43 (4 years ago)
Author:
arango
Comment:

Legend:

Unmodified
Added
Removed
Modified

  • Ticket #870

    • Property Resolution → Done
    • Property Status new → closed

  • Ticket #870 – Description

    initial v1  
    1 Several changes were introduced src:ticket:869 about how to process the surface tracer fluxes that are used as boundary conditions to the vertical diffusion term. A couple of additional changes are made here to remove any ambiguity about the surface freshwater flux needed for ROMS.
     1Several changes were introduced in src:ticket:869 about how to process the surface/bottom tracer fluxes that are used as top/bottom boundary conditions to the vertical diffusion terms. A couple of additional changes are made here to remove any ambiguity about the surface freshwater flux needed for ROMS.
    22
    33Notice that ROMS has the following vertical surface/bottom boundary condition for temperature and salinity:
    … …  
    55[[Image(http://www.myroms.org/trac/stflux.png, center, 600)]]
    66
    7 In particular, the blue square shows what ROMS needs as a freshwater flux.  Usually, '''E-P''' is units of '''kilogram meter-2 second-1'''.  Then, the flux is divided by freshwater density, which converts to '''meter second-1'''.  In routine '''set_vbc.F''' such flux is multiplied by the surface salinity.  Since the salinity does not have [https://www.myroms.org/forum/viewtopic.php?t=294 physical units], the vertical boundary condition is still in '''meter second-1'''.
     7In particular, the blue square shows what ROMS needs as a freshwater flux.  Usually, '''E-P''' is in units of '''kilogram meter-2 second-1'''.  Then, the flux is divided by freshwater density, which converts to '''meter second-1'''.  In routine '''set_vbc.F''' such flux is multiplied by the surface salinity.  Since the salinity does not have [https://www.myroms.org/forum/viewtopic.php?t=294 physical units], the vertical boundary condition is still in '''meter second-1'''.
    88
    99I cleaned the metadata in '''varinfo.dat''', now we have:
    … …  
    3131  'meter second-1'                                 ! computed by NLM ROMS
    3232  'EminusP, scalar, series'                        ! bulk_flux.F or coupling
    33   'ocean_time'                                     ! needed in adjoint-based
    34   'idEmPf'                                         ! applications
     33  'ocean_time'
     34  'idEmPf'
    3535  'r2dvar'
    3636  1.0d0
    3737}}}
    3838
    39 As it is specified in the preamble of '''varinfo.dat''', the '''PSS''' is added on purpose to indicate that the Practical Salinity Scale was use to determine condutivity but it not a physical unit.  We are using the latest UNESCO equation of state.
     39As specified in the preamble of '''varinfo.dat''', the '''PSS''' information is added to indicate that the Practical Salinity Scale was used to determine conductivity, but it is not a physical unit.  We are not using the latest UNESCO equation of state.
    4040
    41 We have two different variables for input freshwater flux: '''swflux''' and '''EminusP'''. Historically, '''swflux''' has been used for available datasets while '''EminusP''' has been used when derived from parameterizations ('''BULK_FLUXES''' and atmospheric coupling).  To avoid any confusion, the units and scale for '''swflux''' were changed to '''meter second-1''' and '''1.0d0''' respectively. It is the standard units.  If such field is in '''centimeter day-1''', I recomend to add the '''scale_factor''' attribute to the NetCDF, so the conversion is done during reading in ROMS:
     41We have two different variables for input freshwater flux: '''swflux''' and '''EminusP'''. Historically, '''swflux''' has been used for available datasets while '''EminusP''' has been used when derived from parameterizations ('''BULK_FLUXES''' and atmospheric coupling).  To avoid any confusion, the units and scale for '''swflux''' were changed to '''meter second-1''' and '''1.0d0''', respectively. It is the standard unit for freshwater flux.  If such field is in '''centimeter day-1''', I recomend to add the '''scale_factor''' attribute to the NetCDF, so the conversion is done during reading in ROMS:
    4242{{{
    4343        float swflux(swf_time, eta_rho, xi_rho) ;
    … …  
    5050                swflux:scale_factor = 1.157407e-07f ;
    5151}}}
     52We could use the Matlab script to adjust the NetCDF file:
     53{{{
     54nc_attadd('my_file.nc','scale_factor',0.01/86400,'swflux')
     55}}}
     56
     57----
    5258
    5359* Corrected '''get_data.F''' and '''set_data.F''' processing of the freshwater flux when the '''SCORRECTION''' is activated.
    5460
    55 * Corrected issue in '''regrid.F'''.  It was missing the '''MASKING''' conditional.
     61* Corrected compiling issue in '''regrid.F'''.  It was missing the '''MASKING''' conditional.