Open dabail10 opened 5 months ago
Is this already in python? @nusbaume @justin-richling
Does the ADF call the CVDP?
Here is some python code. Assuming that you have sea ice files uploaded:
ds_area = (ds.TAREA*ds.aice).where(TLAT>0.).sum(dim=['ni','nj'])
@dabail10 Yeah, we can call the CVDP from the ADF
Aha. I see that the ADF calls the CVDP NCL scripts.
Yeah, I should've mentions that. It runs in the background of the ADF via NCL scripts. I'm working with Adam Phillips to eventually get the major parts of the CVDP and CVDP-LE into python.
@dabail10 the current plan is for a python-only CVDP release to happen in the latter half of 2024. That release will not have all the current CVDP/CVDP-LE capabilities in it. I will speak with you at a later date about what sea ice metrics we should carry over into the new python CVDP package.
Got it. I still think we should work on this one step asap. I will prototype something in python.
@klindsay28 has a number of functions for doing this type of calculation.
This is something the CVDP does. It creates global and hemispheric integrated quantities including:
others?
For the sea ice we need:
This uses tarea from the sea ice model (constant gridcell area from sea ice history files or grid files).
Convert aice to fraction if it is not already.
NH extent = where(aice.ge.0.15,tarea,0.) where tlat > 0. NH area = sum(aicetarea) for tlat > 0. x 1.0e-12 NH ice volume = sum(hitarea) for tlat > 0. x 1.0e-13 NH snow volume = sum(hs*tarea) for tlat > 0. x 1.0e-13