search

ACCESS-NRI / accessdev-Trac-archive

Archive accessdev Trac contents as issues
Apache License 2.0
0 stars 0 forks source link

Create UM land sea mask and surface ancillaries #383

Open penguian opened 3 years ago

penguian commented 3 years ago

type_non-CMIP6 | by mrd599@nci.org.au

UM land sea mask needs to match ocean model for flux conservation.

The strict requirement is that the UM has non-zero ocean fraction intersection with every ocean grid box from the ocean model.

Issue migrated from trac:383 at 2024-01-31 18:36:19 +1100

penguian commented 3 years ago

@martin.dix@anu.edu.au changed status from new to assigned

penguian commented 3 years ago

@martin.dix@anu.edu.au set owner to mrd599

penguian commented 3 years ago

@martin.dix@anu.edu.au set milestone to ACCESS-CM2-025

penguian commented 3 years ago

@martin.dix@anu.edu.au changed status from assigned to accepted

penguian commented 3 years ago

@martin.dix@anu.edu.au commented

Mask file is /g/data/ik11/inputs/access-om2/input_20201102/mom_025deg/ocean_mask.nc. This has a few points different to the mask that was used in the earlier CM2 0.025 degree configuration so coupling files need to be recreated.

Created ~access/access-cm2/input_O025/cpl_n96/oasis3_masks_N96_20201102.nc from oasis3_masks_N96_30102014.nc by replacing the cice.msk with the mask from /g/data/ik11/inputs/access-om2/input_20201102/mom_025deg/ocean_mask.nc.

Ran the coupled model suite u-ch054 with no mapping files so that OASIS generated them. Original files (e.g. ~access/access-cm2/input_O025/cpl_n96/rmp_cice_to_um1t_CONSERV_FRACNNEI.nc) moved to rmp_old subdirectory and new ones copied in place.

To calculate the UM land fraction, used the ESMF regridding weights generation described at https://github.com/aidanheerdegen/access_cm2_025_remap_weights with the conserve option, e.g. ESMF_RegridWeightGen -i -m conserve -s MOM025_scrip.nc -d um1t_scrip.nc -w rmp_mom025_um1t_conserve.nc These weights are unmasked. Apply these weights to the ocean mask to calculate a land fraction on the UM grid. Points with land fraction < 0.01 were set to ocean.

The resulting land fraction file has only 5 points different to that from the old version ~(access/data/ancil/access_cm2_n96e/O025/qrarm.landfrac) corresponding to the number of ocean mask changes, confirming that the procedure is valid. The UM land sea mask does not change.

UM ancillary generating suite is u-cd511. Used the orography source patched with the AUSDEM data over Australia. To use the correct land fraction rather than one created by the suite, copy qrparm.mask and qrparm.landfrac to /cylc-run/u-cd511/share/data/n96e_O025/land_sea_mask/etop01/ after the ancilMask task has run and created the directory.

Only the ancilVegfrac, ancilSmcsnow, ancilOrog tasks from the CAP group are required, as well as the ANTS group.

Programs for creating the CABLE vegetation ancillaries are in the repository https://trac.nci.org.au/svn/access_tools/cable_ancils/. This uses NCAR CMIP5 data from ~wan028/lulcc-ncar/pftlandusedyn.0.5x0.5.simyr1850-2005.c090630 on ruby.

penguian commented 3 years ago

@martin.dix@anu.edu.au changed _comment0 which not transferred by tractive

penguian commented 3 years ago

@martin.dix@anu.edu.au changed _comment1 which not transferred by tractive