Closed eguil closed 5 years ago
@eguil what are the dates of the erroneous file? This is what I am working with (latest)
-rwxrwxr-x. 1 401 2019-03-05 14:50:10.701365017 -0800 obs.EN4.historical.r0i0p0.an.ocn.Omon.density.ver-1.latestX.nc
-rwxrwxr-x. 1 52K 2019-03-05 14:49:27.540365647 -0800 obs.EN4.historical.r0i0p0.mo.ocn.Omon.areacello.ver-1.latestX.nc
-rwxrwxr-x. 1 20K 2019-02-13 17:06:14.827453366 -0800 obs.EN4.historical.r0i0p0.mo.ocn.Omon.so.ver-1.latestX.xml
-rwxrwxr-x. 1 20K 2019-02-13 17:06:17.549453326 -0800 obs.EN4.historical.r0i0p0.mo.ocn.Omon.thetao.ver-1.latestX.xml
-rwxrwxr-x. 1 49K 2019-03-05 14:49:25.276365680 -0800 obs.IPRC.historical.r0i0p0.mo.ocn.Omon.areacello.ver-1.latestX.nc
-rwxrwxr-x. 1 6.2K 2019-02-13 17:06:25.438453211 -0800 obs.IPRC.historical.r0i0p0.mo.ocn.Omon.so.ver-1.latestX.xml
-rwxrwxr-x. 1 6.2K 2019-02-13 17:06:26.390453197 -0800 obs.IPRC.historical.r0i0p0.mo.ocn.Omon.thetao.ver-1.latestX.xml
-rwxrwxr-x. 1 49K 2019-03-05 14:49:27.088365653 -0800 obs.JAMSTEC.historical.r0i0p0.mo.ocn.Omon.areacello.ver-1.latestX.nc
-rwxrwxr-x. 1 8.7K 2019-02-13 17:06:29.983453144 -0800 obs.JAMSTEC.historical.r0i0p0.mo.ocn.Omon.so.ver-1.latestX.xml
-rwxrwxr-x. 1 8.7K 2019-02-13 17:06:32.138453113 -0800 obs.JAMSTEC.historical.r0i0p0.mo.ocn.Omon.thetao.ver-1.latestX.xml
-rwxrwxr-x. 1 49K 2019-03-05 14:49:25.694365674 -0800 obs.SmithAndMurphy2007.historical.r0i0p0.mo.ocn.Omon.areacello.ver-1.latestX.nc
-rwxrwxr-x. 1 12K 2019-02-13 17:06:35.036453071 -0800 obs.SmithAndMurphy2007.historical.r0i0p0.mo.ocn.Omon.so.ver-1.latestX.xml
-rwxrwxr-x. 1 12K 2019-02-13 17:06:37.882453029 -0800 obs.SmithAndMurphy2007.historical.r0i0p0.mo.ocn.Omon.thetao.ver-1.latestX.xml
-rwxrwxr-x. 1 49K 2019-03-05 14:49:25.579365675 -0800 obs.UCSD.historical.r0i0p0.mo.ocn.Omon.areacello.ver-1.latestX.nc
-rwxrwxr-x. 1 8.6K 2019-02-13 17:06:41.139452982 -0800 obs.UCSD.historical.r0i0p0.mo.ocn.Omon.so.ver-1.latestX.xml
-rwxrwxr-x. 1 8.6K 2019-02-13 17:06:43.871452942 -0800 obs.UCSD.historical.r0i0p0.mo.ocn.Omon.thetao.ver-1.latestX.xmlobs.JAMSTEC.historical.r0i0p0.mo.ocn.Omon.thetao.ver-1.latestX.xml
I believed I copied this one: -rwxrwxr-x. 1 52K 2019-03-05 14:49:27.540365647 -0800 obs.EN4.historical.r0i0p0.mo.ocn.Omon.areacello.ver-1.latestX.nc
In [20]: fEN4 = '.../190213_data_density/obs.EN4.historical.r0i0p0.mo.ocn.Omon.areacello.ver-1.latestX.nc'
In [21]: fHEN = cdm.open(fEN4)
In [22]: areacelloEN4 = fHEN('areacello')
In [23]: areacelloEN4.sum()
...:
Out[24]: 1.1297774063073618e+17
This should resemble the WOA18 number of (see https://github.com/eguil/Density_bining/issues/80#issuecomment-488370264)
float basinmask3_area(latitude, longitude) ;
basinmask3_area:_FillValue = 1.e+20f ;
basinmask3_area:missing_value = 1.e+20f ;
basinmask3_area:earthLandAreaM2 = 148940000000000. ;
basinmask3_area:oceanSurfaceAreaM2 = 352342040224260. ; <- 352.3 x 10^12 m^2
basinmask3_area:earthWaterAreaM2 = 361132000000000. ;
basinmask3_area:earthSurfaceAreaM2 = 510100000000000. ; <- 5.1 x 10^14 m^2
basinmask3_area:units = "m^2" ;
so we clearly have a problem here, it seems the maps look right, just the units/magnitude has an issue:
EN4
WOA18
Revised EN4:
There is factor ~300 in the areacello of obs.EN4.historical.r0i0p0.mo.ocn.Omon.areacello.ver-1.latestX.nc