[MITgcm-support] Script for computing CFL conditions
Habib Toye Mahamadou Kele
habib.toyemahamadoukele at kaust.edu.sa
Thu Dec 3 08:37:14 EST 2015
Thanks Ryan for your reply,
Is 'advcfl_uvel_max' the same as the parameter S_a in equation (3.9) at
http://mitgcm.org/public/r2_manual/latest/online_documents/node115.html ?
I would like a script which will compute conditions like the ones given in:
http://mitgcm.org/public/r2_manual/latest/online_documents/node86.html#SECTION003211600000000000000
(equations (2.217) and (2.218))
http://mitgcm.org/public/r2_manual/latest/online_documents/node131.html#SECTION004122100000000000000
(from (3.43) to (3.48))
http://mitgcm.org/public/r2_manual/latest/online_documents/node115.html
(from (3.24) to (3.29))
http://mitgcm.org/public/r2_manual/latest/online_documents/node120.html
(from (3.6) to (3.9))
The result would look like:
(for the values at
http://mitgcm.org/public/r2_manual/latest/online_documents/node115.html )
NUMERICAL STABILITY PARAMETERS COMPUTATION FROM STDOUT.0000:
A_h = 400 ms^(-1)
delta_t = 1200 s
Delta_x = ...
Condition Name | Condition | User | Is user
condition fulfilled? |
numerical stability for inertial oscillations | S_l < 0.3 (S_l
= 4* A_h * delta_t / Delta_x^2) | 0.012 <
0.3 | Yes |
stability parameter to the horizontal Laplacian friction | S_i <
0.5 (S_i = f^2 * delta_t^2) | 0.0144 < 0.5 |
Yes |
...
Regards,
Habib
On Tue, Dec 1, 2015 at 5:18 PM, Ryan Abernathey <ryan.abernathey at gmail.com>
wrote:
> Hi Habib,
>
> MITgcm automatically computes the CFL for you as part of the monitor
> statistics.
>
> $ grep advcfl STDOUT.0000
> (PID.TID 0000.0001) %MON advcfl_uvel_max =
> 1.1796518160724E-01
> (PID.TID 0000.0001) %MON advcfl_vvel_max =
> 7.8377125552173E-02
> (PID.TID 0000.0001) %MON advcfl_wvel_max =
> 3.8263958212481E+00
> (PID.TID 0000.0001) %MON advcfl_W_hf_max =
> 1.8899738930714E+00
> ...
>
> Cheers,
> Ryan
>
> On Tue, Dec 1, 2015 at 8:47 AM, Habib Toye Mahamadou Kele <
> habib.toyemahamadoukele at kaust.edu.sa> wrote:
>
>> Greetings,
>>
>> Is there any script or routine for computing a mitgcmuv run CFL
>> conditions for example from the "output.txt" file obtained by "./mitgcmuv >
>> output.txt" ?
>> If not, could you please provide us with one which output would look like
>> the following table?
>>
>> Condition Name | Condition | User
>> | Is user condition fulfilled? |
>>
>> condition_1 | C1 <= 1.5 |
>> 5.4 <= 1.5 | No |
>>
>> condition_2 | var_2_1*var_2_2 / var_2_3 < 1 | 0.7 < 1
>> | Yes |
>>
>> ...
>>
>>
>> Thanks and Regards,
>> Habib TOYE
>>
>> ------------------------------
>> This message and its contents, including attachments are intended solely
>> for the original recipient. If you are not the intended recipient or have
>> received this message in error, please notify me immediately and delete
>> this message from your computer system. Any unauthorized use or
>> distribution is prohibited. Please consider the environment before printing
>> this email.
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>>
>>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
>
--
------------------------------
This message and its contents, including attachments are intended solely
for the original recipient. If you are not the intended recipient or have
received this message in error, please notify me immediately and delete
this message from your computer system. Any unauthorized use or
distribution is prohibited. Please consider the environment before printing
this email.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20151203/0fde1ef2/attachment.htm>
More information about the MITgcm-support
mailing list