{"action":"created","issue":{"url":"https://api.github.com/repos/NOAA-GFDL/MOM6-examples/issues/77","repository_url":"https://api.github.com/repos/NOAA-GFDL/MOM6-examples","labels_url":"https://api.github.com/repos/NOAA-GFDL/MOM6-examples/issues/77/labels{/name}","comments_url":"https://api.github.com/repos/NOAA-GFDL/MOM6-examples/issues/77/comments","events_url":"https://api.github.com/repos/NOAA-GFDL/MOM6-examples/issues/77/events","html_url":"https://github.com/NOAA-GFDL/MOM6-examples/issues/77","id":144699209,"number":77,"title":"ocean analysis scripts do not work for 1/2 degree models","user":{"login":"nikizadehgfdl","id":6394244,"avatar_url":"https://avatars.githubusercontent.com/u/6394244?v=3","gravatar_id":"","url":"https://api.github.com/users/nikizadehgfdl","html_url":"https://github.com/nikizadehgfdl","followers_url":"https://api.github.com/users/nikizadehgfdl/followers","following_url":"https://api.github.com/users/nikizadehgfdl/following{/other_user}","gists_url":"https://api.github.com/users/nikizadehgfdl/gists{/gist_id}","starred_url":"https://api.github.com/users/nikizadehgfdl/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/nikizadehgfdl/subscriptions","organizations_url":"https://api.github.com/users/nikizadehgfdl/orgs","repos_url":"https://api.github.com/users/nikizadehgfdl/repos","events_url":"https://api.github.com/users/nikizadehgfdl/events{/privacy}","received_events_url":"https://api.github.com/users/nikizadehgfdl/received_events","type":"User","site_admin":false},"labels":[],"state":"closed","locked":false,"assignee":null,"milestone":null,"comments":2,"created_at":"2016-03-30T19:40:18Z","updated_at":"2016-04-08T19:50:33Z","closed_at":"2016-04-08T19:50:33Z","body":"The ocean analysis scripts were made for 1/4 degree ocean models and they are not working for 1/2 degree models.\r\nE.g., see\r\n```\r\n/home/Amy.Langenhorst/ulm_201510_awg_v20160212_mom6_2016.03.22/CM4_c96L32_am4g9_fullaero_2000_t1_OMp5_lmix_H5_nmle_ndiff_meke/gfdl.ncrc3-intel15-prod-openmp/scripts/analysis/ocn_monthly.frepp.0001-0005.printout\r\n\r\n\r\n/nbhome/Amy.Langenhorst/ulm_201510_awg_v20160212_mom6_2016.03.22/CM4_c96L32_am4g9_fullaero_2000_t1_OMp5_lmix_H5_nmle_ndiff_meke/mom6/tools/analysis/MLD_003.py -g /archive/gold/datasets/OM4_025/mosaic.v20140610.unpacked -o /nbhome/Amy.Langenhorst/ulm_201510_awg_v20160212_mom6_2016.03.22/CM4_c96L32_am4g9_fullaero_2000_t1_OMp5_lmix_H5_nmle_ndiff_meke/ocean_0001-0005/Hosoda_MLD -l 0001-0005 /archive/Amy.Langenhorst/ulm_201510_awg_v20160212_mom6_2016.03.22/CM4_c96L32_am4g9_fullaero_2000_t1_OMp5_lmix_H5_nmle_ndiff_meke/gfdl.ncrc3-intel15-prod-openmp/pp/ocean_monthly/ts/monthly/5yr//ocean_monthly.000101-000512.MLD_003.nc\r\n/nbhome/Amy.Langenhorst/ulm_201510_awg_v20160212_mom6_2016.03.22/CM4_c96L32_am4g9_fullaero_2000_t1_OMp5_lmix_H5_nmle_ndiff_meke/mom6/tools/analysis/m6plot.py:561: VisibleDeprecationWarning: boolean index did not match indexed array along dimension 0; dimension is 1080 but corresponding boolean dimension is 576\r\n  if not numpy.ma.getmask(s).any()==numpy.ma.nomask: weight[s.mask] = 0.\r\nTraceback (most recent call last):\r\n  File \"/nbhome/Amy.Langenhorst/ulm_201510_awg_v20160212_mom6_2016.03.22/CM4_c96L32_am4g9_fullaero_2000_t1_OMp5_lmix_H5_nmle_ndiff_meke/mom6/tools/analysis/MLD_003.py\", line 56, in <module>\r\n    save=cmdLineArgs.outdir+'/MLD_003_minimum.png')\r\n  File \"/nbhome/Amy.Langenhorst/ulm_201510_awg_v20160212_mom6_2016.03.22/CM4_c96L32_am4g9_fullaero_2000_t1_OMp5_lmix_H5_nmle_ndiff_meke/mom6/tools/analysis/m6plot.py\", line 67, in xyplot\r\n    sMin, sMax, sMean, sStd, sRMS = myStats(maskedField, area, debug=debug)\r\n  File \"/nbhome/Amy.Langenhorst/ulm_201510_awg_v20160212_mom6_2016.03.22/CM4_c96L32_am4g9_fullaero_2000_t1_OMp5_lmix_H5_nmle_ndiff_meke/mom6/tools/analysis/m6plot.py\", line 566, in myStats\r\n    mean = numpy.ma.sum(weight*s)/sumArea\r\n  File \"/usr/local/python/2.7.3/lib/python2.7/site-packages/numpy-1.11.0b3-py2.7-linux-x86_64.egg/numpy/ma/core.py\", line 3945, in __rmul__\r\n    return multiply(other, self)\r\n  File \"/usr/local/python/2.7.3/lib/python2.7/site-packages/numpy-1.11.0b3-py2.7-linux-x86_64.egg/numpy/ma/core.py\", line 966, in __call__\r\n    result = self.f(da, db, *args, **kwargs)\r\nValueError: operands could not be broadcast together with shapes (1080,1440) (576,720)\r\n\r\n```"},"comment":{"url":"https://api.github.com/repos/NOAA-GFDL/MOM6-examples/issues/comments/207575859","html_url":"https://github.com/NOAA-GFDL/MOM6-examples/issues/77#issuecomment-207575859","issue_url":"https://api.github.com/repos/NOAA-GFDL/MOM6-examples/issues/77","id":207575859,"user":{"login":"nikizadehgfdl","id":6394244,"avatar_url":"https://avatars.githubusercontent.com/u/6394244?v=3","gravatar_id":"","url":"https://api.github.com/users/nikizadehgfdl","html_url":"https://github.com/nikizadehgfdl","followers_url":"https://api.github.com/users/nikizadehgfdl/followers","following_url":"https://api.github.com/users/nikizadehgfdl/following{/other_user}","gists_url":"https://api.github.com/users/nikizadehgfdl/gists{/gist_id}","starred_url":"https://api.github.com/users/nikizadehgfdl/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/nikizadehgfdl/subscriptions","organizations_url":"https://api.github.com/users/nikizadehgfdl/orgs","repos_url":"https://api.github.com/users/nikizadehgfdl/repos","events_url":"https://api.github.com/users/nikizadehgfdl/events{/privacy}","received_events_url":"https://api.github.com/users/nikizadehgfdl/received_events","type":"User","site_admin":false},"created_at":"2016-04-08T19:50:33Z","updated_at":"2016-04-08T19:50:33Z","body":"John is right. They worked when I changed the xml to use OM4_05 examples directory."}}