Daymet Mask NaN issue in the Watershed Workflow

[Zexuan Xu]

HI ATS/Watershed Workflow users,

This is a question for Watershed Workflow but I guess it would be appropriate to post it here

I also submit an issue on https://github.com/environmental-modeling-workflows/watershed-workflow/issues/51

Basically, it seems that when trying to write the daymet data in the format that ATS can read in the get_Daymet.ipynb, in the last cell of watershed_workflow.daymet.daymet_to_daily_averages, the mask NaN does not work on my end, and the output still contains -9999, or other negative values (e.g., -9725.85 for temperature) that makes them a problem for ATS simulations.

Looking into the code I feel this line (https://github.com/environmental-modeling-workflows/watershed-workflow/blob/ac430e2f3de2fe239e1000a02feb708f14d8e817/watershed_workflow/daymet.py#L42) suppose to change all -9999 to np.NaN but it doesn't for some reason. 

Here is an output of temperature in the Coweeta example

[[-9725.85       -9725.85       -9725.85       -9725.85

  -9725.85       -9725.85         266.32999995 -9725.85

  -9725.85      ]

 [-9725.85         262.34000006   264.08999982   265.47999992

    265.40499988   265.06499996   264.8599998  -9725.85

  -9725.85      ]

 [  262.96499982   262.29499998   261.76500025   262.48000016

    263.8599998    264.61999979   264.41500034   266.74000003

  -9725.85      ]

 [  262.62500014   261.83499994   261.69999971   261.83499994

    264.57499995   265.23500004   264.86499991   266.35000005

  -9725.85      ]

 [-9725.85         261.36000004   261.05499973   264.45000019

    265.78         266.11499979   267.27999988   267.79499992

  -9725.85      ]

 [-9725.85         261.47499981   261.7300004    264.7599999

    266.02500024   267.02499988   267.60499986   268.01999995

  -9725.85      ]

 [-9725.85         261.07500019   263.37000003   264.15999999

    266.05000021   267.03000011   267.20999994   267.9

    268.08999988]

 [-9725.85         263.07000008   264.16499987   265.32499971

    264.61000028   266.48000004   266.98000004   267.4999999

    268.05999991]

 [-9725.85         265.25500002   265.48500028   266.36999979

    265.60000029   265.02000012   265.95000007   267.43499979

  -9725.85      ]

 [-9725.85       -9725.85         266.91499999   266.4450002

  -9725.85       -9725.85       -9725.85       -9725.85

  -9725.85      ]]

I was wondering if this is only an issue on my end?

Thanks

Zexuan  

--

Zexuan Xu, Ph.D

Climate & Ecosystem Sciences Division

Lawrence Berkeley National Laboratory

1 Cyclotron Road, MS 74-316C

Berkeley, CA 94720

Website:

http://eesa.lbl.gov/profiles/zexuan-xu/

[Coon, Ethan]

Let’s discuss this on the ticket – definitely a WW bug!

Ethan

 

-- 

-------------------------------------------------------------------------

Ethan Coon

Senior Research Scientist

Oak Ridge National Laboratory

 

865-241-1296

https://www.ornl.gov/staff-profile/ethan-t-coon

-------------------------------------------------------------------------

 

--
You received this message because you are subscribed to the Google Groups "Amanzi-ATS Users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to ats-users+...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/ats-users/CAFie-4yxf-QKWpLAr05NNz_pX0Y4Q3Xuw-suFX495b2Gi9LNYg%40mail.gmail.com.