questions about data structure

[arthur_panch@gmail_com]

Hi,I have two questions about the code of the WACCM model:Firstly, when the model is running, it transforms the space structure (longitude×latitude×vertical) to the structure (column×chunk×vertical). Now, I want to close the radiative heating over Tibet Plateau. Now I use one node and 40 CPUs to compute the model, and the horizontal resolution is 1.9 degree at latitude and 2.5 degree at longitude. I need to know the transformational relation between (chunk, column) and (latitude, longitude). Could you explain it for me? Thank you.Secondly, I want to output the chemical loss or produce rate of the reactants and products of 17 heterogeneous reactions on the surface of PSCs. I know the source file mo_prod_loss.F90 is used to compute the chemical loss and production. But I don’t understand the meaning of each dimension of the variables loss(:, :, :) and prod(:, :, :). Please explain it for me. Thank you.

 

[jedwards]

You can get the column lat and lon from the phys_grid.F90 module in cam physics.   Look at for example get_rlat_p.

 

[santos]

As an alternative, if you have a module with access to the physics state, you can access these as state%lon and state%lat. (These are in radians.)

 

[fvitt]

To answer the question on the dimensions of loss and prod in mo_prod_loss.F90 see mo_exp_sol.F90:    real(r8), dimension(ncol,pver,clscnt1) :: &
         prod, &
         loss, &

ncol is the number of physics columnspver is the number of vertical levelsclscnt1 is the number of species in the explicit solver class 

 

[arthur_panch@gmail_com]

thank for your help

 

[arthur_panch@gmail_com]

thank you.

 

[arthur_panch@gmail_com]

Thank you very much. I have solved this problem.

 

[r01229009@ntu_edu_tw]

Hello, May I ask you about how to do the "close the radiative heating " thing ??I also want to do the similar thing. I try close off the radiative heating in atmosphere, but I have no idea.Cloud you please give me some steps to do such simulation??  It would be very appreciated.