Single Column Simulation with CCLM – in #8: General Questions

in #8: General Questions

Dear all,
I am planning to do single column simulations with CCLM 5.0.
Now I am wondering whether there is a maximum number of grid cells needed.

When I choose my domain for example to be 9×9 grid cells there comes the following error in INT2LM:
ERROR in tautsp(), src_vert_interpol.f90: while processing j-index:

When I extend the domain a little bit, I get the following error in INT2LM:
Caught signal 8 (Floating point exception)
mxm_handle_error()

And when I extend my domain up to more then 20×20 grid cells without changing anything else it works fine.
Is there any namelist variable I have to change when using a smaller domain or is this not possible with CCLM ?

Thank you very much for your help and best regards,
Eva

  @evanowatzki in #cdf57aa

Dear all,
I am planning to do single column simulations with CCLM 5.0.
Now I am wondering whether there is a maximum number of grid cells needed.

When I choose my domain for example to be 9×9 grid cells there comes the following error in INT2LM:
ERROR in tautsp(), src_vert_interpol.f90: while processing j-index:

When I extend the domain a little bit, I get the following error in INT2LM:
Caught signal 8 (Floating point exception)
mxm_handle_error()

And when I extend my domain up to more then 20×20 grid cells without changing anything else it works fine.
Is there any namelist variable I have to change when using a smaller domain or is this not possible with CCLM ?

Thank you very much for your help and best regards,
Eva

Single Column Simulation with CCLM

Dear all,
I am planning to do single column simulations with CCLM 5.0.
Now I am wondering whether there is a maximum number of grid cells needed.

When I choose my domain for example to be 9×9 grid cells there comes the following error in INT2LM:
ERROR in tautsp(), src_vert_interpol.f90: while processing j-index:

When I extend the domain a little bit, I get the following error in INT2LM:
Caught signal 8 (Floating point exception)
mxm_handle_error()

And when I extend my domain up to more then 20×20 grid cells without changing anything else it works fine.
Is there any namelist variable I have to change when using a smaller domain or is this not possible with CCLM ?

Thank you very much for your help and best regards,
Eva

View in channel

I am not aware of a single column version in the CLM -Community. Maybe there is one in the COSMO consortium?

  @burkhardtrockel in #e64448d

I am not aware of a single column version in the CLM -Community. Maybe there is one in the COSMO consortium?

I am not aware of a single column version in the CLM -Community. Maybe there is one in the COSMO consortium?

There is a single column model from Matthias Raschendorfer, but this is not public and I don’t think that you used that. Maybe you configured the COSMO -Model in a way that it runs on a very small domain (9×9 grid points). Or what did you do exactly?

  @ulrichschättler in #a379436

There is a single column model from Matthias Raschendorfer, but this is not public and I don’t think that you used that. Maybe you configured the COSMO -Model in a way that it runs on a very small domain (9×9 grid points). Or what did you do exactly?

There is a single column model from Matthias Raschendorfer, but this is not public and I don’t think that you used that. Maybe you configured the COSMO -Model in a way that it runs on a very small domain (9×9 grid points). Or what did you do exactly?

Yes, that is what I did. I run the model on a very small domain. But it only runs without errors when the domain is bigger than 20×20 grid points.

  @evanowatzki in #084e938

Yes, that is what I did. I run the model on a very small domain. But it only runs without errors when the domain is bigger than 20×20 grid points.

Yes, that is what I did. I run the model on a very small domain. But it only runs without errors when the domain is bigger than 20×20 grid points.

How do you control it? Is it an artificial data run? Or do you have initial and boundary data for it? Could you show us the namelist configuration ( YUSPECIF )?

  @ulrichschättler in #5a2fc18

How do you control it? Is it an artificial data run? Or do you have initial and boundary data for it? Could you show us the namelist configuration ( YUSPECIF )?

How do you control it? Is it an artificial data run? Or do you have initial and boundary data for it? Could you show us the namelist configuration ( YUSPECIF )?

I use ERA -Interim as initial data. Here the YUSPECIF of my simulation. Thank you very much

  @evanowatzki in #d91321f

I use ERA -Interim as initial data. Here the YUSPECIF of my simulation. Thank you very much

I use ERA -Interim as initial data. Here the YUSPECIF of my simulation. Thank you very much

Ok, I see one problem: You have ie_tot=je_tot=21 and nprocx=nprocy=24! I wonder that it works at all. Because this means that you distribute the 21×21=441 grid points to 24×24=576 processors.
You have to use nprocx=nprocy=1 (!), so only 1 CPU . Then it should also work with less grid points. You have nboundlines=3 (what you need for the dynamics), so can try to use ie_tot=je_tot=7. This is one grid point in the interior domain and 3 grid points of boundaries around.

I am sorry that I cannot help further, because this is my last day before I leave for 3 weeks of holidays. Good Luck!

  @ulrichschättler in #cb8361a

Ok, I see one problem: You have ie_tot=je_tot=21 and nprocx=nprocy=24! I wonder that it works at all. Because this means that you distribute the 21×21=441 grid points to 24×24=576 processors.
You have to use nprocx=nprocy=1 (!), so only 1 CPU . Then it should also work with less grid points. You have nboundlines=3 (what you need for the dynamics), so can try to use ie_tot=je_tot=7. This is one grid point in the interior domain and 3 grid points of boundaries around.

I am sorry that I cannot help further, because this is my last day before I leave for 3 weeks of holidays. Good Luck!

Ok, I see one problem: You have ie_tot=je_tot=21 and nprocx=nprocy=24! I wonder that it works at all. Because this means that you distribute the 21×21=441 grid points to 24×24=576 processors.
You have to use nprocx=nprocy=1 (!), so only 1 CPU . Then it should also work with less grid points. You have nboundlines=3 (what you need for the dynamics), so can try to use ie_tot=je_tot=7. This is one grid point in the interior domain and 3 grid points of boundaries around.

I am sorry that I cannot help further, because this is my last day before I leave for 3 weeks of holidays. Good Luck!

Thank you so much for your help. I will try it that way. I wish you happy holidays!

  @evanowatzki in #d31ac61

Thank you so much for your help. I will try it that way. I wish you happy holidays!

Thank you so much for your help. I will try it that way. I wish you happy holidays!