Hello everyone,
I'm trying to obtain the vertical eddy viscosity to calculate turbulent fluxes (w'u' = -Km du/dz or w'v' = -Km dv/dz ), however, I'm failing to obtain Km or the eddy viscosity to calculate the fluxes.
I have been reading the forum and I found some threads related to this topic (such as T1, T2, or T3), but I haven't found a straightforward solution for MYJ in those threads. Looking at the module_bl_myjpbl.F, it is clear that MYJ does not output EXCH_M (only EXCH_H which is for heat).
I have edited the Registry.EM to obtain different variables such as EXCH_M, XKMV, or XKHV, but all of them are zero. The only ones that are not zero are EXCH_H, XKMH, and XKHH (the last two are horizontal eddy viscosities). I think XKHV should be equal to EXCH_H, but I'm still missing EXCH_M.
I know that one solution is to edit the MYJ module to output EXCH_M as it does with EXCH_H. However, I have to create many variables in order to perform this calculation, which for me is strange because I have seen many publications where they use MYJ and they plot the turbulent fluxes. Most of them, don't mention that they edited the codes.
Am I doing something wrong? Is there any solution to obtain Km or EXCH_M with MYJ? This is the namelist that I'm using (My simulation is only until 1km resolution, I understand for values lower than that I have to avoid PBL schemes):
Thank you everyone in advance,
I'm trying to obtain the vertical eddy viscosity to calculate turbulent fluxes (w'u' = -Km du/dz or w'v' = -Km dv/dz ), however, I'm failing to obtain Km or the eddy viscosity to calculate the fluxes.
I have been reading the forum and I found some threads related to this topic (such as T1, T2, or T3), but I haven't found a straightforward solution for MYJ in those threads. Looking at the module_bl_myjpbl.F, it is clear that MYJ does not output EXCH_M (only EXCH_H which is for heat).
I have edited the Registry.EM to obtain different variables such as EXCH_M, XKMV, or XKHV, but all of them are zero. The only ones that are not zero are EXCH_H, XKMH, and XKHH (the last two are horizontal eddy viscosities). I think XKHV should be equal to EXCH_H, but I'm still missing EXCH_M.
I know that one solution is to edit the MYJ module to output EXCH_M as it does with EXCH_H. However, I have to create many variables in order to perform this calculation, which for me is strange because I have seen many publications where they use MYJ and they plot the turbulent fluxes. Most of them, don't mention that they edited the codes.
Am I doing something wrong? Is there any solution to obtain Km or EXCH_M with MYJ? This is the namelist that I'm using (My simulation is only until 1km resolution, I understand for values lower than that I have to avoid PBL schemes):
Code:
&domains
/
&physics
...
sf_sfclay_physics = 2, 2, 2, 2,2,2,
sf_surface_physics = 2, 2, 2, 2, 2,2,
bl_pbl_physics = 2, 2, 2, 2, 2,2,
bldt = 0, 0, 0, 0, 0,0,
...
sf_urban_physics = 0, 0, 0, 0,0,0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
sst_update = 1,
/
&fdda
/
&dynamics
...
/
&bdy_control
...
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
Thank you everyone in advance,
Last edited: