[ecco-support] [EXTERNAL] CORRECTION: How to get the total diffusivity in ECCOv4
Paola Cessi
pcessi at ucsd.edu
Fri Jul 23 20:06:22 EDT 2021
Hi Ou,
I am having the same problem as Christopher, but with release 3: I downloaded files <https://ecco.jpl.nasa.gov/drive/files> / Version4 <https://ecco.jpl.nasa.gov/drive/files/Version4> / Release3 <https://ecco.jpl.nasa.gov/drive/files/Version4/Release3> / input_init <https://ecco.jpl.nasa.gov/drive/files/Version4/Release3/input_init> /xx_diffkr.0000000059.data <https://ecco.jpl.nasa.gov/drive/files/Version4/Release3/input_init/xx_diffkr.0000000059.data>, and the range of values is [-0.6 0.4]: too large to be vertical diffusivities (negative?). Can you please point me to the diffusivity files for release3?
Many thanks,
Paola
--------------------------------------------------------------------------------------------
Paola Cessi Tel: +1 858 534 0622
Scripps Institution of Oceanography Fax: +1 858 534 8045
9500 Gilman Drive #0213 e-mail: pcessi at ucsd.edu
La Jolla, CA 92093-0213 Web: http://pordlabs.ucsd.edu/pcessi
USA
Delivery address (Fed-Ex, DHL, UPS, etc.):
University of California San Diego
Keck Center - Room 366
8851 Shellback Way
La Jolla, CA 92093-0213, USA
> On Jun 4, 2021, at 11:40 AM, Wang, Ou (US 329B) <ou.wang at jpl.nasa.gov> wrote:
>
> Hi Christopher,
>
> One correction to the alternative way I mentioned;
> Alternatively, you can calculate it by summing up the first guess of a constant 1e-5 m^2/s (not total_diffkr_r009bit11.bin)and xx_diffkr.dim.0000000129.data.
>
> See Forget et al. (2015; doi:10.5194/gmd-8-3071-2015), Table 1 for the first guess of the mixing parameters.
>
> Ou
>
>
> From: <ecco-support-bounces at mit.edu <mailto:ecco-support-bounces at mit.edu>> on behalf of "Wang, Ou (US 329B)" <ou.wang at jpl.nasa.gov <mailto:ou.wang at jpl.nasa.gov>>
> Reply-To: "ECCO support list, wider membership" <ecco-support at mit.edu <mailto:ecco-support at mit.edu>>
> Date: Friday, June 4, 2021 at 9:17 AM
> To: "ECCO support list, wider membership" <ecco-support at mit.edu <mailto:ecco-support at mit.edu>>
> Subject: Re: [ecco-support] [EXTERNAL] How to get the total diffusivity in ECCOv4
>
> Hi Christopher,
>
> PO.DAAC is now hosting V4r4’s files. You can get V4r4’s dimensional vertical diffusivity field (first guess plus the control adjustment) from PO.DAAC at https://archive.podaac.earthdata.nasa.gov/podaac-ops-cumulus-protected/ECCO_L4_OCEAN_3D_MIX_COEFFS_LLC0090GRID_V4R4/OCEAN_3D_MIXING_COEFFS_ECCO_V4r4_native_llc0090.nc <https://urldefense.proofpoint.com/v2/url?u=https-3A__urldefense.us_v3_-5F-5Fhttps-3A_archive.podaac.earthdata.nasa.gov_podaac-2Dops-2Dcumulus-2Dprotected_ECCO-5FL4-5FOCEAN-5F3D-5FMIX-5FCOEFFS-5FLLC0090GRID-5FV4R4_OCEAN-5F3D-5FMIXING-5FCOEFFS-5FECCO-5FV4r4-5Fnative-5Fllc0090.nc-5F-5F-3B-21-21PvBDto6Hs4WbVuu7-21bO3-2DZfkkeOEujVIhL2nn3AEb9q1Zl86CHb1ZTrYDdah4HPKuSZ-5FF2v6ZPAv3JA9Giy4-24&d=DwMGaQ&c=-35OiAkTchMrZOngvJPOeA&r=YvEcUc3iPS8v6LXnkXPI-w&m=6Ilpubv4qn3OD41-8yhyk6uu5j1XJQRLMfnf7OSEpRQ&s=FPXRnuyQe_XQBb92q8_lgGGF7xNgSrdEC2S7NliIZnw&e=>. The variable is called “DIFFKR” in the NetCDF file. Alternatively, you can calculate it by summing up total_diffkr_r009bit11.bin and xx_diffkr.dim.0000000129.data that are available on ECCO Drive at https://ecco.jpl.nasa.gov/drive/files/Version4/Release4/input_init/total_diffkr_r009bit11.bin <https://urldefense.proofpoint.com/v2/url?u=https-3A__ecco.jpl.nasa.gov_drive_files_Version4_Release4_input-5Finit_total-5Fdiffkr-5Fr009bit11.bin&d=DwQGaQ&c=-35OiAkTchMrZOngvJPOeA&r=YvEcUc3iPS8v6LXnkXPI-w&m=6Ilpubv4qn3OD41-8yhyk6uu5j1XJQRLMfnf7OSEpRQ&s=At_qtxweGPSK3iWrt0MTdwFZh-Qh5o8DZk6kxGJy3qc&e=> andhttps://ecco.jpl.nasa.gov/drive/files/Version4/Release4/other/adjustments/xx_diffkr.dim.0000000129.data <https://urldefense.proofpoint.com/v2/url?u=https-3A__ecco.jpl.nasa.gov_drive_files_Version4_Release4_other_adjustments_xx-5Fdiffkr.dim.0000000129.data&d=DwMGaQ&c=-35OiAkTchMrZOngvJPOeA&r=YvEcUc3iPS8v6LXnkXPI-w&m=6Ilpubv4qn3OD41-8yhyk6uu5j1XJQRLMfnf7OSEpRQ&s=PoZK-6xokNwmLvGaaEYxY1OlakgA_2PpUzEKmntKQwA&e=>.
>
> The description about total_kapredi_r009bit11.bin in the ECCOv4 Python tutorial is incorrect and will be corrected. The field in total_kapredi_r009bit11.bin is the Redi parameter from Release 1, not the adjusted Redi parameter from Release 4.
>
> Ou
>
> From: <ecco-support-bounces at mit.edu <mailto:ecco-support-bounces at mit.edu>> on behalf of Christopher Wolfe <christopher.wolfe at stonybrook.edu <mailto:christopher.wolfe at stonybrook.edu>>
> Reply-To: "ECCO support list, wider membership" <ecco-support at mit.edu <mailto:ecco-support at mit.edu>>
> Date: Friday, June 4, 2021 at 6:44 AM
> To: "ecco-support at mit.edu <mailto:ecco-support at mit.edu>" <ecco-support at mit.edu <mailto:ecco-support at mit.edu>>
> Subject: [EXTERNAL] [ecco-support] How to get the total diffusivity in ECCOv4
>
> Hi all:
>
> I’m trying to figure out what the vertical diffusivity field in ECCOv4r4 is, but I’m a little confused about which input file actually contains this data.
>
> The ECCOv4 python tutorial (at https://ecco-v4-python-tutorial.readthedocs.io/ECCO_v4_Loading_LLC_compact_binary_files.html#Example-3:-Load-a-3D-'compact'-llc-binary-file-with-3rd-dimension-=-Depth <https://urldefense.proofpoint.com/v2/url?u=https-3A__urldefense.us_v3_-5F-5Fhttps-3A_ecco-2Dv4-2Dpython-2Dtutorial.readthedocs.io_ECCO-5Fv4-5FLoading-5FLLC-5Fcompact-5Fbinary-5Ffiles.html-2AExample-2D3-3A-2DLoad-2Da-2D3D-2D-27compact-27-2Dllc-2Dbinary-2Dfile-2Dwith-2D3rd-2Ddimension-2D-3D-2DDepth-5F-5F-3BIw-21-21PvBDto6Hs4WbVuu7-21eGOlWbFIVdKaXoIQQEcSSqRPb2UW2kabBqJXw8Hdg5AizGFceCu2pSW1rdXfRM0vvIA-24&d=DwMGaQ&c=-35OiAkTchMrZOngvJPOeA&r=YvEcUc3iPS8v6LXnkXPI-w&m=6Ilpubv4qn3OD41-8yhyk6uu5j1XJQRLMfnf7OSEpRQ&s=fxELJ2y1vojHpoBCo-jcZf5rct4ipSYukwbcHp3nRWM&e=>) says that total_kapredi_r009bit11.bin is the "adjusted GM redi parameter (first guess + adjustments)”, which implies that total_diffkr_r009bit11.bin should be the adjusted diffusivity.
>
> On the other hand, the README at https://ecco.jpl.nasa.gov/drive/files/Version4/Release4/input_init/README <https://urldefense.proofpoint.com/v2/url?u=https-3A__ecco.jpl.nasa.gov_drive_files_Version4_Release4_input-5Finit_README&d=DwMGaQ&c=-35OiAkTchMrZOngvJPOeA&r=YvEcUc3iPS8v6LXnkXPI-w&m=6Ilpubv4qn3OD41-8yhyk6uu5j1XJQRLMfnf7OSEpRQ&s=Y5km8rLOIAScTdlOx-H4gujN5dag4QnliKlbo22oRnY&e=> says that total_diffkr_r009bit11.bin is the "vert. diff. of release 1 (this field plus xx is the total)”, which means that I should add xx_diffkr.0000000129.data to this to get the full diffusivity. However, the values in this file are order 1, so adding it to total_diffkr_r009bit11 would produce diffusivities far outside the expected range.
>
>
>
> The gencost documentation of the ctrl package (https://mitgcm.readthedocs.io/en/latest/ocean_state_est/ocean_state_est.html#ctrl-model-parameter-adjustment-capability <https://urldefense.proofpoint.com/v2/url?u=https-3A__urldefense.us_v3_-5F-5Fhttps-3A_mitgcm.readthedocs.io_en_latest_ocean-5Fstate-5Fest_ocean-5Fstate-5Fest.html-2Actrl-2Dmodel-2Dparameter-2Dadjustment-2Dcapability-5F-5F-3BIw-21-21PvBDto6Hs4WbVuu7-21eGOlWbFIVdKaXoIQQEcSSqRPb2UW2kabBqJXw8Hdg5AizGFceCu2pSW1rdXfrjEs71w-24&d=DwMGaQ&c=-35OiAkTchMrZOngvJPOeA&r=YvEcUc3iPS8v6LXnkXPI-w&m=6Ilpubv4qn3OD41-8yhyk6uu5j1XJQRLMfnf7OSEpRQ&s=LhszeqkXgaIfNrgQxtvATYOLrdq9U9ivph736gNqDmA&e=>) says that the xx adjustments are non-dimensional and need to be scaled by one over the square root of the weights, which are in r2.wdiffkrFldv2.data. Using this produces adjustments on the order of 1e-6 to 1e-5, which seems about right. However, the data.ctrl namelist also specifies WC01 smoothing and some thresholding. Do these need to be applied as well?
>
>
>
> Is the best way to get this data simply to run the model for a bit and have it dump diffkr? It would be convenient to be able to calculate these quantities without having to run the model, however.
>
>
>
> Thanks in advance for your advise,
> Christopher
>
>
>
> ----------------------------------------------------------------------------
> Christopher L. Pitt Wolfe
> Associate Professor (Physical Oceanography)
> School of Marine and Atmospheric Sciences
> Stony Brook University
> christopher.wolfe at stonybrook.edu <mailto:christopher.wolfe at stonybrook.edu> 631-632-3152
> ----------------------------------------------------------------------------
>
>
>
> _______________________________________________
> ecco-support mailing list
> ecco-support at mit.edu <mailto:ecco-support at mit.edu>
> https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mit.edu_mailman_listinfo_ecco-2Dsupport&d=DwICAg&c=-35OiAkTchMrZOngvJPOeA&r=YvEcUc3iPS8v6LXnkXPI-w&m=6Ilpubv4qn3OD41-8yhyk6uu5j1XJQRLMfnf7OSEpRQ&s=T3TSbllVDDLo9752ro7_295rrHclwPhf-7vHqty_f4M&e= <https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mit.edu_mailman_listinfo_ecco-2Dsupport&d=DwICAg&c=-35OiAkTchMrZOngvJPOeA&r=YvEcUc3iPS8v6LXnkXPI-w&m=6Ilpubv4qn3OD41-8yhyk6uu5j1XJQRLMfnf7OSEpRQ&s=T3TSbllVDDLo9752ro7_295rrHclwPhf-7vHqty_f4M&e=>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.mit.edu/pipermail/ecco-support/attachments/20210723/9582828f/attachment-0001.html
More information about the ecco-support
mailing list